leaf_gbw_nikolov.h

Go to the documentation of this file.

#ifndef LEAF_GBW_NIKOLOV_H
#define LEAF_GBW_NIKOLOV_H
 
#include "../framework/module.h"
#include "../framework/state_map.h"
#include "conductance_helpers.h"         // for g_to_mass
#include "water_and_air_properties.h"    // for saturation_vapor_pressure, molar_volume
#include "boundary_layer_conductance.h"  // for leaf_boundary_layer_conductance_nikolov
 
namespace standardBML
{
class leaf_gbw_nikolov : public direct_module
{
   public:
    leaf_gbw_nikolov(
        state_map const& input_quantities,
        state_map* output_quantities)
        : direct_module{},
 
          // Get references to input quantities
          windspeed{get_input(input_quantities, "windspeed")},
          leafwidth{get_input(input_quantities, "leafwidth")},
          air_temperature{get_input(input_quantities, "temp")},
          leaf_temperature{get_input(input_quantities, "leaf_temperature")},
          Gs{get_input(input_quantities, "Gs")},
          rh{get_input(input_quantities, "rh")},
          air_pressure{get_input(input_quantities, "air_pressure")},
 
          // Get pointers to output quantities
          gbw_leaf_op{get_op(output_quantities, "gbw_leaf")}
    {
    }
    static string_vector get_inputs();
    static string_vector get_outputs();
    static std::string get_name() { return "leaf_gbw_nikolov"; }
 
   private:
    // References to input quantities
    double const& windspeed;
    double const& leafwidth;
    double const& air_temperature;
    double const& leaf_temperature;
    double const& Gs;
    double const& rh;
    double const& air_pressure;
 
    // Pointers to output quantities
    double* gbw_leaf_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector leaf_gbw_nikolov::get_inputs()
{
    return {
        "windspeed",         // m / s
        "leafwidth",         // m
        "temp",              // degrees C
        "leaf_temperature",  // degrees C
        "Gs",                // mol / m^2 / s
        "rh",                // dimensionless
        "air_pressure"       // Pa
    };
}
 
string_vector leaf_gbw_nikolov::get_outputs()
{
    return {
        "gbw_leaf"  // m / s
    };
}
 
void leaf_gbw_nikolov::do_operation() const
{
    // Determine the partial pressure of water vapor from the air temperature
    // and relative humidity
    const double SWVP = saturation_vapor_pressure(air_temperature);  // Pa
    const double water_vapor_pressure = rh * SWVP;                   // Pa
 
    // Calculate the boundary layer conductance. Here we need to convert
    // stomatal conductance from mol / m^2 / s to m / s.
    const double gbw_leaf = leaf_boundary_layer_conductance_nikolov(
        air_temperature,
        leaf_temperature - air_temperature,
        water_vapor_pressure,
        g_to_mass(air_pressure, Gs, air_temperature),
        leafwidth,
        windspeed,
        air_pressure);  // m / s
 
    // Update the output quantity list
    update(gbw_leaf_op, gbw_leaf);  // m / s
}
 
}  // namespace standardBML
#endif