ld_plantfunctions.h

 
#ifndef  LD_PLANTFUNCTIONS_H_
#define  LD_PLANTFUNCTIONS_H_
 
#include  "ld_modelsconfig.h"
 
#include  "mbe_plant.h"
#include  "state/mbe_state.h"
 
#include  <input/climate/climate.h>
#include  <input/soillayers/soillayers.h>
#include  <input/setup/setup.h>
#include  <input/species/species.h>
#include  "physiology/ld_allometry.h"
 
#include  <containers/lbstack.h>
 
 
namespace ldndc {
 
class LDNDC_API LD_PlantFunctions
{
public:
    LD_PlantFunctions(
                      MoBiLE_State *,
                      cbm::io_kcomm_t *);
    ~LD_PlantFunctions();
 
    lerr_t  initialize_crop( MoBiLE_Plant *, species::crop_properties_t const *);
    lerr_t  initialize_grass( MoBiLE_Plant *, species::grass_properties_t const *);
    lerr_t  initialize_tree( MoBiLE_Plant *, species::wood_properties_t const *,
                             //cbm::string_t, /* branchfraction */
                             cbm::string_t /* crownlength */,
                             bool /* competition */);
 
protected:
 
    climate::input_class_climate_t const &  cl_;
    setup::input_class_setup_t const &  se_;
    soillayers::input_class_soillayers_t const &  sl_;
 
    substate_microclimate_t &  mc_;
    substate_physiology_t &  ph_;
    substate_soilchemistry_t &  sc_;
 
    LD_Allometry m_alm;
    MoBiLE_PlantVegetation *  m_veg;
 
    /* ##############   WOOD   ############## */
    lerr_t  set_wood_defaults_( MoBiLE_Plant *, species::wood_properties_t const *);
    lerr_t  gapfill_wood_defaults_( MoBiLE_Plant *,
                                    species::wood_properties_t const *,
                                    bool /* competition */);
    lerr_t  initialize_vegetationstructure_wood_( MoBiLE_Plant *,
                                                  /* cbm::string_t, // branchfraction*/
                                                  cbm::string_t /* crownlength */);
 
    /* ##############   GRASS  ############## */
    lerr_t  set_grass_defaults_( MoBiLE_Plant *, species::grass_properties_t const *);
    lerr_t  gapfill_grass_defaults_( MoBiLE_Plant *, species::species_properties_t const *);
    lerr_t  reset_grass_properties_( MoBiLE_Plant *);

    lerr_t  initialize_vegetationstructure_grass_( MoBiLE_Plant *);
 
    /* ##############   CROP   ############## */
    lerr_t  set_crop_defaults_( MoBiLE_Plant *, species::crop_properties_t const *);
    lerr_t  reset_crop_properties_( MoBiLE_Plant *);
    
 
 

    lerr_t  initialize_nitrogen_concentrations( MoBiLE_Plant *);

    lerr_t  canopy_microclimate();
 
    
public:
 
    void
    set_rooting_depth( input_class_soillayers_t const *  m_soillayers,
                       substate_soilchemistry_t const *  soilchem,
                       MoBiLE_Plant *,
                       bool consider_texture = false);
    
    double
    get_optimum_foliage_biomass_trees( MoBiLE_Plant * /* plant species */);

    double
    relative_crown_cover_area( MoBiLE_Plant * /* plant species */);
 
    
    double
    get_optimum_foliage_biomass_grass( MoBiLE_Plant * /* plant species */);
 
 
    int
    get_foliage_age(
                    int _day_of_year_start,
                    int _doy,
                    int _nd_age_classes);
 
    double
    get_dvs_mortality_of_age_class(int /* day of start senescence */,
                                   int /* number of days of senescence */,
                                   int /* day of year */);
 
    double
    get_dvs_mortality_contribution_of_age_class(
                                                MoBiLE_Plant *,
                                                int /* age class */,
                                                double /* mortality of age class */);

    lerr_t
    initialize_seasonality(
                           MoBiLE_Plant *,
                           unsigned int /* day of year start */,
                           unsigned int /* current day of year */,
                           unsigned int &  /* bud state */);

    lerr_t
    initialize_new_foliage_biomass(
                                   MoBiLE_Plant *,
                                   double /*initial biomass*/,
                                   unsigned int /* day of year start */,
                                   unsigned int /* current day of year */,
                                   unsigned int /* bud state */);

    lerr_t
    initialize_new_biomass(
                           MoBiLE_Plant *,
                           double /*initial biomass*/);

    double
    update_dvs_flush(
                     MoBiLE_Plant *  vt,
                     int _days_since_emergence);
    
    lerr_t
    initialize_wood_biomass(
                            MoBiLE_Plant *,
                            double /*initial biomass*/);
 
    lerr_t
    initialize_foliage_dependent_fine_root_biomass(
                                                   MoBiLE_Plant *,
                                                   double /*fom*/);
 
    lerr_t
    update_foliage_layer_biomass_and_lai(
                                         MoBiLE_Plant *,
                                         size_t, /*foliage layer count*/
                                         double /*minimum allowed lai*/);

    double
    aboveground_wood_biomass_from_structure(
                                            MoBiLE_Plant * /* species */);
    
    lerr_t
    update_vertical_biomass_distribution(
                                         MoBiLE_Plant *);
    
    double
    leaf_temperature_(
                      lvector_t< double > /* canopy air temperature */,
                      MoBiLE_Plant * /* species */);
 
    double
    nitrogen_uptake_temperature_dependency(
                                           double /* minimum temperature */,
                                           double /* optimum temperature */,
                                           double /* temperature */);

    double
    optimum_nitrogen_content_foliage(
                                     MoBiLE_Plant * /* plant species */);

    double
    optimum_nitrogen_content_buds(
                                  MoBiLE_Plant * /* plant species */);

    double
    get_foliage_age_of_senescence(
                                  MoBiLE_Plant * /* species */,
                                  double /* fraction */);
 
    double
    accumulated_n_uptake();
};
 

void
estimate_single_layer_biomass(
                              double const _ps,
                              double const _h_bottom,
                              double const _h_top,
                              double  _lw,
                              double  &_fract,
                              double  _h_cum_bottom,
                              double  _h_cum_top);
 
void
fineroots_biomass_distribution(
                               size_t _sl_cnt_litter,
                               size_t _sl_cnt_soil,
                               double _ps,
                               double _rooting_depth,
                               double *_f_frt_sl,
                               double const *_h_sl);
 
} /* namespace ldndc */
 
 
#endif  /*  !LD_PLANTFUNCTIONS_H_   */