vocguenther.h

 
#ifndef  LM_PHYSIOLOGY_VOCGUENTHER_H_
#define  LM_PHYSIOLOGY_VOCGUENTHER_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
 
namespace ldndc {
 
class  LDNDC_API  PhysiologyVOCGUENTHER  :  public  MBE_LegacyModel
{
    LMOD_EXPORT_MODULE_INFO(PhysiologyVOCGUENTHER,"physiology:vocguenther","Physiology VOCGUENTHER");
 
    /***  module specific constants  ***/
 
    static double const  TREF;
    static double const  TOPT;
    static double const  ALPHA;
    static double const  BETA;
    static double const  CT1;
    static double const  CT2;
    static double const  CL1;
 
    public:
        PhysiologyVOCGUENTHER(
                MoBiLE_State *,
                                cbm::io_kcomm_t *,
                                timemode_e);
 
        ~PhysiologyVOCGUENTHER();
 
        lerr_t  configure(
                ldndc::config_file_t const *);
 
        lerr_t  initialize();
 
        lerr_t  solve();
 
        lerr_t  finalize() { return  LDNDC_ERR_OK; }
 
        lerr_t  sleep() { return  LDNDC_ERR_OK; }
        lerr_t  wake() { return  LDNDC_ERR_OK; }
 
    private:
        /* required state components */
        substate_microclimate_t *  mc;
        substate_physiology_t *  phys;
        MoBiLE_PlantVegetation *  m_veg;
        input_class_setup_t const *  m_setup;
 
        struct  photosynth_t
        {
            /* aggregated photosynthetic active radiation (umol m-2 s-1) */
            double  par240;
            /* aggregated photosynthetic active radiation (umol m-2 s-1) */
            double  par24;
            /* photosynthetic active radiation (umol m-2 s-1) */
            double  par;
 
            photosynth_t();
        };
        struct  foliage_t
        {
            /* aggregated foliage temperature within a canopy layer (K) */
            double  tempK240;
            /* aggregated foliage temperature within a canopy layer (K) */
            double  tempK24;
            /* foliage temperature within a canopy layer (K) */
            double  tempK;
//            double  ef_isoAct;   // emission factor of isoprene(ug gDW-1 h-1)
//            double  ef_monoAct;  // emission factor of monoterpenes (ug gDW-1 h-1)
//            double  isoprene;    // isoprene emission (ug m-2ground h-1)
//            double  monoterp;    // monoterpene emission (ug m-2ground h-1)
 
            foliage_t();
        };
        struct  enzyme_activity_t
        {
            /* emission factor of isoprene(ug gDW-1 h-1) */
            double  iso;
            /* isoprene emission (ug m-2ground h-1) */
            double  isoprene;
            /* emission factor of monoterpenes (ug gDW-1 h-1) */
            double  mono;
            /* monoterpene emission (ug m-2ground h-1) */
            double  monoterp;
 
            enzyme_activity_t();
        };
        struct  leaf_emission_t
        {
            size_t  foliage_layer;
 
            photosynth_t  pho;
            foliage_t  fol;
            enzyme_activity_t  ef_act;
        };
 
        void  CalcLeafEmission(
                leaf_emission_t *,
                MoBiLE_Plant const &) const;
 
        void  step_init();
};
} /*namespace ldndc*/
 
#endif  /*  !LM_PHYSIOLOGY_VOCGUENTHER_H_  */