fix "sending pbl_utils the wrong definition of 'theta' in clubb_intr.…

…F90 ESCOMP#1222"
PeterHjortLauritzen · Jan 15, 2025 · e229ed5 · e229ed5
1 parent e5ceda5
commit e229ed5
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/physics/cam/clubb_intr.F90 b/src/physics/cam/clubb_intr.F90
@@ -4701,7 +4701,9 @@ subroutine clubb_tend_cam( state,   ptend_all,   pbuf,     hdtime, &
     do i=1,ncol
       do k=1,pver
          !use local exner since state%exner is not a proper exner
-         th(i,k) = state1%t(i,k)*inv_exner_clubb(i,k)
+!+++arh
+         !th(i,k) = state1%t(i,k)*inv_exner_clubb(i,k)
+         th(i,k) = state1%t(i,k)*state1%exner(i,k)
          !thv should have condensate loading to be consistent with earlier def's in this module
          thv(i,k) = th(i,k)*(1.0_r8+zvir*state1%q(i,k,ixq) - state1%q(i,k,ixcldliq))
       enddo