| 1 | title NS94 NB parameters and model
|
|---|
| 2 | c
|
|---|
| 3 | c This version creates a grid of models for inclusion in xspec
|
|---|
| 4 | c
|
|---|
| 5 | c commands controlling continuum =======
|
|---|
| 6 | c Interpolate input spectrum from Jon's for power law
|
|---|
| 7 | c extending down to 0.1 keV
|
|---|
| 8 | interpolate (0.00001 -10)
|
|---|
| 9 | continue ( 7.63227 -7.12709) ( 87.8176 -7.95081)
|
|---|
| 10 | continue ( 89.9845 -7.94449) ( 92.2610 -7.93789)
|
|---|
| 11 | continue ( 94.6559 -7.93101) ( 97.1787 -7.92385)
|
|---|
| 12 | continue ( 99.8390 -7.91642) ( 102.649 -7.90872)
|
|---|
| 13 | continue ( 105.623 -7.90077) ( 108.773 -7.89256)
|
|---|
| 14 | continue ( 112.118 -7.88409) ( 115.674 -7.87538)
|
|---|
| 15 | continue ( 119.463 -7.86649) ( 123.510 -7.85740)
|
|---|
| 16 | continue ( 127.839 -7.84813) ( 132.484 -7.83871)
|
|---|
| 17 | continue ( 137.479 -7.82923) ( 142.865 -7.81972)
|
|---|
| 18 | continue ( 148.690 -7.81020) ( 155.011 -7.80078)
|
|---|
| 19 | continue ( 161.892 -7.79161) ( 169.413 -7.78274)
|
|---|
| 20 | continue ( 177.668 -7.77435) ( 186.767 -7.76669)
|
|---|
| 21 | continue ( 196.849 -7.75997) ( 208.081 -7.75454)
|
|---|
| 22 | continue ( 220.673 -7.75086) ( 234.887 -7.74951)
|
|---|
| 23 | continue ( 251.059 -7.75124) ( 269.621 -7.75718)
|
|---|
| 24 | continue ( 291.148 -7.76871) ( 316.410 -7.78799)
|
|---|
| 25 | continue ( 346.473 -7.81798) ( 382.848 -7.86306)
|
|---|
| 26 | continue ( 427.756 -7.93007) ( 484.601 -8.02977)
|
|---|
| 27 | continue ( 558.870 -8.18014) ( 660.023 -8.41298)
|
|---|
| 28 | continue ( 1276.47 -10.0456)
|
|---|
| 29 | continue (10,000,000 -30)
|
|---|
| 30 | luminosity = 37.53
|
|---|
| 31 | c
|
|---|
| 32 | c add background so that FIR/radio does not have zero flux
|
|---|
| 33 | background
|
|---|
| 34 | iterate
|
|---|
| 35 | c
|
|---|
| 36 | c density
|
|---|
| 37 | c hden 13.25 vary
|
|---|
| 38 | hden 13.25, power=-2.3
|
|---|
| 39 | c vary
|
|---|
| 40 | c grid range from 13.05 to 13.45 in 0.2 size steps
|
|---|
| 41 | c
|
|---|
| 42 | c set the radius range
|
|---|
| 43 | radius 10.5 vary
|
|---|
| 44 | grid range from 10.5 to 11.5 in 1.0 size steps
|
|---|
| 45 | stop thickness 12.7
|
|---|
| 46 | c
|
|---|
| 47 | c column density 23.0 vary
|
|---|
| 48 | c grid range 23.0 to 24.0 in 0.2 size steps
|
|---|
| 49 | c
|
|---|
| 50 | c set abundances ========
|
|---|
| 51 | c element scale factor oxygen 2
|
|---|
| 52 | c element scale factor neon 2
|
|---|
| 53 | element scale factor calcium 2
|
|---|
| 54 | element scale factor scandium 2
|
|---|
| 55 | element scale factor titanium 2
|
|---|
| 56 | element scale factor vanadium 2
|
|---|
| 57 | element scale factor chromium 2
|
|---|
| 58 | element scale factor manganese 2
|
|---|
| 59 | element scale factor iron 2.0
|
|---|
| 60 | c element scale factor cobalt 2.25
|
|---|
| 61 | c element scale factor nickel 2.25
|
|---|
| 62 | c
|
|---|
| 63 | c write stuff ===========
|
|---|
| 64 | set continuum resolution 0.2
|
|---|
| 65 | turbulence = 300
|
|---|
| 66 | punch overview last "ns94nb_full.ovr"
|
|---|
| 67 | punch column density last "ns94nb_full.col"
|
|---|
| 68 | punch continuum last "ns94nb_full.spc"
|
|---|
| 69 | punch heating last "ns94nb_full.het"
|
|---|
| 70 | c punch last xspec atable spectrum "trialspec.fit"
|
|---|
| 71 | c
|
|---|
| 72 | c punch out lots of useful info
|
|---|
| 73 | c
|
|---|
| 74 | punch line optical depths, limit=-2 "ns94nb_depths.dat"
|
|---|
| 75 | punch attenuated incident continuum "ns94nb_atten.dat"
|
|---|
| 76 | c punch last xspec mtable attenuated incident continuum "trialspec.fit"
|
|---|
