arable-dndc.h

 
#ifndef  LM_PHYSIOLOGY_ARABLE_DNDC_H_
#define  LM_PHYSIOLOGY_ARABLE_DNDC_H_
 
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
#include  "ld_shared.h"
 
#include  "physiology/ld_plantfunctions.h"
#include  "eventhandler/graze/graze.h"
#include  "eventhandler/cut/cut.h"
 
#include  <containers/lgrowarray.h>
#include  <utils/cbm_utils.h>
 
namespace ldndc {
 
    struct BaseRootSystemDNDC;
    class  LDNDC_API  PhysiologyArableDNDC  :  public  MBE_LegacyModel
    {
        LMOD_EXPORT_MODULE_INFO(PhysiologyArableDNDC,"physiology:arabledndc","ArableDNDC");
 
        /***  module specific constants  ***/
 
        static double const  GROWTH_RESPIRATION_FRACTION;
 
        static const double  INITIAL_BIOMASS;
 
        static const double  S_TLIMIT;
 
    public:
        PhysiologyArableDNDC( MoBiLE_State *,
                             cbm::io_kcomm_t *, timemode_e);
        ~PhysiologyArableDNDC();
 
        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() { return  LDNDC_ERR_OK; }
        lerr_t  wake() { return  LDNDC_ERR_OK; }
 
    private:
        MoBiLE_State *  m_state;
        cbm::io_kcomm_t *  io_kcomm;
 
        /* required input classes */
        input_class_setup_t const &  m_setup;
        input_class_soillayers_t const &  sl_;
        input_class_siteparameters_t const * m_sipar;
        input_class_species_t const *  m_species;
        input_class_speciesparameters_t const * m_speciesparameters;
 
        /* required state components */
        substate_airchemistry_t const &  ac_;
        substate_microclimate_t &  mc_;
        substate_physiology_t &  ph_;
        substate_soilchemistry_t &  sc_;
        substate_watercycle_t &  wc_;
 
        MoBiLE_PlantVegetation *  m_veg;
 
        SubscribedEvent<LD_EventHandlerQueue>  m_HarvestEvents;
        SubscribedEvent<LD_EventHandlerQueue>  m_PlantEvents;
 
        ldndc::growable_array < BaseRootSystemDNDC*, 1000, 1 >  root_system;
 
        double const *mc_temp;
 
        double const FTS_TOT_;
 
        LD_PlantFunctions  m_pf;
 
        double tot_c_;
        double tot_n_;
 
        lvector_t< double >  n_demand_vt_;
        lvector_t< double >  n_uptake_vt_;
        lvector_t< double >  n_avail_tot_vt_;
        lvector_t< double >  n_plant_vt_;
 
        lvector_t< double >  dw_fol_old_vt;
        lvector_t< double >  dw_frt_old_vt;
        lvector_t< double >  dw_lst_old_vt;
        lvector_t< double >  dw_bud_old_vt;
 
        lvector_t< double >  chill_factor;
 
        lvector_t< double >  chill_units;
 
        lvector_t< unsigned int >  root_q_vt_;
 
        lvector_t< double > temp_sl;
 
        /***  methods  ***/
 
        lerr_t  ArableDNDC_step_init();
        lerr_t  ArableDNDC_step_out();
        lerr_t  ArableDNDC_resize();
 
        void  update_species_related_state_items( MoBiLE_Plant *);
 
        lerr_t  ArableDNDC_event_plant( MoBiLE_Plant *, EventAttributes const &);
 
        lerr_t  event_harvest( MoBiLE_Plant *, EventAttributes const &);
 
        lerr_t  development_stage( MoBiLE_Plant *);
 
        void  growth_respiration( MoBiLE_Plant *, double);
 
        void  maintenance_respiration( MoBiLE_Plant *);
 
        double  nitrogen_demand( MoBiLE_Plant *);
 
        void  nitrogen_uptake();
 
        double  nitrogen_fixation( MoBiLE_Plant *, double, double);
        double  total_fixed_nitrogen( MoBiLE_Plant const *, double) const;
 
        void  biomass_partition( MoBiLE_Plant *, double);
 
        void  senescence( MoBiLE_Plant *);
 
        /* coverage */
        void  update_canopy_coverage( MoBiLE_Plant *);
 
        /* plant height */
        void  update_canopy_height( MoBiLE_Plant *);
 
        /* canopy structure */
        void  update_canopy_items( MoBiLE_Plant *);
 
        /* C & N balance checker */
        void  balance_check( int /*stage*/);
    };
} /*namespace ldndc*/
 
#endif  /*  !LM_PHYSIOLOGY_ARABLE_DNDC_H_  */