1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<Scenery>
<Metric kind = "RezzollaZhidenko">
<Mass unit="sunmass">4.31e6</Mass>
<Epsilon>0.01</Epsilon> Horizon parameter rH=2/(1+eps)
<AParam>0.01 0.01</AParam> RZ14 a-parameters [a0,a1,a2,a3];
are zero if not provided
<BParam>0.01 0.01</BParam> same for RZ14 b-parameters
</Metric>
Compare with Schwarzschild metric:
<!-- <Metric kind = "KerrBL"> -->
<!-- <Mass unit="sunmass">4.31e6</Mass> -->
<!-- <Spin>0.</Spin> -->
<!-- </Metric> -->
<Screen>
<Distance unit="kpc"> 8.33 </Distance>
Time of start of ray tracing
<Time unit="yr"> 100e3 </Time>
<FieldOfView unit="µas"> 300 </FieldOfView>
This inclination is 180° minus the standard disk
inclination (angle between spin vector and line of sight)
<Inclination unit="°"> 95 </Inclination>
<PALN> 3.1415 </PALN>
<Argument> -1.5708 </Argument>
<Resolution>30</Resolution>
Compute a spectrum with 50 frequencies from 10^10 Hz to 10^18 Hz
<Spectrometer kind="freqlog" nsamples="5">10. 18.</Spectrometer>
</Screen>
<Quantities>Spectrum</Quantities>
<Astrobj kind = "PolishDoughnut">
Compute Bremsstrahlung spectrum?
<!--<Bremsstrahlung/>-->
Compute non-thermal synchrotron spectrum
with a fraction delta of non-thermal electrons
and a power-law exponent?
<!--<NonThermalDeltaExpo>0.01 3.</NonThermalDeltaExpo>-->
The default is to compute only the thermal synchrotron spectrum.
Compute synchrotron averaged over angle between
Bfield and direciton of emission?
<AngleAveraged/>
Constant angular momentum -u_phi / u_t, and inner radius
<AngMomRinner> 4. 6.2 </AngMomRinner>
Typical SgrA* enthalpy hc = 9000 cgs; ne [cm^{-3}] ~ 700 * hc [cgs]
<CentralEnthalpyPerUnitVolume unit="erg/cm3">
9000
</CentralEnthalpyPerUnitVolume>
Typical SgrA* central temperature: 5e10K
<CentralTemperature>5.3e10</CentralTemperature>
beta = p_gas / p_mag
<Beta>10</Beta>
<OpticallyThin/>
RMax: caution, make sure the torus is smaller than that!
<RMax>100.</RMax>
</Astrobj>
Computation stops for time less than following value
to prevent infinite integration
<MinimumTime> -1e10 </MinimumTime>
Do we want to use multi-threading?
<NThreads> 1 </NThreads>
</Scenery>
|
