From 982fd9152e7cc1e6e17fb624afe770dea8bb06c4 Mon Sep 17 00:00:00 2001 From: Robert Pincus Date: Tue, 20 Jun 2023 16:00:46 -0400 Subject: [PATCH] Recompute Tv in stratosphere --- examples/all-sky/rrtmgp_allsky.F90 | 2 ++ 1 file changed, 2 insertions(+) diff --git a/examples/all-sky/rrtmgp_allsky.F90 b/examples/all-sky/rrtmgp_allsky.F90 index 2558785d1..f10d1b0d4 100644 --- a/examples/all-sky/rrtmgp_allsky.F90 +++ b/examples/all-sky/rrtmgp_allsky.F90 @@ -445,6 +445,7 @@ subroutine compute_profiles(SST, ncol, nlay, p_lay, t_lay, p_lev, t_lev, q_lay, if (z > z_trop) then q = q_t T = SST - gamma*z_trop/(1. + 0.608*q_0) + Tv = (1. + 0.608*q ) * T p = p0 * (Tv/Tv0)**(g/(Rd*gamma)) * exp( -((g*(z-z_trop))/(Rd*Tv)) ) end if p_lay(icol,ilay) = p @@ -469,6 +470,7 @@ subroutine compute_profiles(SST, ncol, nlay, p_lay, t_lay, p_lev, t_lev, q_lay, if (z > z_trop) then q = q_t T = SST - gamma*z_trop/(1. + 0.608*q_0) + Tv = (1. + 0.608*q ) * T p = p0 * (Tv/Tv0)**(g/(Rd*gamma)) * exp( -((g*(z-z_trop))/(Rd*Tv)) ) end if p_lev(icol,ilay) = p