attempt to address warm high-lat low lev temperatures

Author: rtodling

GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/Process_Library.F90

@@ -405,26 +405,16 @@ function ICE_FRACTION_SC (TEMP,CNV_FRACTION,SRF_TYPE) RESULT(ICEFRCT)
      ! Sigmoidal functions like figure 6b/6c of Hu et al 2010, doi:10.1029/2009JD012384
       if (SRF_TYPE >= 2.0) then
         ! Over snow (SRF_TYPE == 2.0) and ice (SRF_TYPE == 3.0)
-        if (ICE_RADII_PARAM == 1) then
-          ! Jason formula
-          ICEFRCT_M  = 0.00
-          if ( TEMP <= JiT_ICE_ALL ) then
-             ICEFRCT_M = 1.000
-          else if ( (TEMP > JiT_ICE_ALL) .AND. (TEMP <= JiT_ICE_MAX) ) then
-             ICEFRCT_M = 1.00 -  ( TEMP - JiT_ICE_ALL ) / ( JiT_ICE_MAX - JiT_ICE_ALL )
-          end if
-        else
-          ICEFRCT_M  = 0.00
-          if ( TEMP <= iT_ICE_ALL ) then
-             ICEFRCT_M = 1.000
-          else if ( (TEMP > iT_ICE_ALL) .AND. (TEMP <= iT_ICE_MAX) ) then
-             ICEFRCT_M = SIN( 0.5*MAPL_PI*( 1.00 - ( TEMP - iT_ICE_ALL ) / ( iT_ICE_MAX - iT_ICE_ALL ) ) )
-          end if
+        ICEFRCT_M  = 0.00
+        if ( TEMP <= iT_ICE_ALL ) then
+           ICEFRCT_M = 1.000
+        else if ( (TEMP > iT_ICE_ALL) .AND. (TEMP <= iT_ICE_MAX) ) then
+           ICEFRCT_M = SIN( 0.5*MAPL_PI*( 1.00 - ( TEMP - iT_ICE_ALL ) / ( iT_ICE_MAX - iT_ICE_ALL ) ) )
         end if
         ICEFRCT_M = MIN(ICEFRCT_M,1.00)
         ICEFRCT_M = MAX(ICEFRCT_M,0.00)
         ICEFRCT_M = ICEFRCT_M**iICEFRPWR
-      else if (SRF_TYPE > 1.0) then
+      else if (SRF_TYPE == 1.0) then
         ! Over Land
         ICEFRCT_M  = 0.00
         if ( TEMP <= lT_ICE_ALL ) then