Fix error in parameter ensembles example

CryoGrid · Oct 20, 2024 · 89b9fd9 · 89b9fd9
1 parent d2744b2
commit 89b9fd9
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/examples/cglite_parameter_ensembles.jl b/examples/cglite_parameter_ensembles.jl
@@ -35,8 +35,7 @@ upperbc = TemperatureBC(
 )
 ssinit = ThermalSteadyStateInit(T0=-15.0u"°C")
 heatop = Heat.EnthalpyImplicit()
-freezecurve = FreeWater()
-heat = HeatBalance(heatop; freezecurve)
+heat = HeatBalance(heatop)
 soil_layers = map(para -> Ground(para; heat), soilprofile)
 strat = Stratigraphy(
     z_top => Top(upperbc),

diff --git a/src/Physics/Snow/snow_lite.jl b/src/Physics/Snow/snow_lite.jl
@@ -138,6 +138,7 @@ function CryoGrid.initialcondition!(snowpack::LiteSnowpack, ::SnowMassBalance, s
     @. state.ρsn = snowpack.para.ρsn_0
     @. state.swe = 0.0
     @. state.dsn = 0.0
+    @. state.ubc_idx = length(state.T)
 end
 
 function Hydrology.watercontent!(snow::LiteSnowpack, ::WaterBalance, state)