psim.h

 
#ifndef  LM_PHYSIOLOGY_GROWTHPSIM_H_
#define  LM_PHYSIOLOGY_GROWTHPSIM_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
#include  "ld_shared.h"
 
#include  "physiology/ld_nitrogenfixation.h"
#include  "physiology/ld_plantfunctions.h"
#include  "physiology/ld_rootsystem.h"
#include  "physiology/ld_treedyn.h"
#include  "physiology/photofarquhar/berryball.h"
 
#include  "eventhandler/graze/graze.h"
#include  "eventhandler/cut/cut.h"
 
namespace ldndc {
 
extern cbm::logger_t *  PhysiologyGrowthPSIMLogger;

class  LDNDC_API  PhysiologyGrowthPSIM  :  public  MBE_LegacyModel
{
    LMOD_EXPORT_MODULE_INFO(PhysiologyGrowthPSIM,"physiology:psim","Physiology GrowthPSIM");
 
    public:
        static const unsigned int  IMAX_W;
        static const double  C1;
        static const double  CCOMNOX;
        static const double  CSATNH3;
        static const double  CSATNOX;
        static const double  EPOTNOX;
        static const double  FRFRT;
        static const double  FMIN;
        static const double  KMMM;
        static const double  PAMM;
        static const double  PGDD;
        static const double  PMIN;
        static const double  PNIT;
        static const double  PPHLOE;
        static const double  PREDFRT;
        static const double  PREDSAP;
        static const double  PTW;
        static const double  TAU;
        static const double  TGLIM;
        static const double  TRMAX;
        static const double  TRMIN;
        static const double  TROPT;
        static const double  UPOTNH3;
 
    public:
        PhysiologyGrowthPSIM(
                             MoBiLE_State *,
                             ::cbm::io_kcomm_t *,
                             timemode_e);
 
        ~PhysiologyGrowthPSIM();
 
        lerr_t  configure( ::ldndc::config_file_t const *);
        lerr_t  register_ports( cbm::io_kcomm_t *);
        lerr_t  initialize();
 
        lerr_t  solve();
        lerr_t  unregister_ports( cbm::io_kcomm_t *);
        lerr_t  finalize() { return  LDNDC_ERR_OK; }
 
        lerr_t  sleep();
        lerr_t  wake();

        cbm::string_t branchfraction_opt;

        cbm::string_t crownlength_opt;

        bool competition_opt;
        bool rootgrowth_opt;
 
        double clumping_factor;

        cbm::string_t timemode_opt;

        int stand_structure_opt;
        double target_biomass;
        int target_treenumber;
 
        enum canopy_nitrogen_uptake_e
        {
            NOUPTAKE,
            STOMATALUPTAKE,
            DIRECTUPTAKE
        };
 
        canopy_nitrogen_uptake_e n_uptake_canopy_method;
 
        bool have_forest_structure();
 
    protected:
    
        lerr_t PSIM_ResizeVegetationVectors();

        lerr_t PSIM_StepInit();

        lerr_t PSIM_PlantingEvent();

        lerr_t PSIM_HydraulicConductance();

        lerr_t PSIM_Photosynthesis();

        lerr_t  PSIM_Potentialtranspiration();
    
        lerr_t  PSIM_ResetFlushing();
        
        lerr_t  PSIM_ResetPhenology( MoBiLE_Plant *);
 
        lerr_t  PSIM_Reset( MoBiLE_Plant *);
 
        lerr_t  PSIM_ResetBiomassFromLastYear();
 
        lerr_t  PSIM_UpdateNitrogenConc();
 
        lerr_t  PSIM_AgriculturalManagement();
 
        lerr_t  PSIM_BudBurst();
 
        lerr_t  PSIM_Senescence();
 
        lerr_t  PSIM_CarAllocation();
 
        lerr_t  PSIM_SoilCarbonRelease();
 
        lerr_t  PSIM_NitAllocation();
 
        lerr_t  PSIM_Respiration();
 
        lerr_t  PSIM_NitrogenFixation();
 
        lerr_t  PSIM_NitrogenUptake();
 
        lerr_t  PSIM_BiomassUpdate();
 
        lerr_t  PSIM_PhotosynthesisRates();
 
        lerr_t  PSIM_DroughtStress();
 
        lerr_t  PSIM_NitrogenBalance();
 
        lerr_t  PSIM_VegStructure();
 
        lerr_t  PSIM_BalanceCheck(size_t);
 
        bool  subdaily_timemode();
 
        lerr_t PSIM_Harvest();
 
        lerr_t PSIM_StepOut();
 
    private:
 
        MoBiLE_State *  m_state;
        cbm::io_kcomm_t *  io_kcomm;
 
        input_class_climate_t const &  cl_;
        input_class_soillayers_t const & sl_;
        input_class_setup_t const *  m_setup;
        input_class_siteparameters_t const * m_sipar;
        input_class_species_t const *  m_species;
 
        substate_airchemistry_t &  ac_;
        substate_microclimate_t &  mc_;
        substate_physiology_t &  ph_;
        substate_soilchemistry_t &  sc_;
        substate_watercycle_t &  wc_;
 
        MoBiLE_PlantVegetation *  m_veg;
 
        NitrogenFixation *NitrogenFixation_;
 
        SubscribedEvent<LD_EventHandlerQueue>  m_PlantEvents;
        SubscribedEvent<LD_EventHandlerQueue>  m_HarvestEvents;
    
        lvector_t< bool > additional_season;              //set to true if additional growing season within the current phenological year is possible
        lvector_t< unsigned int > phenological_year;      //stores year of start of current growing season
        lvector_t< int > foliage_age_of_start_senescence; // age of foliage in days at which leaves start to shed
 
        lvector_t< double > carbonuptake_vt;
 
        lvector_t< int > days_since_emergence;
        std::vector< std::vector< int > > foliage_age_na;
 
        lvector_t< double > temp_sl;
        lvector_t< double > ts_temp_fl;
 
        lvector_t< double > no_concentration_fl;
        lvector_t< double > no2_concentration_fl;
        lvector_t< double > nh3_concentration_fl;
 
        //lvector_t< double > xylem_resistance_smoothed_vt;
        lvector_t< double > xylem_resistance_smoothed_old_vt;
        lvector_t< double > psidecline_cum_vt;
        lvector_t< double > psidecline_daycum_vt;
        lvector_t< double > rehydrationindex_cum_vt;
 
        lvector_t< double > mSapOld_vt;
        lvector_t< double > mCorOld_vt;
 
        lvector_t< double > qsaparea_vt;
        lvector_t< double > fsub_vt;
 
        lvector_t< double > mSapOpt_vt;
        lvector_t< double > mBudLoss_vt;
        lvector_t< double > mSapLoss_vt;
        lvector_t< double > ncFolOpt;
        lvector_t< double > ncBudOpt;
        lvector_t< double > ts_leaftemp_vt;
        lvector_t< double > nd_leaftemp_vt;
        lvector_t< double > nDem_vt;
        lvector_t< double > nFol_vt;
        lvector_t< double > nFrt_vt;
        lvector_t< double > nSap_vt;
        lvector_t< double > nCor_vt;
        lvector_t< double > nBud_vt;
        lvector_t< double > cNOx_fl;
 
        lvector_t< double > uptNH4Max_vt;
        lvector_t< double > uptNO3Max_vt;
        lvector_t< double > uptDONMax_vt;
        lvector_t< double > uptNWetMax_vt;
        lvector_t< double > uptNWet_vt;
        lvector_t< double > uptNTot_vt;
        std::vector< std::vector< double > > uptNH4Max_vtsl;
        std::vector< std::vector< double > > uptNO3Max_vtsl;
        std::vector< std::vector< double > > uptDONMax_vtsl;
 
        lvector_t< double > rFolOld_vt;
        lvector_t< double > rBudOld_vt;
        lvector_t< double > rSapOld_vt;
        lvector_t< double > rFrtOld_vt;
        lvector_t< double > rTraOld_vt;
        lvector_t< double > exsuLossOld_vt;
 
        lvector_t< double > uptNH4_vt;
        lvector_t< double > uptNO3_vt;
        lvector_t< double > uptNH3_vt;
        lvector_t< double > uptDON_vt;
 
        lvector_t< double > uptNH4Old_vt;
        lvector_t< double > uptNO3Old_vt;
        lvector_t< double > uptNH3Old_vt;
        lvector_t< double > uptDONOld_vt;
 
        lvector_t< double > uptNOxOld_vt;
        lvector_t< double > uptNHyOld_vt;
 
        lvector_t< double > n_bud_to_fol_vt;
 
        lvector_t< double > n_sap_to_cor_vt;
        lvector_t< double > n_bud_to_cor_vt;
        lvector_t< double > n_fol_to_cor_vt;
        lvector_t< double > n_frt_to_cor_vt;
 
        lvector_t< double > n_sap_to_litter_vt;
    
        lvector_t< double > n_sap_to_redistribute_vt;
        lvector_t< double > n_bud_to_redistribute_vt;
        lvector_t< double > n_fol_to_redistribute_vt;
        lvector_t< double > n_frt_to_redistribute_vt;
 
        lvector_t< double > n_sap_gain_vt;
        lvector_t< double > n_bud_gain_vt;
        lvector_t< double > n_fol_gain_vt;
        lvector_t< double > n_frt_gain_vt;
 
        lvector_t< double > nhy_uptake_fl;
        lvector_t< double > nox_uptake_fl;
 
        double mFol_emergence;
        double accumulated_wateruptake_old;
        double accumulated_potentialtranspiration_old;
 
        double plant_waterdeficit_old;
        double psi_deficit_old;
 
        double  tot_c;
        double  tot_n;
 
        bool daytime;
        double SIM_DAY_FRACTION;
 
        size_t  fl_cnt_;
 
        LD_PlantFunctions  m_pf;
        TreeDynamics  m_treedyn;
 
        lvector_t< double > no3_sl;
        lvector_t< double > nh4_sl;
        lvector_t< double > don_sl;
 
        EventHandlerGraze  m_eventgraze;
        EventHandlerCut  m_eventcut;
 
        BerryBall m_photo;

        std::map< int, RootSystemDNDC > root_system;
};
} /*namespace ldndc*/
 
#endif  /*  !LM_PHYSIOLOGY_GROWTHPSIM_H_  */