| 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 = 36.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 |
c hden 13.25, power=-2.3 |
|---|
| 39 |
hden 13.25 |
|---|
| 40 |
c vary |
|---|
| 41 |
c grid range from 13.05 to 13.45 in 0.2 size steps |
|---|
| 42 |
c |
|---|
| 43 |
c set the radius range |
|---|
| 44 |
radius 10.5 vary |
|---|
| 45 |
grid range from 10.5 to 11.5 in 1.0 size steps |
|---|
| 46 |
stop thickness 12.7 |
|---|
| 47 |
c |
|---|
| 48 |
c column density 23.0 vary |
|---|
| 49 |
c grid range 23.0 to 24.0 in 0.2 size steps |
|---|
| 50 |
c |
|---|
| 51 |
c set abundances ======== |
|---|
| 52 |
c element scale factor oxygen 2 |
|---|
| 53 |
c element scale factor neon 2 |
|---|
| 54 |
element scale factor calcium 2 |
|---|
| 55 |
element scale factor scandium 2 |
|---|
| 56 |
element scale factor titanium 2 |
|---|
| 57 |
element scale factor vanadium 2 |
|---|
| 58 |
element scale factor chromium 2 |
|---|
| 59 |
element scale factor manganese 2 |
|---|
| 60 |
element scale factor iron 2.0 |
|---|
| 61 |
c element scale factor cobalt 2.25 |
|---|
| 62 |
c element scale factor nickel 2.25 |
|---|
| 63 |
c |
|---|
| 64 |
c write stuff =========== |
|---|
| 65 |
set continuum resolution 0.2 |
|---|
| 66 |
turbulence = 300 |
|---|
| 67 |
punch overview last "miller.ovr" |
|---|
| 68 |
punch column density last "miller.col" |
|---|
| 69 |
punch continuum last "miller.spc" |
|---|
| 70 |
punch heating last "miller.het" |
|---|
| 71 |
c punch last xspec atable spectrum "miller.fit" |
|---|
| 72 |
c |
|---|
| 73 |
c punch out lots of useful info |
|---|
| 74 |
c |
|---|
| 75 |
c punch line optical depths, limit=-2 "ns94nb_depths.dat" |
|---|
| 76 |
punch attenuated incident continuum "miller.dat" |
|---|
| 77 |
punch last xspec mtable attenuated incident continuum "miller.fit" |
|---|
| 78 |
c |
|---|
| 79 |
c additional punches requested by Ferland |
|---|
| 80 |
c |
|---|
| 81 |
punch transmitted continuum last "miller.trn" |
|---|
| 82 |
punch continuum "miller.con" |
|---|
