BlueEdge_Plotter_V4.cpp

//V1: Currrently reading LINA data into arrays correctly, but if statement for outputting parameters at blue edge not working.
	//Note: if statement worked when using 2D arrays (reading all of LINA data) and inside of NumOfModes loop..however outputted 10 times..did not work outside of loop.
	//Note: From testing problem seems to be with outFile not the if statement? outFile worked when right before external loop test
	//Note: WORKAROUND - works when no header is added to outputfile. Why? Who knows.. Working program!
	
//V2: Goal is to add graphics to program so HR diagram of blue edge can be plotted.
	//Note: Graphics not working, Root isnt finding the points? Fixed this! Problem was in selectors 
	//Note: Root cant have x-axis go from highest to lowest? Maybe we dont use root?
	
//V3: Giving up on root graphics, this is V2 but without them.

//V4: Adding option to get all positive growth rate models so they can be plotted. 
	//Note: Reorganize program so the outside loop doesnt have to be done twice when getting positive growth and blue edge? Using more selectors inside instead of one outside?

# include <iostream>
# include <fstream>
# include <math.h>
# include <iomanip>
# include <cmath>
# include <stdlib.h>
# include <cstdlib>
//# include <fstream.h>
# include <string.h>
# include <string>
//# include <dos.h> //For Sleep() 


using namespace std;
int main(){
	ifstream inFile;
	
	//Controls//
	
	int NumOfModels = 20412; //20412 is total TypeII models
	int NumOfModes = 10;
	int NumOfSpaces = 1; // Number of spaces between each models output 
	
	char inputFileName[30] = "InputDataFULL.dat"; //File with the model parameters
	char outputFileName_ForBlueEdge[30] = "BlueEdge_SetA.dat"; //File with model parameters on the blue edge
	char outputFileName_ForPositiveGrowthModels[30] = "PostiveGrowth_SetD.dat"; //File with model parameters that have positve growth rates
	char LINAFileName[30] = "LINA_period_growth.data"; //Normally should be "LINA_period_growth.data"
	char log_directory_prefix[30] = "Raw_LINA_data/SetD/SetD/LOGS_"; //Prefix to log_directory, suffix is model number. This is where LINA file should be
	
	bool GetPositiveGrowthRateModels = true; // Positive Growth. When true, postive growth rate models will be outputted to a file.
	bool GetFirstNegativeGrowthRateModels = false; // Blue Edge. When true, the first negative growth rate models will be outputted to a file (For blue edge!)
	
	///////////////////////////////
	
	//Arrays//
	
	const int nArray = NumOfModels + 1;
	const int nArray_modes = NumOfModes + 1;
	
	double Mod[nArray];
	double z[nArray];
	double x[nArray];
	double M[nArray];
	double L[nArray];
	double T[nArray];
	
	double L_Blue[nArray];
	double T_Blue[nArray];
	
	double mode[nArray];
	double period[nArray];
	double growth[nArray];
	
	/////////////////////////////////////////
	
	//Reading input parameters from file into arrays//
	
	inFile.open(inputFileName,ios::in);
	
	// Warning if file cant be opened
	if(!inFile.is_open()){ 
		cout << "Error opening file. \n";
		//cout << "Giving Retry... \n";
	}
	if(inFile.is_open()){
		cout<<"Input File was opened successfully"<<endl;
	}
	if(inFile.good()){
		cout<<"Input File is ready for reading"<<endl;
	}
	
	
	//reading file
	if(inFile.is_open()){
	
		// Putting cursor at start of file
		inFile.clear();
		
		//Reading first line - This is needed outside loop otherwise array won't fill.
		inFile >> Mod[0] >> z[0] >> x[0] >> M[0] >> L[0] >> T[0];
		cout<<"Model"<<setw(10)<<"z"<<setw(15)<<"x"<<setw(15)<<"M"<<setw(15)<<"L"<<setw(15)<<"T"<<endl;
		cout<<Mod[0]<<setw(15)<<z[0]<<setw(15)<<x[0]<<setw(15)<<M[0]<<setw(15)<<L[0]<<setw(15)<<T[0]<<endl;
		for (int a = 1; a < NumOfModels; a++){
			
			// Reading rest of file
			inFile >> Mod[a] >> z[a] >> x[a] >> M[a] >> L[a] >> T[a];
			//cout<<Mod[a]<<setw(15)<<z[a]<<setw(15)<<x[a]<<setw(15)<<M[a]<<setw(15)<<L[a]<<setw(15)<<T[a]<<endl;
				
		}
		
		//To show last and first index only, have this line uncommmented and the cout line in loop commented
		cout<<Mod[NumOfModels-1]<<setw(15)<<z[NumOfModels-1]<<setw(15)<<x[NumOfModels-1]<<setw(15)<<M[NumOfModels-1]<<setw(15)<<L[NumOfModels-1]<<setw(15)<<T[NumOfModels-1]<<endl;
		
		// Close the file.
		inFile.close();
	}
			
			
	if(!inFile.is_open()){
		cout<<"Input File closed successfully"<<endl;
	}
	
	///////////////////////////////////////////////////////////////
	
	//BLUE EDGE SELECTOR//
	
	if(GetFirstNegativeGrowthRateModels){
		
		////////////////////////////////////////
		
		cout<<"Creating output file"<<endl;
		ofstream outFile(outputFileName_ForBlueEdge);
		outFile<<setw(5)<<"Model"<<setw(20)<<"z"<<setw(20)<<"x"<<setw(20)<<"M"<<setw(20)<<"L"<<setw(20)<<"T"<<endl; //Header
		//outFile<<""<<endl;
		
		////////////////////////////////////
		
		//Variables to keep track of stuff//
		
		int LINAi = 0;
		int BlueEdgei = 0;
		int NumOfFailedModels = 0;
		int NumOfBlueEdgePoints = 0;
		
		double LowestTemp_Blue = 10000;
		double HighestTemp_Blue = 0;
		double LowestLum_Blue = 10000;
		double HighestLum_Blue = 0;
		
		//Blue Edge Loop//
		for(int i = 0; i < NumOfModels; i++){
			
			/////Reading LINA file////////
			
			//Getting model number into string
			stringstream stream;
				stream<<Mod[i];
				string Model = stream.str();
				//cout<<"Current Model in string form: "<<Model<<endl;
			
			string LINA_header; //Varible for header of LINA file
			string readout; //String thats viewing each line
			
			char LINAFilePath[50]; //log_directory_prefix + Model + "/" + LINAFileName;
			strcpy(LINAFilePath, log_directory_prefix);
			strcat(LINAFilePath, Model.c_str());
			strcat(LINAFilePath, "/");
			strcat(LINAFilePath, LINAFileName);
			cout<<"Path to LINA file: "<<LINAFilePath<<endl;
			
			inFile.open(LINAFilePath,ios::in);
		
			// Warning if file cant be opened
			if(!inFile.is_open()){ 
				cout << "Error opening LINA file. Model did not run. \n";
				NumOfFailedModels = NumOfFailedModels + 1;
			}
			if(inFile.is_open()){
				cout<<"LINA File was opened successfully"<<endl;
			}
			
			
			if(inFile.good()){
				cout<<"LINA File is being read..."<<endl;
			
				//loop that reads LINA file into temp arrays, saving L and T of input file if Blue Edge 
				getline(inFile,LINA_header);
				cout<<LINA_header<<endl;
				inFile>>mode[LINAi]>>period[LINAi]>>growth[LINAi];
				cout<<setw(3)<<mode[LINAi]<<setw(17)<<period[LINAi]<<setw(17)<<growth[LINAi]<<endl;
				
				
				// for(int l = 1; l < NumOfModes; l++){
					
					// inFile>>mode[i]>>period[i]>>growth[i];
					
					// cout<<setw(3)<<mode[i]<<setw(17)<<period[i]<<setw(17)<<growth[i]<<endl;
				
					
				// }
				
				

				//if f0 is negative and last was positive, then output model parameters and assign T and L to blue edge arrays
				if(LINAi != 0){
					
					if(growth[LINAi] < 0. && growth[LINAi-1] > 0.){
						
						outFile<<setw(5)<<Mod[i]<<setw(20)<<z[i]<<setw(20)<<x[i]<<setw(20)<<M[i]<<setw(20)<<L[i]<<setw(20)<<T[i]<<endl;
						
						L_Blue[BlueEdgei] = L[i];
						T_Blue[BlueEdgei] = T[i];
						
						//Selectors to keep track of graphics limits
						if(L_Blue[BlueEdgei] > HighestLum_Blue){
							HighestLum_Blue = L_Blue[BlueEdgei];
						}
						if(L_Blue[BlueEdgei] < LowestLum_Blue){
							LowestLum_Blue = L_Blue[BlueEdgei];
						}
						
						if(T_Blue[BlueEdgei] > HighestTemp_Blue){
							HighestTemp_Blue = T_Blue[BlueEdgei];
						}
						if(T_Blue[BlueEdgei] < LowestTemp_Blue){
							LowestTemp_Blue = T_Blue[BlueEdgei];
						}
						
						BlueEdgei = BlueEdgei + 1;
						NumOfBlueEdgePoints = NumOfBlueEdgePoints + 1;
					}
						
				}

				
				LINAi = LINAi + 1; //Incrementing LINAi
				
				//Space loop - puts space between each model
				// for(int s = 0; s < NumOfSpaces; s++){
					// outFile<<""<<endl;
				// }
				
				inFile.close();
				
				if(!inFile.is_open()){
					cout<<"LINA File closed successfully"<<endl;
				}
		
			}
			
		}
		
	}
	
	///////////////////////////////////////
	
	//POSITVE GROWTH SELECTOR//	
	
	if(GetPositiveGrowthRateModels){
		
		////////////////////////////////////////
		
		cout<<"Creating output file"<<endl;
		ofstream outFile(outputFileName_ForPositiveGrowthModels);
		outFile<<setw(5)<<"Model"<<setw(20)<<"z"<<setw(20)<<"x"<<setw(20)<<"M"<<setw(20)<<"L"<<setw(20)<<"T"<<endl; //Header
		//outFile<<""<<endl;
		
		////////////////////////////////////
		
		//Variables to keep track of stuff
		
		int LINAi = 0;
		int PosGrowthi = 0;
		int NumOfFailedModels = 0;
		int NumOfPosGrowthPoints = 0;
		
		double LowestTemp_Pos = 10000;
		double HighestTemp_Pos = 0;
		double LowestLum_Pos = 10000;
		double HighestLum_Pos = 0;
		
		//Positive Growth Loop//
		for(int i = 0; i < NumOfModels; i++){
			
			/////Reading LINA file////////
			
			//Getting model number into string
			stringstream stream;
				stream<<Mod[i];
				string Model = stream.str();
				//cout<<"Current Model in string form: "<<Model<<endl;
			
			string LINA_header; //Varible for header of LINA file
			string readout; //String thats viewing each line
			
			char LINAFilePath[50]; //log_directory_prefix + Model + "/" + LINAFileName;
			strcpy(LINAFilePath, log_directory_prefix);
			strcat(LINAFilePath, Model.c_str());
			strcat(LINAFilePath, "/");
			strcat(LINAFilePath, LINAFileName);
			cout<<"Path to LINA file: "<<LINAFilePath<<endl;
			
			inFile.open(LINAFilePath,ios::in);
		
			// Warning if file cant be opened
			if(!inFile.is_open()){ 
				cout << "Error opening LINA file. Model did not run. \n";
				NumOfFailedModels = NumOfFailedModels + 1;
			}
			if(inFile.is_open()){
				cout<<"LINA File was opened successfully"<<endl;
			}
			
			
			if(inFile.good()){
				cout<<"LINA File is being read..."<<endl;
			
				//loop that reads LINA file into temp arrays, saving L and T of input file if Blue Edge 
				getline(inFile,LINA_header);
				cout<<LINA_header<<endl;
				inFile>>mode[LINAi]>>period[LINAi]>>growth[LINAi];
				cout<<setw(3)<<mode[LINAi]<<setw(17)<<period[LINAi]<<setw(17)<<growth[LINAi]<<endl;
				
				
				// for(int l = 1; l < NumOfModes; l++){
					
					// inFile>>mode[i]>>period[i]>>growth[i];
					
					// cout<<setw(3)<<mode[i]<<setw(17)<<period[i]<<setw(17)<<growth[i]<<endl;
				
					
				// }
				
				

				//if f0 is positive, add to file
				if(LINAi != 0){
					
					if(growth[LINAi] > 0){
						
						outFile<<setw(5)<<Mod[i]<<setw(20)<<z[i]<<setw(20)<<x[i]<<setw(20)<<M[i]<<setw(20)<<L[i]<<setw(20)<<T[i]<<endl;
						
						L_Blue[PosGrowthi] = L[i];
						T_Blue[PosGrowthi] = T[i];
						
						//Selectors to keep track of graphics limits
						if(L_Blue[PosGrowthi] > HighestLum_Pos){
							HighestLum_Pos = L_Blue[PosGrowthi];
						}
						if(L_Blue[PosGrowthi] < LowestLum_Pos){
							LowestLum_Pos = L_Blue[PosGrowthi];
						}
						
						if(T_Blue[PosGrowthi] > HighestTemp_Pos){
							HighestTemp_Pos = T_Blue[PosGrowthi];
						}
						if(T_Blue[PosGrowthi] < LowestTemp_Pos){
							LowestTemp_Pos = T_Blue[PosGrowthi];
						}
						
						PosGrowthi = PosGrowthi + 1;
						NumOfPosGrowthPoints = NumOfPosGrowthPoints + 1;
					}
						
				}

				
				LINAi = LINAi + 1; //Incrementing LINAi
				
				//Space loop - puts space between each model
				// for(int s = 0; s < NumOfSpaces; s++){
					// outFile<<""<<endl;
				// }
				
				inFile.close();
				
				if(!inFile.is_open()){
					cout<<"LINA File closed successfully"<<endl;
				}
		
			}
			
		}
		
	}
	
	///////////////////////////////////////////////
	
	//Selectors for total number printouts//
	
	if(GetFirstNegativeGrowthRateModels){
		cout<<"NumOfBlueEdgePoints: "<<NumOfBlueEdgePoints<<endl;
	}
	
	
	if(GetPositiveGrowthRateModels){
		cout<<"NumOfPosGrowthPoints: "<<NumOfPosGrowthPoints<<endl;
	}
	
	
	////////////////////////////////////
	
	// cout<<"Creating output file"<<endl;
	// ofstream outFile(outputFileName);
	// outFile<<"Model"<<setw(10)<<"z"<<setw(15)<<"x"<<setw(15)<<"M"<<setw(15)<<"L"<<setw(15)<<"T"<<endl; //Header
	// outFile<<""<<endl;
	
	////////////////////////////////////
	
	///Test loop for LINA data 
	// int LINAi_max = LINAi;
	// for(int l = 0; l < LINAi_max; l++){
		// cout<<setw(3)<<mode[l]<<setw(17)<<period[l]<<setw(17)<<growth[l]<<endl;
		
		// if(l != 0){
			// if(growth[l-1] > 0 && growth[l] < 0){
				// outFile<<Mod[i]<<setw(15)<<z[i]<<setw(15)<<x[i]<<setw(15)<<M[i]<<setw(15)<<L[i]<<setw(15)<<T[i]<<endl;
				
				// L_Blue[i] = L[i];
				// T_Blue[i] = T[i];
			// }
		// }
	// }
	
	////////////////////////////////////////
	
	
	return 0;
}