From ef31f665d516dad6fe6f92af72182fa8ff21feee Mon Sep 17 00:00:00 2001 From: Manda Chasteen Date: Mon, 23 Jan 2023 22:15:39 -0700 Subject: [PATCH] Initialize tend_u_phys as zero I've added code to initialize and reset tend_u_phys as 0._RKIND during each time step. Otherwise, these tendency values accumulate over successive time steps, which impacts the PV budget calculation in pv_diagnostics.F --- src/core_atmosphere/physics/mpas_atmphys_todynamics.F | 1 + 1 file changed, 1 insertion(+) diff --git a/src/core_atmosphere/physics/mpas_atmphys_todynamics.F b/src/core_atmosphere/physics/mpas_atmphys_todynamics.F index 8bf7279e75..29fb2780b5 100755 --- a/src/core_atmosphere/physics/mpas_atmphys_todynamics.F +++ b/src/core_atmosphere/physics/mpas_atmphys_todynamics.F @@ -235,6 +235,7 @@ subroutine physics_get_tend( block, mesh, state, diag, tend, tend_physics, confi do i = 1, nEdgesScalar+1 do k = 1, nVertLevels tend_ru_physics(k,i) = 0._RKIND + tend_u_phys(k,i) = 0._RKIND enddo enddo !$acc loop gang vector collapse(3)