soilchemistry-dndc.h

 
#ifndef  LM_SOILCHEMISTRY_DNDC_H_
#define  LM_SOILCHEMISTRY_DNDC_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
#include  "ld_shared.h"
#include  "ld_isotopes.h"
#include  "soilchemistry/dndc/output-soilchemistry-dndc.h"
 
#include  "soilchemistry/ld_enhanced_efficiency_nitrogen_fertilizers.h"
 
#include  "eventhandler/fertilize/fertilize.h"
#include  "eventhandler/till/till.h"
 
#include  <input/event/events.h>
 
class  solver_tridiagonal_t;
 
namespace ldndc {
class  LDNDC_API  SoilChemistryDNDC  :  public  MBE_LegacyModel
{
    /***  module specific constants  ***/
 
    static const unsigned int  IMAX_S;
    
    static const unsigned int  NN;
    
    static const double  TNORM;
    
    static const double  HLIM;
    
    static const double  FRN2;
 
    LMOD_EXPORT_MODULE_INFO(SoilChemistryDNDC,"soilchemistry:dndc","Soilchemistry DNDC");
    public:
        SoilChemistryDNDC(
                MoBiLE_State *,
                cbm::io_kcomm_t *,
                timemode_e);
 
        ~SoilChemistryDNDC();
 
        lerr_t  configure( ldndc::config_file_t const *);
        lerr_t  initialize();
 
        lerr_t  register_ports( cbm::io_kcomm_t *);
        lerr_t  unregister_ports( cbm::io_kcomm_t *);
 
        lerr_t  solve();
        lerr_t  finalize();
 
        lerr_t  sleep() { return  LDNDC_ERR_OK; }
        lerr_t  wake() { return  LDNDC_ERR_OK; }
 
    private:
        cbm::io_kcomm_t *  io_kcomm;
 
        input_class_groundwater_t const *  gw_;
        input_class_siteparameters_t const &  sipar_;
        input_class_soillayers_t const &  sl_;
 
        substate_surfacebulk_t &  sb_;
        substate_airchemistry_t const &  ac_;
        substate_microclimate_t const &  mc_;
        substate_physiology_t const &  ph_;
        substate_soilchemistry_t &  sc_;
        substate_watercycle_t &  wc_;
 
        MoBiLE_PlantVegetation *  m_veg;
 
        timemode_e const timemode_;
 
        OutputSoilchemistryDNDC  output_writer;
        ldndc_flt64_t *  output_buffer;
 
        EventHandlerFertilize  m_eventfertilize;
        EventHandlerTill  m_eventtill;
 
        CBM_Handle  m_PlantHandle;
        int  m_PlantCounter;
 
    private:
        
        double *  d_soil_eff_;
 
        solver_tridiagonal_t *  tdma_solver_;
 
        double *waterflux_sl;
        double *accumulated_waterflux_old_sl;
 
        lvector_t< double >  crb_l;             // very labile active organic carbon
        lvector_t< double >  crb_r;             // labile active organic carbon
        lvector_t< double >  hdc_l;             // recalcitrant active organic carbon
        lvector_t< double >  hdc_r;             // very recalcitrant active organic carbon
 
        lvector_t< double >  p_o2_sl;           
        lvector_t< double >  fdeni;             
        lvector_t< double >  day_co2;           
        lvector_t< double >  co2_auto_sl;
        lvector_t< double >  anvf_old;          
        lvector_t< double >  nitrify_n2o;       
 
        lvector_t< double >  communicate_sl;
 
        lvector_t< double >  factorA;           
        lvector_t< double >  factorB;           
        lvector_t< double >  flux_chem_no_sl;   
        lvector_t< double >  flux_chem_n2o_sl;  
 
        double  no_0;                           
        double  n2o_0;                          
        double  an_no_0;                        
        double  an_n2o_0;                       
 
        double day_floor_no;                    
        double day_floor_n2o;                   
        double day_floor_an_no;                 
        double day_floor_an_n2o;                
 
        double day_no3_groundwater_access;      
        double year_no3_groundwater_access;     
 
        double n2o_atm;                         
        double no_atm;                          
 
        double day_mineral_litter;
        double day_mineral_humus;
        double day_mineral_microbes;
        double day_mineral_aorg;
        double day_mineral_humads;
 
        lerr_t  SCDNDC_step_initialize_();
        lerr_t  SCDNDC_step_finalize_();
 
 
 
 
        lerr_t
        SCDNDC_write_output();
        lerr_t
        SCDNDC_write_output_daily();
 
 
        void SCDNDC_urea_hydrolysis();
 
        void SCDNDC_autotrophic_respiration();
 
        void SCDNDC_littering();
 
        void SCDNDC_pertubation();
 
         void SCDNDC_distribute_among_litter_layers(
                double  _add_c[], 
                double  _n_diff);
 
        void SCDNDC_nh4_clay_equilibrium();
 
        void SCDNDC_calc_sl_arrays();
    
        void SCDNDC_split_aorg();
    
        void SCDNDC_combine_aorg();
    
        void SCDNDC_litter_decomposition();
        
        void SCDNDC_humus_decomposition();
        
        void SCDNDC_microbial_dynamics();
    
        void SCDNDC_humads_decomposition();
        
        void SCDNDC_nh4_nh3_dynamics();
 
        void SCDNDC_nitrification();
 
        void SCDNDC_denitrification();
 
        void SCDNDC_gas_diffusion();
 
        void SCDNDC_groundwater_access();
 
        void SCDNDC_leaching();
 
        void SCDNDC_o2_diffusion();
 
        void SCDNDC_anvf();
 
        double SCDNDC_temp_mois(
                size_t);
 
        void SCDNDC_decomp_litter(
                size_t sl,
                double rcn,                 
                double decoConstant,        
                double microUse,            
                double temp_mois_fact,      
                double &crb_l,              
                double &rc);                
 
        void SCDNDC_wood_decomposition();
        
        void SCDNDC_nh4_nh3_equilibrium();
 
        void SCDNDC_nh3_volatilization();
 
        void SCDNDC_diffusion(
                size_t,
                double,
                double,
                double &,
                double &);
 
        void SCDNDC_conversion(
                size_t sl,
                double CNEFF,
                double MWN,
                double n_sum,
                double doc_MM,
                double fn,
                double MAINTENANCE,
                double MAXRGR,
                double MAXYLD,
                double phk,
                double f_MM,
                double &deniC,
                lvector_t< double > &  an_doc,
                lvector_t< double > &  day_co2,
                double &chem_from,
                double &chem_to);
 
        void SCDNDC_chemodenitrification(
                size_t sl,
                double fts);
 
        void SCDNDC_emission();
 
        double SCDNDC_get_fact_o2(size_t /* soil layer */);
 
        double SCDNDC_get_wfps(
                               size_t /* soil layer */);
    
        /* balance check switch */
        bool  bc_warn;
        double  bc_tolerance;
 
        lerr_t  SCDNDC_balance_check(
                unsigned int /*stage*/, double * /*carbon*/, double * /*nitrogen*/);
        double SCDNDC_get_c_tot();
        double SCDNDC_get_n_tot();
        
        Isotopes  isotope;
};
 
} /*namespace ldndc*/
 
#endif    /*  LM_SOILCHEMISTRY_DNDC_H_  */