Refactors the configuration of freeze curves in heat conduction and soil types

Since the removal of the old SFCC type in FreezeCurves which included a nonlinear solver with the freeze curve, the solver is specified in the Soil type. This doesn't really make sense, however, since the solver is not directly related to the soil itself and is rather a numerical detail. Furthermore, it is conceivable that other types of materials could potentially also induce implicit constitutive relations that would require such a solver for the inverse enthalpy formulation. The solver is therefore now specified in the operator type when necessary. No solver need be specified for heat conduction formulations that do not require it, e.g. the apparent heat capacity formulation (now referred to as ForwardTemperature in the code).

examples/heat_freeW_samoylov_bucketW.jl

@@ -14,15 +14,14 @@ grid = CryoGrid.Presets.DefaultGrid_5cm
 _, tempprofile = CryoGrid.Presets.SamoylovDefault
 initT = initializer(:T, tempprofile)
 initsat = initializer(:sat, (l,p,state) -> state.sat .= l.para.sat)
-heatop = Heat.Diffusion(:H)
 # @Stratigraphy macro lets us list multiple subsurface layers
 strat = @Stratigraphy(
     -2.0u"m" => Top(TemperatureGradient(tair), Rainfall(pr)),
-    0.0u"m" => :topsoil1 => Soil(HomogeneousMixture(por=0.80,sat=0.7,org=0.75), heat=HeatBalance(heatop), water=WaterBalance(BucketScheme())),
-    0.1u"m" => :topsoil2 => Soil(HomogeneousMixture(por=0.80,sat=0.8,org=0.25), heat=HeatBalance(heatop), water=WaterBalance(BucketScheme())),
-    0.4u"m" => :sediment1 => Soil(HomogeneousMixture(por=0.55,sat=0.9,org=0.25), heat=HeatBalance(heatop), water=WaterBalance(BucketScheme())),
-    3.0u"m" => :sediment2 => Soil(HomogeneousMixture(por=0.50,sat=1.0,org=0.0), heat=HeatBalance(heatop), water=WaterBalance(BucketScheme())),
-    10.0u"m" => :sediment3 => Soil(HomogeneousMixture(por=0.30,sat=1.0,org=0.0), heat=HeatBalance(heatop), water=WaterBalance(BucketScheme())),
+    0.0u"m" => :topsoil1 => Soil(HomogeneousMixture(por=0.80,sat=0.7,org=0.75), heat=HeatBalance(), water=WaterBalance(BucketScheme())),
+    0.1u"m" => :topsoil2 => Soil(HomogeneousMixture(por=0.80,sat=0.8,org=0.25), heat=HeatBalance(), water=WaterBalance(BucketScheme())),
+    0.4u"m" => :sediment1 => Soil(HomogeneousMixture(por=0.55,sat=0.9,org=0.25), heat=HeatBalance(), water=WaterBalance(BucketScheme())),
+    3.0u"m" => :sediment2 => Soil(HomogeneousMixture(por=0.50,sat=1.0,org=0.0), heat=HeatBalance(), water=WaterBalance(BucketScheme())),
+    10.0u"m" => :sediment3 => Soil(HomogeneousMixture(por=0.30,sat=1.0,org=0.0), heat=HeatBalance(), water=WaterBalance(BucketScheme())),
     1000.0u"m" => Bottom(GeothermalHeatFlux(0.053u"W/m^2"))
 );
 tile = Tile(strat, grid, initT, initsat);