/* Physical constant */
#define SECPERRAD	13750.9871	/* seconds per radian of hour angle */
#define RADPERDAY	0.017214	/* radians of Earth orbit per julian day */
#define RADPERDEG	0.01745329	/* radians per degree */
#define MINDECL		-0.4092797      /* minimum declination (radians) */
#define DAYSOFF		11.25           /* julian day offset of winter solstice */
#define SRADDT		600.0           /* timestep for radiation routine (seconds) */
#define MC		12.011e-3       /* (kg mol-1) molecular weight of co2 */
#define MA		28.9644e-3      /* (kg mol-1) molecular weight of air */
#define MW		18.0148e-3      /* (kg mol-1) molecular weight of water */
#define R		8.3143          /* (m3 Pa mol-1 K-1) gas law constant */
#define G_STD		9.80665         /* (m s-2) standard gravitational accel. */
#define P_STD		101325.0        /* (Pa) standard pressure at 0.0 m elevation */
#define T_STD		288.15          /* (K) standard temp at 0.0 m elevation  */
#define CP		1010.0          /* (J kg-1 K-1) specific heat of air */
#define SBC		5.67e-8         /* (W/(m2 K4)) Stefan-Boltzmann constant */
#define LR_STD		0.0065          /* (-K m-1) standard temperature lapse rate */
#define EPS		0.62196351      /* (MW/MA) unitless ratio of molec weights */
#define PI		3.14159265      /* pi */
#define EPAR		4.55            /* (umol/J) PAR photon energy ratio */
#define ppe             2.6		//(mol/mol) photons absorbed by PSII per e- transported
#define rhow		1000		// density of water (kg m-3)
#define grav		9.81		// acceleration of gravity (m s-2)

/* Parameters for the Tair algorithm */
#define TDAYCOEF	0.45  // (dim) daylight air temperature coefficient (dim) // !!!0.45

/* Parameters for the snowpack algorithm */
#define SNOW_TCRIT	-6.0  // (deg C) critical temperature for snowmelt
#define SNOW_TRATE	0.042  // (cm/degC/day) snowmelt rate
#define T_F	0.0  // (deg C) fusion temperature

/* Parameters for the radiation algorithm */
#define TBASE          0.870  // (dim) max inst. trans., 0m, nadir, dry atm
#define ABASE        -6.1e-5  // (1/Pa) vapor pressure effect on transmittance
#define C                1.5  // (dim) radiation parameter
#define B0             0.013  // (dim) radiation parameter
#define B1             0.201  // (dim) radiation parameter
#define B2             0.185  // (dim) radiation parameter
#define RAIN_SCALAR     0.75  // (dim) correction to trans. for rain day
#define DIF_ALB          0.6  // (dim) diffuse albedo for horizon correction
#define SC_INT          1.32  // (MJ/m2/day) snow correction intercept
#define SC_SLOPE       0.096  // (MJ/m2/day/cm) snow correction slope

/* Parameters for the soils water transfer model */
#define mkpa		9.807	// 1 m of water is 9.807 kPa
extern int nslay; // #define nslay		4	// numbers of soil layers // FG
#define nslay_max   5
#define swp_fc		33.0	// soil water potential at field capacity (kPa)
// #define swp_wp		1500.0	// soil water potential at wilting point  (kPa) // FG

/* Various */
#define albedo_sw	0.2	// (DIM) site shortwave albedo
#define sm_init		0.5	// initial soil moisture content

//////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////
//Site Parameters QPFFN
//#define exp_site                2	// 1: Deciduous 2: Evergreen
extern int exp_site;
extern int model_version;
//#define meteo_site	       	4	// 1: complete	2: reduced
#define test		       	1	// 1: no	2: yes
//#define base_elev_cst		144.0	//MeteoFrance ???  base station elev m FG
//#define base_isoh_cst		90.0	//base station isoh cm FG
//#define site_lat_cst		48.38	//Etien et al.,  site lat deg FG
//#define site_elev_cst		144.0	//Etien et al.,  site elev m FG
//#define site_slp_cst		0.0	//site slope deg FG
//#define site_asp_cst		0.0	//site aspect deg FG
//#define site_isoh_cst		90.0	//site isoh cm FG
//#define site_ehoriz_cst		0.0	//site ehoriz deg FG
//#define site_whoriz_cst		0.0	//site whoriz deg FG
//#define tmax_lr_cst		-6.0	//tmax_lr_degC/km FG
//#define tmin_lr_cst		-3.0	//tmin_lr_degC/km FG
//#define CO2_lr_cst		-999.0	//ppmV/km PAD 291008 attente valeurs ???

//biophysical parameters
#define zref			17.5
#define z0			0.3
#define disd			12.5

// GGI constant for allocation module
//#define phenol_sd                    171.0  //ggi: standard deviation (in days) gaussian function phenology
//#define phenol_alpha                 6.0    //ggi: 04072014 exponent gaussian function phenology
//#define percent_fall                 0.33     //ggi: % of leaves that fall every year (evergreen in phenology type exp_site=3)//FG

//Control parameters
extern int ndays; // #define ndays        21535  FG
extern int year0; // #define year0         1941 FG
// #define yearn	      1999 FG
// #define year0_sim	1961 FG
// #define yearn_sim	1999 FG