Merge pull request #21 from GEOS-ESM/feature/pnorris/RRTMGP-min-temp-…

…selectable pmn: use selectable tmin_OK
GEOS-ESM · Aug 1, 2023 · 2f3cafd · 2f3cafd
2 parents be8a4b4 + 02259e7
commit 2f3cafd
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 b/GEOSirrad_GridComp/GEOS_IrradGridComp.F90
@@ -1466,7 +1466,7 @@ subroutine LW_Driver(IM,JM,LM,LATS,LONS,RC)
    real(wp) :: press_ref_min, ptop
    real(wp) ::  temp_ref_min, tmin
    real(wp), parameter :: ptop_increase_OK_fraction = 0.01_wp
-   real(wp), parameter :: tmin_increase_OK_Kelvin   = 10.0_wp
+   real(wp) :: tmin_increase_OK_Kelvin
 
    ! block size for efficient column processing (set from resource file)
    integer :: rrtmgp_blockSize
@@ -1952,6 +1952,9 @@ subroutine LW_Driver(IM,JM,LM,LATS,LONS,RC)
       tmin = minval(t_lay)
       if (temp_ref_min > tmin) then
         ! allow a small increase of tmin
+        call MAPL_GetResource (MAPL, &
+           tmin_increase_OK_Kelvin, 'RRTMGP_LW_TMIN_INC_OK_K:', &
+           DEFAULT = 10., __RC__)
         if (temp_ref_min - tmin <= tmin_increase_OK_Kelvin) then
           where (t_lay < temp_ref_min) t_lay = temp_ref_min
         else

diff --git a/GEOSsolar_GridComp/GEOS_SolarGridComp.F90 b/GEOSsolar_GridComp/GEOS_SolarGridComp.F90
@@ -2144,7 +2144,7 @@ subroutine SORADCORE(IM,JM,LM,include_aerosols,CURRTIME,LoadBalance,RC)
       real(wp) :: press_ref_min, ptop
       real(wp) ::  temp_ref_min, tmin
       real(wp), parameter :: ptop_increase_OK_fraction = 0.01_wp
-      real(wp), parameter :: tmin_increase_OK_Kelvin   = 10.0_wp
+      real(wp) :: tmin_increase_OK_Kelvin
 
       ! block size for efficient column processing (set from resource file)
       integer :: rrtmgp_blockSize
@@ -3014,6 +3014,9 @@ subroutine SORADCORE(IM,JM,LM,include_aerosols,CURRTIME,LoadBalance,RC)
       tmin = minval(t_lay)
       if (temp_ref_min > tmin) then
         ! allow a small increase of tmin
+        call MAPL_GetResource (MAPL, &
+           tmin_increase_OK_Kelvin, 'RRTMGP_SW_TMIN_INC_OK_K:', &
+           DEFAULT = 10., __RC__)
         if (temp_ref_min - tmin <= tmin_increase_OK_Kelvin) then
           where (t_lay < temp_ref_min) t_lay = temp_ref_min
         else