Fixing U/V tendency calculation.

mwaxmonsky · Apr 12, 2024 · 1f11c14 · 1f11c14
1 parent 2cd5f29
commit 1f11c14
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/tj2016/tj2016_sfc_pbl_hs.F90 b/tj2016/tj2016_sfc_pbl_hs.F90
@@ -467,8 +467,8 @@ subroutine tj2016_sfc_pbl_hs_run(ncol, pver, index_surface, gravit, cappa, rairv
 
     do i = 1, ncol
         do k = 1, pver
-            dudt(i, k)                    = UCopy(i, k) - U(i, k) / dtime
-            dvdt(i, k)                    = VCopy(i, k) - V(i, k) / dtime
+            dudt(i, k)                    = (UCopy(i, k) - U(i, k)) / dtime
+            dvdt(i, k)                    = (VCopy(i, k) - V(i, k)) / dtime
             tendency_of_air_enthalpy(i,k) = (T(i,k) - stateT(i,k)) / dtime * cpairv(i,k)
         end do
     end do