Merge commit '3c36b0e92817ce4d282ca813b8991e9fa43c7d53'

Author: ftomei

agrolib/hydrall/hydrall.cpp

@@ -0,0 +1,40 @@
+/*!
+    \name hydrall.cpp
+    \brief 
+    \authors Antonio Volta, Caterina Toscano
+	
+*/
+
+
+#include <stdio.h>
+#include <math.h>
+#include "crit3dDate.h"
+#include "commonConstants.h"
+#include "hydrall.h"
+
+
+bool computeHydrall(Crit3DDate myDate, double myTemperature, double myElevation)
+{
+    getCO2(myDate, myTemperature, myElevation);
+    return true;
+}
+
+double getCO2(Crit3DDate myDate, double myTemperature, double myElevation)
+{
+    double atmCO2 ; // fitting from data of Mauna Loa,Hawaii
+    if (myDate.year < 1990)
+    {
+        atmCO2= 280 * exp(0.0014876*(myDate.year -1840));//exponential change in CO2 concentration (ppm)
+    }
+    else
+    {
+        atmCO2= 350 * exp(0.00630*(myDate.year - 1990));
+    }
+    atmCO2 += 3*cos(2*PI*getDoyFromDate(myDate)/365.0);		     // to consider the seasonal effects
+    return atmCO2*getPressureFromElevation(myTemperature, myElevation)/1000000  ;   // [Pa] in +- ppm/10
+}
+
+double getPressureFromElevation(double myTemperature, double myElevation)
+{
+    return SEA_LEVEL_PRESSURE * exp((- GRAVITY * M_AIR * myElevation) / (R_GAS * myTemperature));
+}

agrolib/hydrall/hydrall.h

@@ -0,0 +1,31 @@
+#ifndef HYDRALL_H
+#define HYDRALL_H
+
+    #ifndef COMMONCONSTANTS_H
+        #include "commonConstants.h"
+    #endif
+    #ifndef CRIT3DDATE_H
+        #include "crit3dDate.h"
+    #endif
+
+    #define UPSCALINGFUNC(z,LAI) ((1.0 - exp(-(z)*(LAI))) / (z))
+
+    // Tree-plant properties
+    #define FORM   0.5          // stem form factor
+    #define RHOF   0.1          // [KgDM m-3] foliage density
+    #define RHOS   750          // [KgDM m-3] wood-stem density
+
+    // Hydraulic properties
+    #define H50     0.4         // height for 50% maturation of xylem cells (m) [not relevant]
+    #define KR      4.0E-7      // root specific conductance (m3 MPa-1 s-1 kg-1) [not relevant]
+    #define KSMAX   2.5E-3      // max. sapwood specific conductivity (m2 MPa-1 s-1) [not relevant]
+    #define PSITHR  -2.5        // water potential threshold for cavitation (MPa) [not relevant]
+
+
+    #define NOT_INITIALIZED_VINE -1
+
+    double getCO2(Crit3DDate myDate, double myTemperature, double myElevation);
+    double getPressureFromElevation(double myTemperature, double myElevation);
+
+
+#endif // HYDRALL_H

agrolib/hydrall/hydrall.pro

@@ -0,0 +1,36 @@
+#---------------------------------------------------
+#
+#   hydrall library
+#   This project is part of CRITERIA3D distribution
+#
+#---------------------------------------------------
+
+QT       -= core gui
+
+TEMPLATE = lib
+CONFIG += staticlib
+
+CONFIG += debug_and_release
+CONFIG += c++11 c++14 c++17
+
+#DEFINES += _CRT_SECURE_NO_WARNINGS
+
+
+unix:{
+    CONFIG(debug, debug|release) {
+        TARGET = debug/hydrall
+    } else {
+        TARGET = release/hydrall
+    }
+}
+win32:{
+    TARGET = hydrall
+}
+
+INCLUDEPATH += ../crit3dDate ../mathFunctions ../soil ../crop
+
+SOURCES += hydrall.cpp 
+
+
+HEADERS += hydrall.h 
+

agrolib/mathFunctions/commonConstants.h

@@ -159,6 +159,8 @@
     #define	MO2		0.032
     // [kg mol-1] mass of molecular nitrogen (N2)
     #define	MN2		0.028
+    // [kg mol-1] mass of air
+    #define M_AIR    0.029
     // [K] zero Celsius
     #define	ZEROCELSIUS	273.15
     // [] ratio molecular weight of water vapour/dry air