mirror of
https://asciireactor.com/otho/industry-website.git
synced 2024-11-22 18:55:05 +00:00
133 lines
3.9 KiB
C++
133 lines
3.9 KiB
C++
#include "agn.hpp"
|
|
#include "spectral_lines.hpp"
|
|
|
|
// When this number of iterations is achieved, the program flags it
|
|
// as a possible failure.
|
|
double MAX_ITERATIONS = 15;
|
|
|
|
|
|
int main(int argc, char const *argv[]) {
|
|
std::ifstream cloudy_result_file;
|
|
cloudy_result_file.open(argv[1]);
|
|
|
|
std::cout
|
|
<< "Reading cloudy grid from "
|
|
<< argv[1]
|
|
<< ".\n";
|
|
agn::cloudy_grid grid = agn::read_cloudy_grid(cloudy_result_file);
|
|
|
|
std::cout
|
|
<< "Reading line list from "
|
|
<< argv[2]
|
|
<< ".\n";
|
|
std::ifstream line_list_file;
|
|
line_list_file.open(argv[2]);
|
|
agn::line_list lines_to_print = agn::read_line_list(line_list_file);
|
|
std::cout
|
|
<< "Compiling table2ds for "
|
|
<< lines_to_print.size()
|
|
<< " emission lines and continuum bins.\n";
|
|
std::list<agn::eqwidth_table> tables =
|
|
agn::compile_eqwidth_tables(grid,lines_to_print,1215.00);
|
|
|
|
|
|
// Remove any tables that are zero.
|
|
std::list<agn::eqwidth_table>::iterator table_it = tables.begin();
|
|
std::ofstream zeroreport;
|
|
zeroreport.open("zero_report");
|
|
zeroreport << "These headers were pulled from tables"
|
|
<< " that returned minimum values."
|
|
<< std::endl;
|
|
int num_zeroes=0;
|
|
while(table_it != tables.end()) {
|
|
if(agn::is_zero(*table_it)) {
|
|
std::string linetype = table_it->header[0].substr(0,4);
|
|
//if (linetype == "nFnu" ||
|
|
// linetype == "nInu" ||
|
|
// linetype == "InwT" ||
|
|
// linetype == "InwC") {
|
|
// table_it++;
|
|
// continue;
|
|
//}
|
|
|
|
zeroreport << " "
|
|
<< table_it->header[0]
|
|
<< std::endl;
|
|
num_zeroes++;
|
|
//table_it = tables.erase(table_it);
|
|
table_it++; // needed since the erase() would have incremented
|
|
continue;
|
|
}
|
|
|
|
table_it++;
|
|
}
|
|
std::cout << "Removed (currently disabled) "
|
|
<< num_zeroes
|
|
<< " tables from the list because"
|
|
<< " they had zero value."
|
|
<< std::endl;
|
|
|
|
// Write the tables to the fortfile stack.
|
|
std::cout
|
|
<< "Printing "
|
|
<< tables.size()
|
|
<< " tables to fortfiles.\n";
|
|
table_it = tables.begin();
|
|
int fortfilenum=11;
|
|
while(table_it != tables.end()) {
|
|
std::ofstream outfile;
|
|
std::stringstream filename;
|
|
filename << "fort.";
|
|
filename << fortfilenum;
|
|
outfile.open(filename.str().c_str());
|
|
outfile << *table_it;
|
|
outfile.close();
|
|
table_it++;
|
|
fortfilenum++;
|
|
}
|
|
|
|
int num_unconverged = 0;
|
|
agn::cloudy_grid::iterator result_it = grid.begin();
|
|
std::ofstream cautionreportfile;
|
|
cautionreportfile.open("cautions");
|
|
cautionreportfile
|
|
<< "The following solutions probably did not converge."
|
|
<< std::endl << std::endl;
|
|
while(result_it != grid.end()) {
|
|
if (result_it->second.iterations >= MAX_ITERATIONS) {
|
|
num_unconverged++;
|
|
cautionreportfile
|
|
<< "hden = "
|
|
<< std::fixed
|
|
<< std::setprecision(3)
|
|
<< result_it->second.hden
|
|
<< ", phi = "
|
|
<< result_it->second.phi
|
|
<< std::endl
|
|
<< "───────────────────────────"
|
|
<< std::endl;
|
|
std::list<std::string>::iterator caution_it =
|
|
result_it->second.cautions.begin();
|
|
while(caution_it != result_it->second.cautions.end()) {
|
|
cautionreportfile
|
|
<< *caution_it
|
|
<< std::endl;
|
|
caution_it++;
|
|
}
|
|
cautionreportfile << std::endl << std::endl;
|
|
}
|
|
result_it++;
|
|
}
|
|
|
|
std::cout
|
|
<< "Saved cautions for "
|
|
<< num_unconverged
|
|
<< " unconverged solutions."
|
|
<< std::endl;
|
|
cautionreportfile.close();
|
|
|
|
std::cout << "Done.\n";
|
|
return 0;
|
|
}
|
|
|