no autoresp change

Author: AlexisRenchon

Project.toml

@@ -4,6 +4,7 @@ version = "1.7.0"
 authors = ["Clima Land Team"]
 
 [deps]
+CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
 ClimaComms = "3a4d1b5c-c61d-41fd-a00a-5873ba7a1b0d"
 ClimaCore = "d414da3d-4745-48bb-8d80-42e94e092884"
 ClimaDiagnostics = "1ecacbb8-0713-4841-9a07-eb5aa8a2d53f"
@@ -25,7 +26,6 @@ Thermodynamics = "b60c26fb-14c3-4610-9d3e-2d17fe7ff00c"
 
 [weakdeps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
-CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
 ClimaAnalysis = "29b5916a-a76c-4e73-9657-3c8fd22e65e6"
 DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
 DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"

src/standalone/Vegetation/autotrophic_respiration.jl

@@ -139,19 +139,12 @@ function compute_autrophic_respiration(
 )
 
     (; ne, ηsl, σl, μr, μs, Rel) = model.parameters
-    FT = typeof(LAI)
-
-    # Root and stem maintenance respiration: LAI-independent (trees keep roots and stems in winter)
-    Ra_root = Rd * μr * RAI
-    Ra_stem = Rd * μs * ηsl * h * SAI / σl
-
-    # Leaf maintenance and growth respiration: canopy-scaled (→ 0 as LAI → 0)
-    Rpm_leaf = Rd * β
-    Rg = plant_respiration_growth(Rel, An, Rpm_leaf)
-    canopy_scale = (1 - exp(-K * LAI * Ω)) / (K * Ω)
-    Ra_canopy = max(FT(0), Rpm_leaf + Rg) * canopy_scale
-
-    return Ra_canopy + Ra_root + Ra_stem
+    Nl, Nr, Ns =
+        nitrogen_content(ne, Vcmax25, LAI, SAI, RAI, ηsl, h, σl, μr, μs)
+    Rpm = plant_respiration_maintenance(Rd, β, Nl, Nr, Ns)
+    Rg = plant_respiration_growth(Rel, An, Rpm)
+    Ra = Rpm + Rg
+    return Ra * (1 - exp(-K * LAI * Ω)) / (K * Ω) # adjust to canopy level
 end
 
 Base.broadcastable(model::AutotrophicRespirationModel) = tuple(model) # this is so that @. does not broadcast on Ref(canopy.autotrophic_respiration)
@@ -211,7 +204,7 @@ function nitrogen_content(
     μr::FT, # Ratio root nitrogen to top leaf nitrogen (-), typical value 1.0
     μs::FT, # Ratio stem nitrogen to top leaf nitrogen (-), typical value 0.1
 ) where {FT}
-    Sc = ηsl * h * SAI
+    Sc = ηsl * h * LAI * ClimaLand.heaviside(SAI)
     Rc = σl * RAI
     nm = Vcmax25 / ne
     Nl = nm * σl * LAI