canopyecm.h

 
#ifndef  LM_MICROCLIMATE_CANOPY_ECM_H_
#define  LM_MICROCLIMATE_CANOPY_ECM_H_
 
#include  "mbe_legacymodel.h"
#include  "state/mbe_state.h"
 
namespace ldndc {
 
 
class  LDNDC_API  MicroClimateCanopyECM  :  public  MBE_LegacyModel
{
    LMOD_EXPORT_MODULE_INFO(MicroClimateCanopyECM,"microclimate:canopyecm","Microclimate CanopyECM");
    
    static const double  QMASSAREA;
    static const double  EXPW;
    static const double  ROUGH;
    static const double  FEXT_W;
 
    static const double  FLITALB;
    
    static const double  HLMIN;
    static const double  HBASE;
    static const double  SDAMP;
 
public:
    MicroClimateCanopyECM(
                          MoBiLE_State *,
                          cbm::io_kcomm_t *,
                          timemode_e);
    
    ~MicroClimateCanopyECM();
    
    
    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:
    input_class_setup_t const &  se_;
    substate_watercycle_t const &  wc_;
    substate_microclimate_t &  mc_;
    substate_physiology_t const &  ph_;
    
    MoBiLE_PlantVegetation *  m_veg;
    
    input_class_climate_t const &  climate_;
    input_class_soillayers_t const &  sl_;
    substate_soilchemistry_t const &  sc_;
    
    substate_surfacebulk_t &  sbl_;
    
    solver_tridiagonal_t  tdma_solver_;
    
    /* pointer to daily or subdaily state items */
    /* radiation per canopy layer (temporarly unit-converted) [W m-2] */
    lvector_t< double >  loc_lw_rad_down_fl;
    lvector_t< double >  loc_lw_rad_up_fl;
    lvector_t< double >  loc_sw_rad_fl;
    lvector_t< double >  loc_sw_rad_refl_fl;
    double loc_sw_rad_refl_a;
    lvector_t< double >  f_area_fl;
    lvector_t< double >  lai_fl;
    lvector_t< double >  loc_temp_fl;
    lvector_t< double >  loc_fsun_fl;
    lvector_t< double >  loc_parshd_fl;
    lvector_t< double >  loc_parsun_fl;
    lvector_t< double >  loc_vpd_fl;
    lvector_t< double >  loc_rh_fl;
    lvector_t< double >  loc_win_fl;
 
    bool al_new;                        // if additional layers above soil occur the first time ala_new is true, else false;
    
    double *  temp_old_sl;              // old soil temperature for soil layers
    double *  temp_old_sbl;             // old soil temperature for bulk surface layers (snow or water)
    
    double *  hcap_sl;                  // volumetric heat capacity of the soil layer (J K-1 m-2)
    double *  lambda_sl;                // heat conductivity (J K-1 m-1 s-1)
    
    double *  h_extended_sl;            // extends the last soil layer for more realistic lower boundary condition of constant annual temp
    double *  som_extended_sl;
    double *  min_extended_sl;
    
    //  depth of soil layer (m)
    double *  dsl;
    // volumetric heat capacity of the soil layer (J K-1 m-2)
    double *  hcap;
    // heat conductivity (J K-1 m-1 s-1)
    double *  lambda;
    // soil temperature to start with (oC)
    double *  tempOld_sl;
    
    bool dirichletboundarycondition;
    double annual_average_temperature;
    
    double domainextension;
    bool energybalance;
    
    cbm::string_t soiltemperature_method;
    
    double accumulated_evapotranspiration_canopy;
    double accumulated_evapotranspiration_soil;
    
    /* outgoing (reflected) short wave radiation [W m-2] */
    double loc_sw_rad_out;
 
    /* outgoing long wave radiation [W m-2] */
    double loc_lw_rad_out;
 
    /* outgoing sensible heat [W m-2] */
    double loc_sensible_heat_out_top;
    double loc_sensible_heat_out_bottom;
    
    /* radiation above litter [W m-2] */
    double loc_rad_a;
    double loc_temp_a;
    double loc_rad_sw;
    double loc_rad_lw;
    double loc_temp;
    double loc_win;
    double loc_vpd;
    double loc_fdif;
    double loc_sinb;
 
    double  lai_sum;
    
    // parameter aggregation
    double  fdrag;
    
    
    lerr_t  CanopyECM_step_init();
    lerr_t  CanopyECM_step_out();
        
    void  CalcTemperature();
 
    lerr_t  CalcSoilTemperature_implicite(
                                 lvector_t< double > & /*soil temperatures*/,
                                 double /*average temperature above canopy*/,
                                 double &/*average temperature above litter*/);
 
    lerr_t  CalcSoilTemperature_implicite_extended(
                                 lvector_t< double > & /*soil temperatures*/,
                                 double /*average temperature above canopy*/,
                                 double &/*average temperature above litter*/);
    
    lerr_t  CalcSoilTemperature_explicite(
                                 lvector_t< double > & /*soil temperatures*/,
                                 double /*average temperature above litter*/);
    
    
    void  CalcShortwaveRadiation();
    
    
    void  CalcHumidity();
    
    void  CalcWindspeed();
    
    void  CalcAlbedo();
 
    double  foliage_layer_heat_capacity( size_t /* foliage layer*/);
    
    double  get_canopy_energy();
 
    double  get_surface_bulk_energy();
    
    double  get_soil_energy( lvector_t< double > &  /* soil temperature */);
    
    double get_lai_fraction( size_t /* foliage layer */,
                             size_t /* no of foliage layer */);
    
    double  get_sky_view( double const & /* lai*/);
    
    double get_lw_radiation_down( double const & /* open area fraction top */,
                                  double const & /* open area fraction */,
                                  double const & /* leaf area index */,
                                  double const & /* incoming longwave radiation */);
};
} /*namespace ldndc*/
 
#endif  /*  !LDNDC_MOD_MICROCLIMATE_CANOPY_ECM_H_  */