Science/cloudy-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/cloudy-agn.git synced 2025-04-05 11:35:46 +00:00
Code Issues Packages Projects Releases Wiki Activity

updated create_fort_files function for c17

Browse Source
This commit is contained in:
caes 2017-06-20 09:28:17 -04:00
parent 528ee08c84
commit 22c0334c0d
2 changed files with 40 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
notes
View File

@ -15,7 +15,7 @@ Run grid with new SED, or at least a point, check values: how to find the approp
✓ faint to equivalent width of 0.01 angstrom ✓ faint to equivalent width of 0.001 angstrom
✓ compare between c13 and c17 -- show diferences in emergent line lists ✓ compare between c13 and c17 -- show diferences in emergent line lists
@ -27,9 +27,9 @@ Run grid with new SED, or at least a point, check values: how to find the approp
□ check if blnds are sums of the expected multiplets □ check if blnds are sums of the expected multiplets
□ add in more levels to the H atom, unresolved, 20, also add some to paschen, maybe others -- check how much time is lost ✓ add in more levels to the H atom, unresolved, 20, also add some to paschen, maybe others -- check how much time is lost
□ punch the paschen and balmer lines up to 18 ✓ punch the paschen and balmer lines up to 18
□ Read Intrinsic Lines □ Read Intrinsic Lines
@ -37,7 +37,7 @@ Run grid with new SED, or at least a point, check values: how to find the approp
□ Contour plots □ Contour plots
□ Run Grids in c17 ✓ Run Grids in c17
□ Email about order of lines □ Email about order of lines

55
src/agn.hpp
View File

@ -52,7 +52,7 @@ typedef std::map<std::string,cloudy_line_data> cloudy_line_output;
struct cloudy_result { struct cloudy_result {
std::string header, footer; std::string header, footer;
std::list<std::string> emergent_line_raw_text; std::list<std::string> intrinsic_line_raw_text;
std::list<std::string> intrinsic_line_raw_text; std::list<std::string> intrinsic_line_raw_text;
std::list<std::string> cautions; std::list<std::string> cautions;
cloudy_line_output emergent_line_intensity; cloudy_line_output emergent_line_intensity;
@ -62,7 +62,7 @@ struct cloudy_result {
cloudy_result(): cloudy_result():
header(""), header(""),
footer(""), footer(""),
emergent_line_raw_text(), intrinsic_line_raw_text(),
cautions(), cautions(),
iterations(0), iterations(0),
phi(0), phi(0),
@ -198,6 +198,7 @@ agn::cloudy_grid agn::read_cloudy_grid(std::ifstream& inputfile) {
<< std::endl; << std::endl;
agn::gridcoordlist coordlist; agn::gridcoordlist coordlist;
std::pair<double,double> xy (0,0); std::pair<double,double> xy (0,0);
// collect parameter coordinates
seek_string="**************************************************"; seek_string="**************************************************";
while(!inputfile.eof()) { while(!inputfile.eof()) {
while (inputline.find(seek_string) == std::string::npos) while (inputline.find(seek_string) == std::string::npos)
@ -234,10 +235,9 @@ agn::cloudy_grid agn::read_cloudy_grid(std::ifstream& inputfile) {
agn::cloudy_result point; agn::cloudy_result point;
point.hden = coords->first; point.hden = coords->first;
point.phi = coords->second; point.phi = coords->second;
seek_string = "Intrinsic line intensities";
if(agn::debug) std::cout if(agn::debug) std::cout
<< "Grabbing header and cautions."; << "Grabbing header and cautions.";
seek_string = "Intrinsic line intensities";
std::string header=""; std::string header="";
std::list<std::string> cautions; std::list<std::string> cautions;
while (inputline.find(seek_string) == std::string::npos) { while (inputline.find(seek_string) == std::string::npos) {
@ -284,19 +284,19 @@ agn::cloudy_grid agn::read_cloudy_grid(std::ifstream& inputfile) {
point.colden = colden; point.colden = colden;
if(agn::debug) std::cout if(agn::debug) std::cout
<< " Grabbing emission lines."; << " Grabbing emission lines.";
agn::seek_to("Emergent line intensities",inputfile); agn::seek_to("Intrinsic line intensities",inputfile);
std::list<std::string> emergent_line_raw_text; std::list<std::string> intrinsic_line_raw_text;
getline(inputfile,inputline); getline(inputfile,inputline);
while (inputline != "") { while (inputline != "") {
if(inputline.find("..") != std::string::npos) { if(inputline.find("..") != std::string::npos) {
getline(inputfile,inputline); getline(inputfile,inputline);
continue; continue;
} }
emergent_line_raw_text.push_back(inputline); intrinsic_line_raw_text.push_back(inputline);
getline(inputfile,inputline); getline(inputfile,inputline);
} }
point.emergent_line_raw_text.swap(emergent_line_raw_text); point.intrinsic_line_raw_text.swap(intrinsic_line_raw_text);
std::list<std::string>::iterator linetext_it=point.emergent_line_raw_text.begin(); std::list<std::string>::iterator linetext_it=point.intrinsic_line_raw_text.begin();
std::list<std::string> duplicate_labels; std::list<std::string> duplicate_labels;
int index=0; int index=0;
// For cloudy 17, the characters reserved for the output quantities: // For cloudy 17, the characters reserved for the output quantities:
@ -308,20 +308,37 @@ agn::cloudy_grid agn::read_cloudy_grid(std::ifstream& inputfile) {
// c17: 0,18; 18,9; 27, 10 // c17: 0,18; 18,9; 27, 10
// c13 settings // c13 settings
int str1_pos = 1; //int labelstr_pos = 1;
int str1_len = 13; //int labelstr_len = 13;
int str2_pos = 14; //int radiatedenergystr_pos = 14;
int str2_len = 8; //int radiatedenergystr_len = 8;
int str3_pos = 22; //int eqwidthstr_pos = 22;
int str3_len = 9; //int eqwidthstr_len = 9;
while(linetext_it != point.emergent_line_raw_text.end()) { // c17 settings
std::string label=(*linetext_it).substr(str1_pos,str1_len); int labelstr_pos = 0;
int labelstr_len = 18;
int radiatedenergystr_pos = 18;
int radiatedenergystr_len = 9;
int eqwidthstr_pos = 27;
int eqwidthstr_len = 10;
while(linetext_it != point.intrinsic_line_raw_text.end()) {
std::string label=(*linetext_it).substr(
labelstr_pos,
labelstr_len
);
agn::cloudy_line_data data; agn::cloudy_line_data data;
data.index = ++index; data.index = ++index;
std::stringstream values; std::stringstream values;
data.radiated_energy = atof((*linetext_it).substr(str2_pos,str2_len).c_str()); data.radiated_energy = atof((*linetext_it).substr(
data.eq_width = atof((*linetext_it).substr(str3_pos,str3_len).c_str()); radiatedenergystr_pos,
radiatedenergystr_len).c_str()
);
data.eq_width = atof((*linetext_it).substr(
eqwidthstr_pos,
eqwidthstr_len).c_str()
);
if(line_debug) { if(line_debug) {
std::cout << *linetext_it std::cout << *linetext_it
<< ": " << ": "