From c19eb7e8f1b1baf4acb69ebd9cd4315d0ff6a27d Mon Sep 17 00:00:00 2001 From: imreddyTeja Date: Thu, 19 Sep 2024 16:43:02 -0700 Subject: [PATCH] Change runs to use spatialllyvarying rooting depth longruns/land.jl, benchmarks/land.jl and longruns/land_regional.jl are changed to use the root depth map from the clm data --- experiments/benchmarks/land.jl | 8 +++++++- experiments/long_runs/land.jl | 8 +++++++- experiments/long_runs/land_region.jl | 8 +++++++- 3 files changed, 21 insertions(+), 3 deletions(-) diff --git a/experiments/benchmarks/land.jl b/experiments/benchmarks/land.jl index 0f67c8e54f..d86776964b 100644 --- a/experiments/benchmarks/land.jl +++ b/experiments/benchmarks/land.jl @@ -402,7 +402,13 @@ function setup_prob(t0, tf, Δt; nelements = (101, 15)) retention_model = Canopy.PlantHydraulics.LinearRetentionCurve{FT}(a) plant_ν = FT(1.44e-4) plant_S_s = FT(1e-2 * 0.0098) # m3/m3/MPa to m3/m3/m - rooting_depth = FT(0.5) # from Natan + rooting_depth = SpaceVaryingInput( + joinpath(clm_artifact_path, "root_map.nc"), + "rooting_depth", + surface_space; + regridder_type, + regridder_kwargs = (; extrapolation_bc,), + ) n_stem = 0 n_leaf = 1 h_stem = FT(0.0) diff --git a/experiments/long_runs/land.jl b/experiments/long_runs/land.jl index e089c71735..9e29a38bb7 100644 --- a/experiments/long_runs/land.jl +++ b/experiments/long_runs/land.jl @@ -410,7 +410,13 @@ function setup_prob(t0, tf, Δt; outdir = outdir, nelements = (101, 15)) retention_model = Canopy.PlantHydraulics.LinearRetentionCurve{FT}(a) plant_ν = FT(1.44e-4) plant_S_s = FT(1e-2 * 0.0098) # m3/m3/MPa to m3/m3/m - rooting_depth = FT(0.5) # from Natan + rooting_depth = SpaceVaryingInput( + joinpath(clm_artifact_path, "root_map.nc"), + "rooting_depth", + surface_space; + regridder_type, + regridder_kwargs = (; extrapolation_bc,), + ) n_stem = 0 n_leaf = 1 h_stem = FT(0.0) diff --git a/experiments/long_runs/land_region.jl b/experiments/long_runs/land_region.jl index f928adcf0f..1780cb9946 100644 --- a/experiments/long_runs/land_region.jl +++ b/experiments/long_runs/land_region.jl @@ -412,7 +412,13 @@ function setup_prob(t0, tf, Δt; outdir = outdir, nelements = (10, 10, 15)) retention_model = Canopy.PlantHydraulics.LinearRetentionCurve{FT}(a) plant_ν = FT(1.44e-4) plant_S_s = FT(1e-2 * 0.0098) # m3/m3/MPa to m3/m3/m - rooting_depth = FT(0.5) # from Natan + rooting_depth = SpaceVaryingInput( + joinpath(clm_artifact_path, "root_map.nc"), + "rooting_depth", + surface_space; + regridder_type, + regridder_kwargs = (; extrapolation_bc,), + ) n_stem = 0 n_leaf = 1 h_stem = FT(0.0)