doc-dev/html/DataClasses_8h_source.html

/*

 *  SNOWPACK stand-alone

 *

 *  Copyright WSL Institute for Snow and Avalanche Research SLF, DAVOS, SWITZERLAND

*/

/*  This file is part of Snowpack.

    Snowpack is free software: you can redistribute it and/or modify

    it under the terms of the GNU General Public License as published by

    the Free Software Foundation, either version 3 of the License, or

    (at your option) any later version.


    Snowpack is distributed in the hope that it will be useful,

    but WITHOUT ANY WARRANTY; without even the implied warranty of

    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

    GNU General Public License for more details.


    You should have received a copy of the GNU General Public License

    along with Snowpack.  If not, see <http://www.gnu.org/licenses/>.

*/

#ifndef DATACLASSES_H

#define DATACLASSES_H


#include <snowpack/SnowpackConfig.h>

#include <snowpack/vanGenuchten.h>

#include <snowpack/snowpackCore/SeaIce.h>


#include <snowpack/Constants.h>

#include <meteoio/MeteoIO.h>


#include <string>

#include <vector>


enum {

        SOLID,

        LIQUID,

        GAS,

        NUMBER_OF_PHASES

};


struct SNOW_OPTIC {

        double ggg;

        double exteff;

        double ssa;

};


struct WL_STRUCT {

        double nm;

        double perc;

};


class ZwischenData {

        public:

                ZwischenData(): hoar24(48, 0.0), drift24(48, 0.0), hn3(144, 0.0), hn24(144, 0.0) {}

                void reset();


                friend std::ostream& operator<<(std::ostream& os, const ZwischenData& data);

                friend std::istream& operator>>(std::istream& is, ZwischenData& data);


                std::vector<double> hoar24;

                std::vector<double> drift24;

                std::vector<double> hn3;

                std::vector<double> hn24;

};


class CurrentMeteo {

        public:

                CurrentMeteo();

                CurrentMeteo(const SnowpackConfig& i_cfg);


                void reset(const SnowpackConfig& i_cfg);

                void setMeasTempParameters(const mio::MeteoData& md);

                size_t getNumberMeasTemperatures() const;

                size_t getNumberFixedRates() const;

                size_t getMaxNumberMeasTemperatures() const;

                void getFixedPositions(std::vector<double>& positions) const;

                size_t getNumberFixedPositions() const;

                void copySnowTemperatures(const mio::MeteoData& md, const unsigned int& current_slope);

                void copySolutes(const mio::MeteoData& md, const size_t& i_number_of_solutes);


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const CurrentMeteo& data);

                friend std::istream& operator>>(std::istream& is, CurrentMeteo& data);


                mio::Date date;

                double ta;

                double rh;

                double rh_avg;

                double vw;

                double vw_avg;

                double vw_max;

                double dw;

                double vw_drift;

                double dw_drift;

                double ustar;

                double z0;

                double psi_s;

                double iswr;

                double rswr;

                double mAlbedo;

                double diff;

                double dir_h;

                double elev;

                double ea;

                double lw_net;

                double tss;

                double tss_a12h;

                double tss_a24h;

                double ts0;

                double psum;

                double psum_ph;

                double psum_tech;

                double hs;

                double hs_a3h;

                double hs_rate;

                double geo_heat;

                double adv_heat;


                std::vector<double> ts;

                std::vector<double> zv_ts;

                std::vector<double> conc;

                double rho_hn;

                double rime_hn;

                double lwc_hn;


                bool poor_ea;


        private:

                size_t getNumberMeasTemperatures(const mio::MeteoData& md);


                std::vector<double> fixedPositions;

                double minDepthSubsurf;

                size_t maxNumberMeasTemperatures;

                size_t numberMeasTemperatures;

                size_t numberFixedRates;

};


enum SN_FIELD{

        TEMPERATURE,

        SEEPAGE,

        SETTLEMENT,

        N_SN_FIELDS

};


enum {

        SOIL,

        ICE,

        WATER,

        WATER_PREF,

        AIR,

        N_COMPONENTS

};


enum SN_SOIL_DATA{

        SOIL_RHO,

        SOIL_K,

        SOIL_C,

        N_SOIL_FIELDS

};


class LayerData {

        public:

                LayerData();


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const LayerData& data);

                friend std::istream& operator>>(std::istream& is, LayerData& data);


                mio::Date depositionDate;

                double hl;

                size_t ne;

                double tl;

                double phiSoil;

                double phiIce;

                double phiIceReservoir;

                double phiIceReservoirCumul;

                double phiWater;

                double phiWaterPref;

                double phiVoids;

                std::vector<double> cSoil;

                std::vector<double> cIce;

                std::vector<double> cWater;

                std::vector<double> cVoids;

                double SoilRho;

                double SoilK;

                double SoilC;

                double rg;

                double sp;

                double dd;

                double rb;

                unsigned short int mk;

                double hr;

                double CDot;

                double metamo;

                double salinity;

                double h;

                double dsm;

};


class SN_SNOWSOIL_DATA {

        public:

                SN_SNOWSOIL_DATA() : meta(), profileDate(), nN(0), Height(0.),

                     nLayers(0), Ldata(), HS_last(0.), Albedo(mio::IOUtils::nodata),

                                                                     SoilAlb(mio::IOUtils::nodata), BareSoil_z0(mio::IOUtils::nodata),

                     Canopy_Height(mio::IOUtils::nodata), Canopy_LAI(mio::IOUtils::nodata),

                                                                                 Canopy_Direct_Throughfall(mio::IOUtils::nodata),

                                                                                 ErosionLevel(static_cast<int>(mio::IOUtils::nodata)), TimeCountDeltaHS(mio::IOUtils::nodata),

                                                                                 Canopy_BasalArea(mio::IOUtils::nodata), Canopy_diameter(mio::IOUtils::nodata),

                                                                                 Canopy_lai_frac_top_default(mio::IOUtils::nodata),Canopy_int_cap_snow(mio::IOUtils::nodata),

                                                                                 Canopy_alb_dry(mio::IOUtils::nodata),Canopy_alb_wet(mio::IOUtils::nodata),

                                                                                 Canopy_alb_snow(mio::IOUtils::nodata),Emissivity_soil(mio::IOUtils::nodata)

                    {}


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const SN_SNOWSOIL_DATA& data);

                friend std::istream& operator>>(std::istream& is, SN_SNOWSOIL_DATA& data);


                mio::StationData meta;

                mio::Date profileDate;

                size_t nN;

                double Height;

                size_t nLayers;

                std::vector<LayerData> Ldata;

                double HS_last;

                double Albedo;

                double SoilAlb;

                double BareSoil_z0;

                double Canopy_Height;

                double Canopy_LAI;

                double Canopy_Direct_Throughfall;

                int    ErosionLevel;

                double TimeCountDeltaHS;

                double Canopy_BasalArea;

                double Canopy_diameter;

                double Canopy_lai_frac_top_default;

                double Canopy_int_cap_snow;

                double Canopy_alb_dry;  // Albedo of dry canopy (calibr: 0.09, Alptal)

                double Canopy_alb_wet;  // Albedo of wet canopy (calibr: 0.09, Alptal)

                double Canopy_alb_snow;  // Albedo of snow covered albedo (calibr: 0.35, Alptal)

                double Emissivity_soil;


};


class ElementData {

        public:

                typedef enum YOUNG_MODULUS {

                            Sigrist,

                            Pow,

                            Exp

                } Young_Modulus;


                ElementData(const unsigned short int& in_ID);

                ElementData(const ElementData& cc); //required to get the correct back-reference in vanGenuchten object

                ElementData& operator=(const ElementData&) = default;


                bool checkVolContent();

                void heatCapacity();

                double coldContent() const;

                void updDensity();

                double extinction() const;

                void opticalEquivalentGrainSize();

                void snowResidualWaterContent();

                static double snowResidualWaterContent(const double& theta_i);

                double soilFieldCapacity() const;

                double soilRelativeHumidity() const;


                double snowElasticity() const;

                double neckStressEnhancement() const;

                double concaveNeckRadius() const;

                double neckLength() const;

                double neck2VolumetricStrain() const;


                void snowType();

                unsigned short int getSnowType() const;

                static unsigned short int snowType(const double& dendricity, const double& sphericity, const double& grain_dia, const unsigned short int& marker,

                        const double& theta_w, const double& res_wat_cont_loc);

                static double getYoungModule(const double& rho_slab, const Young_Modulus& model);


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const ElementData& data);

                friend std::istream& operator>>(std::istream& is, ElementData& data);


                mio::Date depositionDate;

                double L0, L;

                double Te;

                double gradT;

                double meltfreeze_tk;

                std::vector<double> theta;

                double h;

                mio::Array2D<double> conc;

                std::vector<double> k;

                //   Stored in order to visualize constitutive laws

                //   Will be used for creep field hydraulic conductivity in m3 s kg-1

                std::vector<double> c;

                //   Will also be used for creep specific snow water capacity  in m3 J-1

                std::vector<double> soil;

                double Rho;

                double M;

                double sw_abs;

                // Snow Metamorphism Data

                double rg;

                double dd;

                double sp;

                double ogs;

                double rb;

                double N3;

                size_t mk;

                unsigned short int type;

                double metamo;

                double salinity;

                double dth_w;

                double res_wat_cont;

                double Qmf;

                double QIntmf;

                double dEps, Eps, Eps_e, Eps_v;

                double Eps_Dot, Eps_vDot;

                double E;

                double S;

                double C;

                double CDot;

                double ps2rb;

                double s_strength;

                double hard;

                double S_dr;

                double crit_cut_length;

                vanGenuchten VG;

                double lwc_source;

                double PrefFlowArea;

                double theta_w_transfer;

                double theta_i_reservoir;

                double theta_i_reservoir_cumul;

                double SlopeParFlux;

                double Qph_up;

                double Qph_down;

                //NIED (H. Hirashima)

                double dsm;

                double rime;


                unsigned short int ID;

                static const unsigned short int noID;


                double rhov;

                double Qmm;

                double vapTrans_fluxDiff;

                double vapTrans_snowDenChangeRate;

                double vapTrans_cumulativeDenChange;

                double vapTrans_underSaturationDegree;

};


class NodeData {

        public:

                NodeData() : z(0.), u(0.), f(0.), udot(0.), T(0.), S_n(0.), S_s(0.), ssi(6.), hoar(0.),

                             dsm(0.), S_dsm(0.), Sigdsm(0.), rime(0.), water_flux(0.), rhov(0.) {} //HACK: set ssi to max_stability!


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const NodeData& data);

                friend std::istream& operator>>(std::istream& is, NodeData& data);


                double z;

                double u;

                double f;

                double udot;

                double T;

                double S_n;

                double S_s;

                double ssi;

                double hoar;


                //NIED (H. Hirashima)

                double dsm;

                double S_dsm;

                double Sigdsm;

                double rime;


                double water_flux;


                double rhov;

};


class CanopyData {

        public:

                CanopyData() : int_cap_snow(0.), int_cap_rain(0.),  interception_timecoef(0.), can_alb_dry(0.),

                can_alb_wet(0.), can_alb_snow(0.), krnt_lai(0.), can_diameter(0.), biomass_heat_capacity(0.),

                biomass_density(0.), lai_frac_top_default(0.), trunk_frac_height(0.), trunkalb(0.), et(0.),

                canopy_stabilitycorrection(true), roughmom_to_canopyheight_ratio(0.), displ_to_canopyheight_ratio(0.),

                raincrease_snow(0.), canopytemp_maxchange_perhour(0.), roughheat_to_roughmom_ratio(0.),

                can_ch0(0.), can_rs_mult(0.), rsmin(0.), f3_gd(0.), rootdepth(0.), wp_fraction(0.),

                h_wilt(0.), storage(0.), temp(0.), sigf(0.), ec(0.), lai(0.), z0m(0.), z0h(0.), zdispl(0.),

                height(0.), direct_throughfall(0.), ra(0.), rc(0.), rs(0.), rstransp(0.), canopyalb(0.),

                totalalb(0.), wetfraction(0.), intcapacity(0.), rswrac(0.), iswrac(0.), rswrbc(0.), iswrbc(0.),

                ilwrac(0.), rlwrac(0.), ilwrbc(0.), rlwrbc(0.), rsnet(0.), rlnet(0.), sensible(0.), latent(0.),

                latentcorr(0.), transp(0.), intevap(0.), interception(0.), throughfall(0.), snowunload(0.),

                snowfac(0.), rainfac(0.), liquidfraction(0.), sigftrunk(0.), Ttrunk(0.), CondFluxCanop(0.),

                CondFluxTrunks(0.), LWnet_Trunks(0.), SWnet_Trunks(0.), QStrunks(0.), forestfloor_alb(0.),

                BasalArea(0.), HMLeaves(0.), HMTrunks(0.) {}


                void initialize(const SN_SNOWSOIL_DATA& SSdata, const bool useCanopyModel, const bool isAlpine3D);

                void reset(const bool& cumsum_mass);

                void initializeSurfaceExchangeData();

                void multiplyFluxes(const double& factor);


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const CanopyData& data);

                friend std::istream& operator>>(std::istream& is, CanopyData& data);


                double int_cap_snow; //iMax in Gouttevin,2015

                double int_cap_rain;

                double interception_timecoef;


                double can_alb_dry;     // Albedo of dry canopy (calibr: 0.09, Alptal)

                double can_alb_wet;     // Albedo of wet canopy (calibr: 0.09, Alptal)

                double can_alb_snow;    // Albedo of snow covered albedo (calibr: 0.35, Alptal)

                double krnt_lai;        // Radiation transmissivity parameter, in the range 0.4-0.8 if the true LAI is used; higher if optical LAI is used.

                                        // (calibrated on Alptal)

                double can_diameter;    // average canopy (tree) diameter [m], parameter in the new radiation transfer model

                double biomass_heat_capacity;   // from Linroth et al., 2013 (J Kg-1 K-1)

                double biomass_density;         // from Linroth et al., 2013 (Kg m-3)

                double lai_frac_top_default;    // fraction of total LAI that is attributed to the uppermost layer. Here calibrated for Alptal.

                double trunk_frac_height;       // (optional) fraction of total tree height occupied by trunks,

                                                // used to calculate direct solar insolation of trunks.

                double trunkalb;                // trunk albedo

                double et;                      // trunk emissivity

                bool canopy_stabilitycorrection;

                double roughmom_to_canopyheight_ratio;

                double displ_to_canopyheight_ratio;

                double raincrease_snow;


                double canopytemp_maxchange_perhour;

                double roughheat_to_roughmom_ratio;

                double can_ch0;

                double can_rs_mult;

                double rsmin;

                double f3_gd;

                double rootdepth;

                double wp_fraction;

                double h_wilt;


                // State variable

                double storage;

                double temp;

                double sigf;

                double ec;

                // parameters

                double lai;

                double z0m;

                double z0h;

                double zdispl;

                double height;

                double direct_throughfall;

                // aerodynamic resistances

                double ra;

                double rc;

                double rs;

                double rstransp;

                // Averaged variables

                double canopyalb;

                double totalalb;

                double wetfraction;

                double intcapacity;

                // Radiations

                double rswrac;

                double iswrac;

                double rswrbc;

                double iswrbc;

                double ilwrac;

                double rlwrac;

                double ilwrbc;

                double rlwrbc;

                double rsnet;

                double rlnet;

                // Turbulent fluxes

                double sensible;

                double latent;

                double latentcorr;

                // Evap fluxes

                double transp;

                double intevap;

                // Mass fluxes

                double interception;

                double throughfall;

                double snowunload;


                double snowfac;

                double rainfac;

                double liquidfraction;

                double sigftrunk;

                double Ttrunk;

                double CondFluxCanop;

                double CondFluxTrunks;

                double LWnet_Trunks;

                double SWnet_Trunks;

                double QStrunks;

                double forestfloor_alb;

                double BasalArea;

                double HMLeaves;

                double HMTrunks;

};


class SeaIce;   // Foreward-declare sea ice class

class SnowStation {

        public:

                explicit SnowStation(const bool i_useCanopyModel=true, const bool i_useSoilLayers=true,

                         const bool i_isAlpine3D=false, const bool i_useSeaIceModule=false);

                SnowStation(const SnowStation& c);


                ~SnowStation();

                SnowStation& operator=(const SnowStation&);


                void initialize(const SN_SNOWSOIL_DATA& SSdata, const size_t& i_sector);

                void resize(const size_t& number_of_elements);


                void reduceNumberOfElements(const size_t& rnE);

                void combineElements(const size_t& number_top_elements, const bool& reduce_n_elements, const size_t& cond, const double& comb_thresh_l);

                static bool combineCondition(const ElementData& Edata0, const ElementData& Edata1, const double& depth, const bool& reduce_n_elements, const double& comb_thresh_l);

                static void mergeElements(ElementData& Edata0, const ElementData& Edata1, const bool& merge, const bool& topElement);

                void splitElement(const size_t& e);                                                     //Split an element

                void splitElements(const double& max_element_length, const double& comb_thresh_l);      //Check for splitting, calls splitElement(...) for actual splitting


                void compSnowpackMasses();

                void compSnowpackInternalEnergyChange(const double& sn_dt);

                void compSoilInternalEnergyChange(const double& sn_dt);

                double getLiquidWaterIndex() const;

                double getModelledTemperature(const double& z) const;

                double getTotalLateralFlowSnow() const;

                double getTotalLateralFlowSoil() const;

                void resetSlopeParFlux();


                size_t getNumberOfElements() const;

                size_t getNumberOfNodes() const;

                bool isGlacier(const bool& hydro=false) const;

                bool hasSoilLayers() const;

                double findMarkedReferenceLayer() const;


                size_t find_tag(const size_t& tag) const;


                void reset_water_fluxes();


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const SnowStation& data);

                friend std::istream& operator>>(std::istream& is, SnowStation& data);


                mio::StationData meta;

                double cos_sl;

                size_t sector;


                CanopyData Cdata;

                SeaIce* Seaice;

                double pAlbedo;

                double Albedo;

                double SoilAlb;

                double SoilEmissivity;

                double BareSoil_z0;

                size_t SoilNode;

                double Ground;

                double cH;

                double mH;

                double mass_sum;

                double swe;

                double lwc_sum;

                double lwc_sum_soil;

                double swc_sum_soil;

                double hn;

                double rho_hn;

                double rime_hn;

                size_t ErosionLevel;

                double ErosionMass;

                char S_class1;

                char S_class2;

                double S_d;

                double z_S_d;

                double S_n;

                double z_S_n;

                double S_s;

                double z_S_s;

                double S_4;

                double z_S_4;

                double S_5;

                double z_S_5;

                std::vector<NodeData> Ndata;

                std::vector<ElementData> Edata;

                void *Kt;

                double ColdContent;

                double ColdContentSoil;

                double dIntEnergy;

                double dIntEnergySoil;

                double meltFreezeEnergy;

                double meltFreezeEnergySoil;

                double ReSolver_dt;

                bool windward;

                double TimeCountDeltaHS;

                static const double comb_thresh_l_ratio, comb_thresh_ice, comb_thresh_water;

                static const double comb_thresh_dd, comb_thresh_sp, comb_thresh_rg;

                static const double thresh_moist_snow, thresh_moist_soil;

                static const size_t number_top_elements;

                static unsigned short number_of_solutes;


        private:

                size_t nNodes;

                size_t nElems;

                unsigned short int maxElementID;

                bool useCanopyModel, useSoilLayers, isAlpine3D;

                static double flexibleMaxElemLength(const double& depth, const double& comb_thresh_l);

};


class BoundCond {


        public:

                BoundCond() : lw_out(0.), lw_net(0.), qs(0.), ql(0.), qr(0.), qg(Constants::undefined) {}

                const std::string toString() const;

                void reset();


                double lw_out;

                double lw_net;

                double qs;

                double ql;

                double qr;

                double qg;

};


class SurfaceFluxes {

        public:

                enum SN_MASS_CHANGES {

                        MS_TOTALMASS,

                        MS_SWE,

                        MS_WATER,

                        MS_WATER_SOIL,

                        MS_ICE_SOIL,

                        MS_HNW,

                        MS_RAIN,

                        MS_WIND,

                        MS_EVAPORATION,

                        MS_SUBLIMATION,

                        MS_SNOWPACK_RUNOFF,

                        MS_SURFACE_MASS_FLUX,

                        MS_SOIL_RUNOFF,

                        MS_FLOODING,

                        MS_ICEBASE_MELTING_FREEZING,

                        N_MASS_CHANGES

                };


                const std::string toString() const;

                friend std::ostream& operator<<(std::ostream& os, const SurfaceFluxes& data);

                friend std::istream& operator>>(std::istream& is, SurfaceFluxes& data);


                SurfaceFluxes();


                void reset(const bool& cumsum_mass);

                void compSnowSoilHeatFlux(const SnowStation& Xdata);

                void collectSurfaceFluxes(const BoundCond& Bdata, SnowStation& Xdata, const CurrentMeteo& Mdata);

                void multiplyFluxes(const double& factor);


                double lw_in;

                double lw_out;

                double lw_net;

                double qs;

                double ql;

                double hoar;

                double qr;

                double qg;

                double qg0;

                double sw_hor;

                double sw_in;

                double sw_out;

                double qw;

                double sw_dir;

                double sw_diff;

                double pAlbedo;

                double mAlbedo;

                double dIntEnergy;

                double dIntEnergySoil;

                double meltFreezeEnergy;

                double meltFreezeEnergySoil;


                double drift;

                std::vector<double> mass;

                std::vector<double> load;

                double dhs_corr;

                double cRho_hn;

                double mRho_hn;


};


//HACK: could it be moved to plugins? (as well as Aggregate)

class SnowProfileLayer {

        public:

                SnowProfileLayer();


                void average(const double& w1, const double& w2, const SnowProfileLayer& Pdata);

                static std::vector<SnowProfileLayer> generateProfile(const mio::Date& dateOfProfile, const SnowStation& Xdata, const double hoar_density_surf, const double hoar_min_size_surf);


                // Profile meta data

                mio::Date profileDate;

                std::string stationname;

                unsigned char  loc_for_snow;

                unsigned char  loc_for_wind;


                mio::Date depositionDate;

                double height;

                double rho;

                double T;

                double gradT;

                double v_strain_rate;

                double theta_i;

                double theta_w;

                double theta_a;

                double grain_size;

                double bond_size;

                double dendricity;

                double sphericity;

                double ogs;

                double coordin_num;

                unsigned short int marker;

                short unsigned int type;

                double hard;


        private:

                void generateLayer(const ElementData& Edata, const NodeData& Ndata);

                void generateLayer(const ElementData& Edata, const NodeData& Ndata,

                                   const mio::Date& dateOfProfile, const double hoar_density_surf);

};


class RunInfo {

        public:

                RunInfo();

                RunInfo(const RunInfo& orig);

                RunInfo& operator=(const RunInfo&) {return *this;} //everything is static, so we can not change anything


                const std::string version;

                const double version_num;

                const mio::Date computation_date;

                const std::string compilation_date;

                const std::string user;

                const std::string hostname;


        private:

                static double getNumericVersion(std::string version_str);

                static mio::Date getRunDate();

                static std::string getCompilationDate();

};


#endif

Constants.h

SOLID
@ SOLID
Solid.
Definition: DataClasses.h:41

GAS
@ GAS
Gas.
Definition: DataClasses.h:43

LIQUID
@ LIQUID
Liquid.
Definition: DataClasses.h:42

NUMBER_OF_PHASES
@ NUMBER_OF_PHASES
Definition: DataClasses.h:44

SN_SOIL_DATA
SN_SOIL_DATA
The different soil fields (PERMAFROST)
Definition: DataClasses.h:174

SOIL_RHO
@ SOIL_RHO
Bulk density of dry soil component (without accounting for porosity)
Definition: DataClasses.h:175

SOIL_C
@ SOIL_C
Specific heat of dry soil component.
Definition: DataClasses.h:177

SOIL_K
@ SOIL_K
Bulk conductivity of dry soil component.
Definition: DataClasses.h:176

N_SOIL_FIELDS
@ N_SOIL_FIELDS
Definition: DataClasses.h:178

SN_FIELD
SN_FIELD
The 3 mathematical fields that can be solved.
Definition: DataClasses.h:156

SEEPAGE
@ SEEPAGE
Water transport.
Definition: DataClasses.h:158

N_SN_FIELDS
@ N_SN_FIELDS
Definition: DataClasses.h:160

TEMPERATURE
@ TEMPERATURE
Temperature (phase change and metamorphism)
Definition: DataClasses.h:157

SETTLEMENT
@ SETTLEMENT
Creep displacements.
Definition: DataClasses.h:159

AIR
@ AIR
Air.
Definition: DataClasses.h:169

WATER_PREF
@ WATER_PREF
Water in preferential flow.
Definition: DataClasses.h:168

SOIL
@ SOIL
Soil.
Definition: DataClasses.h:165

ICE
@ ICE
Ice.
Definition: DataClasses.h:166

N_COMPONENTS
@ N_COMPONENTS
Definition: DataClasses.h:170

WATER
@ WATER
Water.
Definition: DataClasses.h:167

SeaIce.h

SnowpackConfig.h

BoundCond
BoundCond is used to set Neumann boundary conditions.
Definition: DataClasses.h:711

BoundCond::qs
double qs
sensible heat
Definition: DataClasses.h:720

BoundCond::qg
double qg
geothermal heat flux or heat flux at lower boundary
Definition: DataClasses.h:723

BoundCond::ql
double ql
latent heat
Definition: DataClasses.h:721

BoundCond::BoundCond
BoundCond()
Definition: DataClasses.h:714

BoundCond::lw_out
double lw_out
outgoing longwave radiation
Definition: DataClasses.h:718

BoundCond::toString
const std::string toString() const
Definition: DataClasses.cc:272

BoundCond::qr
double qr
rain energy
Definition: DataClasses.h:722

BoundCond::lw_net
double lw_net
net longwave radiation
Definition: DataClasses.h:719

BoundCond::reset
void reset()
Definition: DataClasses.cc:282

CanopyData
Canopy data used as a pointer in the SnowStation structure.
Definition: DataClasses.h:436

CanopyData::wp_fraction
double wp_fraction
Wilting point, defined as a fraction of water content at field capacity (-)
Definition: DataClasses.h:530

CanopyData::CanopyData
CanopyData()
Definition: DataClasses.h:438

CanopyData::wetfraction
double wetfraction
fraction of canopy covered by interception [-]
Definition: DataClasses.h:555

CanopyData::intcapacity
double intcapacity
maximum interception storage [mm]
Definition: DataClasses.h:556

CanopyData::rc
double rc
from canopy to canopy air
Definition: DataClasses.h:549

CanopyData::f3_gd
double f3_gd
gd (Pa-1) parameter for canopy surface resistance response to vapour pressure:
Definition: DataClasses.h:526

CanopyData::sigftrunk
double sigftrunk
radiation interception cross section for trunk layer ()
Definition: DataClasses.h:583

CanopyData::latentcorr
double latentcorr
Definition: DataClasses.h:571

CanopyData::LWnet_Trunks
double LWnet_Trunks
net LW to trunks (>0 towards trunks)
Definition: DataClasses.h:587

CanopyData::rswrac
double rswrac
upward shortwave above canopy
Definition: DataClasses.h:558

CanopyData::rlnet
double rlnet
net longwave radiation
Definition: DataClasses.h:567

CanopyData::snowfac
double snowfac
snowfall above canopy
Definition: DataClasses.h:580

CanopyData::CondFluxTrunks
double CondFluxTrunks
biomass heat storage flux towards Trunks (if 2L)
Definition: DataClasses.h:586

CanopyData::operator>>
friend std::istream & operator>>(std::istream &is, CanopyData &data)
Definition: DataClasses.cc:920

CanopyData::height
double height
Definition: DataClasses.h:545

CanopyData::CondFluxCanop
double CondFluxCanop
biomass heat storage flux towards Canopy (if 1L) towards Leaves (if 2L). (>0 towards canopy)
Definition: DataClasses.h:585

CanopyData::can_alb_snow
double can_alb_snow
Definition: DataClasses.h:483

CanopyData::canopy_stabilitycorrection
bool canopy_stabilitycorrection
Definition: DataClasses.h:498

CanopyData::krnt_lai
double krnt_lai
Definition: DataClasses.h:484

CanopyData::zdispl
double zdispl
Definition: DataClasses.h:544

CanopyData::totalalb
double totalalb
total albedo above canopy and snow/soil surface [-]
Definition: DataClasses.h:554

CanopyData::iswrbc
double iswrbc
downward shortwave radiation below canopy
Definition: DataClasses.h:561

CanopyData::transp
double transp
Definition: DataClasses.h:573

CanopyData::can_rs_mult
double can_rs_mult
1+CAN_RS_MULT = maximum factor to increase Cdata->rs below canopy
Definition: DataClasses.h:517

CanopyData::rswrbc
double rswrbc
upward shortwave below canopy
Definition: DataClasses.h:560

CanopyData::toString
const std::string toString() const
Definition: DataClasses.cc:1010

CanopyData::rs
double rs
from subsurface to canpopy air
Definition: DataClasses.h:550

CanopyData::rlwrbc
double rlwrbc
upward longwave radiation BELOW canopy
Definition: DataClasses.h:565

CanopyData::multiplyFluxes
void multiplyFluxes(const double &factor)
If multiple fluxes have been summed over multiple time steps, the fluxes then need to be averaged by ...
Definition: DataClasses.cc:801

CanopyData::trunkalb
double trunkalb
Definition: DataClasses.h:494

CanopyData::rsmin
double rsmin
TRANSPIRATION.
Definition: DataClasses.h:520

CanopyData::ra
double ra
from canopy air to reference height
Definition: DataClasses.h:548

CanopyData::SWnet_Trunks
double SWnet_Trunks
net SW to trunks (>0 towards trunks)
Definition: DataClasses.h:588

CanopyData::rstransp
double rstransp
stomatal surface resistance for transpiration
Definition: DataClasses.h:551

CanopyData::reset
void reset(const bool &cumsum_mass)
Definition: DataClasses.cc:761

CanopyData::biomass_heat_capacity
double biomass_heat_capacity
Definition: DataClasses.h:489

CanopyData::HMTrunks
double HMTrunks
Trunks heat mass (J K-1 /m2 ground surface)
Definition: DataClasses.h:593

CanopyData::interception_timecoef
double interception_timecoef
Definition: DataClasses.h:478

CanopyData::z0m
double z0m
Definition: DataClasses.h:542

CanopyData::rsnet
double rsnet
net shortwave radiation
Definition: DataClasses.h:566

CanopyData::interception
double interception
Definition: DataClasses.h:576

CanopyData::lai
double lai
Definition: DataClasses.h:541

CanopyData::can_alb_dry
double can_alb_dry
RADIATION BALANCE.
Definition: DataClasses.h:481

CanopyData::can_ch0
double can_ch0
minimum heat exchange (Wm-2K-1) at zero wind
Definition: DataClasses.h:515

CanopyData::ilwrbc
double ilwrbc
downward longwave radiation BELOW canopy
Definition: DataClasses.h:564

CanopyData::BasalArea
double BasalArea
basal area of trees on the stand
Definition: DataClasses.h:591

CanopyData::direct_throughfall
double direct_throughfall
Definition: DataClasses.h:546

CanopyData::lai_frac_top_default
double lai_frac_top_default
Definition: DataClasses.h:491

CanopyData::biomass_density
double biomass_density
Definition: DataClasses.h:490

CanopyData::latent
double latent
Definition: DataClasses.h:570

CanopyData::initializeSurfaceExchangeData
void initializeSurfaceExchangeData()
Function called to initialize the canopy "Surface" exchange data (to enable accumulated mass and ener...
Definition: DataClasses.cc:1107

CanopyData::HMLeaves
double HMLeaves
Leaves heat mass (J K-1 /m2 ground surface)
Definition: DataClasses.h:592

CanopyData::can_alb_wet
double can_alb_wet
Definition: DataClasses.h:482

CanopyData::initialize
void initialize(const SN_SNOWSOIL_DATA &SSdata, const bool useCanopyModel, const bool isAlpine3D)
Initialize all the CData elements value with values by default or with values read from the SNO file....
Definition: DataClasses.cc:548

CanopyData::int_cap_rain
double int_cap_rain
Specific interception capacity for rain (I_LAI) (mm/LAI)
Definition: DataClasses.h:473

CanopyData::intevap
double intevap
Definition: DataClasses.h:574

CanopyData::iswrac
double iswrac
downward shortwave radiation above canopy
Definition: DataClasses.h:559

CanopyData::roughmom_to_canopyheight_ratio
double roughmom_to_canopyheight_ratio
Ratio between canopy height and roughness length.
Definition: DataClasses.h:500

CanopyData::rainfac
double rainfac
rainfall above canopy
Definition: DataClasses.h:581

CanopyData::int_cap_snow
double int_cap_snow
Definition: DataClasses.h:471

CanopyData::temp
double temp
temperature (K)
Definition: DataClasses.h:537

CanopyData::displ_to_canopyheight_ratio
double displ_to_canopyheight_ratio
As above for displacement height.
Definition: DataClasses.h:502

CanopyData::storage
double storage
intercepted water (mm or kg m-2)
Definition: DataClasses.h:536

CanopyData::h_wilt
double h_wilt
Wilting point pressure head, when using Richards equation for soil.
Definition: DataClasses.h:532

CanopyData::rootdepth
double rootdepth
Root depth, determining the soil layers influenced by root water uptake.
Definition: DataClasses.h:528

CanopyData::sigf
double sigf
radiation transmissivity (1)
Definition: DataClasses.h:538

CanopyData::rlwrac
double rlwrac
upward longwave radiation ABOVE canopy
Definition: DataClasses.h:563

CanopyData::liquidfraction
double liquidfraction
Definition: DataClasses.h:582

CanopyData::snowunload
double snowunload
Definition: DataClasses.h:578

CanopyData::roughheat_to_roughmom_ratio
double roughheat_to_roughmom_ratio
(~=1, but Not allowed to be exactly 1)
Definition: DataClasses.h:513

CanopyData::sensible
double sensible
Definition: DataClasses.h:569

CanopyData::operator<<
friend std::ostream & operator<<(std::ostream &os, const CanopyData &data)
Definition: DataClasses.cc:830

CanopyData::trunk_frac_height
double trunk_frac_height
Definition: DataClasses.h:492

CanopyData::ilwrac
double ilwrac
downward longwave radiation ABOVE canopy
Definition: DataClasses.h:562

CanopyData::Ttrunk
double Ttrunk
trunk temperature (K)
Definition: DataClasses.h:584

CanopyData::forestfloor_alb
double forestfloor_alb
albedo of the forest floor
Definition: DataClasses.h:590

CanopyData::QStrunks
double QStrunks
sensible heat flux from trunks (>0 if heat lost from trunk)
Definition: DataClasses.h:589

CanopyData::z0h
double z0h
Definition: DataClasses.h:543

CanopyData::canopytemp_maxchange_perhour
double canopytemp_maxchange_perhour
Maximum allowed canopy temperature change (K hr-1)
Definition: DataClasses.h:511

CanopyData::ec
double ec
longwave emissivity (1)
Definition: DataClasses.h:539

CanopyData::raincrease_snow
double raincrease_snow
Definition: DataClasses.h:508

CanopyData::can_diameter
double can_diameter
Definition: DataClasses.h:486

CanopyData::canopyalb
double canopyalb
canopy albedo [-]
Definition: DataClasses.h:553

CanopyData::throughfall
double throughfall
Definition: DataClasses.h:577

CanopyData::et
double et
Definition: DataClasses.h:495

CurrentMeteo
CurrentMeteo is the class of interpolated meteo data for the current calculation time step  It contai...
Definition: DataClasses.h:83

CurrentMeteo::conc
std::vector< double > conc
Solute concentrations in precipitation.
Definition: DataClasses.h:138

CurrentMeteo::date
mio::Date date
Date of current meteo data.
Definition: DataClasses.h:102

CurrentMeteo::adv_heat
double adv_heat
Advective heat to inject in the soil (if ADVECTIVE_HEAT and related parameters set to true)
Definition: DataClasses.h:134

CurrentMeteo::getFixedPositions
void getFixedPositions(std::vector< double > &positions) const
Definition: DataClasses.cc:3169

CurrentMeteo::reset
void reset(const SnowpackConfig &i_cfg)
Definition: DataClasses.cc:3091

CurrentMeteo::vw
double vw
Wind velocity at snow station (m s-1)
Definition: DataClasses.h:106

CurrentMeteo::getNumberMeasTemperatures
size_t getNumberMeasTemperatures() const
Returns the number of measured snow/soil temperatures stored in MeteoData.
Definition: DataClasses.cc:3150

CurrentMeteo::CurrentMeteo
CurrentMeteo()
Definition: DataClasses.cc:3065

CurrentMeteo::rime_hn
double rime_hn
riming index of new snow
Definition: DataClasses.h:140

CurrentMeteo::toString
const std::string toString() const
Definition: DataClasses.cc:3352

CurrentMeteo::lwc_hn
double lwc_hn
liquid water content of new snow
Definition: DataClasses.h:141

CurrentMeteo::poor_ea
bool poor_ea
when ilwr has not been measured nor parametrized in good conditions, it could be redone later on
Definition: DataClasses.h:143

CurrentMeteo::z0
double z0
The roughness length computed in SnowDrift and also used later for the MeteoHeat fluxes (m)
Definition: DataClasses.h:113

CurrentMeteo::psum_ph
double psum_ph
precipitation phase for the current timestep (between 0 and 1, 0 is fully solid while 1 is fully liqu...
Definition: DataClasses.h:128

CurrentMeteo::hs_rate
double hs_rate
The rate of change in snow depth (m h-1)
Definition: DataClasses.h:132

CurrentMeteo::vw_max
double vw_max
Maximum wind velocity at snow station (m s-1)
Definition: DataClasses.h:108

CurrentMeteo::iswr
double iswr
Incoming SHORTWAVE radiation (W m-2)
Definition: DataClasses.h:115

CurrentMeteo::getNumberFixedPositions
size_t getNumberFixedPositions() const
Definition: DataClasses.cc:3174

CurrentMeteo::setMeasTempParameters
void setMeasTempParameters(const mio::MeteoData &md)
Definition: DataClasses.cc:3112

CurrentMeteo::rho_hn
double rho_hn
Measured new snow density (kg m-3)
Definition: DataClasses.h:139

CurrentMeteo::operator<<
friend std::ostream & operator<<(std::ostream &os, const CurrentMeteo &data)
Definition: DataClasses.cc:3220

CurrentMeteo::ta
double ta
Air temperature (K)
Definition: DataClasses.h:103

CurrentMeteo::vw_avg
double vw_avg
Running mean of wind velocity at snow station (m s-1)
Definition: DataClasses.h:107

CurrentMeteo::lw_net
double lw_net
Net longwave radiation (W m-2)
Definition: DataClasses.h:122

CurrentMeteo::dw
double dw
Wind direction at snow station (deg)
Definition: DataClasses.h:109

CurrentMeteo::tss_a24h
double tss_a24h
Snow surface temperature averaged over past 24 hours (K)
Definition: DataClasses.h:125

CurrentMeteo::operator>>
friend std::istream & operator>>(std::istream &is, CurrentMeteo &data)
Definition: DataClasses.cc:3284

CurrentMeteo::dw_drift
double dw_drift
Wind direction of blowing and drifting snow (operational: wind ridge station)
Definition: DataClasses.h:111

CurrentMeteo::vw_drift
double vw_drift
Wind velocity for blowing and drifting snow (operational: wind ridge station)
Definition: DataClasses.h:110

CurrentMeteo::psum
double psum
precipitation sum over the current timestep (mm)
Definition: DataClasses.h:127

CurrentMeteo::rswr
double rswr
Reflected SHORTWAVE radiation (W m-2) divide this value by the ALBEDO to get iswr.
Definition: DataClasses.h:116

CurrentMeteo::diff
double diff
Diffuse radiation from the sky (W m-2)
Definition: DataClasses.h:118

CurrentMeteo::ustar
double ustar
The friction velocity (m s-1) computed in mt_MicroMet() and also used later for the MeteoHeat fluxes.
Definition: DataClasses.h:112

CurrentMeteo::zv_ts
std::vector< double > zv_ts
Positions of all measured snow or/and soil temperatures (m)
Definition: DataClasses.h:137

CurrentMeteo::copySnowTemperatures
void copySnowTemperatures(const mio::MeteoData &md, const unsigned int &current_slope)
Definition: DataClasses.cc:3189

CurrentMeteo::psi_s
double psi_s
Stability correction for scalar heat fluxes.
Definition: DataClasses.h:114

CurrentMeteo::psum_tech
double psum_tech
Equivalent precipitation water sum for technical snow over the current timestep (mm)
Definition: DataClasses.h:129

CurrentMeteo::hs_a3h
double hs_a3h
Snow depth averaged over 3 past hours.
Definition: DataClasses.h:131

CurrentMeteo::rh_avg
double rh_avg
Running mean of relative humidity (1)
Definition: DataClasses.h:105

CurrentMeteo::tss_a12h
double tss_a12h
Snow surface temperature averaged over past 12 hours (K)
Definition: DataClasses.h:124

CurrentMeteo::getMaxNumberMeasTemperatures
size_t getMaxNumberMeasTemperatures() const
Definition: DataClasses.cc:3184

CurrentMeteo::tss
double tss
Snow surface temperature (K)
Definition: DataClasses.h:123

CurrentMeteo::ea
double ea
Atmospheric emissivity (1)
Definition: DataClasses.h:121

CurrentMeteo::hs
double hs
The measured height of snow (m)
Definition: DataClasses.h:130

CurrentMeteo::copySolutes
void copySolutes(const mio::MeteoData &md, const size_t &i_number_of_solutes)
Definition: DataClasses.cc:3206

CurrentMeteo::ts0
double ts0
Bottom temperatures of snow/soil pack (K)
Definition: DataClasses.h:126

CurrentMeteo::ts
std::vector< double > ts
Measured snow or/and soil temperatures (K)
Definition: DataClasses.h:136

CurrentMeteo::rh
double rh
Relative humidity (% or 1)
Definition: DataClasses.h:104

CurrentMeteo::geo_heat
double geo_heat
Geo heat flux (W/m^2), for the neumann lower boundary condition in the heat equation.
Definition: DataClasses.h:133

CurrentMeteo::getNumberFixedRates
size_t getNumberFixedRates() const
Definition: DataClasses.cc:3179

CurrentMeteo::elev
double elev
Solar elevation to be used in Canopy.c (rad) => see also.
Definition: DataClasses.h:120

CurrentMeteo::mAlbedo
double mAlbedo
Measured snow albedo.
Definition: DataClasses.h:117

CurrentMeteo::dir_h
double dir_h
Horizontal direct radiation from the sky (W m-2)
Definition: DataClasses.h:119

ElementData
ELEMENT DATA used as a pointer in the SnowStation structure NOTE on M below: this is the mass of an e...
Definition: DataClasses.h:284

ElementData::L
double L
Original and present element thickness (m)
Definition: DataClasses.h:325

ElementData::c
std::vector< double > c
For example, specific heat of TEMPERATURE field (J kg-1 K-1)
Definition: DataClasses.h:335

ElementData::salinity
double salinity
bulk salinity (PSU, which is g/kg)
Definition: DataClasses.h:351

ElementData::Qph_down
double Qph_down
Heat source/sink due to phase changes for the heat equation (W/m^3), at the lower node of the element...
Definition: DataClasses.h:375

ElementData::QIntmf
double QIntmf
Apparent change in internal energy due to phase change (caused by difference in heat capacity of wate...
Definition: DataClasses.h:355

ElementData::neck2VolumetricStrain
double neck2VolumetricStrain() const
Relates the neck strain to the global volumetric strain.
Definition: DataClasses.cc:1655

ElementData::VG
vanGenuchten VG
Van Genuchten Model for water retention.
Definition: DataClasses.h:367

ElementData::soilFieldCapacity
double soilFieldCapacity() const
Field Capacity Soil is dependent on grain properties. This nice formulation is based on some tedious ...
Definition: DataClasses.cc:1541

ElementData::Te
double Te
mean element temperature (K)
Definition: DataClasses.h:326

ElementData::Qmf
double Qmf
Subsurface Melting & Freezing Data: change of energy due to phase changes (melt-freeze)
Definition: DataClasses.h:354

ElementData::CDot
double CDot
Stress rate (Pa s-1), that is the overload change rate.
Definition: DataClasses.h:361

ElementData::Eps_e
double Eps_e
Definition: DataClasses.h:356

ElementData::getSnowType
unsigned short int getSnowType() const
Determine the type of snow  First revisited by Fierz and Bellaire 2006 and 2007 TODO needs to be adap...
Definition: DataClasses.cc:1673

ElementData::S_dr
double S_dr
Stability Index based on deformation rate (Direct Action Avalanching)
Definition: DataClasses.h:365

ElementData::soilRelativeHumidity
double soilRelativeHumidity() const
soilRelativeHumidity
Definition: DataClasses.cc:1580

ElementData::coldContent
double coldContent() const
Computes cold content of an element, taking meltfreeze_tk as reference.
Definition: DataClasses.cc:1443

ElementData::operator=
ElementData & operator=(const ElementData &)=default
Assignement operator.

ElementData::PrefFlowArea
double PrefFlowArea
Preferential flow path relative area (-)
Definition: DataClasses.h:369

ElementData::Eps_v
double Eps_v
Total element strain (GREEN'S strains – TOTAL LAGRANGIAN FORMULATION): Eps_e is elastic and Eps_v is ...
Definition: DataClasses.h:356

ElementData::res_wat_cont
double res_wat_cont
Residual water content.
Definition: DataClasses.h:353

ElementData::extinction
double extinction() const
Density dependent extinction coefficient -> Michi's magic trick... out of his magic hat.
Definition: DataClasses.cc:1482

ElementData::type
unsigned short int type
grain class
Definition: DataClasses.h:349

ElementData::Young_Modulus
Young_Modulus
This enum provides names for possible Young's modulus calculations.
Definition: DataClasses.h:287

ElementData::Pow
@ Pow
another power law
Definition: DataClasses.h:289

ElementData::Exp
@ Exp
exponential law
Definition: DataClasses.h:290

ElementData::Sigrist
@ Sigrist
Sigrist, 2006.
Definition: DataClasses.h:288

ElementData::ElementData
ElementData(const unsigned short int &in_ID)
Definition: DataClasses.cc:1147

ElementData::sw_abs
double sw_abs
total absorbed shortwave radiation by the element (W m-2)
Definition: DataClasses.h:340

ElementData::theta_i_reservoir_cumul
double theta_i_reservoir_cumul
Volumetric ice content in cumulated ice reservoir (1)
Definition: DataClasses.h:372

ElementData::M
double M
the total mass of the element (kg m-2)
Definition: DataClasses.h:339

ElementData::vapTrans_fluxDiff
double vapTrans_fluxDiff
vapor dissusion flux in the case of vapor transport (W/m^2/s)
Definition: DataClasses.h:385

ElementData::vapTrans_underSaturationDegree
double vapTrans_underSaturationDegree
the degree of undersaturation, (rhov-rohv_sat)/rhov_sat (-)
Definition: DataClasses.h:388

ElementData::E
double E
Young's modulus of elasticity (Pa)
Definition: DataClasses.h:358

ElementData::vapTrans_snowDenChangeRate
double vapTrans_snowDenChangeRate
snow density change rate in the case of vapor transport (kg/m^3/s)
Definition: DataClasses.h:386

ElementData::metamo
double metamo
keep track of metamorphism
Definition: DataClasses.h:350

ElementData::lwc_source
double lwc_source
Source/sink term for Richards equation (m^3/m^3 / timestep)
Definition: DataClasses.h:368

ElementData::rb
double rb
grain bond radius (mm)
Definition: DataClasses.h:346

ElementData::operator<<
friend std::ostream & operator<<(std::ostream &os, const ElementData &data)
Definition: DataClasses.cc:1172

ElementData::mk
size_t mk
grain marker (history dependent)
Definition: DataClasses.h:348

ElementData::ID
unsigned short int ID
Element ID used to track elements.
Definition: DataClasses.h:380

ElementData::h
double h
capillary pressure head (m)
Definition: DataClasses.h:330

ElementData::k
std::vector< double > k
For example, heat conductivity of TEMPERATURE field (W m-1 K-1)
Definition: DataClasses.h:332

ElementData::snowType
void snowType()
Definition: DataClasses.cc:1661

ElementData::rg
double rg
grain radius (mm)
Definition: DataClasses.h:342

ElementData::rhov
double rhov
vapor density...(kg/m^3)
Definition: DataClasses.h:383

ElementData::conc
mio::Array2D< double > conc
Concentration for chemical constituents in (kg m-3)
Definition: DataClasses.h:331

ElementData::dsm
double dsm
Dry snow metamorphism factor.
Definition: DataClasses.h:377

ElementData::heatCapacity
void heatCapacity()
Computes heat capacity of an element based on volumetric contents.
Definition: DataClasses.cc:1426

ElementData::hard
double hard
Parameterized hand hardness (1)
Definition: DataClasses.h:364

ElementData::gradT
double gradT
temperature gradient over element (K m-1)
Definition: DataClasses.h:327

ElementData::meltfreeze_tk
double meltfreeze_tk
melt/freeze temperature of layer (principally initialized as 0 degC, but enables possibility for free...
Definition: DataClasses.h:328

ElementData::getYoungModule
static double getYoungModule(const double &rho_slab, const Young_Modulus &model)
Definition: DataClasses.cc:1341

ElementData::theta_w_transfer
double theta_w_transfer
Volumetric content of water transferred from preferential flow to matrix domain (1)
Definition: DataClasses.h:370

ElementData::theta_i_reservoir
double theta_i_reservoir
Volumetric ice content in ice reservoir (1)
Definition: DataClasses.h:371

ElementData::theta
std::vector< double > theta
volumetric contents: SOIL, ICE, WATER, WATER_PREF, AIR (1)
Definition: DataClasses.h:329

ElementData::depositionDate
mio::Date depositionDate
Date of deposition.
Definition: DataClasses.h:324

ElementData::operator>>
friend std::istream & operator>>(std::istream &is, ElementData &data)
Definition: DataClasses.cc:1254

ElementData::s_strength
double s_strength
Parameterized snow shear strength (kPa)
Definition: DataClasses.h:363

ElementData::Qph_up
double Qph_up
Heat source/sink due to phase changes for the heat equation (W/m^3), at the upper node of the element...
Definition: DataClasses.h:374

ElementData::S
double S
Total Element Stress (Pa), S being the energy conjugate stress.
Definition: DataClasses.h:359

ElementData::ps2rb
double ps2rb
proportion of grain bond growth due to pressure sintering (1)
Definition: DataClasses.h:362

ElementData::C
double C
Total Element Stress (Pa), C being the real or the Cauchy stress, which is output.
Definition: DataClasses.h:360

ElementData::Eps_vDot
double Eps_vDot
Total Strain Rate, elastic and viscous, respectively (s-1) (Simply, Eps/sn_dt)
Definition: DataClasses.h:357

ElementData::Eps_Dot
double Eps_Dot
Definition: DataClasses.h:357

ElementData::N3
double N3
grain Coordination number (1)
Definition: DataClasses.h:347

ElementData::vapTrans_cumulativeDenChange
double vapTrans_cumulativeDenChange
cumulative density change in the case of vapor transport (kg/m^3)
Definition: DataClasses.h:387

ElementData::Eps
double Eps
Definition: DataClasses.h:356

ElementData::neckStressEnhancement
double neckStressEnhancement() const
Computes the enhancement of hydrostatically applied stress (overburden) in the necks (or bonds)
Definition: DataClasses.cc:1615

ElementData::dth_w
double dth_w
Subsurface Melting & Freezing Data: change of water content.
Definition: DataClasses.h:352

ElementData::Rho
double Rho
mean element density (or BULK density; kg m-3), that is, rho=M/V=sum( theta(i)*rho(i) )
Definition: DataClasses.h:338

ElementData::SlopeParFlux
double SlopeParFlux
Slope parallel flux (m^3/m^3 * m / timestep)
Definition: DataClasses.h:373

ElementData::neckLength
double neckLength() const
Computes the neck (or bond) length (mm)
Definition: DataClasses.cc:1644

ElementData::updDensity
void updDensity()
Updates element density.
Definition: DataClasses.cc:1452

ElementData::rime
double rime
Rime index.
Definition: DataClasses.h:378

ElementData::L0
double L0
Definition: DataClasses.h:325

ElementData::ogs
double ogs
optical equivalent grain size (mm)
Definition: DataClasses.h:345

ElementData::checkVolContent
bool checkVolContent()
Check volumetric content.
Definition: DataClasses.cc:1369

ElementData::opticalEquivalentGrainSize
void opticalEquivalentGrainSize()
Opical equivalent grain size CROCUS implementation as described in Vionnet et al.,...
Definition: DataClasses.cc:1468

ElementData::noID
static const unsigned short int noID
Definition: DataClasses.h:381

ElementData::snowResidualWaterContent
void snowResidualWaterContent()
Definition: DataClasses.cc:1490

ElementData::soil
std::vector< double > soil
Contains the heat conductivity, capacity and dry density of the soil (solid, non-ice) component phase...
Definition: DataClasses.h:337

ElementData::dd
double dd
snow dendricity: 0 = none, 1 = newsnow
Definition: DataClasses.h:343

ElementData::snowElasticity
double snowElasticity() const
SNOW ELASTICITY : This important routine was programmed by Marc Christen, who took it directly from M...
Definition: DataClasses.cc:1600

ElementData::crit_cut_length
double crit_cut_length
Critical cut length (m)
Definition: DataClasses.h:366

ElementData::toString
const std::string toString() const
Definition: DataClasses.cc:1827

ElementData::sp
double sp
sphericity: 1 = round, 0 = angular
Definition: DataClasses.h:344

ElementData::concaveNeckRadius
double concaveNeckRadius() const
A non-generic function to compute the concave neck radius (mm).  It is assumed that the neck is bound...
Definition: DataClasses.cc:1629

ElementData::Qmm
double Qmm
Heat source/sink due to phase changes in the case of vapor transport (W/m^3)
Definition: DataClasses.h:384

ElementData::dEps
double dEps
Definition: DataClasses.h:356

LayerData
Parameters of the different layers of the snowpack.
Definition: DataClasses.h:186

LayerData::sp
double sp
Micro-structure : Sphericity.
Definition: DataClasses.h:213

LayerData::rg
double rg
Micro-structure : Grainsize in mm.
Definition: DataClasses.h:212

LayerData::phiIceReservoir
double phiIceReservoir
Volumetric ice reservoir content in %.
Definition: DataClasses.h:200

LayerData::salinity
double salinity
bulk salinity (g/kg)
Definition: DataClasses.h:220

LayerData::phiVoids
double phiVoids
Volumetric void content in %.
Definition: DataClasses.h:204

LayerData::phiSoil
double phiSoil
Volumetric soil content in %.
Definition: DataClasses.h:198

LayerData::phiIce
double phiIce
Volumetric ice content in %.
Definition: DataClasses.h:199

LayerData::dsm
double dsm
dry snow metamorphism factor
Definition: DataClasses.h:222

LayerData::h
double h
capillary pressure head (m)
Definition: DataClasses.h:221

LayerData::cVoids
std::vector< double > cVoids
Solute concentrations in Air.
Definition: DataClasses.h:208

LayerData::CDot
double CDot
Stress rate (Pa s-1), that is the LAST overload change rate.
Definition: DataClasses.h:218

LayerData::depositionDate
mio::Date depositionDate
Date of deposition (mainly used for snow layers)
Definition: DataClasses.h:194

LayerData::cWater
std::vector< double > cWater
Solute concentrations in Water.
Definition: DataClasses.h:207

LayerData::phiIceReservoirCumul
double phiIceReservoirCumul
Volumetric cumulated ice reservoir content in %.
Definition: DataClasses.h:201

LayerData::cIce
std::vector< double > cIce
Solute concentrations in Ice.
Definition: DataClasses.h:206

LayerData::cSoil
std::vector< double > cSoil
Solute concentrations in Soil.
Definition: DataClasses.h:205

LayerData::phiWaterPref
double phiWaterPref
Volumetric preferential water content in %.
Definition: DataClasses.h:203

LayerData::metamo
double metamo
keep track of metamorphism
Definition: DataClasses.h:219

LayerData::toString
const std::string toString() const
Definition: DataClasses.cc:3621

LayerData::operator>>
friend std::istream & operator>>(std::istream &is, LayerData &data)
Definition: DataClasses.cc:3569

LayerData::rb
double rb
Micro-structure : Bond Radius in mm.
Definition: DataClasses.h:215

LayerData::tl
double tl
Temperature at the top of the layer in K.
Definition: DataClasses.h:197

LayerData::hl
double hl
The thickness of the layer in m.
Definition: DataClasses.h:195

LayerData::hr
double hr
Surface hoar Mass in kg m-2.
Definition: DataClasses.h:217

LayerData::SoilC
double SoilC
Heat Capacity of soil.
Definition: DataClasses.h:211

LayerData::SoilK
double SoilK
Conductivity of soil.
Definition: DataClasses.h:210

LayerData::operator<<
friend std::ostream & operator<<(std::ostream &os, const LayerData &data)
Definition: DataClasses.cc:3521

LayerData::SoilRho
double SoilRho
Density of soil in kg m-3.
Definition: DataClasses.h:209

LayerData::LayerData
LayerData()
Definition: DataClasses.cc:3513

LayerData::phiWater
double phiWater
Volumetric water content in %.
Definition: DataClasses.h:202

LayerData::dd
double dd
Micro-structure : Dendricity.
Definition: DataClasses.h:214

LayerData::mk
unsigned short int mk
Micro-structure : Marker.
Definition: DataClasses.h:216

LayerData::ne
size_t ne
Number of finite elements in the the layer (hl/ne defines elm. size)
Definition: DataClasses.h:196

NodeData
NODAL DATA used as a pointer in the SnowStation structure.
Definition: DataClasses.h:392

NodeData::udot
double udot
downward creep velocity in m s-1
Definition: DataClasses.h:404

NodeData::rhov
double rhov
nodal vapor density in kg/m^3
Definition: DataClasses.h:419

NodeData::Sigdsm
double Sigdsm
Definition: DataClasses.h:414

NodeData::S_n
double S_n
Stability Index for natural avalanches.
Definition: DataClasses.h:406

NodeData::dsm
double dsm
Dry snow metamorphism factor.
Definition: DataClasses.h:412

NodeData::S_s
double S_s
Stability Index for skier triggered avalanches.
Definition: DataClasses.h:407

NodeData::hoar
double hoar
Mass of surface hoar collected while node was exposed to surface.
Definition: DataClasses.h:409

NodeData::operator<<
friend std::ostream & operator<<(std::ostream &os, const NodeData &data)
Definition: DataClasses.cc:1850

NodeData::NodeData
NodeData()
Definition: DataClasses.h:394

NodeData::f
double f
reaction or unbalanced forces (CREEP)
Definition: DataClasses.h:403

NodeData::u
double u
creep displacements in m
Definition: DataClasses.h:402

NodeData::ssi
double ssi
Structural Stability Index.
Definition: DataClasses.h:408

NodeData::T
double T
nodal temperature in K
Definition: DataClasses.h:405

NodeData::toString
const std::string toString() const
Definition: DataClasses.cc:1890

NodeData::operator>>
friend std::istream & operator>>(std::istream &is, NodeData &data)
Definition: DataClasses.cc:1870

NodeData::water_flux
double water_flux
Water flowing through the node (kg/m2). Positive values denote downward fluxes.
Definition: DataClasses.h:417

NodeData::S_dsm
double S_dsm
Definition: DataClasses.h:413

NodeData::rime
double rime
Definition: DataClasses.h:415

NodeData::z
double z
nodal height from ground in m
Definition: DataClasses.h:401

RunInfo
class to collect the information about the current simulation (version, date)
Definition: DataClasses.h:847

RunInfo::computation_date
const mio::Date computation_date
Date of computation.
Definition: DataClasses.h:855

RunInfo::version_num
const double version_num
SNOWPACK version formatted as a number.
Definition: DataClasses.h:854

RunInfo::version
const std::string version
SNOWPACK version.
Definition: DataClasses.h:853

RunInfo::hostname
const std::string hostname
hostname of the computer running the simulation
Definition: DataClasses.h:858

RunInfo::RunInfo
RunInfo()
Definition: DataClasses.cc:60

RunInfo::compilation_date
const std::string compilation_date
Date of compilation.
Definition: DataClasses.h:856

RunInfo::user
const std::string user
logname of the user running the simulation
Definition: DataClasses.h:857

RunInfo::operator=
RunInfo & operator=(const RunInfo &)
Definition: DataClasses.h:851

SN_SNOWSOIL_DATA
SN_SNOWSOIL_DATA includes all important station parameters as well as LayerData.
Definition: DataClasses.h:230

SN_SNOWSOIL_DATA::toString
const std::string toString() const
Definition: DataClasses.cc:3450

SN_SNOWSOIL_DATA::TimeCountDeltaHS
double TimeCountDeltaHS
Definition: DataClasses.h:263

SN_SNOWSOIL_DATA::Canopy_alb_dry
double Canopy_alb_dry
Definition: DataClasses.h:269

SN_SNOWSOIL_DATA::profileDate
mio::Date profileDate
Date of profile.
Definition: DataClasses.h:249

SN_SNOWSOIL_DATA::Height
double Height
Total height of soil-snow column in m (sum of the layer heights)
Definition: DataClasses.h:251

SN_SNOWSOIL_DATA::nN
size_t nN
Total number of FE nodes after loading.
Definition: DataClasses.h:250

SN_SNOWSOIL_DATA::ErosionLevel
int ErosionLevel
Erosion Level in operational mode (flat field virtual erosion)
Definition: DataClasses.h:262

SN_SNOWSOIL_DATA::Canopy_Direct_Throughfall
double Canopy_Direct_Throughfall
Direct throughfall [fraction of precipitation].
Definition: DataClasses.h:261

SN_SNOWSOIL_DATA::meta
mio::StationData meta
Station meta data.
Definition: DataClasses.h:248

SN_SNOWSOIL_DATA::operator<<
friend std::ostream & operator<<(std::ostream &os, const SN_SNOWSOIL_DATA &data)
Definition: DataClasses.cc:3384

SN_SNOWSOIL_DATA::Canopy_alb_wet
double Canopy_alb_wet
Definition: DataClasses.h:270

SN_SNOWSOIL_DATA::Canopy_BasalArea
double Canopy_BasalArea
OPTIONAL PARAMETERS, a warning will be thrown in CANOPY::Initialize if no value is provided.
Definition: DataClasses.h:265

SN_SNOWSOIL_DATA::Emissivity_soil
double Emissivity_soil
OPTIONAL PARAMETERS, if not provided Constants::emissivity_soil will be used.
Definition: DataClasses.h:273

SN_SNOWSOIL_DATA::nLayers
size_t nLayers
Total number of soil and snow layers at loading.
Definition: DataClasses.h:252

SN_SNOWSOIL_DATA::Canopy_LAI
double Canopy_LAI
Canopy Leaf Area Index in m2 m-2.
Definition: DataClasses.h:260

SN_SNOWSOIL_DATA::Canopy_diameter
double Canopy_diameter
Average canopy (tree) diameter [m], parameter in the new radiation transfer model.
Definition: DataClasses.h:266

SN_SNOWSOIL_DATA::SoilAlb
double SoilAlb
Soil albedo, default 0.2.
Definition: DataClasses.h:257

SN_SNOWSOIL_DATA::Canopy_Height
double Canopy_Height
Canopy Height in m.
Definition: DataClasses.h:259

SN_SNOWSOIL_DATA::SN_SNOWSOIL_DATA
SN_SNOWSOIL_DATA()
Definition: DataClasses.h:232

SN_SNOWSOIL_DATA::Ldata
std::vector< LayerData > Ldata
contains all the information required to construct the Xdata
Definition: DataClasses.h:253

SN_SNOWSOIL_DATA::Canopy_int_cap_snow
double Canopy_int_cap_snow
Specific interception capacity for rain (I_LAI) (mm/LAI)
Definition: DataClasses.h:268

SN_SNOWSOIL_DATA::BareSoil_z0
double BareSoil_z0
Bare soil roughness in m, default 0.02 m.
Definition: DataClasses.h:258

SN_SNOWSOIL_DATA::Canopy_lai_frac_top_default
double Canopy_lai_frac_top_default
fraction of total LAI that is attributed to the uppermost layer. Here calibrated for Alptal.
Definition: DataClasses.h:267

SN_SNOWSOIL_DATA::Canopy_alb_snow
double Canopy_alb_snow
Definition: DataClasses.h:271

SN_SNOWSOIL_DATA::operator>>
friend std::istream & operator>>(std::istream &is, SN_SNOWSOIL_DATA &data)
Definition: DataClasses.cc:3416

SN_SNOWSOIL_DATA::Albedo
double Albedo
REQUIRED PARAMETERS, an error will be thrown at reading (SnowpackIO) if no parameter are provided.
Definition: DataClasses.h:256

SN_SNOWSOIL_DATA::HS_last
double HS_last
Definition: DataClasses.h:254

SeaIce
Definition: SeaIce.h:42

SnowProfileLayer
Defines structure for snow profile layers.
Definition: DataClasses.h:808

SnowProfileLayer::theta_w
double theta_w
0 to 1, volume fraction of water (-)
Definition: DataClasses.h:828

SnowProfileLayer::height
double height
Height of snow or snow depth; 0 to 1000 (cm)
Definition: DataClasses.h:822

SnowProfileLayer::theta_a
double theta_a
0 to 1, volume fraction of air (-)
Definition: DataClasses.h:829

SnowProfileLayer::average
void average(const double &w1, const double &w2, const SnowProfileLayer &Pdata)
Determines the averaged quantities of the current layer with another layer.
Definition: DataClasses.cc:248

SnowProfileLayer::T
double T
-50 to 50, snow temperature at top of layer (degC)
Definition: DataClasses.h:824

SnowProfileLayer::hard
double hard
0. to 5. (1)
Definition: DataClasses.h:838

SnowProfileLayer::profileDate
mio::Date profileDate
Date of profile.
Definition: DataClasses.h:816

SnowProfileLayer::depositionDate
mio::Date depositionDate
Date of deposition (mainly used for snow layers)
Definition: DataClasses.h:821

SnowProfileLayer::marker
unsigned short int marker
0 to 999 (1)
Definition: DataClasses.h:836

SnowProfileLayer::SnowProfileLayer
SnowProfileLayer()
Definition: DataClasses.cc:131

SnowProfileLayer::generateProfile
static std::vector< SnowProfileLayer > generateProfile(const mio::Date &dateOfProfile, const SnowStation &Xdata, const double hoar_density_surf, const double hoar_min_size_surf)
Generates a snow profile from snow station data (1 element = 1 layer)
Definition: DataClasses.cc:203

SnowProfileLayer::v_strain_rate
double v_strain_rate
0 to 1.0e-5, viscous strain rate (s-1)
Definition: DataClasses.h:826

SnowProfileLayer::stationname
std::string stationname
Definition: DataClasses.h:817

SnowProfileLayer::loc_for_wind
unsigned char loc_for_wind
Definition: DataClasses.h:819

SnowProfileLayer::coordin_num
double coordin_num
0 to 10 (1)
Definition: DataClasses.h:835

SnowProfileLayer::loc_for_snow
unsigned char loc_for_snow
Definition: DataClasses.h:818

SnowProfileLayer::ogs
double ogs
0 to 100, optical equivalent grain size (mm)
Definition: DataClasses.h:834

SnowProfileLayer::rho
double rho
0 to 1000 (kg m-3)
Definition: DataClasses.h:823

SnowProfileLayer::type
short unsigned int type
0 to 999 (1)
Definition: DataClasses.h:837

SnowProfileLayer::sphericity
double sphericity
0 to 1 (1)
Definition: DataClasses.h:833

SnowProfileLayer::grain_size
double grain_size
0 to 100 (mm)
Definition: DataClasses.h:830

SnowProfileLayer::gradT
double gradT
-1000 to 1000, temperature gradient across layer (K m-1)
Definition: DataClasses.h:825

SnowProfileLayer::dendricity
double dendricity
0 to 1 (1)
Definition: DataClasses.h:832

SnowProfileLayer::bond_size
double bond_size
0 to 100 (mm)
Definition: DataClasses.h:831

SnowProfileLayer::theta_i
double theta_i
0 to 1, volume fraction of ice (-)
Definition: DataClasses.h:827

SnowStation
Definition: DataClasses.h:603

SnowStation::combineElements
void combineElements(const size_t &number_top_elements, const bool &reduce_n_elements, const size_t &cond, const double &comb_thresh_l)
If more than NUMBER_TOP_ELEMENTS snow elements exist, attempt to reduce their number in the FEM mesh,...
Definition: DataClasses.cc:2238

SnowStation::findMarkedReferenceLayer
double findMarkedReferenceLayer() const
returns the height of a marked reference layer inside the model domain Searches for the layer that is...
Definition: DataClasses.cc:2844

SnowStation::comb_thresh_sp
static const double comb_thresh_sp
Sphericity (1)
Definition: DataClasses.h:695

SnowStation::SoilNode
size_t SoilNode
The top soil node, 0 in case of SNP_SOIL == 0.
Definition: DataClasses.h:656

SnowStation::getNumberOfNodes
size_t getNumberOfNodes() const
Definition: DataClasses.cc:2190

SnowStation::hasSoilLayers
bool hasSoilLayers() const
Definition: DataClasses.cc:2211

SnowStation::z_S_n
double z_S_n
Depth of Minimum S_n.
Definition: DataClasses.h:675

SnowStation::initialize
void initialize(const SN_SNOWSOIL_DATA &SSdata, const size_t &i_sector)
This routine initializes the snow cover structure which contains all information about a station incl...
Definition: DataClasses.cc:2349

SnowStation::getNumberOfElements
size_t getNumberOfElements() const
Definition: DataClasses.cc:2185

SnowStation::S_class1
char S_class1
Stability class based on hand hardness, grain class ...
Definition: DataClasses.h:670

SnowStation::thresh_moist_soil
static const double thresh_moist_soil
Definition: DataClasses.h:696

SnowStation::operator>>
friend std::istream & operator>>(std::istream &is, SnowStation &data)
Definition: DataClasses.cc:2937

SnowStation::S_n
double S_n
Minimum natural stability index.
Definition: DataClasses.h:674

SnowStation::S_s
double S_s
Minimum skier stability index sk38 (Sk38)
Definition: DataClasses.h:676

SnowStation::hn
double hn
Depth of new snow to be used on slopes.
Definition: DataClasses.h:665

SnowStation::meta
mio::StationData meta
Station meta data.
Definition: DataClasses.h:645

SnowStation::find_tag
size_t find_tag(const size_t &tag) const
Find element with corresponding tag or return IOUtils::npos if not found.
Definition: DataClasses.cc:2200

SnowStation::SoilAlb
double SoilAlb
Soil albedo.
Definition: DataClasses.h:653

SnowStation::rime_hn
double rime_hn
rime of new snow to be used on slopes
Definition: DataClasses.h:667

SnowStation::ErosionMass
double ErosionMass
Eroded mass either real or virtually (storage if less than one element)
Definition: DataClasses.h:669

SnowStation::number_top_elements
static const size_t number_top_elements
Number of top elements left untouched by the join functions.
Definition: DataClasses.h:697

SnowStation::Seaice
SeaIce * Seaice
Pointer to sea ice class.
Definition: DataClasses.h:650

SnowStation::S_d
double S_d
Minimum deformation rate stability index.
Definition: DataClasses.h:672

SnowStation::swc_sum_soil
double swc_sum_soil
Total solid water in soil.
Definition: DataClasses.h:664

SnowStation::dIntEnergy
double dIntEnergy
Internal energy change of snowpack (J m-2)
Definition: DataClasses.h:687

SnowStation::Kt
void * Kt
Pointer to pseudo-conductivity and stiffnes matrix.
Definition: DataClasses.h:684

SnowStation::cos_sl
double cos_sl
Cosinus of slope angle, initialized once!
Definition: DataClasses.h:646

SnowStation::resize
void resize(const size_t &number_of_elements)
Reallocate element and node data  Edata and Ndata as well as nElems and nNodes are reallocated or res...
Definition: DataClasses.cc:2165

SnowStation::z_S_4
double z_S_4
Depth of Minimum S_4.
Definition: DataClasses.h:679

SnowStation::rho_hn
double rho_hn
Density of new snow to be used on slopes.
Definition: DataClasses.h:666

SnowStation::~SnowStation
~SnowStation()
Definition: DataClasses.cc:2003

SnowStation::z_S_s
double z_S_s
Depth of Minimum S_s.
Definition: DataClasses.h:677

SnowStation::meltFreezeEnergy
double meltFreezeEnergy
Melt freeze part of internal energy change of snowpack (J m-2)
Definition: DataClasses.h:689

SnowStation::getTotalLateralFlowSoil
double getTotalLateralFlowSoil() const
Computes the total lateral flow in soil (kg / m2)
Definition: DataClasses.cc:2137

SnowStation::lwc_sum_soil
double lwc_sum_soil
Total liquid water in soil.
Definition: DataClasses.h:663

SnowStation::ColdContent
double ColdContent
Cold content of snowpack (J m-2)
Definition: DataClasses.h:685

SnowStation::S_4
double S_4
placeholder - currently Minimum structural stability index (SSI)
Definition: DataClasses.h:678

SnowStation::windward
bool windward
True for windward (luv) slope.
Definition: DataClasses.h:692

SnowStation::S_5
double S_5
placeholder
Definition: DataClasses.h:680

SnowStation::splitElement
void splitElement(const size_t &e)
Split the element provided as first argument.
Definition: DataClasses.cc:2628

SnowStation::ColdContentSoil
double ColdContentSoil
Cold content of soil (J m-2)
Definition: DataClasses.h:686

SnowStation::getModelledTemperature
double getModelledTemperature(const double &z) const
Returns modelled internal snow or/and soil temperature (instantaneous value; degC),...
Definition: DataClasses.cc:2104

SnowStation::mass_sum
double mass_sum
Total mass summing mass of snow elements.
Definition: DataClasses.h:660

SnowStation::swe
double swe
Total mass summing snow water equivalent of elements.
Definition: DataClasses.h:661

SnowStation::comb_thresh_rg
static const double comb_thresh_rg
Grain radius (mm)
Definition: DataClasses.h:695

SnowStation::compSnowpackMasses
void compSnowpackMasses()
Computes the actual total masses of the snowpack (kg m-2)
Definition: DataClasses.cc:2026

SnowStation::reduceNumberOfElements
void reduceNumberOfElements(const size_t &rnE)
Remove the upper "marked" element of two (snow only)  .
Definition: DataClasses.cc:2284

SnowStation::z_S_d
double z_S_d
Depth of Minimum S_d.
Definition: DataClasses.h:673

SnowStation::Ndata
std::vector< NodeData > Ndata
pointer to nodal data array (e.g. T, z, u, etc..)
Definition: DataClasses.h:682

SnowStation::mergeElements
static void mergeElements(ElementData &Edata0, const ElementData &Edata1, const bool &merge, const bool &topElement)
Merging two elements.
Definition: DataClasses.cc:2709

SnowStation::ReSolver_dt
double ReSolver_dt
Last used RE time step in the previous SNOWPACK time step.
Definition: DataClasses.h:691

SnowStation::number_of_solutes
static unsigned short number_of_solutes
The model treats that number of solutes.
Definition: DataClasses.h:698

SnowStation::BareSoil_z0
double BareSoil_z0
Bare soil roughness in m.
Definition: DataClasses.h:655

SnowStation::comb_thresh_dd
static const double comb_thresh_dd
Dendricity (1)
Definition: DataClasses.h:695

SnowStation::reset_water_fluxes
void reset_water_fluxes()
Definition: DataClasses.cc:2216

SnowStation::mH
double mH
The MEASURED height, including soil depth if SNP_SOIL == 1.
Definition: DataClasses.h:659

SnowStation::operator=
SnowStation & operator=(const SnowStation &)
Assignement operator.
Definition: DataClasses.cc:1938

SnowStation::Edata
std::vector< ElementData > Edata
pointer to element data array (e.g. Te, L, Rho, etc..)
Definition: DataClasses.h:683

SnowStation::Albedo
double Albedo
Snow albedo used by the model.
Definition: DataClasses.h:652

SnowStation::dIntEnergySoil
double dIntEnergySoil
Internal energy change of soil (J m-2)
Definition: DataClasses.h:688

SnowStation::z_S_5
double z_S_5
Depth of Minimum S_5.
Definition: DataClasses.h:681

SnowStation::operator<<
friend std::ostream & operator<<(std::ostream &os, const SnowStation &data)
Definition: DataClasses.cc:2857

SnowStation::getTotalLateralFlowSnow
double getTotalLateralFlowSnow() const
Computes the total lateral flow in snow (kg / m2)
Definition: DataClasses.cc:2122

SnowStation::comb_thresh_water
static const double comb_thresh_water
Water content (1)
Definition: DataClasses.h:694

SnowStation::sector
size_t sector
current slope sector of width 360./max(1, nSlopes-1)
Definition: DataClasses.h:647

SnowStation::Ground
double Ground
The ground height – meaning the height of the top soil node.
Definition: DataClasses.h:657

SnowStation::toString
const std::string toString() const
Definition: DataClasses.cc:3022

SnowStation::resetSlopeParFlux
void resetSlopeParFlux()
Reset lateral flow.
Definition: DataClasses.cc:2152

SnowStation::comb_thresh_ice
static const double comb_thresh_ice
Volumetric ice content (1), i.e., about 46 kg m-3.
Definition: DataClasses.h:694

SnowStation::SnowStation
SnowStation(const bool i_useCanopyModel=true, const bool i_useSoilLayers=true, const bool i_isAlpine3D=false, const bool i_useSeaIceModule=false)
Definition: DataClasses.cc:1904

SnowStation::S_class2
char S_class2
Stability class based on hand hardness, grain class ...
Definition: DataClasses.h:671

SnowStation::getLiquidWaterIndex
double getLiquidWaterIndex() const
Computes the liquid water index defined as the ratio of total liquid water content (in mm w....
Definition: DataClasses.cc:2092

SnowStation::compSoilInternalEnergyChange
void compSoilInternalEnergyChange(const double &sn_dt)
Computes the internal energy change of the soil during one computation time step (J m-2)
Definition: DataClasses.cc:2067

SnowStation::isGlacier
bool isGlacier(const bool &hydro=false) const
returns if a snow profile can be considered as a glacier. Practically, the hydrological criteria is t...
Definition: DataClasses.cc:2803

SnowStation::SoilEmissivity
double SoilEmissivity
Soil emissivity.
Definition: DataClasses.h:654

SnowStation::TimeCountDeltaHS
double TimeCountDeltaHS
Time counter tracking erroneous settlement in operational mode.
Definition: DataClasses.h:693

SnowStation::lwc_sum
double lwc_sum
Total liquid water in snowpack.
Definition: DataClasses.h:662

SnowStation::comb_thresh_l_ratio
static const double comb_thresh_l_ratio
The default ratio between height_new_elem and comb_thresh_l, in case comb_thresh_l is not explicitly ...
Definition: DataClasses.h:694

SnowStation::cH
double cH
The CALCULATED height, including soil depth if SNP_SOIL == 1.
Definition: DataClasses.h:658

SnowStation::meltFreezeEnergySoil
double meltFreezeEnergySoil
Melt freeze part of internal energy change of soil (J m-2)
Definition: DataClasses.h:690

SnowStation::pAlbedo
double pAlbedo
Parameterized snow albedo.
Definition: DataClasses.h:651

SnowStation::splitElements
void splitElements(const double &max_element_length, const double &comb_thresh_l)
Split elements when they are near the top of the snowpack, when REDUCE_N_ELEMENTS is used.
Definition: DataClasses.cc:2674

SnowStation::ErosionLevel
size_t ErosionLevel
Element where snow erosion stopped previously for the drift index.
Definition: DataClasses.h:668

SnowStation::combineCondition
static bool combineCondition(const ElementData &Edata0, const ElementData &Edata1, const double &depth, const bool &reduce_n_elements, const double &comb_thresh_l)
Boolean routine to check whether two snow elements can be combined.
Definition: DataClasses.cc:2578

SnowStation::compSnowpackInternalEnergyChange
void compSnowpackInternalEnergyChange(const double &sn_dt)
Computes the internal energy change of the snowpack during one computation time step (J m-2)
Definition: DataClasses.cc:2045

SnowStation::Cdata
CanopyData Cdata
Pointer to canopy data.
Definition: DataClasses.h:649

SnowStation::thresh_moist_snow
static const double thresh_moist_snow
Snow elements with a LWC above this threshold are considered at least to be moist.
Definition: DataClasses.h:696

SnowpackConfig
Definition: SnowpackConfig.h:28

SurfaceFluxes
Definition: DataClasses.h:731

SurfaceFluxes::SN_MASS_CHANGES
SN_MASS_CHANGES
The different types of mass fluxes: Mass fluxes in kg m-2  Rates in kg m-2 h-1 (MS_HNW,...
Definition: DataClasses.h:738

SurfaceFluxes::MS_WATER
@ MS_WATER
The total amount of water in the snowpack at the present time.
Definition: DataClasses.h:741

SurfaceFluxes::MS_SUBLIMATION
@ MS_SUBLIMATION
The mass loss or gain of the top element due to snow (ice) sublimating.
Definition: DataClasses.h:748

SurfaceFluxes::MS_TOTALMASS
@ MS_TOTALMASS
This of course is the total mass of the snowpack at the present time.
Definition: DataClasses.h:739

SurfaceFluxes::MS_HNW
@ MS_HNW
Solid precipitation rate.
Definition: DataClasses.h:744

SurfaceFluxes::MS_WIND
@ MS_WIND
Mass loss rate due to wind erosion.
Definition: DataClasses.h:746

SurfaceFluxes::MS_FLOODING
@ MS_FLOODING
The mass gain due to adding ocean water to snow- seaice by flodding process.
Definition: DataClasses.h:752

SurfaceFluxes::MS_ICE_SOIL
@ MS_ICE_SOIL
The total amount of ice in the soil at the present time.
Definition: DataClasses.h:743

SurfaceFluxes::MS_SWE
@ MS_SWE
This too, of course, but summing rho*L.
Definition: DataClasses.h:740

SurfaceFluxes::MS_SOIL_RUNOFF
@ MS_SOIL_RUNOFF
Equivalent to MS_SNOWPACK_RUNOFF but at bottom soil node.
Definition: DataClasses.h:751

SurfaceFluxes::MS_SNOWPACK_RUNOFF
@ MS_SNOWPACK_RUNOFF
The mass loss of snowpack from snow melt due to water transport (virtual lysimeter)
Definition: DataClasses.h:749

SurfaceFluxes::MS_WATER_SOIL
@ MS_WATER_SOIL
The total amount of water in the soil at the present time.
Definition: DataClasses.h:742

SurfaceFluxes::MS_EVAPORATION
@ MS_EVAPORATION
The mass loss or gain of the top element due to water evaporating.
Definition: DataClasses.h:747

SurfaceFluxes::N_MASS_CHANGES
@ N_MASS_CHANGES
Total number of different mass change types.
Definition: DataClasses.h:754

SurfaceFluxes::MS_RAIN
@ MS_RAIN
Rain rate.
Definition: DataClasses.h:745

SurfaceFluxes::MS_SURFACE_MASS_FLUX
@ MS_SURFACE_MASS_FLUX
The total mass loss of snowpack due to water transport (virtual lysimeter)
Definition: DataClasses.h:750

SurfaceFluxes::MS_ICEBASE_MELTING_FREEZING
@ MS_ICEBASE_MELTING_FREEZING
The mass gain/loss of the ice base due to melting-freezing.
Definition: DataClasses.h:753

SurfaceFluxes::load
std::vector< double > load
Total load (kg m-2) in water runoff from solutes like nitrate.
Definition: DataClasses.h:798

SurfaceFluxes::operator<<
friend std::ostream & operator<<(std::ostream &os, const SurfaceFluxes &data)
Definition: DataClasses.cc:454

SurfaceFluxes::lw_in
double lw_in
Energy fluxes: Energy change of snowpack in kJ m-2 (dIntEnergy) Fluxes in W m-2.
Definition: DataClasses.h:773

SurfaceFluxes::sw_dir
double sw_dir
incoming direct shortwave radiation; on slopes: projected
Definition: DataClasses.h:786

SurfaceFluxes::sw_hor
double sw_hor
incoming global shortwave radiation on horizontal surface
Definition: DataClasses.h:782

SurfaceFluxes::sw_diff
double sw_diff
incoming diffuse shortwave radiation
Definition: DataClasses.h:787

SurfaceFluxes::mAlbedo
double mAlbedo
measured Albedo (used for OUTPUT only)
Definition: DataClasses.h:789

SurfaceFluxes::cRho_hn
double cRho_hn
Computed new snow density (kg m-3)
Definition: DataClasses.h:800

SurfaceFluxes::SurfaceFluxes
SurfaceFluxes()
Definition: DataClasses.cc:292

SurfaceFluxes::operator>>
friend std::istream & operator>>(std::istream &is, SurfaceFluxes &data)
Definition: DataClasses.cc:494

SurfaceFluxes::collectSurfaceFluxes
void collectSurfaceFluxes(const BoundCond &Bdata, SnowStation &Xdata, const CurrentMeteo &Mdata)
Assign surface data from SnowStation and BoundCond to SurfaceFluxes.
Definition: DataClasses.cc:361

SurfaceFluxes::qr
double qr
rain energy
Definition: DataClasses.h:779

SurfaceFluxes::mass
std::vector< double > mass
Total mass of snowpack PLUS different amounts of total mass change, sublimation, runoff,...
Definition: DataClasses.h:797

SurfaceFluxes::lw_net
double lw_net
net longwave radiation
Definition: DataClasses.h:775

SurfaceFluxes::qg
double qg
geothermal heat flux or heat flux at lower boundary
Definition: DataClasses.h:780

SurfaceFluxes::sw_in
double sw_in
incoming global shortwave radiation; on slopes: projected
Definition: DataClasses.h:783

SurfaceFluxes::reset
void reset(const bool &cumsum_mass)
Definition: DataClasses.cc:297

SurfaceFluxes::drift
double drift
Other surface data:
Definition: DataClasses.h:796

SurfaceFluxes::pAlbedo
double pAlbedo
parameterized Albedo (used for OUTPUT only)
Definition: DataClasses.h:788

SurfaceFluxes::lw_out
double lw_out
outgoing longwave radiation
Definition: DataClasses.h:774

SurfaceFluxes::meltFreezeEnergySoil
double meltFreezeEnergySoil
Melt freeze part of internal energy change in J m-2 in soil (used for OUTPUT only)
Definition: DataClasses.h:793

SurfaceFluxes::hoar
double hoar
mass of surface hoar formed or sublimated
Definition: DataClasses.h:778

SurfaceFluxes::dIntEnergySoil
double dIntEnergySoil
Internal energy change in J m-2 in soil (used for OUTPUT only)
Definition: DataClasses.h:791

SurfaceFluxes::dIntEnergy
double dIntEnergy
Internal energy change in J m-2 in snowpack (used for OUTPUT only)
Definition: DataClasses.h:790

SurfaceFluxes::sw_out
double sw_out
reflected shortwave radiation
Definition: DataClasses.h:784

SurfaceFluxes::qw
double qw
net shortwave radiation at the surface (absorbed within the snowpack)
Definition: DataClasses.h:785

SurfaceFluxes::toString
const std::string toString() const
Definition: DataClasses.cc:3484

SurfaceFluxes::multiplyFluxes
void multiplyFluxes(const double &factor)
If multiple surface fluxes have been summed over multiple time steps, the fluxes then need to be aver...
Definition: DataClasses.cc:432

SurfaceFluxes::compSnowSoilHeatFlux
void compSnowSoilHeatFlux(const SnowStation &Xdata)
Compute ground heat flux at soil/snow boundary.
Definition: DataClasses.cc:331

SurfaceFluxes::dhs_corr
double dhs_corr
operational mode only: snow depth correction in case of squeezing or blow-up (m)
Definition: DataClasses.h:799

SurfaceFluxes::qg0
double qg0
ground heat flux at soil-snow interface
Definition: DataClasses.h:781

SurfaceFluxes::qs
double qs
sensible heat
Definition: DataClasses.h:776

SurfaceFluxes::mRho_hn
double mRho_hn
Measured new snow density (kg m-3)
Definition: DataClasses.h:801

SurfaceFluxes::meltFreezeEnergy
double meltFreezeEnergy
Melt freeze part of internal energy change in J m-2 in snowpack (used for OUTPUT only)
Definition: DataClasses.h:792

SurfaceFluxes::ql
double ql
latent heat
Definition: DataClasses.h:777

ZwischenData
ZwischenData contains "memory" information mainly for operational use It is used to prepare some para...
Definition: DataClasses.h:65

ZwischenData::ZwischenData
ZwischenData()
Definition: DataClasses.h:67

ZwischenData::reset
void reset()
Sets all the values in the vectors to 0.0.
Definition: DataClasses.cc:92

ZwischenData::hn24
std::vector< double > hn24
Twenty-four hour snow heights every half-hour over three days 144.
Definition: DataClasses.h:76

ZwischenData::hn3
std::vector< double > hn3
Three hour new snow heights every half-hour over three days 144.
Definition: DataClasses.h:75

ZwischenData::drift24
std::vector< double > drift24
Twenty-four hour hoar index every half-hour over one day 48.
Definition: DataClasses.h:74

ZwischenData::operator<<
friend std::ostream & operator<<(std::ostream &os, const ZwischenData &data)
Definition: DataClasses.cc:100

ZwischenData::hoar24
std::vector< double > hoar24
Twenty-four hour hoar index every half-hour over one day 48.
Definition: DataClasses.h:73

ZwischenData::operator>>
friend std::istream & operator>>(std::istream &is, ZwischenData &data)
Definition: DataClasses.cc:114

vanGenuchten
This module contains the van Genuchten model for the water retention curve.
Definition: vanGenuchten.h:36

Constants
Definition: Constants.h:43

Constants::undefined
const double undefined
This is the snowpack undefined value.
Definition: Constants.h:44

SNOW_OPTIC
Optical Properties of snow.
Definition: DataClasses.h:48

SNOW_OPTIC::ssa
double ssa
Single Scattering Albedo.
Definition: DataClasses.h:51

SNOW_OPTIC::exteff
double exteff
Extinction Efficiency.
Definition: DataClasses.h:50

SNOW_OPTIC::ggg
double ggg
Asymmetry Parameter.
Definition: DataClasses.h:49

WL_STRUCT
Spectrum of incoming solar radiation.
Definition: DataClasses.h:55

WL_STRUCT::perc
double perc
Percentage of Energy.
Definition: DataClasses.h:57

WL_STRUCT::nm
double nm
Wavelength.
Definition: DataClasses.h:56

vanGenuchten.h