added subroutine to identify nFnu lines and add nFnu, nInu, inwT, and inwC in such cases.

This commit is contained in:
caes 2016-05-12 22:48:04 -04:00
parent e9cdd548d1
commit 77cd7e0333
3 changed files with 55 additions and 64 deletions

View File

@ -189,119 +189,55 @@ FeKa 1.78000A
#Diffuse continua:
nFnu 4.11370m
nInu 4.11370m
inwT 4.11370m
inwC 4.11370m
# need to add in remaining inwT and inwC
nFnu 2.32440m
nInu 2.32440m
nFnu 1.64313m
nInu 1.64313m
nFnu 1.48022m
nInu 1.48022m
nFnu 1.43649m
nInu 1.43649m
nFnu 1.11867m
nInu 1.11867m
nFnu 8329.68A
nInu 8329.68A
nFnu 8078.61A
nInu 8078.61A
nFnu 6209.66A
nInu 6209.66A
nFnu 4885.36A
nInu 4885.36A
nFnu 3704.34A
nInu 3704.34A
nFnu 3587.67A
nInu 3587.67A
nFnu 2430.05A
nInu 2430.05A
nFnu 2392.03A
nInu 2392.03A
nFnu 2071.34A
nInu 2071.34A
nFnu 2052.68A
nInu 2052.68A
nFnu 1793.44A
nInu 1793.44A
nFnu 1585.12A
nInu 1585.12A
nFnu 1458.33A
nInu 1458.33A
nFnu 1356.96A
nInu 1356.96A
nFnu 1297.29A
nInu 1297.29A
nFnu 1262.79A
nInu 1262.79A
nInu 1218.08A
nFnu 1196.33A
nInu 1196.33A
nFnu 1175.06A
nInu 1175.06A
nFnu 1143.64A
nInu 1143.64A
nFnu 1113.22A
nInu 1113.22A
nFnu 1103.23A
nInu 1103.23A
nFnu 1073.86A
nInu 1073.86A
nFnu 1064.34A
nInu 1064.34A
nFnu 1035.92A
nInu 1035.92A
nFnu 1000.00A
nInu 1000.00A
nFnu 925.144A
nInu 925.144A
nFnu 897.800A
nInu 897.800A
nInu 760.023A
nInu 701.514A
nInu 608.160A
nInu 575.146A
nInu 511.948A
nInu 496.874A
nInu 399.153A
nFnu 910.000A
nInu 910.000A
nFnu 913.000A
nInu 913.000A
nFnu 1650.00A
nInu 1650.00A
nFnu 1900.00A
nInu 1900.00A
nFnu 3644.00A # changed this wavelength bin
nInu 3644.00A
nFnu 3645.00A
nInu 3645.00A
nFnu 3648.00A
nInu 3648.00A
nFnu 5100.00A
nInu 5100.00A
nFnu 6500.00A
nInu 6500.00A
nFnu 7000.00A
nInu 7000.00A
nFnu 7500.00A
nInu 7500.00A
nFnu 8000.00A
nInu 8000.00A
nFnu 8204.00A
nInu 8204.00A
nFnu 8207.00A
nInu 8207.00A
nFnu 8500.00A
nInu 8500.00A
nFnu 9000.00A
nInu 9000.00A
nFnu 9500.00A
nInu 9500.00A
nFnu 10000.0A
nInu 10000.0A
nFnu 1.05000m
nInu 1.05000m
nFnu 1.10000m
nInu 1.10000m

View File

@ -0,0 +1,47 @@
how many did we have that don't converge
optical depth
thermal structure
dr. korista to talk about new version of cloudy availability
AGN - other spectrum
want 10^22 10^23 10^24 with 18 15 15 models
get these results together ^
half decades would be good
numerical recipes
try 45-degree interpolation
abundances?
envelope trend?
quasars are episode
seyfert --
inward fraction of lines
define nfnu, ninu, inwC, inwT
SED go out to 200KeV

View File

@ -189,6 +189,14 @@ agn::line_list agn::read_line_list(std::ifstream& inputfile) {
<< "Adding label "
<< label
<< " to capture list.\n";
if (label.substr(0,4).compare(std::string("nFnu")) == 0) {
std::string applabel = label.substr(3);
list.push_back(std::string("nFnu").append(applabel));
list.push_back(std::string("nInu").append(applabel));
list.push_back(std::string("InwT").append(applabel));
list.push_back(std::string("InwC").append(applabel));
continue;
}
list.push_back(label);
}
return list;