Canopy.h

Go to the documentation of this file.

/*
 *  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 CANOPY_H
#define CANOPY_H

#include <snowpack/DataClasses.h>

#include <string>
#include <fstream>

class Canopy {

        public:
                Canopy(const SnowpackConfig& i_cfg);

                static void DumpCanopyData(std::ofstream &fout, const CanopyData *Cdata, const SurfaceFluxes *Sdata, const double cos_sl);
                bool runCanopyModel(CurrentMeteo &Mdata, SnowStation &Xdata, const double& roughness_length,
                                    const double& height_of_wind_val, const bool& adjust_VW_height=true);
                static void writeTimeSeriesAdd2LCanopy(std::ofstream &fout, const CanopyData *Cdata);
                static const double can_alb_dry, can_alb_wet, can_alb_snow, krnt_lai; //public constants

        private:
                static double get_f1(const double& ris);
                static double RootFraction(const double& zupper, const double& zlower);
                void SoilWaterUptake(const size_t& SoilNode, const double& transpiration, ElementData* EMS) const;
                static double get_f4(const double& tempC);
                static double get_f2f4(const size_t& SoilNode, ElementData* EMS);
                static double get_f3(const double& vpd);
                double IntCapacity(const CurrentMeteo& Mdata, const SnowStation& Xdata, const bool& force_rain=false) const;
                static double IntUnload(const double& capacity, const double& storage);
                static double IntRate(const double& capacity, const double& storage, const double& prec,
                                  const double& direct);

                static double CanopyAlbedo(const double& tair, const double& wetfrac);
                static double TotalAlbedo(double CanAlb, double sigf, double SurfAlb, double DirectThroughfall,
                                      double CanopyClosureDirect, double RadFracDirect, double sigfdirect);

                static double CanopyShadeSoilCover(const double& HEIGHT, const double& COVER, const double& ELEV);
                static double CanopyWetFraction(const double& capacity, const double& storage);
                static double CanopyTransmissivity(const double& lai, const double& elev);

                void LineariseNetRadiation(const CurrentMeteo& Mdata,const CanopyData& Cdata, const SnowStation& Xdata,
                                              double& iswrac, double& rsnet, double& ilwrac, double& r0,double& r1,
                                              const double& canopyalb, double& CanopyClosureDirect, double& RadFracDirect,
                                              const double& sigfdirect, double& r1p) const;
                void LineariseNetRadiation2L(const CurrentMeteo& Mdata, const CanopyData& Cdata, const SnowStation& Xdata,
                                      double& iswrac, double& rsnet, double& ilwrac, double& r0,double& r1, double& r2,
                                      double& rt0, double& rt1, double& rt2, const double& canopyalb, double& CanopyClosureDirect, double& RadFracDirect,
                                      const double& sigfdirect, const double& sigftrunkdirect, double& r1p, double& r2p) const;
                static void LineariseSensibleHeatFlux(const double& ch_canopy, const double& tair, double& h0, double& h1, double scalingfactor);

                static double DSaturationPressureDT(const double& L, const double& T);
                static void LineariseLatentHeatFlux(const double& ce_canopy, const double& tc_old, const double& vpair,
                                                double& le0, double& le1, double scalingfactor);
                static void CalculateHeatMass(const double& height, const double& BasalArea, double& lai ,double& HMLeaves,  double& HMTrunks);

                void LineariseConductiveHeatFlux(const double& tc_old, const double& HM, double& HM0, double& HM1,  const double& DT, const double& scalingfactor) const;

                static void CanopyEnergyBalance(const double& h0, const double& h1, const double& le0,
                                                         const double& le1, const double& HM0,  const double& HM1,
                                                         const double& ce_canopy,
                                                         const double& ce_condensation,
                                                         double& r0, double& r1, double& TCANOPY, double& RNCANOPY,
                                                         double& HCANOPY, double& LECANOPY);

                static void CanopyEnergyBalance2L(double& h0, double& h1, double& le0,
                                                         double& le1, double& HM0, double& HM1, double& TT0, double& TT1,
                                                         const double& ce_canopy,
                                                         const double& ce_condensation,
                                                         double& r0, double& r1, double& r2, double& TCANOPY, double& Ttrunk, double& RNCANOPY,
                                                         double& HCANOPY, double& LECANOPY);

                static void CanopyEvaporationComponents(const double& ce_canopy,
                                      const double& ce_transpiration, double& LECANOPY,
                                      const double& ta, const double& I, const double DT,
                                      double& CanopyEvaporation,
                                      double& INTEVAP, double& TRANSPIRATION,
                                      double& RNCANOPY, double& HCANOPY,double& TCANOPY,
                                      const double& r0, const double& r1, const double& h0, const double& h1,
                                      double& LECANOPYCORR,
                                      const double& wetfraction, const double& HM0, const double& HM1);

                static void CanopyEvaporationComponents2L(const double& ce_canopy,
                                      const double& ce_transpiration, double& LECANOPY,
                                      const double& ta, const double& I, const double DT,
                                      double& CanopyEvaporation,
                                      double& INTEVAP, double& TRANSPIRATION,
                                      double& RNCANOPY, double& HCANOPY,double& TCANOPY, double& Ttrunk,
                                      const double& TT0, const double& TT1,
                                      const double& r0, const double& r1, const double& r2, const double& h0, const double& h1,
                                      double& LECANOPYCORR,
                                      const double& wetfraction,
                                      const double& HM0, const double& HM1);
                static double get_psim(const double& xi);
                static double get_psih(const double& xi);
                static double RichardsonToAeta(double za, double TempAir, double DiffTemp, double Windspeed, double zom, double zoh, int maxitt);

                void CanopyTurbulentExchange(const CurrentMeteo& Mdata, const double& refheight, const double& zomg,
                                                                  const double& wetfraction, SnowStation& Xdata, double& ch_canopy,
                                                                  double& ce_canopy, double& ce_transpiration,
                                                                  double& ce_interception, double& ce_condensation) const;

                void CanopyRadiationOutput(SnowStation& Xdata, const CurrentMeteo& Mdata, double ac,
                                                                double &iswrac, double &rswrac,
                                                                double &iswrbc, double &rswrbc, double &ilwrac,
                                                                double &rlwrac, double &ilwrbc, double &rlwrbc,
                                                                double CanopyClosureDirect, double RadFracDirect, double sigfdirect, double sigftrunkdirect) const;

                static const double int_cap_snow, int_cap_rain, interception_timecoef;
                static const bool canopy_stabilitycorrection;
                static const double can_diameter, roughmom_to_canopyheight_ratio, displ_to_canopyheight_ratio, raincrease_snow;
                static const double canopytemp_maxchange_perhour, roughheat_to_roughmom_ratio, can_ch0, can_rs_mult, rsmin;
                static const double f3_gd, rootdepth, wp_fraction;

                std::string hn_density, hn_density_parameterization, variant, watertransportmodel_soil;
                double hn_density_fixedValue, calculation_step_length;
                bool useSoilLayers;
                // variables for canopy heat mass and 2-layer canopy
                bool CanopyHeatMass;
                bool Twolayercanopy, Twolayercanopy_user;
                bool canopytransmission;
                bool forestfloor_alb;
                static const double biomass_heat_capacity, biomass_density, lai_frac_top_default, trunk_frac_height, trunkalb, et;
};

#endif