Now sets rootdepth for crops to 0 outside growing season.

See #1616 (comment)
ESCOMP · Mar 11, 2022 · 3aa7d65 · 3aa7d65
1 parent 8517646
commit 3aa7d65
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/biogeochem/CNRootDynMod.F90 b/src/biogeochem/CNRootDynMod.F90
@@ -135,7 +135,9 @@ subroutine CNRootDyn(bounds, num_soilc, filter_soilc, num_soilp, filter_soilp, &
             c = pcolumn(p)
             if (ivt(p) /= noveg) then
                 if((ivt(p)) >= npcropmin)then !skip generic crop types
-                    if(huigrain(p) > 0._r8)then
+                    if (.not. croplive(p)) then
+                        root_depth(p) = 0._r8
+                    else if(huigrain(p) > 0._r8)then
                         root_depth(p) = max(zi(c,2), min(hui(p)/huigrain(p)* root_dmx(ivt(p)), root_dmx(ivt(p))))
                     end if
                 else