vocmegan.h

 
#ifndef  LM_PHYSIOLOGY_VOCMEGAN_H_
#define  LM_PHYSIOLOGY_VOCMEGAN_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
 
namespace ldndc {
 
class  LDNDC_API  PhysiologyVOCMEGAN  :  public  MBE_LegacyModel
{
    LMOD_EXPORT_MODULE_INFO(PhysiologyVOCMEGAN,"physiology:vocmegan","Physiology VOCMEGAN");
 
    /***  module specific constants  ***/
 
    static double const  TREF;
    static double const  BETA;
    static double const  TOPT;
    static double const  CT1;
    static double const  CT2;
 
    public:
        PhysiologyVOCMEGAN(
                MoBiLE_State *,
                                cbm::io_kcomm_t *,
                                timemode_e);
 
        ~PhysiologyVOCMEGAN();
 
        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_airchemistry_t *  ac;
        substate_microclimate_t *  mc;
        substate_physiology_t *  phys;
        substate_watercycle_t *  wc;
 
        input_class_setup_t const *  m_setup;
        MoBiLE_PlantVegetation *  m_veg;
 
        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;
 
            foliage_t();
        };
 
        struct  leaf_emission_t
        {
            size_t  foliage_layer;
            
            double  isoprene;       // leaf level isoprene emission (ug gDW-1 h-1)
            double  monoterpene;    // leaf level monoterpene emission (ug gDW-1 h-1) */
            double  monoterpene_s;  // leaf level monoterpene emission from storages (ug gDW-1 h-1) */
            double  ovoc;           // leaf level ovoc emission (only from storages) (ug gDW-1 h-1) */
 
            photosynth_t  pho;
            foliage_t  fol;
        };
 
        lerr_t  CalcLeafEmission(
                leaf_emission_t *,
                MoBiLE_Plant const & /*species*/);
 
        void  step_init();
};
} /*namespace ldndc*/
 
#endif  /*  !LM_PHYSIOLOGY_VOCMEGAN_H_  */