soilchemistry-deconit.h

 
 
#ifndef  LM_SOILCHEMISTRY_DECONIT_H_
#define  LM_SOILCHEMISTRY_DECONIT_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
 
#include  "eventhandler/fertilize/fertilize.h"
#include  "eventhandler/till/till.h"
 
#include  <input/event/events.h>
 
namespace ldndc {
 
class  LDNDC_API  SoilChemistryDECONIT  :  public  MBE_LegacyModel
{
    LMOD_EXPORT_MODULE_INFO(SoilChemistryDECONIT,"soilchemistry:deconit","Soilchemistry DECONIT");
    public:
        SoilChemistryDECONIT(
                MoBiLE_State *,
                cbm::io_kcomm_t *,
                timemode_e);
 
        ~SoilChemistryDECONIT();
 
        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:
        input_class_siteparameters_t const &  sipar_;
        input_class_soillayers_t const &  sl_;
 
        substate_surfacebulk_t &  sb_;
        substate_airchemistry_t &  ac_;
        substate_microclimate_t &  mc_;
        substate_physiology_t &  ph_;
        substate_soilchemistry_t &  sc_;
        substate_watercycle_t &  wc_;
 
        MoBiLE_PlantVegetation *  m_veg;
 
        timemode_e const timemode_;
    
        EventHandlerFertilize  m_eventfertilize;
 
        EventQueue  m_TillEvents;
        CBM_Handle  m_TillHandle;
 
        lvector_t< double > addC_lit1_Frt_sl;
        lvector_t< double > addC_lit2_Frt_sl;
        lvector_t< double > addC_lit3_Frt_sl;
 
        lvector_t< double > addC_lit1_Fol_sl;
        lvector_t< double > addC_lit2_Fol_sl;
        lvector_t< double > addC_lit3_Fol_sl;
 
        lvector_t< double > fact_m;
        lvector_t< double > fact_t;
        lvector_t< double > fact_tm;
 
        lvector_t< double > fact_carbonQuality;
 
        lvector_t< double > fdeniFROST_sl;
        lvector_t< double > fdeni;
        lvector_t< double > sumHeterotrophs;
 
        lvector_t< double > const &mc_temp_sl;
 
        double *waterflux_sl;
        double *accumulated_waterflux_old_sl;
 
        double tot_c_balance_1;
        double tot_n_balance_1;
 
        double K_LITTER_TYPE;                                                        // Decomposition reduction factor for coniferous.    [-]
        double ts_scale;
 
        double DECONIT_Fact_T(size_t);
        double DECONIT_Fact_M(size_t);
 
        double DECONIT_afgen(
                     const double ,
                     const double ,
                     const double ,
                     const double );
 
 
        void DECONIT_GasSpecDiffusion(
                              const int sl,          // Soil layer counter                      (-)
                              double &chem_sl,       // Concentration in the current soil layer (kgNm-2)
                              double &chem_last_sl);   // Concentration in a deeper soil layer    (kgNm-2)
 
        lerr_t DECONIT_step_init();
 
        void DECONIT_AnaerobicVolumeFraction();
 
        void DECONIT_EmpiricalAnaerobicVolumeFraction();
 
        lerr_t DECONIT_fertilize();
        lerr_t DECONIT_manure();
        lerr_t DECONIT_till();
        lerr_t DECONIT_urea_hydrolysis();
 
        void DECONIT_ArahVintenAnaerobicVolumeFraction();
 
        void DECONIT_DenitrifierFraction();
 
        void DECONIT_MicrobialGrowth();
 
        void DECONIT_Denitrification();
 
        void DECONIT_UpwardDiffusion();
 
        void DECONIT_Leaching();
 
        void DECONIT_MicrobialDecay();
 
        void DECONIT_Chemodenitrification();
 
        void DECONIT_CalcRootResp();
 
        void DECONIT_AllocLitter();
 
        void DECONIT_Decomposition();
 
        void DECONIT_Fragmentation();
 
        void DECONIT_UpdateTempMoist();
 
        double DECONIT_get_wfps(
                                size_t /* soil layer */);
 
        lerr_t
        DECONIT_balance_check(
                              unsigned int);
 
        double
        DECONIT_get_c_tot();
 
        double
        DECONIT_get_n_tot();
 
        lerr_t DECONIT_CHECK_FOR_NEGATIVE_VALUE(
                                                char const *  /* position */);
};
 
} /*namespace ldndc*/
 
#endif    /*  LM_SOILCHEMISTRY_DECONIT_H_  */