ld_icesnow.h

 
#ifndef  LD_ICESNOW_H_
#define  LD_ICESNOW_H_
 
#include  "ld_modelsconfig.h"
#include  <input/siteparameters/siteparameters.h>
#include  <containers/cbm_vector.h>
#include  "state/mbe_state.h"
 
namespace ldndc {
 
class  LDNDC_API  WaterCycleSnowDNDC
{
public:
 
    static const double  TLIMIT;
 
    static const double  RMCOEFF;
 
    static const double  SNALBEDO;
 
public:
    WaterCycleSnowDNDC( cbm::io_kcomm_t *);
    ~WaterCycleSnowDNDC();
 
    double dt;
 
public:
 
    struct SnowPackStateIn
    {
        double  air_temperature_above_canopy;
        double  air_temperature_below_canopy;
 
        double  throughfall_rain_snow;
 
        double  canopy_snow;
        double  surface_snow;
 
        double  canopy_water;
        double  surface_water;
 
        double  shortwave_radiation_above_canopy;
        double  shortwave_radiation_below_canopy;
    };
 
    struct SnowPackStateOut
    {
        double  rain_throughfall;
 
        double  snow_throughfall;
 
        double  canopy_snow_melt;
        double  surface_snow_melt;
 
        double  canopy_snow;
        double  surface_snow;
 
        double  canopy_water;
        double  surface_water;
    };
    lerr_t  SnowPack( SnowPackStateIn const &, SnowPackStateOut &);
 
    struct IceContentStateIn
    {
        /*in*/
        lvector_t<double>  h_sl;
        lvector_t<double>  bulkdensity_sl;
        lvector_t<double>  fcorg_sl;
 
        /*in/out*/
        lvector_t<double>  wc_sl;
        lvector_t<double>  ice_sl;
        lvector_t<double>  soil_temperature_sl;
    };
 
    struct IceContentStateOut
    {
        /*in/out*/
        lvector_t<double>  wc_sl;
        lvector_t<double>  ice_sl;
        lvector_t<double>  soil_temperature_sl;
    };
 
    lerr_t  IceContent( IceContentStateIn const &, IceContentStateOut &);
 
private:
        static double snow_age_top_layer; // for SnowDNDCGetSnowAlbedo but not in use yet
    
    siteparameters::input_class_siteparameters_t const *  m_param;
 
    double
    SnowDNDCGetIceTemperature(
                              double /* water mass [kg] */,
                              double /* ice mass [kg] */);
 
    double
    SnowDNDCGetSoilHeatCapacity(
                                double /* organic carbon fraction [-] */,
                                double /* soil mass [kg] */,
                                double /* water mass [kg] */,
                                double /* ice mass [kg] */);
 
    double
    SnowDNDCGetSnowAlbedo(
            double /* age of the snow in days */,
            double /* amount of snow in m */,
            double /* amount of surface water in m */);
 
    double
    SnowDNDCGetSnowAgeTopLayer(
            double /* amount of snow in the previous timestep [m] */,
            double /* amount of snow now [m] */,
            double /* age of the snow top layer in days */);
};
 
} /* namespace ldndc */
 
#endif  /*  LD_ICESNOW_H_  */