Fix issues with example 6

CryoGrid · Aug 29, 2023 · 09ecd90 · 09ecd90
1 parent 15845a1
commit 09ecd90
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/examples/06_heat_freeW_seb_snow_bucketW_samoylov.jl b/examples/06_heat_freeW_seb_snow_bucketW_samoylov.jl
@@ -19,7 +19,7 @@ soilprofile = SoilProfile(
 ## mid-winter temperature profile
 tempprofile = CryoGrid.Presets.SamoylovDefault.tempprofile
 forcings = loadforcings(CryoGrid.Presets.Forcings.Samoylov_ERA_obs_fitted_1979_2014_spinup_extended_2044);
-soilprofile, tempprofile = CryoGrid.Presets.SamoylovDefault
+tempprofile = CryoGrid.Presets.SamoylovDefault.tempprofile
 initT = initializer(:T, tempprofile)
 z = 2.0u"m"; # height [m] for which the forcing variables (Temp, humidity, wind, pressure) are provided
 seb = SurfaceEnergyBalance(forcings.Tair, forcings.pressure, forcings.q, forcings.wind, forcings.Lin, forcings.Sin, z)

diff --git a/src/Physics/Hydrology/water_ET.jl b/src/Physics/Hydrology/water_ET.jl
@@ -40,7 +40,7 @@ function evapotranspiration!(
         state.w_tr[i] = Δz[i]*exp(-z[i] / et.d_tr)
         let θwi = state.θwi[i],
             θfc = minwater(sub, water, state, i);
-            state.αᶿ[i] = ifelse(θwi < θfc, 0.25(1-cos(π*θwi/θfc))^2, one(θwi))
+            state.αᶿ[i] = ifelse(θwi < θfc, 0.25*(1-cos(π*θwi/θfc))^2, one(θwi))
         end
     end
     @inbounds @. state.f_et = et.f_tr*(state.αᶿ*state.w_tr) + (1-et.f_tr)*(state.αᶿ*state.w_ev)

diff --git a/src/Physics/Soils/soil_water.jl b/src/Physics/Soils/soil_water.jl
@@ -1,7 +1,7 @@
 """
     RREqForm
 
-Base type for different formulations of Richard's equation.
+Base type for different formulations of Richardson-Richard's equation.
 """
 abstract type RREqForm end
 """