explain lakedepth corruption safeguards

SamuelTrahanNOAA · Dec 18, 2023 · b1bb75a · b1bb75a
1 parent 929b716
commit b1bb75a
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/physics/clm_lake.f90 b/physics/clm_lake.f90
@@ -239,6 +239,9 @@ subroutine calculate_z_dz_lake(i,input_lakedepth,clm_lakedepth,z_lake,dz_lake)
       real(kind_lake) :: depthratio
 
       if (input_lakedepth(i) == spval .or. input_lakedepth(i) < 0.1) then
+        ! This is a safeguard against:
+        ! 1. missing in the lakedepth database (== spval)
+        ! 2. errors in model cycling or unexpected changes in the orography database (< 0.1)
         clm_lakedepth(i) = zlak(nlevlake) + 0.5_kind_lake*dzlak(nlevlake)
         z_lake(1:nlevlake) = zlak(1:nlevlake)
         dz_lake(1:nlevlake) = dzlak(1:nlevlake)