fix a bug in snowwater_lake when dealing with snow melt above unfroze…

…n top lake layer
bnuweinan · Jan 21, 2025 · 2e14a75 · 2e14a75
1 parent 673269e
commit 2e14a75
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/main/MOD_Lake.F90 b/main/MOD_Lake.F90
@@ -1799,7 +1799,7 @@ SUBROUTINE snowwater_lake ( USE_Dynamic_Lake, &
          ! all snow melt
          IF (c>=a+b)THEN
             t_lake(1) = (cpliq*(denh2o*dz_lake(1)*t_lake(1) + (sumsnowice+sumsnowliq)*tfrz) - a - b) / &
-                        (cpliq*(denh2o*dz_lake(1) + sumsnowice+ sumsnowice))
+                        (cpliq*(denh2o*dz_lake(1) + sumsnowice + sumsnowliq))
             sm = sm + scv/deltim
             qout_snowb = qout_snowb + scv/deltim
             scv = 0.