global run of soil/canopy model

.buildkite/pipeline.yml

@@ -109,6 +109,12 @@ steps:
         agents:
           slurm_mem: 16G
 
+      - label: "Global Run CPU"
+        command: "julia --color=yes --project=.buildkite experiments/integrated/global/global_soil_canopy.jl"
+        artifact_paths: "experiments/integrated/global/plots/*png"
+        agents:
+          slurm_mem: 16G
+
   - group: "Experiments on GPU"
     steps:
       - label: "Richards Runoff GPU"

.buildkite/target/pipeline.yml

@@ -45,3 +45,13 @@ steps:
         agents:
           slurm_mem: 8GB
           slurm_gpus: 1
+
+      - label: "Soil/Canopy"
+        command: "julia --color=yes --project=.buildkite experiments/benchmarks/land.jl"
+        artifact_paths:
+          - "land_benchmark_gpu/*html"
+        env:
+          CLIMACOMMS_DEVICE: CUDA
+        agents:
+          slurm_mem: 8GB
+          slurm_gpus: 1

docs/Manifest.toml

@@ -2338,9 +2338,9 @@ version = "0.6.12"
 
 [[deps.RecursiveArrayTools]]
 deps = ["Adapt", "ArrayInterface", "DocStringExtensions", "GPUArraysCore", "IteratorInterfaceExtensions", "LinearAlgebra", "RecipesBase", "SparseArrays", "StaticArraysCore", "Statistics", "SymbolicIndexingInterface", "Tables"]
-git-tree-sha1 = "3400ce27995422fb88ffcd3af9945565aad947f0"
+git-tree-sha1 = "980aabbeac7aee70d0e452a72b0c68b5b266cc7b"
 uuid = "731186ca-8d62-57ce-b412-fbd966d074cd"
-version = "3.23.1"
+version = "3.24.0"
 
     [deps.RecursiveArrayTools.extensions]
     RecursiveArrayToolsFastBroadcastExt = "FastBroadcast"
@@ -2745,9 +2745,9 @@ weakdeps = ["ClimaParams"]
 
 [[deps.SymbolicIndexingInterface]]
 deps = ["Accessors", "ArrayInterface", "RuntimeGeneratedFunctions", "StaticArraysCore"]
-git-tree-sha1 = "a5f6f138b740c9d93d76f0feddd3092e6ef002b7"
+git-tree-sha1 = "2d1b8941f58d0c878b3c71d9b17188a2eb649a8b"
 uuid = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
-version = "0.3.22"
+version = "0.3.24"
 
 [[deps.TOML]]
 deps = ["Dates"]

docs/tutorials/integrated/soil_canopy_tutorial.jl

@@ -339,14 +339,14 @@ Y.soil.ϑ_l = FT(0.4)
 Y.soil.θ_i = FT(0.0)
 T_0 = FT(288.7)
 ρc_s =
-    volumetric_heat_capacity.(Y.soil.ϑ_l, Y.soil.θ_i, Ref(land.soil.parameters))
-Y.soil.ρe_int =
-    volumetric_internal_energy.(
+    volumetric_heat_capacity.(
+        Y.soil.ϑ_l,
         Y.soil.θ_i,
-        ρc_s,
-        T_0,
-        Ref(land.soil.parameters),
+        land.soil.parameters.ρc_ds,
+        earth_param_set,
     )
+Y.soil.ρe_int =
+    volumetric_internal_energy.(Y.soil.θ_i, ρc_s, T_0, earth_param_set)
 
 Y.soilco2.C .= FT(0.000412) # set to atmospheric co2, mol co2 per mol air
 

docs/tutorials/standalone/Soil/evaporation.jl

@@ -158,9 +158,14 @@ function init_soil!(Y, z, params)
     Y.soil.ϑ_l .= hydrostatic_equilibrium.(z, FT(-0.001), params)
     Y.soil.θ_i .= 0
     T = FT(296.15)
-    ρc_s = @. Soil.volumetric_heat_capacity(Y.soil.ϑ_l, FT(0), params)
+    ρc_s = @. Soil.volumetric_heat_capacity(
+        Y.soil.ϑ_l,
+        FT(0),
+        params.ρc_ds,
+        params.earth_param_set,
+    )
     Y.soil.ρe_int =
-        Soil.volumetric_internal_energy.(FT(0), ρc_s, T, Ref(params))
+        Soil.volumetric_internal_energy.(FT(0), ρc_s, T, params.earth_param_set)
 end
 init_soil!(Y, z, soil.parameters);
 

docs/tutorials/standalone/Soil/evaporation_gilat_loess.jl

@@ -146,9 +146,19 @@ function init_soil!(Y, z, params)
     Y.soil.ϑ_l .= estimated_ic.(z)
     Y.soil.θ_i .= 0
     T = FT(294.15)
-    ρc_s = @. Soil.volumetric_heat_capacity(Y.soil.ϑ_l, Y.soil.θ_i, params)
+    ρc_s = @. Soil.volumetric_heat_capacity(
+        Y.soil.ϑ_l,
+        Y.soil.θ_i,
+        params.ρc_ds,
+        params.earth_param_set,
+    )
     Y.soil.ρe_int =
-        Soil.volumetric_internal_energy.(Y.soil.θ_i, ρc_s, T, Ref(params))
+        Soil.volumetric_internal_energy.(
+            Y.soil.θ_i,
+            ρc_s,
+            T,
+            params.earth_param_set,
+        )
 end
 
 init_soil!(Y, z, soil.parameters)

docs/tutorials/standalone/Soil/freezing_front.jl

@@ -183,9 +183,19 @@ function init_soil!(Ysoil, z, params)
     Ysoil.soil.ϑ_l .= FT(0.33)
     Ysoil.soil.θ_i .= FT(0.0)
     T = FT(279.85)
-    ρc_s = Soil.volumetric_heat_capacity(FT(0.33), FT(0.0), params)
+    ρc_s = Soil.volumetric_heat_capacity(
+        FT(0.33),
+        FT(0.0),
+        params.ρc_ds,
+        params.earth_param_set,
+    )
     Ysoil.soil.ρe_int .=
-        Soil.volumetric_internal_energy.(FT(0.0), ρc_s, T, Ref(params))
+        Soil.volumetric_internal_energy.(
+            FT(0.0),
+            ρc_s,
+            T,
+            params.earth_param_set,
+        )
 end
 
 init_soil!(Y, coords.subsurface.z, soil.parameters);

docs/tutorials/standalone/Soil/soil_energy_hydrology.jl

@@ -221,9 +221,20 @@ function init_soil!(Y, z, params)
     T = @.(T_min + (T_max - T_min) * exp(-(z - zmax) / (zmin - zmax) * c))
 
     θ_l = Soil.volumetric_liquid_fraction.(Y.soil.ϑ_l, ν, θ_r)
-    ρc_s = Soil.volumetric_heat_capacity.(θ_l, Y.soil.θ_i, Ref(params))
+    ρc_s =
+        Soil.volumetric_heat_capacity.(
+            θ_l,
+            Y.soil.θ_i,
+            params.ρc_ds,
+            params.earth_param_set,
+        )
     Y.soil.ρe_int .=
-        Soil.volumetric_internal_energy.(Y.soil.θ_i, ρc_s, T, Ref(params))
+        Soil.volumetric_internal_energy.(
+            Y.soil.θ_i,
+            ρc_s,
+            T,
+            params.earth_param_set,
+        )
 end
 
 init_soil!(Y, coords.subsurface.z, soil.parameters);

docs/tutorials/standalone/Soil/sublimation.jl

@@ -144,9 +144,14 @@ function init_soil!(Y, z, params)
     Y.soil.ϑ_l .= hydrostatic_equilibrium.(z, FT(-0.1), params)
     Y.soil.θ_i .= 0
     T = FT(275.0)
-    ρc_s = @. Soil.volumetric_heat_capacity(Y.soil.ϑ_l, FT(0), params)
+    ρc_s = @. Soil.volumetric_heat_capacity(
+        Y.soil.ϑ_l,
+        FT(0),
+        params.ρc_ds,
+        params.earth_param_set,
+    )
     Y.soil.ρe_int =
-        Soil.volumetric_internal_energy.(FT(0), ρc_s, T, Ref(params))
+        Soil.volumetric_internal_energy.(FT(0), ρc_s, T, params.earth_param_set)
 end
 init_soil!(Y, z, soil.parameters);