16#ifndef __GenericClimateModel_HPP
17#define __GenericClimateModel_HPP
43 const double b = 17.2694;
44 const double T2 = 35.86;
45 return e0*exp(b*T/(T+
T0-T2));
49 const double b = 17.2694;
50 const double T2 = 35.86;
51 return e0*exp(b*T/(T+
T0-T2))*(b/(T+
T0-T2) - b*T/pow(T+
T0-T2,2));
65 const double eps = 0.622;
71 const double eps = 0.622;
74 return eps*e_dT/(
P-e) +
eps*e*e_dT/pow(
P-e,2);
106 virtual MoFEMErrorCode
set(
double t = 0) = 0;
constexpr double t
plate stiffness
double calculateVapourPressure_dT(double T)
double calculateAbsoluteVirtualTempertaure_dT(double T, double Td, double P)
virtual MoFEMErrorCode set(double t=0)=0
double calculateVapourPressureClausiusClapeyron(double T)
Clausius-Clapeyron equation.
double calculateAbsoluteVirtualTempertaure(double T, double Td, double P)
double calculateVapourPressure(double T)
double calculateMixingRatio(double T, double P)
double calculateVapourPressureTetenFormula_dT(double T)
double calculateVapourPressureTetenFormula(double T)
double calculateMixingRatio_dT(double T, double P)