Merge pull request #1151 from GEOS-ESM/feature/sdrabenh/polynomialChange

Attempt to mitigate low-level warm temperatures in the high-latitudes consistent with FPP

Author: sdrabenh

GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/Process_Library.F90

@@ -413,21 +413,11 @@ function ICE_FRACTION_SC (TEMP,CNV_FRACTION,SRF_TYPE) RESULT(ICEFRCT)
       select case (nint(SRF_TYPE))
       case (SRF_TYPE_SNOW, SRF_TYPE_ICE)
         ! 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)