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>

class SeaIce;   //forward declaration to prevent include loop


#include <snowpack/Constants.h>

#include <meteoio/MeteoIO.h>


#include <string>

#include <vector>

#include <functional>


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.), rta(0.), 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 rta;

                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, const bool& VapourTransport = false);

                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

                bool needDeflateInflate(const double &sn_dt);

                void deflateInflate(const CurrentMeteo& Mdata, double& dhs_corr, double& mass_corr, const bool &prn_check);

                double forcedErosion(const double hs);


                void compSnowpackMasses();

                void compSnowpackInternalEnergyChange(const double& sn_dt);

                void compSoilInternalEnergyChange(const double& sn_dt);

                double averageFromTop(const double& averagingDepth, const std::function<double(const ElementData&)>& Elementproperty);


                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;

                const std::string history;


        private:

                std::string setHistory();

                static int hexToDecimal(const std::string& hex);

                static double getNumericVersion(std::string version_str);

                static mio::Date getRunDate();

                static std::string getCompilationDate();

};


class Slope {


        public:

                Slope(const mio::Config& cfg);


                double prevailing_wind_dir;

                unsigned int nSlopes;

                unsigned int mainStation;

                unsigned int sector;

                unsigned int first;

                unsigned int luv;

                unsigned int lee;

                bool north, south;

                bool snow_erosion, mainStationDriftIndex;

                bool snow_redistribution, luvDriftIndex;


                unsigned int getSectorDir(const double& dir_or_expo) const;

                void setSlope(const unsigned int slope_sequence, std::vector<SnowStation>& vecXdata, double& wind_dir);


        private:

                double sector_width;

};


class MainControl {

        public:

                MainControl(const SnowpackConfig& cfg);


                void getOutputControl(const mio::Date& step, const mio::Date& sno_step);

                void reset();


                // New members for output control

                double tsstart;

                double tsdaysbetween;

                double profstart;

                double profdaysbetween;

                double first_backup;

                double backup_days_between;

                double calculation_step_length;

                double sn_dt;


                size_t nStep;

                size_t nAvg;

                size_t HzStep;

                bool   TsDump;

                bool   HzDump;

                bool   PrDump;

                bool   XdataDump;

                bool   sdbDump;

                bool   resFirstDump;

};


class Cumsum {


        public:


                Cumsum(const unsigned int nSlopes) : precip(0.),

          drift(0.), snow(0.), runoff(0.), rain(0.), dhs_corr(0.), mass_corr(0.),

          erosion(nSlopes, 0.) {}


                double precip;

                double drift, snow, runoff, rain;

                double dhs_corr, mass_corr; // inflate/deflate variables

                std::vector<double> erosion; // Cumulated eroded mass; dumped to file as rate

};


#endif

Constants.h

SOLID
@ SOLID
Solid.
Definition DataClasses.h:42

GAS
@ GAS
Gas.
Definition DataClasses.h:44

LIQUID
@ LIQUID
Liquid.
Definition DataClasses.h:43

NUMBER_OF_PHASES
@ NUMBER_OF_PHASES
Definition DataClasses.h:45

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

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

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

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

N_SOIL_FIELDS
@ N_SOIL_FIELDS
Definition DataClasses.h:181

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

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

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

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

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

N_SN_FIELDS
@ N_SN_FIELDS
Definition DataClasses.h:163

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

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

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

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

AIR
@ AIR
Air.
Definition DataClasses.h:172

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

SOIL
@ SOIL
Soil.
Definition DataClasses.h:168

ICE
@ ICE
Ice.
Definition DataClasses.h:169

N_COMPONENTS
@ N_COMPONENTS
Definition DataClasses.h:173

WATER
@ WATER
Water.
Definition DataClasses.h:170

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

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

SnowpackConfig.h

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:844

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:591

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Cumsum
Definition DataClasses.h:939

Cumsum::precip
double precip
Definition DataClasses.h:946

Cumsum::rain
double rain
Definition DataClasses.h:947

Cumsum::dhs_corr
double dhs_corr
Definition DataClasses.h:948

Cumsum::Cumsum
Cumsum(const unsigned int nSlopes)
Definition DataClasses.h:942

Cumsum::snow
double snow
Definition DataClasses.h:947

Cumsum::mass_corr
double mass_corr
Definition DataClasses.h:948

Cumsum::drift
double drift
Definition DataClasses.h:947

Cumsum::runoff
double runoff
Definition DataClasses.h:947

Cumsum::erosion
std::vector< double > erosion
Definition DataClasses.h:949

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

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

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

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:137

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

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

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

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

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

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

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

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

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:146

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

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:131

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:115

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:136

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

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

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

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

ElementData
The physical properties of an individual element used in snow layers (one or multiple elements are co...
Definition DataClasses.h:288

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

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

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

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:379

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

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

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

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

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

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

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

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

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:1718

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

ElementData::soilRelativeHumidity
double soilRelativeHumidity() const
Relative humidity in soil.
Definition DataClasses.cc:1625

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

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

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

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:360

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:378

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

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

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

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

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

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

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

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

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

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

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:342

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

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

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

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

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

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

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

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

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

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

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:341

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

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

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

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

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

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:1674

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

MainControl
Definition DataClasses.h:907

MainControl::backup_days_between
double backup_days_between
Interval (days) between Xdata backups.
Definition DataClasses.h:920

MainControl::tsdaysbetween
double tsdaysbetween
Interval (days) between time series dumps.
Definition DataClasses.h:916

MainControl::sdbDump
bool sdbDump
Dump to data base if required in operational mode.
Definition DataClasses.h:931

MainControl::HzDump
bool HzDump
Calculation of hazard information will be performed.
Definition DataClasses.h:928

MainControl::reset
void reset()
Reset all fields not read from the config file to zero or false.
Definition DataClasses.cc:4054

MainControl::resFirstDump
bool resFirstDump
Flag to dump initial state of snowpack.
Definition DataClasses.h:932

MainControl::profdaysbetween
double profdaysbetween
Interval (days) between profile dumps.
Definition DataClasses.h:918

MainControl::calculation_step_length
double calculation_step_length
Length of each calculation step (minutes)
Definition DataClasses.h:921

MainControl::PrDump
bool PrDump
Flag for profile dump.
Definition DataClasses.h:929

MainControl::sn_dt
double sn_dt
calculation step length in seconds
Definition DataClasses.h:922

MainControl::XdataDump
bool XdataDump
Backup of Xdata will be performed.
Definition DataClasses.h:930

MainControl::nStep
size_t nStep
Time step number.
Definition DataClasses.h:924

MainControl::TsDump
bool TsDump
Flag for time series dump.
Definition DataClasses.h:927

MainControl::tsstart
double tsstart
Start time for time series dump (Julian)
Definition DataClasses.h:915

MainControl::nAvg
size_t nAvg
Number of calculation time steps to average fluxes etc.
Definition DataClasses.h:925

MainControl::first_backup
double first_backup
Start time for first Xdata backup (Julian)
Definition DataClasses.h:919

MainControl::getOutputControl
void getOutputControl(const mio::Date &step, const mio::Date &sno_step)
Determine which outputs need to be done for the current time step.
Definition DataClasses.cc:4027

MainControl::profstart
double profstart
Start time for profile dump (Julian)
Definition DataClasses.h:917

MainControl::HzStep
size_t HzStep
Hazard step number (should be half of nStep in operational mode)
Definition DataClasses.h:926

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

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

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

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

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

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

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

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

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

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

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

NodeData::rta
double rta
Relative Threshold Sum apprach.
Definition DataClasses.h:413

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

RunInfo::history
const std::string history
History string, like in ACDD, combining several of the information below.
Definition DataClasses.h:871

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

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

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

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

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

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

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

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

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

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

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

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

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:268

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

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

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

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

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

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

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

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

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

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

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:270

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

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

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

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

SeaIce
Definition SeaIce.h:41

Slope
Slope a C. Fierz class ;-) to handle the multiple expositions in Snowpack standalone.
Definition DataClasses.h:884

Slope::luv
unsigned int luv
Definition DataClasses.h:894

Slope::south
bool south
Definition DataClasses.h:896

Slope::snow_redistribution
bool snow_redistribution
Definition DataClasses.h:898

Slope::mainStation
unsigned int mainStation
main station, flat field or slope
Definition DataClasses.h:891

Slope::getSectorDir
unsigned int getSectorDir(const double &dir_or_expo) const
Determine either direction of blowing wind or slope exposition. NOTE that station slope....
Definition DataClasses.cc:3958

Slope::snow_erosion
bool snow_erosion
Definition DataClasses.h:897

Slope::luvDriftIndex
bool luvDriftIndex
Definition DataClasses.h:898

Slope::setSlope
void setSlope(const unsigned int slope_sequence, std::vector< SnowStation > &vecXdata, double &wind_dir)
Set slope variables.
Definition DataClasses.cc:3974

Slope::nSlopes
unsigned int nSlopes
Definition DataClasses.h:890

Slope::north
bool north
Definition DataClasses.h:896

Slope::lee
unsigned int lee
Definition DataClasses.h:895

Slope::first
unsigned int first
first virtual slope station in computing sequence
Definition DataClasses.h:893

Slope::mainStationDriftIndex
bool mainStationDriftIndex
Definition DataClasses.h:897

Slope::prevailing_wind_dir
double prevailing_wind_dir
Definition DataClasses.h:889

Slope::sector
unsigned int sector
main station (0) or current slope sector (1:nSlopes)
Definition DataClasses.h:892

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

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

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

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

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:291

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

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

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

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

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

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

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:246

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

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

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

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

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

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

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

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

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

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

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

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

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

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

SnowStation
Station data including all information on snowpack layers (elements and nodes) and on canopy.
Definition DataClasses.h:609

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:2326

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:3138

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

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

SnowStation::deflateInflate
void deflateInflate(const CurrentMeteo &Mdata, double &dhs_corr, double &mass_corr, const bool &prn_check)
Deflates or inflates the snowpack for warning service purposes.
Definition DataClasses.cc:2964

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

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

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

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:2439

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

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

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

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

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

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

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

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

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

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:2288

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

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

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

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

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

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

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

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

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

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

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:2253

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

SnowStation::averageFromTop
double averageFromTop(const double &averagingDepth, const std::function< double(const ElementData &)> &Elementproperty)
Extracts a weighted average of a snow bulk property from the top down to a given depth.
Definition DataClasses.cc:2154

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

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

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

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

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

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

SnowStation::needDeflateInflate
bool needDeflateInflate(const double &sn_dt)
Check if deflateInfate should be called.
Definition DataClasses.cc:2912

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

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:2180

SnowStation::forcedErosion
double forcedErosion(const double hs)
Forced erosion of a missed event.
Definition DataClasses.cc:3072

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:2114

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:3097

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

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

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

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:705

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

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

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

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:2770

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

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:2671

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:2092

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

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:707

SnowpackConfig
Definition SnowpackConfig.h:28

SurfaceFluxes
The surface fluxes data as computed by Snowpack.
Definition DataClasses.h:743

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:750

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:785

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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:475

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

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

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

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

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

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

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

VapourTransport
This module contains water vapour transport routines for the 1d snowpack model.
Definition VapourTransport.h:50

ZwischenData
ZwischenData contains "memory" information mainly for operational outputs.
Definition DataClasses.h:67

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

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

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

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

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

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

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

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

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

Constants
Definition Constants.h:43

vanGenuchten.h