Recompute Tv in stratosphere

earth-system-radiation · Jun 20, 2023 · 982fd91 · 982fd91
1 parent bfffbd3
commit 982fd91
Showing 1 changed file with 2 additions and 0 deletions.
diff --git 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