Class responsable to keep the logbook of the Carbon Balance. More...

#include <Carbon.h>

Inheritance diagram for Carbon:

Collaboration diagram for Carbon:

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

	~Carbon ()

double	getStock (const int &regId, const int &stock_type) const
	Returns the current stock of carbon [Mt CO2]. More...

double	getCumSavedEmissions (const int &regId, const int &em_type) const
	Returns the current cumulative saved emissions by type [Mt CO2]. More...

void	registerHarvesting (const double &value, const int &regId, const string &fType)
	Registers the harvesting of trees increasing the value of cumEmittedForOper. More...

void	registerDeathBiomass (const double &value, const int &regId, const string &fType)
	Registers the "death" of a given amount of biomass, storing it in the deathBiomass map. More...

void	registerProducts (const double &value, const int &regId, const string &productName)
	Registers the production of a given amount of products, storing it in the products maps. Also increase material substitution. More...

void	registerTransports (const double &distQ, const int &regId)
	Registers the quantities emitted by transport of wood FROM a given region. More...

void	initialiseDeathBiomassStocks (const vector< double > &deathByFt, const int &regId)
	Initialises the stocks of death biomass for the avgLive_* years before the simulation starts. More...

void	initialiseProductsStocks (const vector< double > &qByProduct, const int &regId)
	Initialises the stocks of products for the avgLive_* years before the simulation starts. More...

void	initialiseEmissionCounters ()
	Initialises the emission counters to zero. More...

void	HWP_eol2energy ()
	Computes the energy substitution for the quota of HWP that reachs end of life and doesn't go to landfill. More...

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

Private Member Functions
void	addSavedEmissions (const double &value, const int &regId, const int &em_type)
	Increases the value to the saved emissions for a given type and region. More...

double	getRemainingStock (const double &initialValue, const double &halfLife, const double &years) const
	Apply a single exponential decay model to retrieve the remining stock given the initial stock, the half life and the time passed from stock formation. More...

Private Attributes
map< iiskey, double >	deathBiomassInventory
	Map that register the death of biomass by year, l2_region and forest type (inventoried)[Mm^3 wood]. More...

map< iiskey, double >	deathBiomassExtra
	Map that register the death of biomass by year, l2_region and forest type (non-inventoried biomass: branches, roots..) [Mm^3 wood]. More...

map< iiskey, double >	products
	Map that register the production of a given product by year, l2_region and product [Mm^3 wood]. More...

map< int, double >	cumSubstitutedEnergy
	Map that store the cumulative CO2 substituted for energy consumption, by l2_region [Mt CO2]. More...

map< int, double >	cumSubstitutedMaterial
	Map that store the cumulative CO2 substituted using less energivory material, by l2_region [Mt CO2]. More...

map< int, double >	cumEmittedForOper
	Map that store emissions for forest operations, including transport, by l2_region [Mt CO2]. More...

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

Detailed Description

Class responsable to keep the logbook of the Carbon Balance.

Author: Antonello Lobianco

A single istance of this class exists and is available trought the global MTHREAD->CBAL pointer.

It consits of functions to track a carbon-related event and store the information in STL maps that either register the events (for the stocks) or contain the cumulated carbon (for emission flows).

Carbon pools are stored as Mm^3 wood while and emission cumulated counters are directly in Mt CO2.

getStock() and getCumSavedEmissions() are then used to report the current levels of carbon in the stock or emitted/substituted.

Definition at line 50 of file Carbon.h.

Constructor & Destructor Documentation

Carbon ( ThreadManager * MTHREAD_h )

Constructor.

Definition at line 32 of file Carbon.cpp.

                                       {
   MTHREAD=MTHREAD_h;
 }

~Carbon ( )

Definition at line 36 of file Carbon.cpp.

36 {

37 }

Member Function Documentation

void addSavedEmissions	(	const double &	value,
		const int &	regId,
		const int &	em_type
	)

private

Increases the value to the saved emissions for a given type and region.

Definition at line 324 of file Carbon.cpp.

Referenced by HWP_eol2energy(), registerHarvesting(), registerProducts(), and registerTransports().

                                                                                      {
   switch (em_type){
     case EM_ENSUB:
       incrMapValue(cumSubstitutedEnergy, regId, value);
       break;
     case EM_MATSUB:
       incrMapValue(cumSubstitutedMaterial, regId, value);
       break;
     case EM_FOROP:
       incrMapValue(cumEmittedForOper, regId, -value);
       break;
     default:
       msgOut(MSG_CRITICAL_ERROR,"Unexpected em_type in function getCumSavedEmissions");
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

double getCumSavedEmissions	(	const int &	regId,
		const int &	em_type
	)		const

Returns the current cumulative saved emissions by type [Mt CO2].

Definition at line 138 of file Carbon.cpp.

Referenced by Output::printCarbonBalance().

                                                                         {
   switch (em_type){
     case EM_ENSUB:
       return findMap(cumSubstitutedEnergy, regId);
       break;
     case EM_MATSUB:
       return findMap(cumSubstitutedMaterial, regId);
       break;
     case EM_FOROP:
       return -findMap(cumEmittedForOper, regId);
       break;
     default:
       msgOut(MSG_CRITICAL_ERROR,"Unexpected em_type in function getCumSavedEmissions");
   }
   return 0.0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

double getRemainingStock	(	const double &	initialValue,
		const double &	halfLife,
		const double &	years
	)		const

private

Apply a single exponential decay model to retrieve the remining stock given the initial stock, the half life and the time passed from stock formation.

Definition at line 341 of file Carbon.cpp.

Referenced by getStock(), and HWP_eol2energy().

                                                                                                          {
   // // TODO: remove this test
   //if(years>0) return 0.0;
   //return initialValue;
 
   double k  = log(2)/halfLife;
   return initialValue*exp(-k*years);
 }

Here is the caller graph for this function:

double getStock	(	const int &	regId,
		const int &	stock_type
	)		const

Returns the current stock of carbon [Mt CO2].

Parameters

reg
stock_type

Returns: the Carbon stocked in a given sink

For product sink:

for primary products it includes the primary products exported out of the country, but not those exported to other regions or used in the region as these are assumed to be totally transformed to secondary products;
for secondary products it includes those produced in the region from locally or regionally imported primary product plus those secondary products imported from other regions, less those exported to other regions. It doesn't include the secondary products imported from abroad the country.

Definition at line 53 of file Carbon.cpp.

Referenced by Output::printCarbonBalance().

                                                                {
   double toReturn = 0.0;
   int currentYear = MTHREAD->SCD->getYear();
   int initialYear = MTHREAD->MD->getIntSetting("initialYear");
   switch (stock_type){
     case STOCK_PRODUCTS: {
       vector <string> priProducts = MTHREAD->MD->getStringVectorSetting("priProducts");
       vector <string> secProducts = MTHREAD->MD->getStringVectorSetting("secProducts");
       vector <string> allProducts = priProducts;
       allProducts.insert( allProducts.end(), secProducts.begin(), secProducts.end() );
       for(uint i=0;i<allProducts.size();i++){
         double coeff = MTHREAD->MD->getProdData("co2content_products",regId,allProducts[i],DATA_NOW,""); //  [kg CO2/m^3 wood]
         double life  = MTHREAD->MD->getProdData("avgLife_products",regId,allProducts[i],DATA_NOW,"");
         //for(int y=currentYear;y>currentYear-life;y--){ // ok
         //  iiskey key(y,regId,allProducts[i]);
         //  toReturn += findMap(products,key,MSG_NO_MSG,0.0)*coeff/1000;
         //}
         for(int y=(initialYear-100);y<=currentYear;y++){
           iiskey key(y,regId,allProducts[i]);
           double originalStock = findMap(products,key,MSG_NO_MSG,0.0);
           double remainingStock = getRemainingStock(originalStock,life,currentYear-y);
           toReturn += remainingStock*coeff/1000;
         }
       }
       break;
     }
     case STOCK_INV:{
       vector <string> fTypes = MTHREAD->MD->getForTypeIds();
       for(uint i=0;i<fTypes.size();i++){
         // units:
         // co2content_inventory: [Kg CO2 / m^3 wood]
         // co2content_extra:     [Kg CO2 / m^3 inventaried wood]
         double coeff = MTHREAD->MD->getForData("co2content_inventory",regId,fTypes[i],"",DATA_NOW); //  [kg CO2/m^3 wood]
         double life  = MTHREAD->MD->getForData("avgLive_deathBiomass_inventory",regId,fTypes[i],"",DATA_NOW);
         // PART A: from death biomass..
         //for(int y=currentYear;y>currentYear-life;y--){ // ok
         //  iiskey key(y,regId,fTypes[i]);
         //  toReturn += findMap(deathBiomassInventory,key,MSG_NO_MSG)*coeff/1000;
         //}
         for(int y=(initialYear-100);y<=currentYear;y++){
           iiskey key(y,regId,fTypes[i]);
           double originalStock = findMap(deathBiomassInventory,key,MSG_NO_MSG,0.0);
           double remainingStock = getRemainingStock(originalStock,life,currentYear-y);
           toReturn += remainingStock*coeff/1000;
         }
 
         // PART B: from inventory volumes
         toReturn += MTHREAD->MD->getForData("vol",regId,fTypes[i],DIAM_ALL,DATA_NOW)*coeff/1000;
       }
       break;
 
     }
     case STOCK_EXTRA:{
       vector <string> fTypes = MTHREAD->MD->getForTypeIds();
       for(uint i=0;i<fTypes.size();i++){
         // units:
         // co2content_inventory: [Kg CO2 / m^3 wood]
         // co2content_extra:     [Kg CO2 / m^3 inventaried wood]
         double coeff = MTHREAD->MD->getForData("co2content_extra",regId,fTypes[i],"",DATA_NOW); //  [kg CO2/m^3 wood]
         double life  = MTHREAD->MD->getForData("avgLive_deathBiomass_extra",regId,fTypes[i],"",DATA_NOW);
         // PART A: from death biomass..
         //for(int y=currentYear;y>currentYear-life;y--){ // ok
         //  iiskey key(y,regId,fTypes[i]);
         //  toReturn += findMap(deathBiomassExtra,key,MSG_NO_MSG),0.0*coeff/1000;
         //}
         for(int y=(initialYear-100);y<=currentYear;y++){
           iiskey key(y,regId,fTypes[i]);
           double originalStock = findMap(deathBiomassExtra,key,MSG_NO_MSG,0.0);
           double remainingStock = getRemainingStock(originalStock,life,currentYear-y);
           toReturn += remainingStock*coeff/1000;
         }
         // PART B: from inventory volumes
         double extraBiomass_ratio = MTHREAD->MD->getForData("extraBiomass_ratio",regId,fTypes[i],"",DATA_NOW);
         toReturn += MTHREAD->MD->getForData("vol",regId,fTypes[i],DIAM_ALL,DATA_NOW)*extraBiomass_ratio*coeff/1000;
       }
       break;
     }
     default:
       msgOut(MSG_CRITICAL_ERROR,"Unexpected stock_type in function getStock");
   }
   return toReturn;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void HWP_eol2energy ( )

Computes the energy substitution for the quota of HWP that reachs end of life and doesn't go to landfill.

Definition at line 288 of file Carbon.cpp.

Referenced by ModelCoreSpatial::registerCarbonEvents().

                       {
 
   int currentYear = MTHREAD->SCD->getYear();
   int initialYear = MTHREAD->MD->getIntSetting("initialYear");
   vector <string> priProducts = MTHREAD->MD->getStringVectorSetting("priProducts");
   vector <string> secProducts = MTHREAD->MD->getStringVectorSetting("secProducts");
   vector <string> allProducts = priProducts;
   allProducts.insert( allProducts.end(), secProducts.begin(), secProducts.end() );
 
   vector<int> regIds = MTHREAD->MD->getRegionIds(2);
   for (uint r=0;r<regIds.size();r++){
     double regId = regIds[r];
     for(uint i=0;i<allProducts.size();i++){
       string pr = allProducts[i];
       double life  = MTHREAD->MD->getProdData("avgLife_products",regId,pr,DATA_NOW,"");
       double eol2e_share      = MTHREAD->MD->getProdData("eol2e_share",regId,pr,DATA_NOW,"");
       double subEnergyCoeff   = MTHREAD->MD->getProdData("co2sub_energy",regId,pr,DATA_NOW,"");
       if(eol2e_share > 0 && subEnergyCoeff>0){
         for(int y=(initialYear-100);y<currentYear;y++){ // notice the minor operator and not minor equal: energy substitution for products produced this year assigned to the following year, otherwise double counring in the process of making dicrete the exponential function
           iiskey key(y,regId,pr);
           double originalStock = findMap(products,key,MSG_NO_MSG,0.0);
           double remainingStockLastYear = getRemainingStock(originalStock,life,currentYear-y-1);
           double remainingStockThisYear = getRemainingStock(originalStock,life,currentYear-y);
           double eofThisYear = remainingStockLastYear-remainingStockThisYear;
           addSavedEmissions(subEnergyCoeff*eofThisYear*eol2e_share/1000,regId,EM_ENSUB);
         }
       }
     }
   }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initialiseDeathBiomassStocks	(	const vector< double > &	deathByFt,
		const int &	regId
	)

Initialises the stocks of death biomass for the avgLive_* years before the simulation starts.

Definition at line 169 of file Carbon.cpp.

Referenced by ModelCoreSpatial::initialiseCarbonModule().

                                                                                        {
   // it must initialize in the past the death biomass taking the value of the first year
   vector <string> fTypes = MTHREAD->MD->getForTypeIds();
   if(fTypes.size() != deathByFt.size()) {msgOut(MSG_CRITICAL_ERROR,"deathByFt and fTypes have different lenght!");}
   int currentYear = MTHREAD->SCD->getYear();
   //int initialYear = MD->getIntSetting("initialYear");
 
   for(uint i=0;i<fTypes.size();i++){
 //    double life_inventory     = MTHREAD->MD->getForData("avgLive_deathBiomass_inventory",regId,fTypes[i],"",DATA_NOW);
 //    double life_extra         = MTHREAD->MD->getForData("avgLive_deathBiomass_extra",regId,fTypes[i],"",DATA_NOW);
       double extraBiomass_ratio = MTHREAD->MD->getForData("extraBiomass_ratio",regId,fTypes[i],"",DATA_NOW);
 
 //    for(int y=currentYear;y>currentYear-life_inventory;y--){
 //        iiskey key(y,regId,fTypes[i]);
 //        pair<iiskey,double> mypair(key,deathByFt.at(i));
 //        deathBiomassInventory.insert(mypair);
 //    }
 //    for(int y=currentYear;y>currentYear-life_extra;y--){
 //        iiskey key(y,regId,fTypes[i]);
 //        pair<iiskey,double> mypair(key,deathByFt.at(i)*extraBiomass_ratio);
 //        deathBiomassExtra.insert(mypair);
 //    }
 
     for(int y=currentYear;y>currentYear-100;y--){
         iiskey key(y,regId,fTypes[i]);
         pair<iiskey,double> mypairInventory(key,deathByFt.at(i));
         pair<iiskey,double> mypairExtra(key,deathByFt.at(i)*extraBiomass_ratio);
         deathBiomassInventory.insert(mypairInventory);
         deathBiomassExtra.insert(mypairExtra);
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initialiseEmissionCounters ( )

Initialises the emission counters to zero.

Definition at line 158 of file Carbon.cpp.

Referenced by ModelCoreSpatial::initialiseCarbonModule().

                                   {
   vector<int> regIds = MTHREAD->MD->getRegionIds(2);
   for (uint i=0;i<regIds.size();i++){
     pair<int,double> mypair(regIds[i],0.0);
     cumSubstitutedEnergy.insert(mypair);
     cumSubstitutedMaterial.insert(mypair);
     cumEmittedForOper.insert(mypair);
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void initialiseProductsStocks	(	const vector< double > &	qByProduct,
		const int &	regId
	)

Initialises the stocks of products for the avgLive_* years before the simulation starts.

Definition at line 203 of file Carbon.cpp.

Referenced by ModelCoreSpatial::initialiseCarbonModule().

                                                                                     {
   // it must initialize in the past the products taking the value of the first year
   vector <string> priProducts = MTHREAD->MD->getStringVectorSetting("priProducts");
   vector <string> secProducts = MTHREAD->MD->getStringVectorSetting("secProducts");
   vector <string> allProducts = priProducts;
   allProducts.insert( allProducts.end(), secProducts.begin(), secProducts.end() );
   if(allProducts.size() != qByProduct.size()) {msgOut(MSG_CRITICAL_ERROR,"allProducts and qByProduct have different lenght!");}
   int currentYear = MTHREAD->SCD->getYear();
   for(uint i=0;i<allProducts.size();i++){
     double life  = MTHREAD->MD->getProdData("avgLife_products",regId,allProducts[i],DATA_NOW);
     //for(int y=currentYear;y>currentYear-life;y--){
     for(int y=currentYear;y>currentYear-100;y--){
         iiskey key(y,regId,allProducts[i]);
         pair<iiskey,double> mypair(key,qByProduct.at(i));
         products.insert(mypair);
     }
   }
   //cout << "" << endl;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void registerDeathBiomass	(	const double &	value,
		const int &	regId,
		const string &	fType
	)

Registers the "death" of a given amount of biomass, storing it in the deathBiomass map.

Definition at line 243 of file Carbon.cpp.

Referenced by ModelCoreSpatial::registerCarbonEvents().

                                                                                         {
   int year = MTHREAD->SCD->getYear();
   iiskey key(year,regId,fType);
   double extraBiomass_ratio = MTHREAD->MD->getForData("extraBiomass_ratio",regId,fType,"",DATA_NOW);
   //pair<iiskey,double> mypairInventory(key,value);
   //pair<iiskey,double> mypairExtra(key,value*extraBiomass_ratio);
   incrOrAddMapValue(deathBiomassInventory, key, value);
   incrOrAddMapValue(deathBiomassExtra, key, value*extraBiomass_ratio);
   //deathBiomassInventory.insert(mypairInventory);
   //deathBiomassExtra.insert(mypairExtra);
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void registerHarvesting	(	const double &	value,
		const int &	regId,
		const string &	fType
	)

Registers the harvesting of trees increasing the value of cumEmittedForOper.

Definition at line 225 of file Carbon.cpp.

Referenced by ModelCoreSpatial::registerCarbonEvents().

                                                                                        {
   double convCoeff = MTHREAD->MD->getForData("forOperEmissions",regId,fType,""); //  Kg of CO2 emitted per cubic meter of forest operations
   // units:
   // value: Mm^3
   // convCoeff: Kg CO2/m^3 wood
   // desidered output: Mt CO2
   // ==> I must divide by 1000
   addSavedEmissions(-convCoeff*value/1000,regId,EM_FOROP);
   // Add the extraBiomass associated to the harvested volumes to the deathBiomassExtra pool
   int    year               = MTHREAD->SCD->getYear();
   double extraBiomass_ratio = MTHREAD->MD->getForData("extraBiomass_ratio",regId,fType,"",DATA_NOW);
   double newDeathBiomass = value*extraBiomass_ratio;
   iiskey key(year,regId,fType);
   incrOrAddMapValue(deathBiomassExtra, key, newDeathBiomass);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void registerProducts	(	const double &	value,
		const int &	regId,
		const string &	productName
	)

Registers the production of a given amount of products, storing it in the products maps. Also increase material substitution.

Definition at line 257 of file Carbon.cpp.

Referenced by ModelCoreSpatial::registerCarbonEvents().

                                                                                           {
   // Registering the CO2 stock embedded in the product...
   int year = MTHREAD->SCD->getYear();
   iiskey key(year,regId,productName);
   pair<iiskey,double> mypair(key,value);
   products.insert(mypair);
   // registering the substituted CO2 for energy and material..
   double subEnergyCoeff   = MTHREAD->MD->getProdData("co2sub_energy",regId,productName,DATA_NOW,"");
   double subMaterialCoeff = MTHREAD->MD->getProdData("co2sub_material",regId,productName,DATA_NOW,"");
   // units:
   // value: Mm^3
   // subEnergyCoeff and subMaterialCoeff:  [kgCO2/m^3 wood]
   // desidered output: Mt CO2
   // ==> I must divide by 1000
   addSavedEmissions(subMaterialCoeff*value/1000,regId,EM_MATSUB); // EM_ENSUB are now included in eol2energy, instantaneous for fuelwoods
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void registerTransports	(	const double &	distQ,
		const int &	regId
	)

Registers the quantities emitted by transport of wood FROM a given region.

Definition at line 277 of file Carbon.cpp.

Referenced by ModelCoreSpatial::registerCarbonEvents().

                                                                 {
   // units:
   // distQ: km*Mm^3
   // transportEmissionsCoeff:  [Kg CO2 / (km*m^3) ]
   // desidered output: Mt CO2
   // ==> I must divide by 1000
   double transportEmissionsCoeff = MTHREAD->MD->getDoubleSetting("transportEmissionsCoeff",DATA_NOW);
   addSavedEmissions(-transportEmissionsCoeff*distQ/1000,regId,EM_FOROP);
 }