DRTypeFinder_V10.cpp

// Program to read from file into arrays, and to output values for D/R type analysis
// V1: Copy of program ReadingFromFileIntoArrays_V4, fitted to do stuff with the arrays, output results
	// Note: Temp and temp zone findings not working right? Keeps getting either 0 or last num.
// V2: Goal is to fix problems with V1 and make working program. -Still not fixed, arrays are being filled with last number in file.
// V3: Reworking reading from files to reading from one file-> this works
	// Note: It seems array is being completely filled from 0 to nArray with each entry from file each iteration.
// V4: Fixed the arrays! Seems array needed at least one index filled before being filled.
// V5: Program is now working.
	// Note: Truncated files will result in truncated results!! 
	// Note: sqrt function giving different numbers here then in excel? -No, error in Co and Vo indices.
// V6: Optimized version of V5, allowing multiple files do be done at once with the same relative zones.
	// Note: Sometimes program will fail to operate if the first file is not found. Restart program.
// V7: Added loops to determine D/R type and outputted to file.
	// Note: bools here are outputting as 1 or 0 for true or false repectively
// V8: Added neutral and ionized temp used to output file that was found by program followed by user inputted estimates. Added condition for inputted temps being lower then Temp[0]. Added control for percision. Added Ti_zone to output.
	// Note: Adjustments made to setw's in order for file to open in excel nicely, this should be monitored if percision is changed
	// Note: Files will be overwritten if exsiting file is asked to be made. If file is open in another program, this program will not detect this and no changes will be made to the file or program output. (ie you won't know until you check)
// V9: Added check for user input from "Is this correct" question to avoid accidental infinite loop from input not being 1 or 0.
// V10: Made it so file can be opened and multiple To's can be found and used before closing. Input for To's all at the begining. 
	// Note: This is not completely tested, different values for NumOfTi other than 1 or other senarios could break it.
	

# include <iostream>
# include <fstream>
# include <math.h>
# include <iomanip>
# include <cmath>
# include <stdlib.h>
# include <cstdlib>

using namespace std;

int main(){
	
	int version = 8;
	
	cout<<"This program takes files with columns | Mdis | rho | pressure | vel_km_per_s | temperature | that are placed in its directory."<<endl;

	//////// Controls ////////
	
	const int NumOfZones = 300;
	const int nArray = NumOfZones + 1;
	int NumOfFiles = 1;
	int NumOfTo = 1;
	int NumOfTi = 1;
	int percision = 17; // The number of charcters outputted to the file. Default = 17
	
	cout<<"Enter the number of files to read: "<<endl;
	cin>>NumOfFiles;
	
	cout<<"Enter the number of temp values for the neutral side (ie To1 = 6000K, To2 = 8000K -> 2 values)"<<endl;
	cin>>NumOfTo;
	
	cout<<"Enter the number of temp values for the ionized side (Usually 1 at Ti = 10000K)"<<endl;
	cin>>NumOfTi;
	
	const int NumOfToArray = NumOfTo + 1;
	const int NumOfTiArray = NumOfTi + 1;
	
	// Note: Program picks first number before specified values for temp
	double To[NumOfToArray];
	double Ti[NumOfTiArray];
	for(int n = 0; n < NumOfTo; n++){
		cout<<"Enter temp "<<n+1<<" for the neutral side:"<<endl;
		cin>>To[n];
	}
	
	for(int i = 0; i < NumOfTi; i++){
		cout<<"Enter temp "<<i+1<<" for the ionized side:"<<endl;
		cin>>Ti[i];
	}
	
	///////////////////////////////
	
	ifstream inFile;
	double Mdis[nArray];
	double rho[nArray];
	double pressure[nArray];
	double vel_km_per_s[nArray];
	double Temp[nArray];
	
	for(int f = 0; f < NumOfFiles; f++){
		
		Space(2);
		
		////////// Reading Files into Arrays ////////////
			
			///// Finding File //////
			int TryAgain = 1;
			while(TryAgain == 1){
				char InputFileName[30];
				cout<<"Enter Input File "<<(f+1)<<"'s Name (Max Char=30): "<<endl;
				cin>>InputFileName;
				
				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<<"File was opened successfully"<<endl;
					TryAgain = 0;
				}
				if(inFile.good()){
					cout<<"File is ready for reading"<<endl;
				}
			}
				
			//reading file
			if(inFile.is_open()){
			
				// Putting cursor at start of file
				inFile.clear();
				
				
				inFile >> Mdis[0] >> rho[0] >> pressure[0] >> vel_km_per_s[0] >> Temp[0];
				cout<<"zone"<<setw(10)<<"Mdis"<<setw(15)<<"rho"<<setw(15)<<"pressure"<<setw(15)<<"vel_km_per_s"<<setw(15)<<"Temp"<<endl;
				cout<<0<<setw(10)<<Mdis[0]<<setw(15)<<rho[0]<<setw(15)<<pressure[0]<<setw(15)<<vel_km_per_s[0]<<setw(15)<<Temp[0]<<endl;
				for (int a = 1; a < NumOfZones; a++){
					
					inFile >> Mdis[a] >> rho[a] >> pressure[a] >> vel_km_per_s[a] >> Temp[a];
					cout<<(a+1)<<setw(10)<<Mdis[a]<<setw(15)<<rho[a]<<setw(15)<<pressure[a]<<setw(15)<<vel_km_per_s[a]<<setw(15)<<Temp[a]<<endl;
						
				}
					
				// Close the file.
				inFile.close();
			}
			
			
			if(!inFile.is_open()){
				cout<<"File closed successfully"<<endl;
			}
			
			//test loop 
			// cout<<"Lookup table test..."<<endl;
			// for(int a = 0; a < NumOfZones; a++){
				// cout<<(a+1)<<setw(10)<<Mdis[a]<<setw(15)<<rho[a]<<setw(15)<<pressure[a]<<setw(15)<<vel_km_per_s[a]<<setw(15)<<Temp[a]<<endl;
			// }
			
			Space(2);
		
		///////////////////////////////////////////////////////////////////////
		
		//// Making Output arrays ////
		
		double To_found[NumOfToArray];
		int To_zone[NumOfToArray];
		double Ti_found[NumOfTiArray];
		int Ti_zone[NumOfTiArray];
		
		// Temp Finder Loop 1 - Neutral side
		for(int n = 0; n < NumOfTo; n++){
			for(int a = 0; a < nArray && Temp[a] < To[n]; a++){
			
				To_found[n] = Temp[a];
				To_zone[n] = a + 1;
				//cout<<"To_found: "<<To_found<<endl;
				//cout<<"To_zone: "<<To_zone<<endl;
				
			}
			if(Temp[0] > To[n]){
				To_found[n] = Temp[0];
				To_zone[n] = 1;
			}
			cout<<"Final To_found for To = "<<To[n]<< ": "<<To_found[n]<<endl;
			cout<<"Final To_zone for To = "<<To[n]<< ": "<<To_zone[n]<<endl;
			Space(1);
		}
		
		// Temp Finder Loop 2 - Ionized side
		for(int i = 0; i < NumOfTi; i++){
			for(int a = 0; a < nArray && Temp[a] < Ti[i]; a++){
				
				Ti_found[i] = Temp[a];
				Ti_zone[i] = a + 1;
				//cout<<"Ti_found: "<<Ti_found<<endl;
				//cout<<"Ti_zone: "<<Ti_zone<<endl;
				
			}
			if(Temp[0] > Ti){
				Ti_found[i] = Temp[0];
				Ti_zone[i] = 1;
			}
			cout<<"Final Ti_found for Ti = "<<Ti[i]<< ": "<<Ti_found[i]<<endl;
			cout<<"Final Ti_zone for Ti = "<<Ti[i]<< ": "<<Ti_zone[i]<<endl;
			Space(1);
		}	
			
			Space(1);
		
		for(int i = 0; i < NumOfTi; i++){
			for(int n = 0; n < NumOfTo; n++){
				int Response1 = -1;
				while(Response1 != 1){
					
					if(Response1 == 0){
						Response1 = -1;
					}
					
					cout<<"Enter desired output file "<<f+1<<"'s name for To = "<<To[n]<<", Ti = "<<Ti[i]<<", input file = "<<InputFileName<<" (Max Char=30): "<<endl;
					char OutputFileName[30];
					cin>>OutputFileName;
					Space(1);
					
					while(Response1 !=1 && Response1 !=0){
						cout<<"Is this correct? Enter 1 for yes and 0 for no."<<endl;
						cout<<"File name: "<<OutputFileName<<endl;
						cin>>Response1;
					}
					Space(1);
				}
				
				ofstream outputfile1(OutputFileName);
				
				if(outputfile1){
					cout<<"File created successfully"<<endl;
				}
				if(outputfile1.good()){
					cout<<OutputFileName<<" is ready for writting"<<endl;
					Space(1);
				}
				

				//const int nArray2 = nArray - Ti_zone;
				
				double Co;
				double Vo;
				double Co[nArray];
				double Ci[nArray];
				double Vo[nArray];
				double Vd[nArray];
				double Vr[nArray];
				
				// double To = To[n];
				// double Ti = Ti[n];
				// int To_zone = To_zone[n];
				// int Ti_zone = Ti_zone[n];
				
				
				////////// Calculations //////////////
				
				double Co = sqrt(pressure[To_zone[n]-1]/rho[To_zone[n]-1]);
				double Vo = vel_km_per_s[To_zone[n]-1];
				
				outputfile1<<setw(2)<<"N_zone"<<setw(20)<<"Co"<<setw(25)<<"Vo"<<setw(25)<<"To"<<endl;
				outputfile1<<setw(2)<<To_zone[n]<<setprecision(percision)<<setw(25)<<Co<<setw(25)<<setw(25)<<Vo<<setw(25)<<To[n]<<setw(25)<<endl;
				outputfile1<<setw(2)<<"I_zone"<<setw(25)<<"To_found"<<setw(25)<<"Ti_found"<<setw(20)<<"Ti"<<endl;
				outputfile1<<setw(2)<<Ti_zone[i]<<setprecision(percision)<<setw(25)<<To_found[n]<<setw(25)<<Ti_found[i]<<setw(25)<<Ti[i]<<endl;
				
				// Calculation loop, starting at Ti_zone that was found
				outputfile1<<setw(2)<<"zone"<<setw(25)<<"Ci"<<setw(25)<<"Vd"<<setw(25)<<"Vr"<<setw(10)<<"Vr>Co"<<setw(10)<<"Vo>=Vr"<<setw(10)<<"Vd<Co"<<setw(10)<<"Vo<=Vd"<<setw(10)<<"Type"<<endl;
				for(int c = Ti_zone[i]-1; c < NumOfZones; c++){
					//Co[c] = Co;
					//Vo[c] = Vo;
					Ci[c] = sqrt(pressure[c]/rho[c]);
					Vd[c] = Ci[c] - sqrt((Ci[c]*Ci[c])-(Co*Co)); 
					Vr[c] = Ci[c] + sqrt((Ci[c]*Ci[c])-(Co*Co)); 
					
					bool Rcond1 = (Vr[c] > Co);
					bool Rcond2 = (Vo >= Vr[c]);
					bool Dcond1 = (Vd[c] < Co);
					bool Dcond2 = (Vo <= Vd[c]);
					
					if(Rcond1 && Rcond2){
						outputfile1<<setw(2)<<(c+1)<<setprecision(percision)<<setw(25)<<Ci[c]<<setw(25)<<Vd[c]<<setw(25)<<Vr[c]<<setw(10)<<Rcond1<<setw(10)<<Rcond2<<setw(10)<<Dcond1<<setw(10)<<Dcond2<<setw(10)<<"R"<<endl;
					}
					if(Dcond1 && Dcond2){
						outputfile1<<setw(2)<<(c+1)<<setprecision(percision)<<setw(25)<<Ci[c]<<setw(25)<<Vd[c]<<setw(25)<<Vr[c]<<setw(10)<<Rcond1<<setw(10)<<Rcond2<<setw(10)<<Dcond1<<setw(10)<<Dcond2<<setw(10)<<"D"<<endl;
					}
					if(!(Rcond1 && Rcond2) && !(Dcond1 && Dcond2)){
						outputfile1<<setw(2)<<(c+1)<<setprecision(percision)<<setw(25)<<Ci[c]<<setw(25)<<Vd[c]<<setw(25)<<Vr[c]<<setw(10)<<Rcond1<<setw(10)<<Rcond2<<setw(10)<<Dcond1<<setw(10)<<Dcond2<<setw(10)<<"Could not be found"<<endl;
					}
					
					
				}
				
				outputfile1.close();
				if(!outputfile1.is_open()){
					cout<<"File was filled and closed successfully"<<endl;
				}
			}
		}
		
	}
	
	return 0;
	
}



void Space(int Size){
	int i = 0;
	while (i < Size){
		cout<<" "<<endl;
		i = i + 1;
	}
}