Output Class Reference

Output methods More...

#include <Output.h>

Inheritance diagram for Output:

Collaboration diagram for Output:

Public Member Functions
	Output (ThreadManager *MTHREAD_h)
	Constructor.
	~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.
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.
void	msgOut (const int &msgCode_h, const int &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log.
void	msgOut (const int &msgCode_h, const double &msg_h, const bool &refreshGUI_h=true) const
	Overloaded function to print the output log.
int	s2i (const string &string_h) const
	string to integer conversion
double	s2d (const string &string_h) const
	string to double conversion
double	s2d (const string &string_h, const bool &replaceComma) const
	string to double conversion
bool	s2b (const string &string_h) const
	string to bool conversion
string	i2s (const int &int_h) const
	integer to string conversion
string	d2s (const double &double_h) const
	double to string conversion
string	b2s (const bool &bool_h) const
	bool to string conversion
vector< int >	s2i (const vector< string > &string_h) const
	string to integer conversion (vector)
vector< double >	s2d (const vector< string > &string_h, const bool &replaceComma=false) const
	string to double conversion (vector)
vector< bool >	s2b (const vector< string > &string_h) const
	string to bool conversion (vector)
vector< string >	i2s (const vector< int > &int_h) const
	integer to string conversion (vector)
vector< string >	d2s (const vector< double > &double_h) const
	double to string conversion (vector)
vector< string >	b2s (const vector< bool > &bool_h) const
	bool to string conversion (vector)
int	getType (const string &type_h) const
	Return a type according to enum TYPE_* from a string (eg: "string" -> TYPE_STRING (2))
void	refreshGUI () const
	Ping to periodically return the control to the GUI.
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)
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.
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.
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.
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.
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.
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.
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)
template<typename V >
void	debugMap (const map< iisskey, V > &mymap)
	Debug a map.
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)
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)
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)
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)
template<typename K >
double	getAvg (const vector< K > &v)
	Returns the average of the elements in the vector.
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).
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.
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
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
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.

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()

Output

(

ThreadManager *

MTHREAD_h

)

Constructor.

Definition at line 37 of file Output.cpp.

                                      {
  MTHREAD=MTHREAD_h;
}

~Output()

~Output

(

)

Definition at line 41 of file Output.cpp.

               {
}

Member Function Documentation

cleanScenario()

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.

                                                                 {
  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";
  }
}

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

Here is the call graph for this function:

Here is the caller graph for this function:

commonInit()

void commonInit

(

)

Definition at line 61 of file Output.cpp.

                  {
  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);
    }
  }
 
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

getOutputFieldDelimiter()

char getOutputFieldDelimiter

(

)

Definition at line 802 of file Output.cpp.

                               {
  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 ',';
  }
}

Referenced by commonInit().

Here is the call graph for this function:

Here is the caller graph for this function:

initCarbonBalance()

void initCarbonBalance

(

)

Definition at line 370 of file Output.cpp.

                         {
  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();
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initDebugOutput()

void initDebugOutput

(

)

Definition at line 170 of file Output.cpp.

                       {
    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();
    }
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initDebugPixelValues()

void initDebugPixelValues

(

)

Definition at line 207 of file Output.cpp.

                            {
    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;
    */
 
 
 
 
 
 
 
 
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initOptimisationLog()

void initOptimisationLog

(

)

Definition at line 131 of file Output.cpp.

                           {
  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();
  }
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initOutput()

void initOutput

(

)

Definition at line 48 of file Output.cpp.

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

Referenced by Init::setInitLevel3().

Here is the call graph for this function:

Here is the caller graph for this function:

initOutputForestData()

void initOutputForestData

(

)

Definition at line 296 of file Output.cpp.

                            {
  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();
 
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initOutputMaps()

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.

                      {
  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();
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

initOutputProductData()

void initOutputProductData

(

)

Definition at line 338 of file Output.cpp.

                             {
  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();
}

Referenced by initOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

print()

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.

                            {
  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);
}

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

Here is the call graph for this function:

Here is the caller graph for this function:

printCarbonBalance()

void printCarbonBalance

(

)

Definition at line 727 of file Output.cpp.

                          {
 
  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();
}

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

printDebugOutput()

void printDebugOutput

(

)

Definition at line 843 of file Output.cpp.

                        {
  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
}

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

printDebugPixelValues()

void printDebugPixelValues

(

)

Definition at line 880 of file Output.cpp.

                             {
 
  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

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

printDetailedHV()

void printDetailedHV

(

map< tr1::array< string, 4 >, double >

hVol_byPrd

)

Definition at line 1173 of file Output.cpp.

                                                                  {
   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();
 
 }

Referenced by ModelCoreSpatial::sumRegionalForData().

Here is the call graph for this function:

Here is the caller graph for this function:

printFinalOutput()

void printFinalOutput

(

)

Definition at line 462 of file Output.cpp.

                        {
  // 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);
  }
}

Referenced by Init::setInitLevel6().

Here is the call graph for this function:

Here is the caller graph for this function:

printForestData()

void printForestData

(

bool

finalFlush

)

Definition at line 474 of file Output.cpp.

                                      {
 
  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" ||  outForVariables[v]=="usableDeadTimber") {
      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();
}

Referenced by print(), and printFinalOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

printMaps()

void printMaps

(

)

Definition at line 451 of file Output.cpp.

                 {
  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();
  }
}

Referenced by print().

Here is the call graph for this function:

Here is the caller graph for this function:

printOptLog()

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

Definition at line 827 of file Output.cpp.

                                                              {
  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();
 
}

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

Here is the call graph for this function:

Here is the caller graph for this function:

printProductData()

void printProductData

(

bool

finalFlush

)

Definition at line 589 of file Output.cpp.

                                       {
 
  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();
}

Referenced by print(), and printFinalOutput().

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

allPr

vector<string> allPr

private

Definition at line 103 of file Output.h.

Referenced by commonInit(), and printProductData().

baseDir

string baseDir

private

Definition at line 83 of file Output.h.

Referenced by commonInit(), initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), initOutputMaps(), initOutputProductData(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

d

char d

private

Definition at line 80 of file Output.h.

Referenced by cleanScenario(), commonInit(), initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), initOutputProductData(), printCarbonBalance(), printDebugOutput(), printDebugPixelValues(), printDetailedHV(), printForestData(), printOptLog(), and printProductData().

dClasses

vector<string> dClasses

private

Definition at line 107 of file Output.h.

Referenced by commonInit(), initDebugPixelValues(), printDebugOutput(), printDebugPixelValues(), printDetailedHV(), and printForestData().

debugFilename

string debugFilename

private

Definition at line 114 of file Output.h.

Referenced by initDebugOutput(), and printDebugOutput().

debugPxValuesFilename

string debugPxValuesFilename

private

Definition at line 115 of file Output.h.

Referenced by initDebugPixelValues(), and printDebugPixelValues().

expReturnsDebug

vector<vector < vector <vector <vector <double> > > > > expReturnsDebug

l2_region, for type, d.c., pr prod, variable name

Definition at line 75 of file Output.h.

Referenced by printDebugOutput(), and ModelCore::runManagementModule().

expReturnsDebugVariables

vector<string> expReturnsDebugVariables

Definition at line 76 of file Output.h.

Referenced by initDebugOutput(), and printDebugOutput().

forestDiamDetailedOutput

bool forestDiamDetailedOutput

private

Definition at line 100 of file Output.h.

Referenced by commonInit(), and printForestData().

fTypes

vector<string> fTypes

private

Definition at line 106 of file Output.h.

Referenced by commonInit(), initDebugPixelValues(), printDebugOutput(), printDebugPixelValues(), printDetailedHV(), and printForestData().

inYear

int inYear

private

Definition at line 81 of file Output.h.

Referenced by commonInit(), initOutputForestData(), initOutputProductData(), printCarbonBalance(), printForestData(), and printProductData().

l1regIds

vector<int> l1regIds

private

Definition at line 104 of file Output.h.

Referenced by commonInit(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

l2r

vector< vector <int> > l2r

private

Definition at line 105 of file Output.h.

Referenced by commonInit(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

logFilename

string logFilename

private

Definition at line 113 of file Output.h.

Referenced by initOptimisationLog(), and printOptLog().

mapsOYears

vector<int> mapsOYears

private

Definition at line 90 of file Output.h.

Referenced by commonInit(), and printMaps().

nAllPr

int nAllPr

private

Definition at line 111 of file Output.h.

Referenced by commonInit().

nL2r

int nL2r

private

Definition at line 112 of file Output.h.

Referenced by commonInit().

nPriPr

int nPriPr

private

Definition at line 109 of file Output.h.

Referenced by commonInit().

nSecPr

int nSecPr

private

Definition at line 110 of file Output.h.

Referenced by commonInit().

nYears

int nYears

private

Definition at line 82 of file Output.h.

Referenced by commonInit(), initOutputForestData(), initOutputProductData(), print(), printForestData(), and printProductData().

oCarbonData

bool oCarbonData

private

Definition at line 98 of file Output.h.

Referenced by commonInit(), initCarbonBalance(), and printCarbonBalance().

oDir

string oDir

private

Definition at line 84 of file Output.h.

Referenced by commonInit(), initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), initOutputMaps(), initOutputProductData(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

oFileExt

string oFileExt

private

Definition at line 86 of file Output.h.

Referenced by commonInit(), initCarbonBalance(), initOutputForestData(), initOutputProductData(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

oForData

bool oForData

private

Definition at line 96 of file Output.h.

Referenced by commonInit(), initOutputForestData(), and printForestData().

oHRedeable

bool oHRedeable

private

Definition at line 87 of file Output.h.

Referenced by commonInit(), initOutputForestData(), initOutputProductData(), printFinalOutput(), printForestData(), and printProductData().

oLevel

int oLevel

private

Definition at line 79 of file Output.h.

Referenced by commonInit(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputMaps(), printDebugOutput(), printDebugPixelValues(), printMaps(), and printOptLog().

oProdData

bool oProdData

private

Definition at line 97 of file Output.h.

Referenced by commonInit(), initOutputProductData(), and printProductData().

oSingleFile

bool oSingleFile

private

Definition at line 88 of file Output.h.

Referenced by commonInit(), initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), initOutputProductData(), printCarbonBalance(), printForestData(), and printProductData().

outFileName

string outFileName

private

Definition at line 93 of file Output.h.

Referenced by initCarbonBalance(), initOutputForestData(), initOutputProductData(), printCarbonBalance(), printDetailedHV(), printForestData(), and printProductData().

outForVariables

vector<string> outForVariables

private

Definition at line 94 of file Output.h.

Referenced by commonInit(), and printForestData().

outProdVariables

vector<string> outProdVariables

private

Definition at line 95 of file Output.h.

Referenced by commonInit(), and printProductData().

outStepRange

int outStepRange

private

Definition at line 99 of file Output.h.

Referenced by commonInit(), and print().

oYears

vector<int> oYears

private

Definition at line 89 of file Output.h.

Referenced by commonInit(), and print().

pDClasses

vector<string> pDClasses

private

includes an empty string for variables without diameter attribute

production diameter classes: exclude the fist diameter class below 15 cm

Definition at line 108 of file Output.h.

Referenced by commonInit().

priPr

vector<string> priPr

private

Definition at line 101 of file Output.h.

Referenced by commonInit(), printDebugOutput(), and printDetailedHV().

scenarioName

string scenarioName

private

Definition at line 85 of file Output.h.

Referenced by cleanScenario(), commonInit(), initCarbonBalance(), initDebugOutput(), initDebugPixelValues(), initOptimisationLog(), initOutputForestData(), initOutputMaps(), initOutputProductData(), printCarbonBalance(), printDebugOutput(), printDebugPixelValues(), printDetailedHV(), printForestData(), printOptLog(), and printProductData().

secPr

vector<string> secPr

private

Definition at line 102 of file Output.h.

Referenced by commonInit().

spMode

bool spMode

private

Definition at line 116 of file Output.h.

Referenced by commonInit(), initDebugOutput(), initDebugPixelValues(), printDebugOutput(), and printDebugPixelValues().

wRegId_l1

int wRegId_l1

private

Definition at line 91 of file Output.h.

Referenced by commonInit().

wRegId_l2

int wRegId_l2

private

Definition at line 92 of file Output.h.

Referenced by commonInit().

---

The documentation for this class was generated from the following files:

• /exports/ffsm/src_latest/src/Output.h
• /exports/ffsm/src_latest/src/Output.cpp