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 DumpCanopyHeader(std::ofstream &fout);
                static void DumpCanopyUnits(std::ofstream &fout);
                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, const double rootdepth);
 
                void SoilWaterUptake(const size_t& SoilNode, const double& transpiration,
                        ElementData* EMS, const double wp_fraction,
                        const double rootdepth, const double h_wilt) const;
 
                static double get_f4(const double& tempC);
 
                static double get_f2f4(const size_t& SoilNode, ElementData* EMS,
                          const double wp_fraction,const double rootdepth);
 
                static double get_f3(const double& vpd, const double f3_gd);
 
                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, const double interception_timecoef);
 
                static double CanopyAlbedo(const double& tair, const double& wetfrac, const SnowStation& Xdata);
 
                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, const double& can_diameter);
 
                static double CanopyWetFraction(const double& capacity, const double& storage);
 
                static double CanopyTransmissivity(const double& lai, const double& elev, const double krnt_lai);
 
                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, const double biomass_density, const double biomass_heat_capacity);
 
                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;
 
                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;
};
 
#endif