Make warning for non-equilibrium.

mvdebolskiy · Jul 23, 2024 · 88adda4 · 88adda4
1 parent 811659f
commit 88adda4
Showing 1 changed file with 11 additions and 0 deletions.
diff --git a/src/biogeophys/WaterStateType.F90 b/src/biogeophys/WaterStateType.F90
@@ -724,6 +724,17 @@ subroutine Restart(this, bounds, ncid, flag, &
     ! Restart excess ice vars
     if (.not. use_excess_ice) then
        ! no need to even define the restart vars
+       call RestartExcessIceIssue( &
+            ncid = ncid, &
+            flag = flag, &
+            excess_ice_on_restart = excess_ice_on_restart)
+       if( excess_ice_on_restart ) then
+          if (masterproc) then
+             write(iulog,*) '--WARNING-- Starting from initial conditions with excess ice present.'
+             write(iulog,*) 'But use_excess_ice=.false.'
+             write(iulog,*) 'This will cause soil moisture and temperature not being in equilibrium'
+          endif
+       endif
        this%excess_ice_col(bounds%begc:bounds%endc,-nlevsno+1:nlevmaxurbgrnd)=0.0_r8
     else
        call RestartExcessIceIssue( &