diff --git a/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 b/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 index 6695e10..67e6329 100644 --- 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 index d9817d2..cd41e82 100644 --- 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