Update loop to apply ice adjusted only at non-land points

(slmskl is not '1'). Fixes #719.
NCAR · Aug 23, 2021 · fb7cfa8 · fb7cfa8
1 parent cabe68f
commit fb7cfa8
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/physics/sfcsub.F b/physics/sfcsub.F
@@ -2019,12 +2019,14 @@ subroutine sfccycle(lugb,len,lsoil,sig1t,deltsfc                    &
 !    *,' tsffcs=',tsffcs(iprnt),' slianl=',slianl(iprnt)
 
       do i=1,len
+       if (nint(slmskl(i)) /= 1) then
         if (sicanl(i) >= min_ice(i)) then
           slianl(i) = 2.0_kind_io8
         else
           slianl(i) = zero
           sicanl(i) = zero
         endif
+       endif
       enddo
 
       if (fh-deltsfc > -0.001 ) then