Skip to content
GitLab
Select Git revision
  • zonalstats_histogram default
Find file BlameHistoryPermalink
FG.cpp 22.84 KiB
1 
/*================================================*/
2 
/*            Plants Functionals Groups           */
3 
/*                 Definition Class               */
4 
/*================================================*/
5
6
7 
#include <cmath>
8 
#include <iostream>
9 
#include <cstring>
10 
#include <fstream>
11 
#include <cstdio>
12 
#include "FG.h"
13 
#include "FGUtils.h"
14 
#include "Params.h"
15
16 
using namespace std;
17
18 
  /* Note : New version of FG constructor using Matt T. parameters handler utilities (Params.h) */
19
20
21 
/*-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-*/
22 
/* Constructors                                                                                    */
23 
/*-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-*/
24
25 
FG::FG() : m_Name(""), m_M(0), m_L(1), m_MaxA(ANone), m_ImmSize(0.0), m_StrataMax(0), m_Strata(0,1000), /* Life history*/
26 
m_PoolL(PTcount,0), m_InnateDorm(false), m_PotentialFecundity(100), /* Propagule biology */
27 
m_ActiveGerm(Rcount, PC100), m_Tolerance(LScount, vector<bool>(Rcount, true)), /* Light response */
28 
m_Dispersed(false), m_disp50(0.0), m_disp99(0.0), m_dispLD(0.0), /* Dispersal module */
29 
m_SoilContrib(0.0), m_SoilLow(0.0), m_SoilHigh(0.0), /* Soil response */
30 
 m_SoilActiveGerm(Rcount, PC100), m_SoilTolerance(LScount, vector<Fract>(Rcount, PC100)), /* Soil response */
31 
m_DistResponse(FGresponse()), /* Disturbance response */
32 
m_FireResponse(FGresponse()), m_Flamm(0.0), /* Fire response */
33 
m_DroughtResponse(FGresponse()), m_DroughtSD(2,0.0), m_CountModToSev(0), m_CountSevMort(0), m_DroughtRecovery(0), /* Drought response */
34 
m_IsAlien(false) /* Alien module */
35 
{
36 
  	/* Nothing to do */
37 
}
38
39 
/*-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-*/
40
41 
void FG::getSuccParams(const GSP& glob_params, const string& PFG_LifeHistoryFile)
42 
{
43 
	/* 1. check parameter file existence */
44 
	testFileExist("--PFG_LIFE_HISTORY_PARAMS--", PFG_LifeHistoryFile, false);
45
46 
	/* 2. read succession parameters */
47 
	par::Params SuccParams(PFG_LifeHistoryFile.c_str(), " = \"", "#"); /* opening PFG life history traits parameters file */
48
49 
	cout << endl;
50 
	cout << "*********************************************" << endl;
51 
	cout << "** PFG : " << SuccParams.get_val<string>("NAME")[0] << endl;
52 
	cout << "*********************************************" << endl;
53 
	cout << endl;
54 
	cout << "> Succession files opened" << endl;
55
56 
	/* 3. fill FG object according to given parameters */
57
58 
	/* PFG Life History parameters filling =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- */
59 
	m_Name = SuccParams.get_val<string>("NAME")[0];
60 
	m_M = SuccParams.get_val<int>("MATURITY")[0];
61 
	m_L = SuccParams.get_val<int>("LONGEVITY")[0];
62 
	if (m_M >= m_L)
63 
	{
64 
		cerr << "!!! MATURITY is superior or equal to LONGEVITY. Please check!" << endl;
65 
		terminate();
66 
	}
67
68 
	m_MaxA = Abund(SuccParams.get_val<int>("MAX_ABUNDANCE")[0]);
69 
	m_ImmSize = FractToDouble(Fract(SuccParams.get_val<int>("IMM_SIZE")[0]));
70 
	//m_StrataMax = SuccParams.get_val<int>("STRATA")[0];