diff --git a/examples/cglite_parameter_ensembles.jl b/examples/cglite_parameter_ensembles.jl
index 7573f27680dc846e4bb30547dd15174c375f3a7c..1eec62a5633f21bb7778954e0b0198a2451eaf43 100644
--- 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
index adc4ee933af5b9feab9e269bff0927b258849b63..f6c28d30ff6780a3de36262ef7c6484b58af8579 100644
--- 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)