Output methods More...

#include <Output.h>

Inheritance diagram for Output:

Collaboration diagram for Output:

Public Member Functions
	Output (ThreadManager *MTHREAD_h)
	Constructor. More...

	~Output ()

void	initOutput ()

void	commonInit ()

void	initOutputMaps ()

void	initOutputForestData ()

void	initOutputProductData ()

void	initOptimisationLog ()

void	initDebugOutput ()

void	initDebugPixelValues ()

void	initCarbonBalance ()

void	print (bool earlyPrint=true)
	Print output. If earlyPrinting it doesn't print some stuff for which we don't yet have values. More...

void	printMaps ()

void	printForestData (bool finalFlush)

void	printProductData (bool finalFlush)

void	printCarbonBalance ()

void	printFinalOutput ()

void	printDebugOutput ()

void	printDetailedHV (map< tr1::array< string, 4 >, double > hVol_byPrd)

void	printDebugPixelValues ()

void	printOptLog (bool optimal, int &nIterations, double &obj)

char	getOutputFieldDelimiter ()

void	cleanScenario (string fileName, string scenarioName, char d)

Public Member Functions inherited from BaseClass
	BaseClass ()

	~BaseClass ()

void	msgOut (const int &msgCode_h, const string &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log. More...

void	msgOut (const int &msgCode_h, const int &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log. More...

void	msgOut (const int &msgCode_h, const double &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log. More...

int	s2i (const string &string_h) const
	string to integer conversion More...

double	s2d (const string &string_h) const
	string to double conversion More...

double	s2d (const string &string_h, const bool &replaceComma) const
	string to double conversion More...

bool	s2b (const string &string_h) const
	string to bool conversion More...

string	i2s (const int &int_h) const
	integer to string conversion More...

string	d2s (const double &double_h) const
	double to string conversion More...

string	b2s (const bool &bool_h) const
	bool to string conversion More...

vector< int >	s2i (const vector< string > &string_h) const
	string to integer conversion (vector) More...

vector< double >	s2d (const vector< string > &string_h, const bool &replaceComma=false) const
	string to double conversion (vector) More...

vector< bool >	s2b (const vector< string > &string_h) const
	string to bool conversion (vector) More...

vector< string >	i2s (const vector< int > &int_h) const
	integer to string conversion (vector) More...

vector< string >	d2s (const vector< double > &double_h) const
	double to string conversion (vector) More...

vector< string >	b2s (const vector< bool > &bool_h) const
	bool to string conversion (vector) More...

int	getType (const string &type_h) const
	Return a type according to enum TYPE_* from a string (eg: "string" -> TYPE_STRING (2)) More...

void	refreshGUI () const
	Ping to periodically return the control to the GUI. More...

template<typename T >
string	toString (const T &x) const

template<typename T >
T	stringTo (const std::string &s) const

int	vSum (const vector< int > &vector_h) const

double	vSum (const vector< double > &vector_h) const

int	vSum (const vector< vector< int > > &vector_h) const

double	vSum (const vector< vector< double > > &vector_h) const

void	tokenize (const string &str, vector< string > &tokens, const string &delimiter=" ") const
	Tokenize a string using a delimiter (default is space) More...

void	untokenize (string &str, vector< string > &tokens, const string &delimiter=" ") const

template<typename K , typename V >
V	findMap (const map< K, V > &mymap, const K &key, const int &error_level=MSG_CRITICAL_ERROR, const V &notFoundValue=numeric_limits< V >::min()) const
	Lookup a map for a value. Return the value starting from the key. More...

template<typename K , typename V >
void	changeMapValue (map< K, V > &mymap, const K &key, const V &value, const int &error_level=MSG_CRITICAL_ERROR)
	Change the value stored in a map given the key and the new value. More...

template<typename K , typename V >
void	incrMapValue (map< K, V > &mymap, const K &key, const V &value, const int &error_level=MSG_CRITICAL_ERROR)
	Increments a value stored in a map of the specified value, given the key. More...

template<typename K , typename V >
void	incrOrAddMapValue (map< K, V > &mymap, const K &key, const V &value)
	Increments a value stored in a map of the specified value, given the key. More...

template<typename K , typename V >
void	resetMapValues (map< K, V > &mymap, const V &value)
	Reset all values stored in a map to the specified one. More...

template<typename K , typename V >
map< K, V >	vectorToMap (const vector< K > &keys, const V &value=0.0)
	Returns a map built using the given vector and the given (scalar) value as keys/values pairs. More...

template<typename T >
vector< T >	positionsToContent (const vector< T > &vector_h, const vector< int > &positions)
	Return a vector of content from a vector and a vector of positions (int) More...

template<typename V >
void	debugMap (const map< iisskey, V > &mymap)
	Debug a map. More...

template<typename K , typename V >
void	debugMap (const map< K, V > &mymap, const K &key)

template<typename K >
int	getMaxPos (const vector< K > &v)
	Returns the position of the maximum element in the vector (the last one in case of multiple equivalent maxima) More...

template<typename K >
int	getMinPos (const vector< K > &v)
	Returns the position of the minimum element in the vector (the first one in case of multiple equivalent minima) More...

template<typename K >
K	getMax (const vector< K > &v)
	Returns the value of the maximum element in the vector (the last one in case of multiple equivalent maxima) More...

template<typename K >
K	getMin (const vector< K > &v)
	Returns the value of the minimum element in the vector (the first one in case of multiple equivalent minima) More...

template<typename K >
double	getAvg (const vector< K > &v)
	Returns the average of the elements in the vector. More...

template<typename K >
double	getSd (const vector< K > &v, bool sample=true)

template<typename K >
int	getPos (const K &element, const vector< K > &v, const int &msgCode_h=MSG_CRITICAL_ERROR)

template<typename K >
bool	inVector (const K &element, const vector< K > &v)

double	normSample (const double &avg, const double &stdev, const double &minval=NULL, const double &maxval=NULL) const
	Sample from a normal distribution with bounds. Slower (double time, but still you see the diff only after milion of loops). More...

template<typename K >
K	normSample (normal_distribution< K > &d, std::mt19937 &gen, const K &minval=NULL, const K &maxval=NULL) const
	Sample from a normal distribution with bounds. Faster (half time) as the normal_distribution is made only once. More...

template<typename T >
std::string	toString (const T &x) const

Public Attributes
vector< vector< vector< vector< vector< double > > > > >	expReturnsDebug
	l2_region, for type, d.c., pr prod, variable name More...

vector< string >	expReturnsDebugVariables

Private Attributes
int	oLevel

char	d

int	inYear

int	nYears

string	baseDir

string	oDir

string	scenarioName

string	oFileExt

bool	oHRedeable

bool	oSingleFile

vector< int >	oYears

vector< int >	mapsOYears

int	wRegId_l1

int	wRegId_l2

string	outFileName

vector< string >	outForVariables

vector< string >	outProdVariables

bool	oForData

bool	oProdData

bool	oCarbonData

int	outStepRange

bool	forestDiamDetailedOutput

vector< string >	priPr

vector< string >	secPr

vector< string >	allPr

vector< int >	l1regIds

vector< vector< int > >	l2r

vector< string >	fTypes

vector< string >	dClasses

vector< string >	pDClasses
	includes an empty string for variables without diameter attribute More...

int	nPriPr

int	nSecPr

int	nAllPr

int	nL2r

string	logFilename

string	debugFilename

string	debugPxValuesFilename

bool	spMode

Additional Inherited Members
Protected Attributes inherited from BaseClass
ThreadManager *	MTHREAD
	Pointer to the Thread manager. More...

Detailed Description

Output methods

Class responsable to output the data, both as all kind of log as well as georeferenciated one.

Author: Antonello Lobianco

Definition at line 48 of file Output.h.

Constructor & Destructor Documentation

Output ( ThreadManager * MTHREAD_h )

Constructor.

Definition at line 37 of file Output.cpp.

                                       {
   MTHREAD=MTHREAD_h;
 }

~Output ( )

Definition at line 41 of file Output.cpp.

41 {

42 }

Member Function Documentation

void cleanScenario	(	string	fileName,
		string	scenarioName,
		char	d
	)

This routine clean the output scenario from previous outputs of the defined scenario. Other scenarios are untouched. The scenarioName must be in the first row.

Parameters

filename	Filename of the output file to clean
scenarioName	Name of the scenario we are replacing
d	Field delimiter. It must not be changed in the meantime (between the various scenarios)

Definition at line 1149 of file Output.cpp.

Referenced by initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), and initOutputProductData().

                                                                  {
   string dStr(&d,1);
   vector <string> rows;
   string tempRow;
   ifstream inFile (fileName.c_str(), ios::in);
   if (!inFile){
     msgOut(MSG_ERROR,"Error in opening the file "+fileName+" for reading.");
     return;
   }
   while( getline (inFile,tempRow) ){
     vector<string> tokens;
     tokenize(tempRow,tokens,dStr);
     if(tokens[0] != scenarioName)
       rows.push_back( tempRow );
   }
   inFile.close();
   ofstream out(fileName.c_str(), ios::out);
   for(uint i=0;i<rows.size();i++){
     out << rows[i];
     out << "\n";
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void commonInit ( )

Definition at line 61 of file Output.cpp.

Referenced by initOutput().

                   {
   oLevel       = MTHREAD->MD->getIntSetting("outputLevel");
   d            = getOutputFieldDelimiter();
   inYear       = MTHREAD->MD->getIntSetting("initialYear");
   nYears       = MTHREAD->MD->getIntSetting("simulationYears");
   baseDir      = MTHREAD->MD->getBaseDirectory();
   oDir         = MTHREAD->MD->getOutputDirectory();
 //    bool initSeed = MTHREAD->MD->getBoolSetting("newRandomSeed");
 //    if (initSeed){
 //        uniform_int_distribution<> d(1, 1000000);
 //        int random = d(*MTHREAD->gen);
 //        scenarioName = MTHREAD->getScenarioName()+"_"+i2s(random);
 //    } else {
 //        scenarioName = MTHREAD->getScenarioName();
 //    }
   if (MTHREAD->MD->getStringSetting("overridenScenarioName") == "none"){
     scenarioName = MTHREAD->getScenarioName();
   } else {
     scenarioName = MTHREAD->MD->getStringSetting("overridenScenarioName");
   }
   oFileExt     = MTHREAD->MD->getStringSetting("outputFileExtension");
   oHRedeable   = MTHREAD->MD->getBoolSetting("outputHumanReadable");
   oSingleFile  = MTHREAD->MD->getBoolSetting("outputSingleFile");
   oYears       = MTHREAD->MD->getIntVectorSetting("outYears");
   mapsOYears   = MTHREAD->MD->getIntVectorSetting("mapsOutYears");
   wRegId_l1    = MTHREAD->MD->getIntSetting("worldCodeLev1");
   wRegId_l2    = MTHREAD->MD->getIntSetting("worldCodeLev2");
   outForVariables          = MTHREAD->MD->getStringVectorSetting("outForVariables");
   outProdVariables         = MTHREAD->MD->getStringVectorSetting("outProdVariables");
   dClasses                 = MTHREAD->MD->getStringVectorSetting("dClasses");
   pDClasses.insert(pDClasses.end(), dClasses.begin()+1, dClasses.end() ); // production diameter classes
   dClasses.push_back(""); // needed for reporting of variables without diameter attribute
   outStepRange             = MTHREAD->MD->getIntSetting("outStepRange");
   forestDiamDetailedOutput = MTHREAD->MD->getBoolSetting("forestDiamDetailedOutput");
   fTypes                   = MTHREAD->MD->getForTypeIds();
 
   oForData    = (oLevel>=OUTVL_DETAILED && outForVariables.size() > 0)?true:false;
   oProdData   = (oLevel>=OUTVL_DETAILED && outProdVariables.size() > 0)?true:false;
   oCarbonData = (oLevel>=OUTVL_DETAILED)?true:false;
 
   priPr    = MTHREAD->MD->getStringVectorSetting("priProducts");
   secPr    = MTHREAD->MD->getStringVectorSetting("secProducts");
   allPr    = priPr;
   allPr.insert( allPr.end(), secPr.begin(), secPr.end() );
   nPriPr   = priPr.size();
   nSecPr   = secPr.size();
   nAllPr   = allPr.size();
   l1regIds = MTHREAD->MD->getRegionIds(1, true);
   nL2r     = MTHREAD->MD->getRegionIds(2, true).size();
   spMode   = MTHREAD->MD->getBoolSetting("usePixelData");
     //if(spMode) {
     //   pxIds =   getXyNPixels();
     //}
 
 
   for(uint i=0;i<l1regIds.size();i++){
     std::vector<int> l2ChildrenIds;
     ModelRegion* l1Region = MTHREAD->MD->getRegion(l1regIds[i]);
     std::vector<ModelRegion*> l2Childrens = l1Region->getChildren(true);
     for(uint j=0;j<l2Childrens.size();j++){
       l2ChildrenIds.push_back(l2Childrens[j]->getRegId());
     }
     if(l2ChildrenIds.size()){
       l2r.push_back(l2ChildrenIds);
     }
   }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

char getOutputFieldDelimiter ( )

Definition at line 802 of file Output.cpp.

Referenced by commonInit().

                                {
   int delimiterID = MTHREAD->MD->getIntSetting("outputFieldDelimiter");
   switch (delimiterID) {
     case 1:
       return ',';
       break;
     case 2:
       return ';';
       break;
     case 3:
       return ':';
       break;
     case 4:
       return '\t';
       break;
     case 5:
       return ' ';
       break;
     default:
       msgOut(MSG_ERROR, "You have specified an unknow output file field delimiter. Using \";\".");
       return ',';
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initCarbonBalance ( )

Definition at line 370 of file Output.cpp.

Referenced by initOutput().

                          {
   if(!oCarbonData) return;
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/carbonBalance"+oFileExt;
     ifstream in(outFileName.c_str(), ios::in);
     if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
       in.close();
       cleanScenario(outFileName, scenarioName, d);
       return;
     } else {
       in.close();
     }
   } else {
     outFileName = baseDir+oDir+"results/carbonBalance_"+scenarioName+oFileExt;
   }
 
   ofstream out(outFileName.c_str(), ios::out);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
     out << "scen" << d << "country" << d << "region" << d << "balItem" << d;
   //if(oHRedeable){
   //  for(int i=0;i<nYears;i++){
   //    out << i+inYear << d;
   //  }
   //} else {
     out << "year" << d << "value" << d;
   //}
   out << "\n";
   out.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initDebugOutput ( )

Definition at line 170 of file Output.cpp.

Referenced by initOutput().

                        {
     if(oLevel<OUTVL_ALL) return;
 
     // init debugging the expected returns...
     if(spMode) return;
     expReturnsDebugVariables.push_back("hVol_byUPp");
     expReturnsDebugVariables.push_back("hV_byFT");
     expReturnsDebugVariables.push_back("finalHarvestFlag");
     expReturnsDebugVariables.push_back("pondCoeff");
     expReturnsDebugVariables.push_back("pW");
     expReturnsDebugVariables.push_back("cumTp");
     expReturnsDebugVariables.push_back("vHa");
     expReturnsDebugVariables.push_back("expectedReturns");
     expReturnsDebugVariables.push_back("weightedAvgCompModeFlag");
 
     if (oSingleFile){
          debugFilename = baseDir+oDir+"debugs/debugOut.csv";
     } else {
          debugFilename = baseDir+oDir+"debugs/debugOut_"+scenarioName+".csv";
     }
 
     ifstream in(debugFilename.c_str(), ios::in);
     if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
         in.close();
         cleanScenario(debugFilename, scenarioName, d);
         return;
     } else { // file doesn't exist
         in.close();
         ofstream out(debugFilename.c_str(), ios::out);
         if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+debugFilename+" for writing.");}
         out << "scenario" << d << "year" << d << "region or pixel" << d << "forType" << d << "freeDim" << d << "prod" << d << "parName" << d << "value" << d <<"\n";
         out.close();
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initDebugPixelValues ( )

Definition at line 207 of file Output.cpp.

Referenced by initOutput().

                             {
     if(oLevel<OUTVL_ALL) return;
 
     // init debugging the expected returns...
     if(!spMode) return;
 
     if (oSingleFile){
          debugPxValuesFilename = baseDir+oDir+"debugs/debugPxValues.csv";
     } else {
          debugPxValuesFilename = baseDir+oDir+"debugs/debugPxValues_"+scenarioName+".csv";
     }
 
     ifstream in(debugPxValuesFilename.c_str(), ios::in);
     if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
         in.close();
         cleanScenario(debugPxValuesFilename, scenarioName, d);
         return;
     } else { // file doesn't exist
         in.close();
         ofstream out(debugPxValuesFilename.c_str(), ios::out);
         if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+debugPxValuesFilename+" for writing.");}
         out << "scenario" << d << "year" << d << "region" << d << "pxId" << d << "pxX" << d << "pxY" << d ;
         /*for(uint f=0;f<fTypes.size();f++){
           string ft = fTypes[f];
           string header = "tp_multiplier_"+ft;
           out << header <<d;
         }
         for(uint f=0;f<fTypes.size();f++){
           string ft = fTypes[f];
           string header = "mortCoef_multiplier_"+ft;
           out << header <<d;
         }*/
         out << "var" << d ;
 
         for(uint f=0;f<fTypes.size();f++){
           string ft = fTypes[f];
           for (uint u=0;u<dClasses.size();u++){
             string dc=dClasses[u];
             string header = ft+"_"+dc;
             out << header <<d;
           }
         }
         for(uint f=0;f<fTypes.size()+1;f++){
           string ft1 = (f == fTypes.size()) ? "NonFor" : fTypes[f];
           for(uint f2=0;f2<fTypes.size()+1;f2++){
             string ft2 = (f2 == fTypes.size()) ? "NonFor" : fTypes[f2];
             string header = ft1+"_2_"+ft2;
             out << header <<d;
           }
         }
         out << "\n";
 
 
         out.close();
     }
 
 
 
 
     /*
     if(oSingleFile){
       outFileName = baseDir+oDir+"results/forestData"+oFileExt;
       ifstream in(outFileName.c_str(), ios::in);
       if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
         in.close();
         cleanScenario(outFileName, scenarioName, d);
         return;
       } else {
         in.close();
       }
     } else {
       outFileName = baseDir+oDir+"results/forestData_"+scenarioName+oFileExt;
     }
 
     ofstream out(outFileName.c_str(), ios::out);
     if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
       out << "scen" << d << "parName" << d << "country" << d << "region" << d << "forType" << d << "freeDim" << d;
     */
 
 
 
 
 
 
 
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initOptimisationLog ( )

Definition at line 131 of file Output.cpp.

Referenced by initOutput().

                            {
   if(oLevel<OUTVL_AGGREGATED) return;
 
     if (oSingleFile){
          logFilename = baseDir+oDir+"optimisationLogs/optimisationLogs.txt";
 
     } else {
          logFilename = baseDir+oDir+"optimisationLogs/"+scenarioName+".txt";
     }
 
 
   ifstream in(logFilename.c_str(), ios::in);
   if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
     in.close();
     cleanScenario(logFilename, scenarioName, d);
     ofstream out(logFilename.c_str(), ios::app);
     if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+logFilename+" for writing.");}
     time_t now;
     time(&now);
     struct tm *current = localtime(&now);
     string timemessage = i2s(current->tm_hour)+":"+i2s(current->tm_min)+":"+ i2s(current->tm_sec);
     out << scenarioName << d << "0000" << d << timemessage << d << d << d <<"\n";
     out.close();
     return;
   } else { // file doesn't exist
     in.close();
     ofstream out(logFilename.c_str(), ios::out);
     if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+logFilename+" for writing.");}
     out << "scenario" << d << "year" << d << "time" << d << "opt flag" << d << "iterations" << d <<"\n";
     time_t now;
     time(&now);
     struct tm *current = localtime(&now);
     string timemessage = i2s(current->tm_hour)+":"+i2s(current->tm_min)+":"+ i2s(current->tm_sec);
     out << scenarioName << d << "0000" << d << timemessage << d << d << d <<"\n";
     out.close();
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initOutput ( )

Definition at line 48 of file Output.cpp.

Referenced by Init::setInitLevel3().

                   {
   commonInit();
   initOutputMaps();
   initDebugOutput();
   initDebugPixelValues();
   initOutputForestData();
   initOutputProductData();
   initOptimisationLog();
   initCarbonBalance();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initOutputForestData ( )

Definition at line 296 of file Output.cpp.

Referenced by initOutput().

                             {
   if(!oForData) return;
 
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/forestData"+oFileExt;
     ifstream in(outFileName.c_str(), ios::in);
     if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
       in.close();
       cleanScenario(outFileName, scenarioName, d);
       return;
     } else {
       in.close();
     }
   } else {
     outFileName = baseDir+oDir+"results/forestData_"+scenarioName+oFileExt;
   }
 
   ofstream out(outFileName.c_str(), ios::out);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
     out << "scen" << d << "parName" << d << "country" << d << "region" << d << "forType" << d << "freeDim" << d;
   if(oHRedeable){
     for(int i=0;i<nYears;i++){
       out << i+inYear << d;
     }
   } else {
     out << "year" << d << "value" << d;
   }
   out << "\n";
   out.close();
 
   // gonna do the same for detailed hV..
   if(!MTHREAD->MD->getBoolSetting("outDetailedHv",DATA_NOW)) return;
   outFileName = baseDir+oDir+"results/detailedHV_"+scenarioName+oFileExt;
   ofstream outHV(outFileName.c_str(), ios::out);
   if (!outHV){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
   outHV << "scen" << d << "parName" << d << "country" << d << "region" << d << "forType" << d << "freeDim" << d << "prod" << d << "year" << d << "value" << d;
   outHV << "\n";
   outHV.close();
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initOutputMaps ( )

Resetting the list of printed layers and the scenario name..
Printing scenario name for post-processing scripts

Definition at line 406 of file Output.cpp.

Referenced by initOutput().

                       {
   if(oLevel<OUTVL_MAPS) return;
   string mapBaseDirectory = baseDir+oDir+"maps/";
     string filenameToSaveScenarioName = mapBaseDirectory+"scenarioNames/"+scenarioName;
     string filenameListIntLayers = mapBaseDirectory+"integerListLayers/"+scenarioName;
     string filenameListFloatLayers = mapBaseDirectory+"floatListLayers/"+scenarioName;
 
   // printing the scenario name in the "scenarioName file"...
   ofstream outSN(filenameToSaveScenarioName.c_str(), ios::out);
   if (!outSN){ msgOut(MSG_ERROR,"Error in opening the file "+filenameToSaveScenarioName+".");}
   outSN << scenarioName << "\n";
   outSN.close();
   // cleaning the "integerListLayers" and "floatListLayers" file...
   ofstream outi(filenameListIntLayers.c_str(), ios::out);
   outi.close();
   ofstream outf(filenameListFloatLayers.c_str(), ios::out);
   outf.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initOutputProductData ( )

Definition at line 338 of file Output.cpp.

Referenced by initOutput().

                              {
   if(!oProdData) return;
 
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/productData"+oFileExt;
     ifstream in(outFileName.c_str(), ios::in);
     if(in.is_open()) { // file exist, no need to initializate it, but we are gonna clean it of previous data of the same scenario if present...
       in.close();
       cleanScenario(outFileName, scenarioName, d);
       return;
     } else {
       in.close();
     }
   } else {
     outFileName = baseDir+oDir+"results/productData_"+scenarioName+oFileExt;
   }
 
   ofstream out(outFileName.c_str(), ios::out);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
   out << "scen" << d << "parName" << d << "country" << d << "region" << d << "prod" << d << "freeDim" << d;
   if(oHRedeable){
     for(int i=0;i<nYears;i++){
       out << i+inYear << d;
     }
   } else {
     out << "year" << d << "value" << d;
   }
   out << "\n";
   out.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void print ( bool earlyPrint = true )

Print output. If earlyPrinting it doesn't print some stuff for which we don't yet have values.

Definition at line 426 of file Output.cpp.

Referenced by ModelCore::runInitPeriod(), ModelCoreSpatial::runInitPeriod(), ModelCore::runSimulationYear(), and ModelCoreSpatial::runSimulationYear().

                             {
   int cYear = MTHREAD->SCD->getYear();
   int initialSimulationYear = MTHREAD->MD->getIntSetting("initialOptYear");
 
   if (outStepRange != -1 && (cYear-initialSimulationYear)%outStepRange != 0 && cYear != (initialSimulationYear+nYears)-1 ) {
     cout << cYear << " not printed" << endl;
     return;
   }
   bool printThisYear = false;
   for(uint i=0;i<oYears.size();i++){
     if (outStepRange == -1 && oYears[i] == cYear) printThisYear = true;
   }
   if(outStepRange == -1 && !printThisYear) return;
 
   printMaps();
   MTHREAD->MD->setErrorLevel(MSG_NO_MSG);
   printForestData(false);
   printProductData(false);
   printCarbonBalance();
   printDebugOutput();
   if(!earlyPrint) printDebugPixelValues();
   MTHREAD->MD->setErrorLevel(MSG_ERROR);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printCarbonBalance ( )

Definition at line 727 of file Output.cpp.

Referenced by print().

                           {
 
   if(!oCarbonData) return;
   int currentYear = MTHREAD->SCD->getYear();
   if (currentYear == inYear) {return;} // don't print carbon balance on first year, carbon balance containers has not yet been initialised
 
   msgOut(MSG_INFO, "Printing forest data..");
 
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/carbonBalance"+oFileExt;
   } else {
     outFileName = baseDir+oDir+"results/carbonBalance_"+scenarioName+oFileExt;
   }
   ofstream out (outFileName.c_str(), ios::app);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for writing.");}
   double outvalue=0;
 
   vector<int> balItems {STOCK_INV,STOCK_EXTRA,STOCK_PRODUCTS,EM_ENSUB,EM_MATSUB,EM_FOROP};
 
   for (uint r1=0;r1<l2r.size();r1++){
     for (uint r2=0;r2<l2r[r1].size();r2++){
       int regId = l2r[r1][r2];
       for (uint b=0;b<balItems.size();b++){
         out << scenarioName << d;
         out << MTHREAD->MD->regId2RegSName(l1regIds.at(r1)) << d;
         out << MTHREAD->MD->regId2RegSName(l2r[r1][r2]) << d;
         string balItemString;
         switch(balItems[b]){
           case STOCK_INV: {
             balItemString = "STOCK_INV";
             outvalue = MTHREAD->CBAL->getStock(regId, balItems[b]);
             break;
           }
           case STOCK_EXTRA: {
             balItemString = "STOCK_EXTRA";
             outvalue = MTHREAD->CBAL->getStock(regId, balItems[b]);
             break;
           }
           case STOCK_PRODUCTS: {
             balItemString = "STOCK_PRODUCTS";
             outvalue = MTHREAD->CBAL->getStock(regId, balItems[b]);
             break;
           }
           case EM_ENSUB: {
             balItemString = "EM_ENSUB";
             outvalue = MTHREAD->CBAL->getCumSavedEmissions(regId, balItems[b]);
             break;
           }
           case EM_MATSUB: {
             balItemString = "EM_MATSUB";
             outvalue = MTHREAD->CBAL->getCumSavedEmissions(regId, balItems[b]);
             break;
           }
           case EM_FOROP: {
             balItemString = "EM_FOROP";
             outvalue = MTHREAD->CBAL->getCumSavedEmissions(regId, balItems[b]);
             break;
           }
         default:
             msgOut(MSG_CRITICAL_ERROR,"Unexpected balance item type in function printCarbonBalance");
         }
         out << balItemString << d;
         out << currentYear << d;
         out << outvalue << d;
         out << "\n";
 
 
       } // end bal items
     } // end r2
   } // end r1
   out.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printDebugOutput ( )

Definition at line 843 of file Output.cpp.

Referenced by print().

                         {
   if(oLevel<OUTVL_ALL) return;
 
   // print debugging the expected returns...
 
   if (!spMode && !expReturnsDebug.empty()){
     ofstream out (debugFilename.c_str(), ios::app);
     if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+debugFilename+" for writing.");}
     int currentYear = MTHREAD->SCD->getYear();
     vector <int> regIds2 = MTHREAD->MD->getRegionIds(2);
 
         for (uint r2=0;r2<regIds2.size();r2++){
             for(uint ft=0;ft<fTypes.size();ft++){
                 for(uint dc=0;dc<(dClasses.size()-1);dc++){
                     for(uint pp=0;pp<priPr.size();pp++){
                         for(uint dv=0;dv<expReturnsDebugVariables.size();dv++){
                           // vector <vector < vector <vector <vector <double> > > > > expReturnsDebug;
                           double outputValue = expReturnsDebug.at(r2).at(ft).at(dc).at(pp).at(dv);
                           out << scenarioName << d;
                           out << currentYear << d;
                           out << MTHREAD->MD->regId2RegSName(regIds2[r2]) << d;
                           out << fTypes[ft] << d;
                           out << dClasses[dc] << d;
                           out << priPr[pp] << d;
                           out << expReturnsDebugVariables[dv] << d;
                           out << outputValue << d;
                           out << "\n";
                         }
                     }
                 }
             }
         }
 
   } // end initial condition checks
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printDebugPixelValues ( )

Definition at line 880 of file Output.cpp.

Referenced by print().

                              {
 
   if(oLevel<OUTVL_ALL) return;
 
   bool filter;
   filter = true; //use this to  filter output
   if(filter && spMode){
     ofstream out (debugPxValuesFilename.c_str(), ios::app);
     if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+debugPxValuesFilename+" for writing.");}
     int currentYear = MTHREAD->SCD->getYear();
     vector <int> regIds2 = MTHREAD->MD->getRegionIds(2);
     for (uint r=0;r<regIds2.size();r++){
       int rId = regIds2[r];
       //if(rId != 11002) continue;
       ModelRegion* REG = MTHREAD->MD->getRegion(rId);
       vector<Pixel*> regPx = REG->getMyPixels();
       for (uint p=0;p<regPx.size();p++){
         Pixel* px = regPx[p];
         int pxID = px->getID();
         int pxX = px->getX();
         int pxY = px->getY();
         string common = scenarioName + d + i2s(currentYear) + d + i2s(rId) + d+ i2s(pxID) +d +i2s(pxX)+d+i2s(pxY)+d;
 
 
 //        for(uint f=0;f<fTypes.size();f++){
 //          double tp_m = px->getMultiplier("tp_multiplier",fTypes[f]);
 //          common += d2s(tp_m)+d;
 //        }
 //        for(uint f=0;f<fTypes.size();f++){
 //          double m_m = px->getMultiplier("mortCoef_multiplier",fTypes[f]);
 //          common += d2s(m_m)+d;
 //        }
 
 
 //        // First vars by only ft...
 //        // expectedReturns
 //        out << common << "expectedReturns" << d;
 //        for(uint f=0;f<fTypes.size();f++){
 //          for(uint u=0;u<dClasses.size()-1;u++){
 //            out << d;
 //          }
 //          out << px->expectedReturns[f] << d;
 //          //out << 0.0 << d;
 //        }
 //        out << "\n";
 
         //----
 
         vector <string> outVars   = MTHREAD->MD->getStringVectorSetting("debugVariablesToPrint");
 
         if (std::find(outVars.begin(), outVars.end(), "vol") != outVars.end()){
             out << common <<"vol" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->vol[f][u]<< d;
               }
               out << vSum(px->vol[f]) << d;
             }
             out << "\n";
        }
 
         if (std::find(outVars.begin(), outVars.end(), "area") != outVars.end()){
             out << common <<"area" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->area[f][u]<< d;
               }
               out << vSum(px->area[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "hArea") != outVars.end()){
             out << common <<"hArea" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->hArea[f][u]<< d;
               }
               out << vSum(px->hArea[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "hVol") != outVars.end()){
             out << common <<"hVol" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->hVol[f][u]<< d;
               }
               out << vSum(px->hVol[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "tp") != outVars.end()){
             out << common <<"tp" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->tp[f][u]<< d;
               }
               out << 0 << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "mort") != outVars.end()){
             out << common <<"mort" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->mort[f][u]<< d;
               }
               out << 0 << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "tp_multiplier") != outVars.end()){
             out << common <<"tp_multiplier" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->getMultiplier("tp_multiplier",fTypes[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "mortCoef_multiplier") != outVars.end()){
             out << common <<"mortCoef_multiplier" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->getMultiplier("mortCoef_multiplier",fTypes[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "vMort") != outVars.end()){
             out << common <<"vMort" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->vMort[f][u]<< d;
               }
               out << vSum(px->vMort[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "vMortAdd") != outVars.end()){
             out << common <<"vMortAdd" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << px->vMortAdd[f][u]<< d;
               }
               out << vSum(px->vMortAdd[f]) << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "optDc") != outVars.end()){
             out << common <<"optDc" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->optDc[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "expectedAnnualIncome_timber") != outVars.end()){
             out << common <<"expectedAnnualIncome_timber" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                   out << d;
               }
               out << px->expectedAnnualIncome_timber[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "expectedAnnualIncome_carbon") != outVars.end()){
             out << common <<"expectedAnnualIncome_carbon" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                   out << d;
               }
               out << px->expectedAnnualIncome_carbon[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "expectedReturnsNotCorrByRa") != outVars.end()){
             out << common <<"expectedReturnsNotCorrByRa" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->expectedReturnsNotCorrByRa[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "expectedReturns") != outVars.end()){
             out << common <<"expectedReturns" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->expectedReturns[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "optFtChosen") != outVars.end()){
             out << common <<"optFtChosen" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->optFtChosen[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "optDcChosen") != outVars.end()){
             out << common <<"optDcChosen" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << px->optDcChosen[f] << d;
             }
             out << "\n";
         }
 
         if (std::find(outVars.begin(), outVars.end(), "deltaAreas") != outVars.end()){
             out << common <<"deltaAreas" << d;
             for(uint f=0;f<fTypes.size();f++){
               for(uint u=0;u<dClasses.size()-1;u++){
                 out << d;
               }
               out << d;
             }
             for(uint f=0;f<fTypes.size()+1;f++){
               for(uint f2=0;f2<fTypes.size()+1;f2++){ // +1 because there is also noforest
                 out << px->deltaArea[f][f2]<<d;
               }
             }
             out << "\n";
         }
 
 
 
       } // end for each pixel
     } // end for each region
   } // end filter
 } // end function printDebugPixelValues

Here is the call graph for this function:

Here is the caller graph for this function:

void printDetailedHV ( map< tr1::array< string, 4 >, double > hVol_byPrd )

Definition at line 1173 of file Output.cpp.

Referenced by ModelCoreSpatial::sumRegionalForData().

                                                                   {
    if(!MTHREAD->MD->getBoolSetting("outDetailedHv",DATA_NOW)) return;
    outFileName = baseDir+oDir+"results/detailedHV_"+scenarioName+oFileExt;
    int currentYear = MTHREAD->SCD->getYear();
    ofstream outHV(outFileName.c_str(), ios::app);
    if (!outHV){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for reading.");}
    //outHV << "scen" << d << "parName" << d << "country" << d << "region" << d << "forType" << d << "freeDim" << d << "prod" << d << "year" << d << "value" << d;
    for (uint r1=0;r1<l2r.size();r1++){
      string reg1Id = i2s(l1regIds.at(r1));
      string reg1Name = MTHREAD->MD->regId2RegSName(l1regIds.at(r1));
      for (uint r2=0;r2<l2r[r1].size();r2++){
        string reg2Id = i2s(l2r[r1][r2]);
        string reg2Name = MTHREAD->MD->regId2RegSName(l2r[r1][r2]);
        for(uint j=0;j<fTypes.size();j++){
            string ft = fTypes[j];
            for(uint u=0;u<dClasses.size();u++){
               string dc = dClasses[u];
               if(dc == "") continue;
               for(uint p=0;p<priPr.size();p++){
                 string prd = priPr[p];
                 tr1::array<string, 4> hVKey = {reg2Id, ft, dc, prd};
                 double value = findMap(hVol_byPrd,hVKey);
                 outHV << scenarioName << d;
                 outHV << "hVbyPrd" << d;
                 outHV << reg1Name << d;
                 outHV << reg2Name << d;
                 outHV << ft << d;
                 outHV << dc << d;
                 outHV << prd << d;
                 outHV << currentYear << d;
                 outHV << value << d;
                 outHV << "\n";
               }
            } //end dc
         } //end ft
      } // end r2
    } // end r1
 
 
    outHV.close();
 
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void printFinalOutput ( )

Definition at line 462 of file Output.cpp.

Referenced by Init::setInitLevel6().

                         {
   // we do this only if we choosed the outputHumanReadable settings, as we flush the data all in ones at the end.
   // oterwise we flush data every year
   if(oHRedeable){
     MTHREAD->MD->setErrorLevel(MSG_NO_MSG);
     printForestData(true);
     printProductData(true);
     MTHREAD->MD->setErrorLevel(MSG_ERROR);
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printForestData ( bool finalFlush )

Definition at line 474 of file Output.cpp.

Referenced by ModelData::getCachedInitialYear(), print(), and printFinalOutput().

                                       {
 
   if(!oForData) return;
   if(oHRedeable && !finalFlush) return;
 
   msgOut(MSG_INFO, "Printing forest data..");
   int currentYear = MTHREAD->SCD->getYear();
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/forestData"+oFileExt;
   } else {
     outFileName = baseDir+oDir+"results/forestData_"+scenarioName+oFileExt;
   }
   ofstream out (outFileName.c_str(), ios::app);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for writing.");}
   double outvalue;
   for(uint v=0;v<outForVariables.size();v++){
     vector<string>fTypes_temp = fTypes;
     if( outForVariables[v]=="expReturns" || outForVariables[v]=="sumExpReturns" ||  outForVariables[v]=="totalShareInvadedArea") {
       fTypes_temp.push_back(""); // adding an empty forest type to report for variables that doesn't have a forestType dimension
       vector<string> ftParents = MTHREAD->MD->getForTypeParents();
       fTypes_temp.insert(fTypes_temp.end(),ftParents.begin(),ftParents.end()); // also inserting forest type "parents" for expected returns
     }
     for (uint r1=0;r1<l2r.size();r1++){
       for (uint r2=0;r2<l2r[r1].size();r2++){
         for(uint ft=0;ft<fTypes_temp.size();ft++){
           if(forestDiamDetailedOutput){
             for(uint dc=0;dc<dClasses.size();dc++){ // an empty "" dc has been already added to the vector
               out << scenarioName << d;
               out << outForVariables[v] << d;
               out << MTHREAD->MD->regId2RegSName(l1regIds.at(r1)) << d;
               out << MTHREAD->MD->regId2RegSName(l2r[r1][r2]) << d;
               out << fTypes_temp[ft] << d;
               out << dClasses[dc] << d;
               if (oHRedeable){
                 for(int y=0;y<nYears;y++){
                   outvalue = MTHREAD->MD->getForData(outForVariables[v],l2r[r1][r2],fTypes_temp[ft],dClasses[dc],y+inYear);
                   out << outvalue << d;
                 }
                 out << "\n";
               } else {
                 outvalue = MTHREAD->MD->getForData(outForVariables[v],l2r[r1][r2],fTypes_temp[ft],dClasses[dc]);
                 out << currentYear << d;
                 out << outvalue << d;
                 out << "\n";
               }
             }
           } else {
             out << scenarioName << d;
             out << outForVariables[v] << d;
             out << MTHREAD->MD->regId2RegSName(l1regIds.at(r1)) << d;
             out << MTHREAD->MD->regId2RegSName(l2r[r1][r2]) << d;
             out << fTypes_temp[ft] << d;
             out << d;
             if (oHRedeable){
               for(int y=0;y<nYears;y++){
                 outvalue = MTHREAD->MD->getForData(outForVariables[v],l2r[r1][r2],fTypes_temp[ft],DIAM_ALL,y+inYear);
                 out << outvalue << d;
               }
               out << "\n";
             } else {
               outvalue = MTHREAD->MD->getForData(outForVariables[v],l2r[r1][r2],fTypes_temp[ft],DIAM_ALL);
               out << currentYear << d;
               out << outvalue << d;
               out << "\n";
             }
           }
         }
       }
     }
   }
   /*
   DataMap::const_iterator i;
   string key;
   vector <double> values;
   string parName;
   int regId;
   string forType;
   string diamClass;
   for(i=MTHREAD->MD->forDataMap.begin();i!=MTHREAD->MD->forDataMap.end();i++){
     key = i->first;
     values = i->second;
     MTHREAD->MD->unpackKeyForData(key, parName, regId, forType, diamClass);
     ModelRegion* REG = MTHREAD->MD->getRegion(regId);
     // we don't want to output data from residual region unless it's the world region we are speaking of
     if(REG->getIsResidual() && !(regId==wRegId_l1 || regId==wRegId_l2)) continue;
     out << scenarioName << d;
     out << parName << d;
     if (REG->getRegLevel()==2){
       ModelRegion* pREG = MTHREAD->MD->getRegion(REG->getParRegId());
       out << pREG->getRegSName() << d;
       out << REG->getRegSName() << d;
     } else if (REG->getRegLevel()==1){
       out << REG->getRegSName() << d;
       out << d;
     } else {
       out << d << d;
     }
     out << forType << d;
     out << diamClass << d;
     if (oHRedeable){
       for(int y=0;y<nYears;y++){
         out << MTHREAD->MD->getTimedData(values,y+inYear) << d;
       }
       out << "\n";
     } else {
       out << currentYear << d;
       out << MTHREAD->MD->getTimedData(values,currentYear) << d;
       out << "\n";
     }
   }
   */
   out.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printMaps ( )

Definition at line 451 of file Output.cpp.

Referenced by print().

                  {
   if(oLevel<OUTVL_MAPS) return;
   int cYear = MTHREAD->SCD->getYear();
   if ( find(mapsOYears.begin(), mapsOYears.end(), cYear) != mapsOYears.end() ){
     MTHREAD->GIS->printLayers();
     if(oLevel<OUTVL_BINMAPS) return;
     MTHREAD->GIS->printBinMaps();
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printOptLog	(	bool	optimal,
		int &	nIterations,
		double &	obj
	)

Definition at line 827 of file Output.cpp.

Referenced by ModelCore::runMarketModule(), and ModelCoreSpatial::runMarketModule().

                                                               {
   if(oLevel<OUTVL_AGGREGATED) return;
 
   ofstream out(logFilename.c_str(), ios::app);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+logFilename+" for writing.");}
   time_t now;
   time(&now);
   struct tm *current = localtime(&now);
   string timemessage = i2s(current->tm_hour)+":"+i2s(current->tm_min)+":"+ i2s(current->tm_sec);
   out << scenarioName << d << MTHREAD->SCD->getYear() << d << timemessage << d << optimal;
   out << d << nIterations << d << obj << "\n";
   out.close();
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void printProductData ( bool finalFlush )

Definition at line 589 of file Output.cpp.

Referenced by ModelData::getCachedInitialYear(), print(), and printFinalOutput().

                                        {
 
   if(!oProdData) return;
   if(oHRedeable && !finalFlush) return;
 
   msgOut(MSG_INFO, "Printing market data..");
   int currentYear = MTHREAD->SCD->getYear();
 
   if(oSingleFile){
     outFileName = baseDir+oDir+"results/productData"+oFileExt;
   } else {
     outFileName = baseDir+oDir+"results/productData_"+scenarioName+oFileExt;
   }
   ofstream out (outFileName.c_str(), ios::app);
   if (!out){ msgOut(MSG_CRITICAL_ERROR,"Error in opening the file "+outFileName+" for writing.");}
 
 
   //11042  hardWSawnW  11083  0.00230651
   //11042  hardWSawnW  11082  0.0390874
 
   //if(MTHREAD->SCD->getYear() == 2007){
 //  double test  = MTHREAD->MD->getProdData("rt",11042,"hardWSawnW",DATA_NOW);
 //  double test2 = MTHREAD->MD->getProdData("rt",11042,"hardWSawnW",DATA_NOW,"11083");
 //  double test3 = MTHREAD->MD->getProdData("rt",11042,"hardWSawnW",DATA_NOW,"11082");
 //  cout << test << '\t' << test2 << '\t' << test3 << endl;
 //  exit(0);
 //  }
 
   double outvalue;
   for(uint v=0;v<outProdVariables.size();v++){
     for (uint r1=0;r1<l2r.size();r1++){
       for (uint r2=0;r2<l2r[r1].size();r2++){
         for(uint p=0;p<allPr.size();p++){
 
           if(outProdVariables[v]=="rt"){
             for(uint r2b=0;r2b<l2r[r1].size();r2b++){
               out << scenarioName << d;
               out << outProdVariables[v] << d;
               out << MTHREAD->MD->regId2RegSName(l1regIds.at(r1)) << d;
               out << MTHREAD->MD->regId2RegSName(l2r[r1][r2]) << d;
               out << allPr[p] << d;
               out << l2r[r1][r2b] << d;
               if (oHRedeable){
                 for(int y=0;y<nYears;y++){
                   outvalue = MTHREAD->MD->getProdData(outProdVariables[v],l2r[r1][r2],allPr[p],y+inYear,i2s(l2r[r1][r2b]));
                   out << outvalue << d;
                 }
                 out << "\n";
               } else {
 //                if(MTHREAD->SCD->getYear() == 2007 && l2r[r1][r2] == 11042 && allPr[p] == "hardWSawnW" && (l2r[r1][r2b]== 11083 || l2r[r1][r2b]== 11082 )){
 //                outvalue = MTHREAD->MD->getProdData(outProdVariables[v],l2r[r1][r2],allPr[p],currentYear,i2s(l2r[r1][r2b]));
 //                cout << outvalue << endl;
 //                }
                 outvalue = MTHREAD->MD->getProdData(outProdVariables[v],l2r[r1][r2],allPr[p],currentYear,i2s(l2r[r1][r2b]));
                 out << currentYear << d;
                 out << outvalue << d;
                 out << "\n";
               }
             }
           } else {
             out << scenarioName << d;
             out << outProdVariables[v] << d;
             out << MTHREAD->MD->regId2RegSName(l1regIds.at(r1)) << d;
             out << MTHREAD->MD->regId2RegSName(l2r[r1][r2]) << d;
             out << allPr[p] << d;
             out << d;
             if (oHRedeable){
               for(int y=0;y<nYears;y++){
                 outvalue = MTHREAD->MD->getProdData(outProdVariables[v],l2r[r1][r2],allPr[p],y+inYear);
                 out << outvalue << d;
               }
               out << "\n";
             } else {
               outvalue = MTHREAD->MD->getProdData(outProdVariables[v],l2r[r1][r2],allPr[p]);
               out << currentYear << d;
               out << outvalue << d;
               out << "\n";
             }
 
           }
         }
       }
     }
   }
 
 
 
 
 /*
   DataMap::const_iterator i;
   string key;
   vector <double> values;
   string parName;
   int regId;
   string prod;
   string freeDim;
   for(i=MTHREAD->MD->prodDataMap.begin();i!=MTHREAD->MD->prodDataMap.end();i++){
     key = i->first;
     values = i->second;
     MTHREAD->MD->unpackKeyProdData(key, parName, regId, prod, freeDim);
     ModelRegion* REG = MTHREAD->MD->getRegion(regId);
     // we don't want to output data from residual region unless it's the world region we are speaking of
     if(REG->getIsResidual() && !(regId==wRegId_l1 || regId==wRegId_l2)) continue;
     out << scenarioName << d;
     out << parName << d;
     if (REG->getRegLevel()==2){
       ModelRegion* pREG = MTHREAD->MD->getRegion(REG->getParRegId());
       out << pREG->getRegSName() << d;
       out << REG->getRegSName() << d;
     } else if (REG->getRegLevel()==1){
       out << REG->getRegSName() << d;
       out << d;
     } else {
       out << d << d;
     }
     out << prod << d;
     out << freeDim << d;
     if (oHRedeable){
       for(int y=0;y<nYears;y++){
         out << MTHREAD->MD->getTimedData(values,y+inYear) << d;
       }
       out << "\n";
     } else {
       out << currentYear << d;
       out << MTHREAD->MD->getTimedData(values,currentYear) << d;
       out << "\n";
     }
   }
 
 */
   out.close();
 }