diff --git a/reference/linelist.c17 b/reference/linelist.c17 index a7d7d2d..12b91bf 100644 --- a/reference/linelist.c17 +++ b/reference/linelist.c17 @@ -1,5 +1,4 @@ -#Proposed 'emission line list' for dumping fluxes to fort.xxx files. -# Lines appear in approximate order they appear in Cloudy standard output file. +# Proposed 'emission line list' for dumping fluxes to fort.xxx files. # Program should ignore the comments after the line label # Program knows to look for inward fractions as well # Program searches intrinsic emission line fluxes @@ -8,10 +7,9 @@ Inci 1215.00A incident continuum nu*f_nu near Ly-alpha, at illuminated face of cloud Inci 4860.00A incident continuum nu*f_nu near H-beta, at illuminated face of cloud -#Some informative quantities -# Bac 3646.00A residual flux at head of Balmer continuum, nuFnu +#Integrated H free-bound diffuse continua emitted. Ba C 0 integrated Balmer continuum emission -PA C 0 Paschen continuum emission +PA C 0 integrated Paschen continuum emission # Hydrogen transitions @@ -108,8 +106,8 @@ He 2 3923.37A H-like, 7 107, 4^2S - 15^2P #He I Lines: # First, the triplets -TOTL 1.08303m total emission in He-like lines, use average of three line wavelengths -He 1 1.08303m He-like, 2 4, 2^3S - 2^3P +#TOTL 1.08303m total emission in He-like lines, use average of three line wavelengths +He 1 1.08303m He-like, 2 4, 2^3S - 2^3P; same as above. # He 1 7065.22A He-like, 6 8, 2^3P_2 - 3^3S He 1 4713.03A He-like, 6 14, 2^3P_2 - 4^3S @@ -144,29 +142,35 @@ He 1 3964.73A He-like, 3 21, 2^1S - 4^1P_1 #Heavy Element Lines: Blnd 1657.00A Blend: "C 1 1657.01A"+"C 1R 1656.00A" +# no inwd #C 1 1657.01A Stout, 3 9 #C 1 1656.27A Stout, 2 9 #C 1 1656.93A Stout, 1 8 Blnd 2326.00A Blend: "C 2 2324.69A"+"C 2 2323.50A"+"C 2 2328.12A"+"C 2 2326.93A"+"C 2 2325.40A"+"C 2H 2326.00A" -#C 2 2323.50A Stout, 1 4 -#C 2 2324.69A Stout, 1 3 -#C 2 2325.40A Stout, 2 5 -#C 2 2328.12A Stout, 2 3 -#C 2 2326.93A Stout, 2 4 +# no inwd +C 2 2323.50A Stout, 1 4 +C 2 2324.69A Stout, 1 3 +C 2 2325.40A Stout, 2 5 +C 2 2328.12A Stout, 2 3 +C 2 2326.93A Stout, 2 4 Blnd 1335.00A Blend: "C 2 1334.53A"+"C 2 1335.66A"+"C 2 1335.71A"+"C 2R 1335.00A" +# no inwd #C 2 1334.53A Stout, 1 7 #C 2 1335.66A Stout, 2 7 #C 2 1335.71A Stout, 2 6 C 3 977.000A recombination line +# no inwd Blnd 1909.00A Blend: "C 3 1908.73A"+"C 3 1906.68A"+"C 3R 1909.00A"+"C 3H 1909.00A" +# no inwd C 3 1908.73A Stout, 1 3 -C 3 1906.68A Stout, 1 4 +#C 3 1906.68A Stout, 1 4 Blnd 1176.00A Blend: "C 3 1174.61A"+"C 3 1174.93A"+"C 3 1175.26A"+"C 3 1175.59A"+"C 3 1175.71A"+"C 3 1175.99A"+"C 3 1176.37A"+"C 3 1176.77A"+"C 3R 1175.00A" +# no inwd #C 3 1175.26A Stout, 2 7 #C 3 1174.61A Stout, 2 8 #C 3 1175.99A Stout, 3 6 @@ -181,10 +185,12 @@ C 4 1550.78A C 4 1548.19A Blnd 2141.00A Blend: "N 2 2139.01A"+"N 2 2142.78A" +# no inwd #N 2 2139.01A #N 2 2142.78A Blnd 1085.00A Blend: "N 2 1083.99A"+"N 2 1084.56A"+"N 2 1084.58A"+"N 2 1085.53A"+"N 2 1085.55A"+"N 2 1085.70A"+"N 2R 1085.00A" +# no inwd #N 2 1083.99A #N 2 1084.56A #N 2 1084.58A @@ -194,6 +200,7 @@ Blnd 1085.00A Blend: "N 2 1083.99A"+"N 2 1084.56A"+"N 2 #N 2R 1085.00A Blnd 1750.00A Blend: "N 3 1746.82A"+"N 3 1748.65A"+"N 3 1749.67A"+"N 3 1752.16A"+"N 3 1753.99A" +# no inwd N 3 1746.82A N 3 1748.65A N 3 1749.67A @@ -201,17 +208,20 @@ N 3 1752.16A N 3 1753.99A Blnd 990.000A Blend: "N 3 989.799A"+"N 3 991.511A"+"N 3 991.577A"+"N 3R 990.000A" +# no inwd #N 3 989.799A #N 3 991.511A #N 3 991.577A #N 3R 990.000A Blnd 1486.00A Blend: "N 4 1483.32A"+"N 4 1486.50A" -N 4 1483.32A +# no inwd +#N 4 1483.32A N 4 1486.50A Blnd 765.000A Blend: "N 4 765.147A"+"N 4R 765.000A" -#N 4 765.147A +# no inwd +N 4 765.147A #N 4R 765.000A Blnd 1240.00A Blend: "N 5 1238.82A"+"N 5 1242.80A" @@ -219,22 +229,26 @@ N 5 1238.82A N 5 1242.80A Blnd 1304.00A Blend: "O 1 1302.17A"+"O 1 1304.86A"+"O 1 1306.03A" -#O 1 1302.17A -#O 1 1304.86A -#O 1 1306.03A +# no inwd +O 1 1302.17A +O 1 1304.86A +O 1 1306.03A Blnd 8446.00A Blend: "O 1 8446.25A"+"O 1 8446.36A"+"O 1 8446.76A" +# no inwd #O 1 8446.25A #O 1 8446.36A #O 1 8446.76A Blnd 1666.00A Blend: "O 3 1666.15A"+"O 3 1660.81A"+"O 3H 1665.00A"+"O 3A 1665.00A" +# no inwd O 3 1666.15A O 3 1660.81A #O 3H 1665.00A #O 3A 1665.00A Blnd 835.000A Blend: "O 3 832.929A"+"O 3 833.715A"+"O 3 833.749A"+"O 3 835.059A"+"O 3 835.092A"+"O 3 835.289A"+"O 3R 835.000A" +# no inwd #O 3 832.929A #O 3 833.715A #O 3 833.749A @@ -244,6 +258,7 @@ Blnd 835.000A Blend: "O 3 832.929A"+"O 3 833.715A" #O 3R 835.000A Blnd 1402.00A Blend: "O 4 1397.20A"+"O 4 1399.77A"+"O 4 1401.16A"+"O 4 1404.78A"+"O 4 1407.38A"+"O 4H 1401.00A" +# no inwd O 4 1397.20A O 4 1399.77A O 4 1401.16A @@ -252,14 +267,17 @@ O 4 1407.38A #O 4H 1401.00A Blnd 630.000A Blend: "O 5 629.732A"+"O 5R 630.000A" -#O 5 629.732A +# no inwd +O 5 629.732A #O 5R 630.000A Blnd 1218.00A Blend: "O 5 1213.81A"+"O 5 1218.34A" -O 5 1213.81A +# no inwd +#O 5 1213.81A O 5 1218.34A Blnd 1035.00A Blend: "O 6 1031.91A"+"O 6 1037.62A" +# ignore inwd -- same fluxes as Blnd O 6 1031.91A O 6 1037.62A @@ -267,32 +285,39 @@ Ne 5 1136.51A Chianti, 2 6 Ne 5 1145.60A Chianti, 3 6 Blnd 774.000A Blend: "Ne 8 770.410A"+"Ne 8 780.325A" +# no inwd Ne 8 770.410A Ne 8 780.325A Blnd 5892.00A Blend: "Na 1 5889.95A"+"Na 1 5895.92A" +# no inwd Na 1 5889.95A Na 1 5895.92A Mg 1 2852.13A Stout, 1 5 Blnd 2798.00A Blend: "Mg 2 2795.53A"+"Mg 2 2802.71A" +# ignore inwd -- spurious values Mg 2 2795.53A Mg 2 2802.71A Blnd 615.000A Blend: "Mg10 609.794A"+"Mg10 624.943A" +# no inwd Mg10 609.794A Mg10 624.943A Blnd 2665.00A Blend: "Al 2 2669.15A"+"Al 2 2660.35A" +# no inwd Al 2 2669.15A -Al 2 2660.35A +#Al 2 2660.35A Blnd 1860.00A Blend: "Al 3 1854.72A"+"Al 3 1862.79A" +# no inwd Al 3 1854.72A Al 3 1862.79A Blnd 2335.00A Blend: "Si 2 2334.41A"+"Si 2 2328.52A"+"Si 2 2350.17A"+"Si 2 2344.20A"+"Si 2 2334.60A" +# no inwd Si 2 2334.41A Si 2 2328.52A Si 2 2350.17A @@ -305,43 +330,48 @@ Si 2 1533.43A Stout, 2 8 Si 2 1309.28A Stout, 2 9 Si 2 1265.00A Stout, 2 10 Si 2 1264.74A Stout, 2 11 -Si 2 1236.03A Stout, 2 12 -Si 2 1235.12A Stout, 2 13 Blnd 1888.00A Blend: "Si 3 1882.71A"+"Si 3 1892.03A" -Si 3 1882.71A +# no inwd +#Si 3 1882.71A Si 3 1892.03A Blnd 1397.00A Blend: "Si 4 1393.75A"+"Si 4 1402.77A" +# no inwd Si 4 1393.75A Si 4 1402.77A Blnd 1256.00A Blend: "S 2 1250.58A"+"S 2 1253.81A"+"S 2 1259.52A" +# no inwd S 2 1250.58A S 2 1253.81A S 2 1259.52A Blnd 1720.00A Blend: "S 3 1704.39A"+"S 3 1713.11A"+"S 3 1728.94A" -S 3 1704.39A +# no inwd +#S 3 1704.39A S 3 1713.11A S 3 1728.94A Blnd 1198.00A Blend: "S 3 1190.20A"+"S 3 1194.06A"+"S 3 1194.45A"+"S 3 1200.97A"+"S 3 1201.73A"+"S 3 1202.12A" +# no inwd S 3 1190.20A S 3 1194.06A S 3 1194.45A -S 3 1200.97 +#S 3 1200.97 S 3 1201.73A S 3 1202.12A Blnd 1406.00A Blend: "S 4 1404.81A"+"S 4 1398.04A"+"S 4 1423.84A"+"S 4 1416.89A"+"S 4 1406.02A" +# no inwd S 4 1404.81A S 4 1398.04A S 4 1423.84A -S 4 1416.89 +#S 4 1416.89 S 4 1406.02A Blnd 1199.00A Blend: "S 5 1199.14A"+"S 5 1188.28A" +# no inwd S 5 1199.14A S 5 1188.28A @@ -349,11 +379,13 @@ S 6 944.524A Chianti, 1 2 S 6 933.380A Chianti, 1 3 Blnd 8579.00A Blend: "Ca 2 8662.14A"+"Ca 2 8542.09A"+"Ca 2 8498.02A" +# no inwd Ca 2 8662.14A Ca 2 8542.09A Ca 2 8498.02A Blnd 3933.00A Blend: "Ca 2 3933.66A"+"Ca 2 3968.47A" +# no inwd Ca 2 3933.66A Ca 2 3968.47A #Ca2R 3933.00A recombination contribution to CaII emission @@ -378,8 +410,8 @@ FeKa 1.78000A total intensity of K-alpha line #nFnu 6917.22m total continuum at selected energy points # nFnu 122.400m total continuum at selected energy points -nFnu 100.000m supplemental -nFnu 80.0000m supplemental +#nFnu 100.000m supplemental +#nFnu 80.0000m supplemental nFnu 60.8322m total continuum at selected energy points nFnu 41.2152m total continuum at selected energy points nFnu 30.8695m total continuum at selected energy points @@ -422,6 +454,7 @@ nFnu 8250.00A supplemental nFnu 8210.00A supplemental nFnu 8200.00A supplemental nFnu 8150.00A supplemental +nFnu 8100.00A supplemental nFnu 8078.61A total continuum at selected energy points nFnu 7800.00A supplemental nFnu 7600.00A supplemental @@ -493,17 +526,17 @@ nFnu 950.000A supplemental nFnu 925.144A total continuum at selected energy points nFnu 913.000A supplemental nFnu 897.800A total continuum at selected energy points -nFnu 875.000A supplemental +#nFnu 875.000A supplemental nFnu 850.000A supplemental nFnu 800.000A supplemental nFnu 760.023A total continuum at selected energy points nFnu 701.514A total continuum at selected energy points -nFnu 650.000A supplemental +#nFnu 650.000A supplemental nFnu 608.160A total continuum at selected energy points nFnu 575.146A total continuum at selected energy points nFnu 511.948A total continuum at selected energy points nFnu 496.874A total continuum at selected energy points -nFnu 450.000A supplemental +#nFnu 450.000A supplemental nFnu 399.153A total continuum at selected energy points -# 125 nFnu continuum wavelength bins +# approximately 125 nFnu continuum wavelength bins # end of line list \ No newline at end of file diff --git a/scripts/cloudy/magdziarz.mpi.in b/scripts/cloudy/magdziarz.mpi.in deleted file mode 100644 index ac6e5dd..0000000 --- a/scripts/cloudy/magdziarz.mpi.in +++ /dev/null @@ -1,348 +0,0 @@ -////// CLOUDY 17.00 one-shot using magdziarz 1997 ngc5548 SED -// -title magzdiarz 1997 ngc5548 one-shot -// -// single run copied from the standard magdziarz grid -// -////// Output continua -// -save last incident continuum ".incident_continuum" -// -////// Formatting -// -//print line column -// print line column commented for testing, -// will print lines in a single column if uncommented -print line inward -print line iso collapsed -print line faint -3 -print last -print line precision 6 -print line iso collapsed on -Database print -save line labels ".line_labels" -// -////// Resolved energy states for atomic models -// -// Defaults: 10, 10, 6 -// -database H-like element hydrogen resolved levels 18 -database H-like element helium resolved levels 15 -database He-like element helium resolved levels 15 -//database H-like element hydrogen collapsed levels 30 -//database H-like element helium collapsed levels 30 -//database He-like element helium collapsed levels 40 -// -////// low-metal abundance for AGN testing ~2017 -// metals 0.5 -// element scale factor carbon 2 -// -////// Prevent cloudy from stopping prematurely -// -set nend 5000 -// -////// Print diffuse continuum information in cloudy output -// and include additional spectral resolution around Lyman, Briggs, Etc. -// -// Continuum labels by default refer to the following fields: -// nFnu = Total Diffuse (in & out) + Reflected Inward -// nInu = not of interest, described in hazy 1 -// InwT = Inward Diffuse + Reflected Inward -// InwC = Reflected Inward -// -// derived quantities added in Otho's program ~2017 -// InwD = InwT - InwC = Inward Diffuse -// nFnD = nFnu - InwC = Total Diffuse -// -// "print continuum" redundant with set nFnu in c17 -// -print continuum -set nFnu add 800 angstrom -set nFnu add 840 angstrom -set nFnu add 850 angstrom -set nFnu add 880 angstrom -set nFnu add 900 angstrom -set nFnu add 910 angstrom -set nFnu add 913 angstrom -set nFnu add 950 angstrom -set nFnu add 960 angstrom -set nFnu add 1000 angstrom -set nFnu add 1050 angstrom -set nFnu add 1100 angstrom -set nFnu add 1150 angstrom -set nFnu add 1200 angstrom -set nFnu add 1240 angstrom -set nFnu add 1250 angstrom -set nFnu add 1300 angstrom -set nFnu add 1324 angstrom -set nFnu add 1350 angstrom -set nFnu add 1400 angstrom -set nFnu add 1410 angstrom -set nFnu add 1450 angstrom -set nFnu add 1500 angstrom -set nFnu add 1550 angstrom -set nFnu add 1600 angstrom -set nFnu add 1650 angstrom -set nFnu add 1700 angstrom -set nFnu add 1750 angstrom -set nFnu add 1800 angstrom -set nFnu add 1850 angstrom -set nFnu add 1900 angstrom -set nFnu add 1950 angstrom -set nFnu add 2000 angstrom -set nFnu add 2050 angstrom -set nFnu add 2100 angstrom -set nFnu add 2150 angstrom -set nFnu add 2200 angstrom -set nFnu add 2250 angstrom -set nFnu add 2300 angstrom -set nFnu add 2350 angstrom -set nFnu add 2400 angstrom -set nFnu add 2450 angstrom -set nFnu add 2500 angstrom -set nFnu add 2550 angstrom -set nFnu add 2600 angstrom -set nFnu add 2650 angstrom -set nFnu add 2700 angstrom -set nFnu add 2750 angstrom -set nFnu add 2800 angstrom -set nFnu add 2850 angstrom -set nFnu add 2900 angstrom -set nFnu add 2950 angstrom -set nFnu add 3000 angstrom -set nFnu add 3050 angstrom -set nFnu add 3100 angstrom -set nFnu add 3150 angstrom -set nFnu add 3200 angstrom -set nFnu add 3250 angstrom -set nFnu add 3300 angstrom -set nFnu add 3350 angstrom -set nFnu add 3400 angstrom -set nFnu add 3450 angstrom -set nFnu add 3500 angstrom -set nFnu add 3550 angstrom -set nFnu add 3600 angstrom -set nFnu add 3644 angstrom -set nFnu add 3645 angstrom -set nFnu add 3648 angstrom -set nFnu add 3650 angstrom -set nFnu add 3700 angstrom -set nFnu add 3750 angstrom -set nFnu add 3800 angstrom -set nFnu add 3850 angstrom -set nFnu add 3900 angstrom -set nFnu add 3950 angstrom -set nFnu add 4000 angstrom -set nFnu add 4050 angstrom -set nFnu add 4100 angstrom -set nFnu add 4150 angstrom -set nFnu add 4200 angstrom -set nFnu add 4250 angstrom -set nFnu add 4300 angstrom -set nFnu add 4350 angstrom -set nFnu add 4400 angstrom -set nFnu add 4450 angstrom -set nFnu add 4500 angstrom -set nFnu add 4550 angstrom -set nFnu add 4600 angstrom -set nFnu add 4650 angstrom -set nFnu add 4700 angstrom -set nFnu add 4750 angstrom -set nFnu add 4800 angstrom -set nFnu add 4850 angstrom -set nFnu add 4900 angstrom -set nFnu add 4950 angstrom -set nFnu add 5000 angstrom -set nFnu add 5050 angstrom -set nFnu add 5100 angstrom -set nFnu add 5150 angstrom -set nFnu add 5200 angstrom -set nFnu add 5250 angstrom -set nFnu add 5300 angstrom -set nFnu add 5350 angstrom -set nFnu add 5400 angstrom -set nFnu add 5450 angstrom -set nFnu add 5500 angstrom -set nFnu add 5550 angstrom -set nFnu add 5600 angstrom -set nFnu add 5650 angstrom -set nFnu add 5700 angstrom -set nFnu add 5750 angstrom -set nFnu add 5800 angstrom -set nFnu add 5850 angstrom -set nFnu add 5900 angstrom -set nFnu add 5950 angstrom -set nFnu add 6000 angstrom -set nFnu add 6050 angstrom -set nFnu add 6100 angstrom -set nFnu add 6150 angstrom -set nFnu add 6200 angstrom -set nFnu add 6250 angstrom -set nFnu add 6300 angstrom -set nFnu add 6350 angstrom -set nFnu add 6400 angstrom -set nFnu add 6450 angstrom -set nFnu add 6500 angstrom -set nFnu add 6550 angstrom -set nFnu add 6600 angstrom -set nFnu add 6650 angstrom -set nFnu add 6700 angstrom -set nFnu add 6750 angstrom -set nFnu add 6800 angstrom -set nFnu add 6850 angstrom -set nFnu add 6900 angstrom -set nFnu add 6950 angstrom -set nFnu add 7000 angstrom -set nFnu add 7050 angstrom -set nFnu add 7100 angstrom -set nFnu add 7150 angstrom -set nFnu add 7200 angstrom -set nFnu add 7250 angstrom -set nFnu add 7300 angstrom -set nFnu add 7350 angstrom -set nFnu add 7400 angstrom -set nFnu add 7450 angstrom -set nFnu add 7500 angstrom -set nFnu add 7550 angstrom -set nFnu add 7600 angstrom -set nFnu add 7650 angstrom -set nFnu add 7700 angstrom -set nFnu add 7750 angstrom -set nFnu add 7800 angstrom -set nFnu add 7850 angstrom -set nFnu add 7900 angstrom -set nFnu add 7950 angstrom -set nFnu add 8000 angstrom -set nFnu add 8050 angstrom -set nFnu add 8100 angstrom -set nFnu add 8150 angstrom -set nFnu add 8200 angstrom -set nFnu add 8204 angstrom -set nFnu add 8205 angstrom -set nFnu add 8207 angstrom -set nFnu add 8210 angstrom -set nFnu add 8250 angstrom -set nFnu add 8300 angstrom -set nFnu add 8350 angstrom -set nFnu add 8400 angstrom -set nFnu add 8450 angstrom -set nFnu add 8500 angstrom -set nFnu add 8550 angstrom -set nFnu add 8600 angstrom -set nFnu add 8650 angstrom -set nFnu add 8700 angstrom -set nFnu add 8750 angstrom -set nFnu add 8800 angstrom -set nFnu add 8850 angstrom -set nFnu add 8900 angstrom -set nFnu add 8950 angstrom -set nFnu add 9000 angstrom -set nFnu add 9050 angstrom -set nFnu add 9100 angstrom -set nFnu add 9150 angstrom -set nFnu add 9200 angstrom -set nFnu add 9250 angstrom -set nFnu add 9300 angstrom -set nFnu add 9350 angstrom -set nFnu add 9400 angstrom -set nFnu add 9450 angstrom -set nFnu add 9500 angstrom -set nFnu add 9550 angstrom -set nFnu add 9600 angstrom -set nFnu add 9650 angstrom -set nFnu add 9700 angstrom -set nFnu add 9750 angstrom -set nFnu add 9800 angstrom -set nFnu add 9850 angstrom -set nFnu add 9900 angstrom -set nFnu add 9950 angstrom -set nFnu add 10000 angstrom -set nFnu add 10050 angstrom -set nFnu add 10100 angstrom -set nFnu add 10150 angstrom -set nFnu add 10200 angstrom -set nFnu add 10250 angstrom -set nFnu add 10300 angstrom -set nFnu add 10350 angstrom -set nFnu add 10400 angstrom -set nFnu add 10450 angstrom -set nFnu add 10500 angstrom -set nFnu add 10550 angstrom -set nFnu add 10600 angstrom -set nFnu add 10650 angstrom -set nFnu add 10700 angstrom -set nFnu add 10750 angstrom -set nFnu add 10800 angstrom -set nFnu add 10850 angstrom -set nFnu add 10900 angstrom -set nFnu add 10950 angstrom -set nFnu add 11000 angstrom -set nFnu add 11050 angstrom -set nFnu add 11100 angstrom -set nFnu add 11150 angstrom -set nFnu add 11200 angstrom -set nFnu add 11250 angstrom -set nFnu add 11300 angstrom -set nFnu add 11350 angstrom -set nFnu add 11400 angstrom -set nFnu add 11450 angstrom -set nFnu add 11500 angstrom -set nFnu add 11550 angstrom -set nFnu add 11600 angstrom -set nFnu add 11650 angstrom -set nFnu add 11700 angstrom -set nFnu add 11750 angstrom -set nFnu add 11800 angstrom -set nFnu add 11850 angstrom -set nFnu add 11900 angstrom -set nFnu add 11950 angstrom -set nFnu add 12000 angstrom -set nFnu add 12050 angstrom -set nFnu add 12100 angstrom -set nFnu add 12150 angstrom -set nFnu add 12200 angstrom -set nFnu add 12250 angstrom -set nFnu add 12300 angstrom -set nFnu add 12350 angstrom -set nFnu add 12400 angstrom -set nFnu add 12450 angstrom -set nFnu add 12500 angstrom -set nFnu add 12550 angstrom -set nFnu add 12600 angstrom -set nFnu add 12650 angstrom -set nFnu add 12700 angstrom -set nFnu add 12750 angstrom -set nFnu add 12800 angstrom -set nFnu add 12850 angstrom -set nFnu add 12900 angstrom -set nFnu add 12950 angstrom -set nFnu add 13000 angstrom -// -////// Print equivalent widths -// -normalize to "Inci" 1215 scale to 1215 -// -////// Background radiation -// -background z = 0.0172 -// -////// Hydrogen density of cloud is varied -// -hden 11.5 -//hden 11 vary -//grid range from 7 to 14 with .25 dex -// -////// Central AGN source shape generated by external program -// -table SED "magdziarz1997.tab" -// -////// Flux is varied -// -phi(h) 21.5 -// phi(h) 23 vary -//grid range from 17 to 24 with .25 dex -// -////// Stopping condition and run -// -stop total column density 23 -iterate to convergence max=40 diff --git a/scripts/cloudy/magdziarz.mpi_grid.in b/scripts/cloudy/magdziarz.mpi_grid.in index 40d9a03..437fc05 100644 --- a/scripts/cloudy/magdziarz.mpi_grid.in +++ b/scripts/cloudy/magdziarz.mpi_grid.in @@ -2,25 +2,22 @@ // title magdziarz 1997 hden_vs_phi .25 dex // -// single run copied from the standard magdziarz grid -// ////// Output continua and included lines // save last incident continuum "incident_continuum" save line labels "line_labels" // -////// Formatting +////// Report almost everything // print line column print line inward -print line iso collapsed +print line iso collapsed on print line faint -3 print last print line precision 6 -print line iso collapsed on Database print // -////// Resolved energy states for atomic models +////// Resolved energy states for atomic models // // Defaults: 10, 10, 6 // @@ -34,8 +31,9 @@ database He-like element helium resolved levels 15 ////// low-metal abundance for AGN testing ~2017 // metals 0.5 // element scale factor carbon 2 +// element scale factor nitrogen 2 // -////// Prevent cloudy from stopping prematurely +////// Prevent cloudy from stopping prematurely // set nend 5000 // @@ -55,6 +53,9 @@ set nend 5000 // "print continuum" redundant with set nFnu in c17 // print continuum +set nFnu add 450 angstrom +set nFnu add 650 angstrom +set nFnu add 875 angstrom set nFnu add 800 angstrom set nFnu add 840 angstrom set nFnu add 850 angstrom @@ -315,21 +316,23 @@ set nFnu add 12850 angstrom set nFnu add 12900 angstrom set nFnu add 12950 angstrom set nFnu add 13000 angstrom +set nFnu add 800000 angstrom +set nFnu add 1000000 angstrom // -////// Print equivalent widths +////// Print equivalent widths // normalize to "Inci" 1215 scale to 1215 // -////// Background radiation +////// Background radiation // background z = 0.0172 // -////// Hydrogen density of cloud is varied +////// Hydrogen density of cloud is varied // hden 11 vary grid range from 7 to 14 with .25 dex // -////// magdizard et al. 1998 NGC 5548 SED +////// magdizard et al. 1998 NGC 5548 SED // table SED "magdziarz1997.tab" // @@ -339,7 +342,7 @@ table SED "magdziarz1997.tab" phi(h) 23 vary grid range from 17 to 24 with .25 dex // -////// Stopping condition and run +////// Stopping condition and run // stop total column density 23 -iterate to convergence max=40 +iterate to convergence max=40 \ No newline at end of file diff --git a/scripts/cloudy/mehdipour.mpi.in b/scripts/cloudy/mehdipour.mpi.in deleted file mode 100644 index 2bff60b..0000000 --- a/scripts/cloudy/mehdipour.mpi.in +++ /dev/null @@ -1,549 +0,0 @@ -////// CLOUDY 17.00 one-shot using mehdipour ngc5548 SED -// -title mehdipour ngc5548 one-shot -// -// single run copied from the standard mehdipour grid -// -////// Output continua -// -save last incident continuum "incident_continuum" -// -////// Formatting -// -//print line column -// print line column commented for testing, -// will print lines in a single column if uncommented -print line inward -print line iso collapsed -print line faint -3 -print last -print line precision 6 -print line iso collapsed on -Database print -save line labels ".line_labels" -// -////// Resolved energy states for atomic models -// -// Defaults: 10, 10, 6 -// -database H-like element hydrogen resolved levels 18 -database H-like element helium resolved levels 15 -database He-like element helium resolved levels 15 -//database H-like element hydrogen collapsed levels 30 -//database H-like element helium collapsed levels 30 -//database He-like element helium collapsed levels 40 -// -////// low-metal abundance for AGN testing ~2017 -// metals 0.5 -// element scale factor carbon 2 -// -////// Prevent cloudy from stopping prematurely -// -set nend 5000 -// -////// Print diffuse continuum information in cloudy output -// and include additional spectral resolution around Lyman, Briggs, Etc. -// -// Continuum labels by default refer to the following fields: -// nFnu = Total Diffuse (in & out) + Reflected Inward -// nInu = not of interest, described in hazy 1 -// InwT = Inward Diffuse + Reflected Inward -// InwC = Reflected Inward -// -// derived quantities added in Otho's program ~2017 -// InwD = InwT - InwC = Inward Diffuse -// nFnD = nFnu - InwC = Total Diffuse -// -// "print continuum" redundant with set nFnu in c17 -// -print continuum -set nFnu add 800 angstrom -set nFnu add 840 angstrom -set nFnu add 850 angstrom -set nFnu add 880 angstrom -set nFnu add 900 angstrom -set nFnu add 910 angstrom -set nFnu add 913 angstrom -set nFnu add 950 angstrom -set nFnu add 960 angstrom -set nFnu add 1000 angstrom -set nFnu add 1050 angstrom -set nFnu add 1100 angstrom -set nFnu add 1150 angstrom -set nFnu add 1200 angstrom -set nFnu add 1240 angstrom -set nFnu add 1250 angstrom -set nFnu add 1300 angstrom -set nFnu add 1324 angstrom -set nFnu add 1350 angstrom -set nFnu add 1400 angstrom -set nFnu add 1410 angstrom -set nFnu add 1450 angstrom -set nFnu add 1500 angstrom -set nFnu add 1550 angstrom -set nFnu add 1600 angstrom -set nFnu add 1650 angstrom -set nFnu add 1700 angstrom -set nFnu add 1750 angstrom -set nFnu add 1800 angstrom -set nFnu add 1850 angstrom -set nFnu add 1900 angstrom -set nFnu add 1950 angstrom -set nFnu add 2000 angstrom -set nFnu add 2050 angstrom -set nFnu add 2100 angstrom -set nFnu add 2150 angstrom -set nFnu add 2200 angstrom -set nFnu add 2250 angstrom -set nFnu add 2300 angstrom -set nFnu add 2350 angstrom -set nFnu add 2400 angstrom -set nFnu add 2450 angstrom -set nFnu add 2500 angstrom -set nFnu add 2550 angstrom -set nFnu add 2600 angstrom -set nFnu add 2650 angstrom -set nFnu add 2700 angstrom -set nFnu add 2750 angstrom -set nFnu add 2800 angstrom -set nFnu add 2850 angstrom -set nFnu add 2900 angstrom -set nFnu add 2950 angstrom -set nFnu add 3000 angstrom -set nFnu add 3050 angstrom -set nFnu add 3100 angstrom -set nFnu add 3150 angstrom -set nFnu add 3200 angstrom -set nFnu add 3250 angstrom -set nFnu add 3300 angstrom -set nFnu add 3350 angstrom -set nFnu add 3400 angstrom -set nFnu add 3450 angstrom -set nFnu add 3500 angstrom -set nFnu add 3550 angstrom -set nFnu add 3600 angstrom -set nFnu add 3644 angstrom -set nFnu add 3645 angstrom -set nFnu add 3648 angstrom -set nFnu add 3650 angstrom -set nFnu add 3700 angstrom -set nFnu add 3750 angstrom -set nFnu add 3800 angstrom -set nFnu add 3850 angstrom -set nFnu add 3900 angstrom -set nFnu add 3950 angstrom -set nFnu add 4000 angstrom -set nFnu add 4050 angstrom -set nFnu add 4100 angstrom -set nFnu add 4150 angstrom -set nFnu add 4200 angstrom -set nFnu add 4250 angstrom -set nFnu add 4300 angstrom -set nFnu add 4350 angstrom -set nFnu add 4400 angstrom -set nFnu add 4450 angstrom -set nFnu add 4500 angstrom -set nFnu add 4550 angstrom -set nFnu add 4600 angstrom -set nFnu add 4650 angstrom -set nFnu add 4700 angstrom -set nFnu add 4750 angstrom -set nFnu add 4800 angstrom -set nFnu add 4850 angstrom -set nFnu add 4900 angstrom -set nFnu add 4950 angstrom -set nFnu add 5000 angstrom -set nFnu add 5050 angstrom -set nFnu add 5100 angstrom -set nFnu add 5150 angstrom -set nFnu add 5200 angstrom -set nFnu add 5250 angstrom -set nFnu add 5300 angstrom -set nFnu add 5350 angstrom -set nFnu add 5400 angstrom -set nFnu add 5450 angstrom -set nFnu add 5500 angstrom -set nFnu add 5550 angstrom -set nFnu add 5600 angstrom -set nFnu add 5650 angstrom -set nFnu add 5700 angstrom -set nFnu add 5750 angstrom -set nFnu add 5800 angstrom -set nFnu add 5850 angstrom -set nFnu add 5900 angstrom -set nFnu add 5950 angstrom -set nFnu add 6000 angstrom -set nFnu add 6050 angstrom -set nFnu add 6100 angstrom -set nFnu add 6150 angstrom -set nFnu add 6200 angstrom -set nFnu add 6250 angstrom -set nFnu add 6300 angstrom -set nFnu add 6350 angstrom -set nFnu add 6400 angstrom -set nFnu add 6450 angstrom -set nFnu add 6500 angstrom -set nFnu add 6550 angstrom -set nFnu add 6600 angstrom -set nFnu add 6650 angstrom -set nFnu add 6700 angstrom -set nFnu add 6750 angstrom -set nFnu add 6800 angstrom -set nFnu add 6850 angstrom -set nFnu add 6900 angstrom -set nFnu add 6950 angstrom -set nFnu add 7000 angstrom -set nFnu add 7050 angstrom -set nFnu add 7100 angstrom -set nFnu add 7150 angstrom -set nFnu add 7200 angstrom -set nFnu add 7250 angstrom -set nFnu add 7300 angstrom -set nFnu add 7350 angstrom -set nFnu add 7400 angstrom -set nFnu add 7450 angstrom -set nFnu add 7500 angstrom -set nFnu add 7550 angstrom -set nFnu add 7600 angstrom -set nFnu add 7650 angstrom -set nFnu add 7700 angstrom -set nFnu add 7750 angstrom -set nFnu add 7800 angstrom -set nFnu add 7850 angstrom -set nFnu add 7900 angstrom -set nFnu add 7950 angstrom -set nFnu add 8000 angstrom -set nFnu add 8050 angstrom -set nFnu add 8100 angstrom -set nFnu add 8150 angstrom -set nFnu add 8200 angstrom -set nFnu add 8204 angstrom -set nFnu add 8205 angstrom -set nFnu add 8207 angstrom -set nFnu add 8210 angstrom -set nFnu add 8250 angstrom -set nFnu add 8300 angstrom -set nFnu add 8350 angstrom -set nFnu add 8400 angstrom -set nFnu add 8450 angstrom -set nFnu add 8500 angstrom -set nFnu add 8550 angstrom -set nFnu add 8600 angstrom -set nFnu add 8650 angstrom -set nFnu add 8700 angstrom -set nFnu add 8750 angstrom -set nFnu add 8800 angstrom -set nFnu add 8850 angstrom -set nFnu add 8900 angstrom -set nFnu add 8950 angstrom -set nFnu add 9000 angstrom -set nFnu add 9050 angstrom -set nFnu add 9100 angstrom -set nFnu add 9150 angstrom -set nFnu add 9200 angstrom -set nFnu add 9250 angstrom -set nFnu add 9300 angstrom -set nFnu add 9350 angstrom -set nFnu add 9400 angstrom -set nFnu add 9450 angstrom -set nFnu add 9500 angstrom -set nFnu add 9550 angstrom -set nFnu add 9600 angstrom -set nFnu add 9650 angstrom -set nFnu add 9700 angstrom -set nFnu add 9750 angstrom -set nFnu add 9800 angstrom -set nFnu add 9850 angstrom -set nFnu add 9900 angstrom -set nFnu add 9950 angstrom -set nFnu add 10000 angstrom -set nFnu add 10050 angstrom -set nFnu add 10100 angstrom -set nFnu add 10150 angstrom -set nFnu add 10200 angstrom -set nFnu add 10250 angstrom -set nFnu add 10300 angstrom -set nFnu add 10350 angstrom -set nFnu add 10400 angstrom -set nFnu add 10450 angstrom -set nFnu add 10500 angstrom -set nFnu add 10550 angstrom -set nFnu add 10600 angstrom -set nFnu add 10650 angstrom -set nFnu add 10700 angstrom -set nFnu add 10750 angstrom -set nFnu add 10800 angstrom -set nFnu add 10850 angstrom -set nFnu add 10900 angstrom -set nFnu add 10950 angstrom -set nFnu add 11000 angstrom -set nFnu add 11050 angstrom -set nFnu add 11100 angstrom -set nFnu add 11150 angstrom -set nFnu add 11200 angstrom -set nFnu add 11250 angstrom -set nFnu add 11300 angstrom -set nFnu add 11350 angstrom -set nFnu add 11400 angstrom -set nFnu add 11450 angstrom -set nFnu add 11500 angstrom -set nFnu add 11550 angstrom -set nFnu add 11600 angstrom -set nFnu add 11650 angstrom -set nFnu add 11700 angstrom -set nFnu add 11750 angstrom -set nFnu add 11800 angstrom -set nFnu add 11850 angstrom -set nFnu add 11900 angstrom -set nFnu add 11950 angstrom -set nFnu add 12000 angstrom -set nFnu add 12050 angstrom -set nFnu add 12100 angstrom -set nFnu add 12150 angstrom -set nFnu add 12200 angstrom -set nFnu add 12250 angstrom -set nFnu add 12300 angstrom -set nFnu add 12350 angstrom -set nFnu add 12400 angstrom -set nFnu add 12450 angstrom -set nFnu add 12500 angstrom -set nFnu add 12550 angstrom -set nFnu add 12600 angstrom -set nFnu add 12650 angstrom -set nFnu add 12700 angstrom -set nFnu add 12750 angstrom -set nFnu add 12800 angstrom -set nFnu add 12850 angstrom -set nFnu add 12900 angstrom -set nFnu add 12950 angstrom -set nFnu add 13000 angstrom -// -////// Print equivalent widths -// -normalize to "Inci" 1215 scale to 1215 -// -////// Background radiation -// -background z = 0.0172 -// -////// Hydrogen density of cloud is varied -// -hden 10.5 -//hden 11 vary -//grid range from 7 to 14 with .25 dex -// -////// Central AGN source shape generated by external program -// -interpolate (1.001e-09 1e-36) -continue (0.000734986 -33) (0.000746929 -33.007) (0.000759066 -33.014) (0.000771399 -33.021) (0.000783934 -33.028) -continue (0.000796672 -33.035) (0.000809616 -33.042) (0.000822772 -33.049) (0.000836141 -33.056) (0.000849727 -33.063) -continue (0.000863534 -33.07) (0.000877565 -33.077) (0.000891824 -33.084) (0.000906315 -33.091) (0.000921042 -33.098) -continue (0.000936008 -33.105) (0.000951216 -33.112) (0.000966672 -33.119) (0.00098238 -33.126) (0.000998342 -33.133) -continue (0.00101456 -33.14) (0.00103105 -33.147) (0.0010478 -33.154) (0.00106483 -33.161) (0.00108213 -33.168) -continue (0.00109971 -33.175) (0.00111758 -33.182) (0.00113574 -33.189) (0.0011542 -33.196) (0.00117295 -33.203) -continue (0.00119201 -33.21) (0.00121138 -33.2169) (0.00123106 -33.2236) (0.00125106 -33.2294) (0.00127139 -33.2321) -continue (0.00129205 -33.2245) (0.00131304 -33.1847) (0.00133438 -33.0632) (0.00135606 -32.8038) (0.0013781 -32.4153) -continue (0.00140049 -31.9597) (0.00142324 -31.483) (0.00144637 -31.0045) (0.00146987 -30.5306) (0.00149375 -30.0633) -continue (0.00151803 -29.6033) (0.00154269 -29.1506) (0.00156776 -28.7052) (0.00159323 -28.267) (0.00161912 -27.8359) -continue (0.00164543 -27.4117) (0.00167217 -26.9944) (0.00169934 -26.5838) (0.00172695 -26.1799) (0.00175501 -25.7825) -continue (0.00178353 -25.3915) (0.00181251 -25.0068) (0.00184196 -24.6284) (0.00187189 -24.2561) (0.0019023 -23.8899) -continue (0.00193321 -23.5295) (0.00196462 -23.175) (0.00199655 -22.8263) (0.00202899 -22.4832) (0.00206196 -22.1456) -continue (0.00209546 -21.8135) (0.00212951 -21.4869) (0.00216411 -21.1655) (0.00219927 -20.8493) (0.00223501 -20.5383) -continue (0.00227133 -20.2323) (0.00230823 -19.9312) (0.00234574 -19.6351) (0.00238385 -19.3438) (0.00242259 -19.0572) -continue (0.00246195 -18.7752) (0.00250195 -18.4979) (0.00254261 -18.225) (0.00258392 -17.9566) (0.00262591 -17.6925) -continue (0.00266857 -17.4328) (0.00271194 -17.1773) (0.002756 -16.9259) (0.00280078 -16.6787) (0.00284629 -16.4354) -continue (0.00289254 -16.1961) (0.00293954 -15.9608) (0.0029873 -15.7292) (0.00303584 -15.5015) (0.00308517 -15.2774) -continue (0.0031353 -15.0571) (0.00318625 -14.8403) (0.00323802 -14.627) (0.00329063 -14.4173) (0.0033441 -14.2109) -continue (0.00339844 -14.008) (0.00345366 -13.8084) (0.00350978 -13.612) (0.0035668 -13.4188) (0.00362476 -13.2288) -continue (0.00368366 -13.042) (0.00374351 -12.8581) (0.00380434 -12.6773) (0.00386616 -12.4995) (0.00392898 -12.3246) -continue (0.00399282 -12.1525) (0.00405769 -11.9833) (0.00412363 -11.8169) (0.00419063 -11.6532) (0.00425872 -11.4922) -continue (0.00432792 -11.3338) (0.00439824 -11.1781) (0.00446971 -11.0249) (0.00454234 -10.8742) (0.00461614 -10.726) -continue (0.00469115 -10.5803) (0.00476738 -10.437) (0.00484484 -10.296) (0.00492356 -10.1574) (0.00500356 -10.021) -continue (0.00508486 -9.88693) (0.00516749 -9.75506) (0.00525145 -9.62537) (0.00533678 -9.49783) (0.0054235 -9.37241) -continue (0.00551162 -9.24907) (0.00560118 -9.12777) (0.00569219 -9.00849) (0.00578468 -8.89119) (0.00587867 -8.77584) -continue (0.00597419 -8.66241) (0.00607127 -8.55087) (0.00616992 -8.44118) (0.00627017 -8.33333) (0.00637205 -8.22728) -continue (0.00647559 -8.12299) (0.00658081 -8.02045) (0.00668774 -7.91962) (0.00679641 -7.82048) (0.00690684 -7.723) -continue (0.00701907 -7.62716) (0.00713312 -7.53292) (0.00724902 -7.44026) (0.00736681 -7.34916) (0.00748651 -7.25959) -continue (0.00760815 -7.17153) (0.00773178 -7.08495) (0.00785741 -6.99983) (0.00798508 -6.91614) (0.00811483 -6.83387) -continue (0.00824668 -6.75299) (0.00838068 -6.67348) (0.00851686 -6.59531) (0.00865524 -6.51847) (0.00879588 -6.44293) -continue (0.0089388 -6.36867) (0.00908405 -6.29568) (0.00923165 -6.22392) (0.00938165 -6.15339) (0.00953409 -6.08405) -continue (0.00968901 -6.0159) (0.00984644 -5.94892) (0.0100064 -5.88307) (0.010169 -5.81835) (0.0103343 -5.75474) -continue (0.0105022 -5.69221) (0.0106728 -5.63075) (0.0108462 -5.57034) (0.0110225 -5.51097) (0.0112016 -5.45261) -continue (0.0113836 -5.39524) (0.0115686 -5.33885) (0.0117565 -5.28342) (0.0119476 -5.22894) (0.0121417 -5.17537) -continue (0.012339 -5.12271) (0.0125395 -5.07094) (0.0127432 -5.02003) (0.0129503 -4.96997) (0.0131607 -4.92074) -continue (0.0133746 -4.87231) (0.0135919 -4.82466) (0.0138127 -4.77778) (0.0140372 -4.73164) (0.0142652 -4.68621) -continue (0.014497 -4.64148) (0.0147326 -4.59742) (0.014972 -4.554) (0.0152153 -4.5112) (0.0154625 -4.46899) -continue (0.0157137 -4.42733) (0.0159691 -4.38621) (0.0162285 -4.34558) (0.0164922 -4.30543) (0.0167602 -4.26571) -continue (0.0170325 -4.22639) (0.0173093 -4.18744) (0.0175905 -4.14882) (0.0178764 -4.1105) (0.0181668 -4.07245) -continue (0.018462 -4.03462) (0.018762 -3.99698) (0.0190669 -3.95951) (0.0193767 -3.92216) (0.0196915 -3.8849) -continue (0.0200115 -3.8477) (0.0203366 -3.81053) (0.0206671 -3.77337) (0.0210029 -3.73619) (0.0213442 -3.69898) -continue (0.021691 -3.6617) (0.0220434 -3.62434) (0.0224016 -3.5869) (0.0227656 -3.54937) (0.0231355 -3.51173) -continue (0.0235114 -3.474) (0.0238935 -3.43617) (0.0242817 -3.39824) (0.0246763 -3.36023) (0.0250772 -3.32215) -continue (0.0254847 -3.28402) (0.0258988 -3.24585) (0.0263196 -3.20766) (0.0267473 -3.16948) (0.0271819 -3.13132) -continue (0.0276235 -3.09322) (0.0280724 -3.05519) (0.0285285 -3.01727) (0.0289921 -2.97947) (0.0294632 -2.94183) -continue (0.0299419 -2.90438) (0.0304284 -2.86712) (0.0309228 -2.8301) (0.0314253 -2.79334) (0.0319359 -2.75685) -continue (0.0324548 -2.72065) (0.0329822 -2.68478) (0.0335181 -2.64924) (0.0340627 -2.61406) (0.0346162 -2.57925) -continue (0.0351787 -2.54483) (0.0357503 -2.5108) (0.0363312 -2.47719) (0.0369215 -2.444) (0.0375214 -2.41124) -continue (0.0381311 -2.37893) (0.0387507 -2.34707) (0.0393803 -2.31566) (0.0400202 -2.28472) (0.0406705 -2.25425) -continue (0.0413313 -2.22425) (0.0420029 -2.19472) (0.0426854 -2.16567) (0.043379 -2.1371) (0.0440838 -2.10901) -continue (0.0448001 -2.0814) (0.0455281 -2.05427) (0.0462679 -2.02763) (0.0470196 -2.00146) (0.0477837 -1.97577) -continue (0.0485601 -1.95055) (0.0493491 -1.92581) (0.050151 -1.90154) (0.0509659 -1.87773) (0.051794 -1.85439) -continue (0.0526356 -1.83151) (0.0534908 -1.80909) (0.05436 -1.78713) (0.0552433 -1.76561) (0.0561409 -1.74454) -continue (0.0570531 -1.7239) (0.0579802 -1.70371) (0.0589223 -1.68394) (0.0598797 -1.66461) (0.0608526 -1.6457) -continue (0.0618414 -1.6272) (0.0628462 -1.60912) (0.0638674 -1.59144) (0.0649052 -1.57417) (0.0659598 -1.5573) -continue (0.0670316 -1.54082) (0.0681207 -1.52472) (0.0692276 -1.50902) (0.0703525 -1.49369) (0.0714956 -1.47873) -continue (0.0726573 -1.46414) (0.0738379 -1.44991) (0.0750377 -1.43605) (0.0762569 -1.42253) (0.077496 -1.40937) -continue (0.0787552 -1.39655) (0.0800349 -1.38406) (0.0813354 -1.37192) (0.0826569 -1.3601) (0.084 -1.3486) -continue (0.0853649 -1.33743) (0.086752 -1.32657) (0.0881616 -1.31602) (0.0895941 -1.30578) (0.0910499 -1.29584) -continue (0.0925293 -1.2862) (0.0940328 -1.27685) (0.0955607 -1.26779) (0.0971134 -1.25901) (0.0986914 -1.25052) -continue (0.100295 -1.2423) (0.101925 -1.23435) (0.103581 -1.22667) (0.105264 -1.21925) (0.106974 -1.2121) -continue (0.108712 -1.2052) (0.110479 -1.19855) (0.112274 -1.19215) (0.114098 -1.18599) (0.115952 -1.18008) -continue (0.117836 -1.1744) (0.119751 -1.16895) (0.121697 -1.16373) (0.123674 -1.15874) (0.125684 -1.15398) -continue (0.127726 -1.14943) (0.129801 -1.1451) (0.13191 -1.14098) (0.134054 -1.13707) (0.136232 -1.13336) -continue (0.138446 -1.12986) (0.140695 -1.12656) (0.142981 -1.12346) (0.145305 -1.12054) (0.147666 -1.11782) -continue (0.150065 -1.11529) (0.152503 -1.11294) (0.154981 -1.11077) (0.1575 -1.10879) (0.160059 -1.10697) -continue (0.162659 -1.10534) (0.165302 -1.10387) (0.167988 -1.10257) (0.170718 -1.10143) (0.173492 -1.10046) -continue (0.176311 -1.09965) (0.179176 -1.099) (0.182087 -1.0985) (0.185046 -1.09815) (0.188053 -1.09796) -continue (0.191108 -1.09791) (0.194213 -1.09801) (0.197369 -1.09825) (0.200576 -1.09863) (0.203835 -1.09916) -continue (0.207147 -1.09982) (0.210513 -1.10061) (0.213934 -1.10154) (0.21741 -1.10259) (0.220943 -1.10378) -continue (0.224533 -1.10509) (0.228181 -1.10652) (0.231889 -1.10808) (0.235656 -1.10976) (0.239486 -1.11156) -continue (0.243377 -1.11347) (0.247331 -1.1155) (0.25135 -1.11765) (0.255434 -1.1199) (0.259585 -1.12226) -continue (0.263803 -1.12473) (0.268089 -1.12731) (0.272445 -1.12999) (0.276872 -1.13278) (0.281371 -1.13566) -continue (0.285943 -1.13865) (0.290589 -1.14173) (0.295311 -1.14491) (0.300109 -1.14819) (0.304986 -1.15156) -continue (0.309941 -1.15502) (0.314977 -1.15857) (0.320095 -1.16221) (0.325296 -1.16594) (0.330582 -1.16975) -continue (0.335954 -1.17365) (0.341412 -1.17764) (0.34696 -1.1817) (0.352598 -1.18585) (0.358327 -1.19008) -continue (0.364149 -1.19439) (0.370066 -1.19877) (0.376079 -1.20323) (0.38219 -1.20776) (0.3884 -1.21237) -continue (0.394711 -1.21705) (0.401125 -1.22181) (0.407642 -1.22663) (0.414266 -1.23152) (0.420997 -1.23648) -continue (0.427838 -1.24151) (0.43479 -1.24661) (0.441855 -1.25177) (0.449034 -1.25699) (0.45633 -1.26228) -continue (0.463745 -1.26763) (0.47128 -1.27304) (0.478938 -1.27851) (0.48672 -1.28404) (0.494629 -1.28963) -continue (0.502666 -1.29527) (0.510833 -1.30098) (0.519134 -1.30674) (0.527569 -1.31255) (0.536141 -1.31842) -continue (0.544853 -1.32434) (0.553706 -1.33031) (0.562703 -1.33634) (0.571846 -1.34241) (0.581138 -1.34854) -continue (0.590581 -1.35472) (0.600177 -1.36094) (0.609929 -1.36721) (0.61984 -1.37353) (0.629911 -1.3799) -continue (0.640146 -1.38631) (0.650548 -1.39277) (0.661119 -1.39927) (0.671861 -1.40582) (0.682778 -1.41241) -continue (0.693872 -1.41904) (0.705146 -1.42571) (0.716604 -1.43242) (0.728248 -1.43918) (0.740081 -1.44597) -continue (0.752107 -1.45281) (0.764327 -1.45968) (0.776747 -1.46659) (0.789368 -1.47354) (0.802194 -1.48053) -continue (0.815228 -1.48755) (0.828475 -1.49461) (0.841937 -1.50171) (0.855617 -1.50884) (0.86952 -1.51601) -continue (0.883648 -1.52321) (0.898006 -1.53044) (0.912598 -1.53771) (0.927426 -1.54501) (0.942496 -1.55234) -continue (0.95781 -1.5597) (0.973373 -1.5671) (0.989189 -1.57453) (1.00526 -1.58198) (1.0216 -1.58947) -continue (1.0382 -1.59699) (1.05507 -1.60453) (1.07221 -1.61211) (1.08963 -1.61971) (1.10734 -1.62735) -continue (1.12533 -1.63501) (1.14361 -1.6427) (1.1622 -1.65041) (1.18108 -1.65816) (1.20027 -1.66593) -continue (1.21977 -1.67372) (1.23959 -1.68154) (1.25974 -1.68939) (1.2802 -1.69726) (1.30101 -1.70516) -continue (1.32215 -1.71308) (1.34363 -1.72103) (1.36546 -1.729) (1.38765 -1.737) (1.4102 -1.74501) -continue (1.43311 -1.75306) (1.4564 -1.76112) (1.48006 -1.76921) (1.50411 -1.77732) (1.52855 -1.78545) -continue (1.55339 -1.79361) (1.57863 -1.80178) (1.60428 -1.80998) (1.63034 -1.8182) (1.65684 -1.82644) -continue (1.68376 -1.83471) (1.71112 -1.84299) (1.73892 -1.85129) (1.76717 -1.85962) (1.79589 -1.86796) -continue (1.82507 -1.87633) (1.85472 -1.88471) (1.88486 -1.89312) (1.91549 -1.90154) (1.94661 -1.90999) -continue (1.97824 -1.91845) (2.01039 -1.92693) (2.04305 -1.93543) (2.07625 -1.94395) (2.10998 -1.95249) -continue (2.14427 -1.96105) (2.17911 -1.96963) (2.21452 -1.97822) (2.2505 -1.98683) (2.28707 -1.99546) -continue (2.32423 -2.00411) (2.362 -2.01278) (2.40038 -2.02146) (2.43938 -2.03017) (2.47902 -2.03889) -continue (2.5193 -2.04763) (2.56023 -2.05638) (2.60183 -2.06515) (2.64411 -2.07394) (2.68707 -2.08275) -continue (2.73073 -2.09158) (2.7751 -2.10042) (2.8202 -2.10928) (2.86602 -2.11815) (2.91259 -2.12705) -continue (2.95992 -2.13596) (3.00801 -2.14489) (3.05689 -2.15383) (3.10656 -2.16279) (3.15703 -2.17177) -continue (3.20833 -2.18076) (3.26046 -2.18977) (3.31344 -2.1988) (3.36728 -2.20785) (3.42199 -2.21691) -continue (3.4776 -2.22599) (3.5341 -2.23508) (3.59153 -2.24419) (3.64989 -2.25332) (3.70919 -2.26247) -continue (3.76946 -2.27163) (3.83071 -2.28081) (3.89295 -2.29) (3.95621 -2.29922) (4.02049 -2.30845) -continue (4.08582 -2.31769) (4.15221 -2.32695) (4.21968 -2.33623) (4.28824 -2.34553) (4.35792 -2.35484) -continue (4.42873 -2.36417) (4.50069 -2.37352) (4.57382 -2.38288) (4.64814 -2.39227) (4.72367 -2.40166) -continue (4.80042 -2.41108) (4.87842 -2.42051) (4.95769 -2.42996) (5.03825 -2.43943) (5.12011 -2.44892) -continue (5.20331 -2.45842) (5.28785 -2.46794) (5.37377 -2.47747) (5.46109 -2.48703) (5.54983 -2.4966) -continue (5.64 -2.50619) (5.73165 -2.5158) (5.82478 -2.52543) (5.91942 -2.53507) (6.01561 -2.54473) -continue (6.11335 -2.55441) (6.21269 -2.56411) (6.31363 -2.57383) (6.41622 -2.58356) (6.52048 -2.59331) -continue (6.62643 -2.60308) (6.7341 -2.61287) (6.84352 -2.62268) (6.95472 -2.63251) (7.06772 -2.64236) -continue (7.18256 -2.65222) (7.29927 -2.6621) (7.41787 -2.67201) (7.5384 -2.68193) (7.66089 -2.69187) -continue (7.78537 -2.70183) (7.91187 -2.71181) (8.04043 -2.7218) (8.17108 -2.73182) (8.30385 -2.74186) -continue (8.43877 -2.75191) (8.57589 -2.76199) (8.71524 -2.77208) (8.85685 -2.7822) (9.00076 -2.79233) -continue (9.14701 -2.80248) (9.29564 -2.81266) (9.44668 -2.82285) (9.60018 -2.83306) (9.75617 -2.84329) -continue (9.91469 -2.85354) (10.0758 -2.86381) (10.2395 -2.8741) (10.4059 -2.88442) (10.575 -2.89475) -continue (10.7468 -2.90509) (10.9214 -2.91546) (11.0989 -2.92585) (11.2792 -2.93626) (11.4625 -2.94669) -continue (11.6488 -2.95714) (11.838 -2.9676) (12.0304 -2.97809) (12.2259 -2.98859) (12.4245 -2.99911) -continue (12.6264 -3.00966) (12.8316 -3.02022) (13.0401 -3.03079) (13.2519 -3.04139) (13.4673 -3.05201) -continue (13.6861 -3.06264) (13.9085 -3.07329) (14.1345 -3.08396) (14.3641 -3.09464) (14.5975 -3.10535) -continue (14.8347 -3.11606) (15.0758 -3.1268) (15.3207 -3.13755) (15.5697 -3.14832) (15.8227 -3.1591) -continue (16.0798 -3.16989) (16.341 -3.1807) (16.6065 -3.19153) (16.8764 -3.20237) (17.1506 -3.21322) -continue (17.4293 -3.22408) (17.7125 -3.23495) (18.0003 -3.24584) (18.2928 -3.25674) (18.59 -3.26764) -continue (18.8921 -3.27856) (19.199 -3.28948) (19.511 -3.30042) (19.828 -3.31135) (20.1502 -3.3223) -continue (20.4776 -3.33325) (20.8103 -3.34421) (21.1485 -3.35516) (21.4921 -3.36612) (21.8413 -3.37709) -continue (22.1962 -3.38805) (22.5569 -3.39901) (22.9234 -3.40997) (23.2959 -3.42093) (23.6744 -3.43189) -continue (24.0591 -3.44284) (24.45 -3.45378) (24.8473 -3.46472) (25.251 -3.47564) (25.6613 -3.48656) -continue (26.0783 -3.49746) (26.502 -3.50835) (26.9327 -3.51923) (27.3703 -3.53009) (27.815 -3.54093) -continue (28.267 -3.55175) (28.7263 -3.56255) (29.193 -3.57333) (29.6674 -3.58408) (30.1494 -3.59481) -continue (30.6393 -3.60551) (31.1372 -3.61617) (31.6431 -3.62681) (32.1573 -3.63741) (32.6798 -3.64798) -continue (33.2108 -3.65851) (33.7504 -3.669) (34.2988 -3.67945) (34.8561 -3.68985) (35.4225 -3.70021) -continue (35.9981 -3.71053) (36.583 -3.72079) (37.1774 -3.731) (37.7815 -3.74116) (38.3954 -3.75127) -continue (39.0193 -3.76131) (39.6533 -3.7713) (40.2976 -3.78123) (40.9524 -3.79109) (41.6178 -3.80089) -continue (42.2941 -3.81062) (42.9813 -3.82028) (43.6797 -3.82988) (44.3894 -3.8394) (45.1107 -3.84884) -continue (45.8437 -3.85822) (46.5886 -3.86751) (47.3456 -3.87673) (48.1149 -3.88586) (48.8967 -3.89492) -continue (49.6912 -3.90389) (50.4986 -3.91278) (51.3191 -3.92159) (52.153 -3.93031) (53.0004 -3.93894) -continue (53.8616 -3.94749) (54.7368 -3.95595) (55.6262 -3.96431) (56.53 -3.97259) (57.4486 -3.98078) -continue (58.3821 -3.98888) (59.3307 -3.99689) (60.2947 -4.00481) (61.2744 -4.01264) (62.2701 -4.02038) -continue (63.2819 -4.02803) (64.3101 -4.03558) (65.3551 -4.04305) (66.417 -4.05043) (67.4962 -4.05773) -continue (68.5929 -4.06493) (69.7075 -4.07205) (70.8401 -4.07908) (71.9912 -4.08603) (73.161 -4.09289) -continue (74.3497 -4.09968) (75.5578 -4.10638) (76.7855 -4.113) (78.0332 -4.11954) (79.3011 -4.12601) -continue (80.5897 -4.1324) (81.8991 -4.13872) (83.2299 -4.14496) (84.5823 -4.15114) (85.9566 -4.15725) -continue (87.3533 -4.16329) (88.7727 -4.16926) (90.2151 -4.17517) (91.681 -4.18103) (93.1707 -4.18682) -continue (94.6846 -4.19255) (96.2231 -4.19823) (97.7866 -4.20386) (99.3755 -4.20943) (100.99 -4.21496) -continue (102.631 -4.22044) (104.299 -4.22587) (105.994 -4.23125) (107.716 -4.2366) (109.466 -4.2419) -continue (111.245 -4.24717) (113.052 -4.2524) (114.889 -4.25759) (116.756 -4.26275) (118.653 -4.26788) -continue (120.581 -4.27298) (122.54 -4.27805) (124.532 -4.28309) (126.555 -4.28811) (128.611 -4.29311) -continue (130.701 -4.29808) (132.825 -4.30303) (134.983 -4.30795) (137.176 -4.31286) (139.405 -4.31776) -continue (141.67 -4.32263) (143.972 -4.32749) (146.312 -4.33234) (148.689 -4.33717) (151.105 -4.34199) -continue (153.56 -4.3468) (156.056 -4.3516) (158.591 -4.35638) (161.168 -4.36116) (163.787 -4.36593) -continue (166.448 -4.37069) (169.153 -4.37545) (171.901 -4.38019) (174.695 -4.38493) (177.533 -4.38967) -continue (180.418 -4.3944) (183.349 -4.39913) (186.329 -4.40385) (189.356 -4.40857) (192.433 -4.41328) -continue (195.56 -4.418) (198.737 -4.4227) (201.967 -4.42741) (205.248 -4.43212) (208.583 -4.43682) -continue (211.972 -4.44152) (215.417 -4.44622) (218.917 -4.45091) (222.474 -4.45561) (226.089 -4.4603) -continue (229.763 -4.465) (233.496 -4.46969) (237.29 -4.47438) (241.146 -4.47908) (245.064 -4.48377) -continue (249.046 -4.48846) (253.093 -4.49315) (257.205 -4.49784) (261.384 -4.50253) (265.631 -4.50722) -continue (269.948 -4.5119) (274.334 -4.51659) (278.791 -4.52128) (283.321 -4.52597) (287.925 -4.53066) -continue (292.603 -4.53535) (297.358 -4.54004) (302.189 -4.54472) (307.1 -4.54941) (312.09 -4.5541) -continue (317.161 -4.55879) (322.314 -4.56348) (327.551 -4.56816) (332.874 -4.57285) (338.282 -4.57754) -continue (343.779 -4.58223) (349.365 -4.58692) (355.042 -4.5916) (360.811 -4.59629) (366.673 -4.60098) -continue (372.631 -4.60567) (378.686 -4.61036) (384.839 -4.61504) (391.092 -4.61973) (397.447 -4.62442) -continue (403.905 -4.62911) (410.468 -4.6338) (417.138 -4.63849) (423.916 -4.64317) (430.804 -4.64786) -continue (437.804 -4.65255) (444.917 -4.65724) (452.147 -4.66193) (459.493 -4.66662) (466.96 -4.6713) -continue (474.547 -4.67599) (482.258 -4.68068) (490.094 -4.68537) (498.057 -4.69006) (506.15 -4.69475) -continue (514.374 -4.69944) (522.732 -4.70412) (531.226 -4.70881) (539.858 -4.7135) (548.63 -4.71819) -continue (557.544 -4.72288) (566.604 -4.72757) (575.81 -4.73226) (585.166 -4.73695) (594.675 -4.74163) -continue (604.337 -4.74632) (614.157 -4.75101) (624.136 -4.7557) (634.278 -4.76039) (644.584 -4.76508) -continue (655.057 -4.76977) (665.701 -4.77446) (676.518 -4.77915) (687.511 -4.78383) (698.682 -4.78852) -continue (710.034 -4.79321) (721.571 -4.7979) (733.296 -4.80259) (745.211 -4.80728) (757.32 -4.81197) -continue (769.625 -4.81666) (782.131 -4.82135) (794.839 -4.82604) (807.754 -4.83073) (820.879 -4.83541) -continue (834.218 -4.8401) (847.773 -4.84479) (861.548 -4.84948) (875.547 -4.85417) (889.773 -4.85886) -continue (904.231 -4.86355) (918.923 -4.86824) (933.855 -4.87293) (949.029 -4.87762) (964.449 -4.88231) -continue (980.12 -4.887) (996.046 -4.89169) (1012.23 -4.89637) (1028.68 -4.90106) (1045.39 -4.90575) -continue (1062.38 -4.91044) (1079.64 -4.91513) (1097.18 -4.91982) (1115.01 -4.92451) (1133.13 -4.9292) -continue (1151.54 -4.93389) (1170.25 -4.93858) (1189.27 -4.94327) (1208.59 -4.94796) (1228.23 -4.95265) -continue (1248.19 -4.95734) (1268.47 -4.96203) (1289.08 -4.96672) (1310.02 -4.9714) (1331.31 -4.97609) -continue (1352.94 -4.98078) (1374.93 -4.98547) (1397.27 -4.99016) (1419.97 -4.99485) (1443.04 -4.99954) -continue (1466.49 -5.00423) (1490.32 -5.00892) (1514.54 -5.01361) (1539.14 -5.0183) (1564.15 -5.02299) -continue (1589.57 -5.02768) (1615.4 -5.03237) (1641.65 -5.03706) (1668.32 -5.04175) (1695.43 -5.04644) -continue (1722.98 -5.05113) (1750.97 -5.05582) (1779.42 -5.06051) (1808.34 -5.0652) (1837.72 -5.06989) -continue (1867.58 -5.07457) (1897.93 -5.07926) (1928.77 -5.08395) (1960.11 -5.08864) (1991.95 -5.09333) -continue (2024.32 -5.09802) (2057.21 -5.10271) (2090.64 -5.1074) (2124.61 -5.11209) (2159.13 -5.11678) -continue (2194.22 -5.12147) (2229.87 -5.12616) (2266.1 -5.13085) (2302.92 -5.13554) (2340.34 -5.14023) -continue (2378.37 -5.14492) (2417.02 -5.14961) (2456.29 -5.1543) (2496.2 -5.15899) (2536.76 -5.16368) -continue (2577.98 -5.16837) (2619.87 -5.17306) (2662.44 -5.17775) (2705.7 -5.18244) (2749.66 -5.18713) -continue (2794.34 -5.19182) (2839.75 -5.19651) (2885.89 -5.2012) (2932.78 -5.20589) (2980.43 -5.21058) -continue (3028.86 -5.21527) (3078.08 -5.21996) (3128.09 -5.22465) (3178.92 -5.22934) (3230.57 -5.23402) -continue (3283.06 -5.23871) (3336.41 -5.2434) (3390.62 -5.24809) (3445.71 -5.25278) (3501.7 -5.25747) -continue (3558.6 -5.26216) (3616.42 -5.26685) (3675.19 -5.27154) (3734.9 -5.27623) (3795.59 -5.28092) -continue (3857.26 -5.28561) (3919.94 -5.2903) (3983.63 -5.29499) (4048.36 -5.29968) (4114.14 -5.30437) -continue (4180.99 -5.30906) (4248.93 -5.31375) (4317.97 -5.31844) (4388.13 -5.32313) (4459.43 -5.32782) -continue (4531.89 -5.33251) (4605.53 -5.3372) (4680.36 -5.34189) (4756.41 -5.34658) (4833.7 -5.35127) -continue (4912.24 -5.35596) (4992.05 -5.36065) (5073.17 -5.36534) (5155.6 -5.37003) (5239.37 -5.37472) -continue (5324.51 -5.37941) (5411.02 -5.3841) (5498.94 -5.38879) (5588.3 -5.39348) (5679.1 -5.39817) -continue (5771.38 -5.40286) (5865.15 -5.40755) (5960.45 -5.41224) (6057.3 -5.41693) (6155.73 -5.42162) -continue (6255.75 -5.42631) (6357.4 -5.431) (6460.7 -5.43569) (6565.67 -5.44038) (6672.36 -5.44507) -continue (6780.78 -5.44976) (6890.95 -5.45445) (7002.92 -5.45914) (7116.71 -5.46383) (7232.35 -5.46852) -continue (7349.86 -5.47321)(73498.6 -7.47321)(734986 -9.47321)(7.34986e+06 -11.4732)(7.34986e+07 -13.4732) -// -////// Flux is varied -// -phi(h) 20.5 -//phi(h) 23 vary -//grid range from 17 to 24 with .25 dex -// -////// Stopping condition and run -// -stop total column density 23 -iterate to convergence max=40 diff --git a/scripts/cloudy/mehdipour.mpi_grid.in b/scripts/cloudy/mehdipour.mpi_grid.in index 518bb2c..b72d77a 100644 --- a/scripts/cloudy/mehdipour.mpi_grid.in +++ b/scripts/cloudy/mehdipour.mpi_grid.in @@ -1,23 +1,20 @@ -////// CLOUDY 17.00 4thdex grid using mehdipour ngc5548 SED +////// Cloudy 17.00 4thdex grid using mehdipour ngc5548 SED // -title mehdipour ngc5548 hden_vs_phi 4thdex -// -// single run copied from the standard mehdipour grid +title mehdipour ngc5548 4thdex grid // ////// Output continua and included lines // save last incident continuum "incident_continuum" save line labels "line_labels" // -////// Formatting +////// Report almost everything // print line column print line inward -print line iso collapsed +print line iso collapsed on print line faint -3 print last print line precision 6 -print line iso collapsed on Database print // ////// Resolved energy states for atomic models @@ -34,6 +31,7 @@ database He-like element helium resolved levels 15 ////// low-metal abundance for AGN testing ~2017 // metals 0.5 // element scale factor carbon 2 +// element scale factor nitrogen 2 // ////// Prevent cloudy from stopping prematurely // @@ -55,6 +53,9 @@ set nend 5000 // "print continuum" redundant with set nFnu in c17 // print continuum +set nFnu add 450 angstrom +set nFnu add 650 angstrom +set nFnu add 875 angstrom set nFnu add 800 angstrom set nFnu add 840 angstrom set nFnu add 850 angstrom @@ -315,6 +316,8 @@ set nFnu add 12850 angstrom set nFnu add 12900 angstrom set nFnu add 12950 angstrom set nFnu add 13000 angstrom +set nFnu add 800000 angstrom +set nFnu add 1000000 angstrom // ////// Print equivalent widths // diff --git a/src/create_fort_files.cpp b/src/create_fort_files.cpp index 76f6fcd..a3bb95c 100644 --- a/src/create_fort_files.cpp +++ b/src/create_fort_files.cpp @@ -22,7 +22,7 @@ int main(int argc, char const *argv[]) { std::cout << "Compiling table2ds for " << lines_to_print.size() - << " emission lines.\n"; + << " emission lines and continuum bins.\n"; std::list tables = agn::compile_eqwidth_tables(grid,lines_to_print,1215.00); diff --git a/src/spectral_lines.hpp b/src/spectral_lines.hpp index 0f7ce81..908c389 100644 --- a/src/spectral_lines.hpp +++ b/src/spectral_lines.hpp @@ -191,6 +191,96 @@ std::ostream& agn::operator<< ( std::ostream& outstream, } +// Returns a table as the difference between two tables +//agn::eqwidth_table main(int argc, char const *argv[]) +//{ // Copied from fortfile subtractor +// std::ifstream table1; +// table1.open(argv[1]); +// std::ifstream table2; +// table2.open(argv[2]); +// +// std::string line; +// +// int table1_numlines = 0; +// while (std::getline(table1, line)) +// ++table1_numlines; +// table1.clear(); +// table1.seekg(0); +// +// int table_entries = table1_numlines - 2; +// +//// int table2_numlines = 0; +//// while (std::getline(table1, line)) +//// ++table2_numlines; +//// +//// int table2_entries = table2_numlines - 2; +//// table2.clear(); +//// table2.seekg(0); +// +// double hden1[table_entries]; +// double phi1[table_entries]; +// double magnitude1[table_entries]; +// double resultant_magnitude[table_entries]; +// +// std::string first_line, second_line, scratch; +// getline(table1,first_line); +// getline(table1,second_line); +// getline(table2,scratch); +// getline(table2,scratch); +// +// first_line[0] = argv[3][0]; +// first_line[1] = argv[3][1]; +// first_line[2] = argv[3][2]; +// first_line[3] = argv[3][3]; +// +// for (int i=0; i> hden1[i] >> phi1[i] >> magnitude1 [i]; +// } +// +// +// int i=0; +// double hden2, phi2, magnitude2; +// for (int i=0; i> hden2 >> phi2 >> magnitude2; +// resultant_magnitude[i] = magnitude1[i] - magnitude2; +// } +// if ( resultant_magnitude[i] < 1e-10 ) +// resultant_magnitude[i]=1e-10; +// } +// +// +// +// std::cout +// << first_line +// << std::endl +// //<< "HDen Phi(H) ScaInten" +// << second_line +// << std::endl; +// +// for (int i=0; i agn::compile_eqwidth_tables(agn::cloudy_grid grid, << std::endl; std::stringstream header0; header0 - << label - << " relative to Inci 1215.00A"; + << label; new_table.header[0] = header0.str(); new_table.header[1] = "Hden Phi(H) F(line)/Inci(1215) "; agn::cloudy_grid::iterator result_it = grid.begin(); @@ -293,6 +382,48 @@ std::list agn::compile_eqwidth_tables(agn::cloudy_grid grid, label.substr(agn::LABELSTR_POS,3) != "nFn" && label.substr(agn::LABELSTR_POS,3) != "Inw" && label.substr(LABELSTR_LEN - 14) != " 0 ") { + // These lines are known to not have inward fractions, so + // should be skipped. + if (label == "Blnd 1657.00A" || + label == "Blnd 2326.00A" || + label == "Blnd 2798.00A" || + label == "Blnd 1335.00A" || + label == "C 3 977.000A" || + label == "Blnd 1909.00A" || + label == "Blnd 1176.00A" || + label == "Blnd 2141.00A" || + label == "Blnd 1085.00A" || + label == "Blnd 1750.00A" || + label == "Blnd 990.000A" || + label == "Blnd 1486.00A" || + label == "Blnd 765.000A" || + label == "Blnd 1304.00A" || + label == "Blnd 8446.00A" || + label == "Blnd 1666.00A" || + label == "Blnd 835.000A" || + label == "Blnd 1402.00A" || + label == "Blnd 630.000A" || + label == "Blnd 1218.00A" || + label == "Blnd 774.000A" || + label == "Blnd 5892.00A" || + label == "Blnd 615.000A" || + label == "Blnd 2665.00A" || + label == "Blnd 1860.00A" || + label == "Blnd 2335.00A" || + label == "Blnd 1888.00A" || + label == "Blnd 1397.00A" || + label == "Blnd 1256.00A" || + label == "Blnd 1720.00A" || + label == "Blnd 1198.00A" || + label == "Blnd 1406.00A" || + label == "Blnd 1199.00A" || + label == "Blnd 8579.00A" || + label == "Blnd 3933.00A" + ) { + line_label_it++; + continue; + } + if(agn::debug) std::cout << "Processing Inwd for label " << label @@ -301,8 +432,7 @@ std::list agn::compile_eqwidth_tables(agn::cloudy_grid grid, agn::eqwidth_table inwd_table; inwdheader0 << "Inwd " - << label.substr(5) - << " relative to Inci 1215.00A"; + << label.substr(5); inwd_table.header[0] = inwdheader0.str(); inwd_table.header[1] = "Hden Phi(H) F(line)/Inci(1215)"; result_it = grid.begin();