From 158ca814e72519539058b4118d233842a1658966 Mon Sep 17 00:00:00 2001 From: Manda Chasteen Date: Tue, 24 Jan 2023 12:06:30 -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 | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/core_atmosphere/physics/mpas_atmphys_todynamics.F b/src/core_atmosphere/physics/mpas_atmphys_todynamics.F index c7baeb213a..d9778fbfb5 100644 --- a/src/core_atmosphere/physics/mpas_atmphys_todynamics.F +++ b/src/core_atmosphere/physics/mpas_atmphys_todynamics.F @@ -181,6 +181,8 @@ subroutine physics_get_tend( block, mesh, state, diag, tend, tend_physics, confi tend_rtheta_physics(:,:) = 0._RKIND tend_rho_physics(:,:) = 0._RKIND ! NB: rho tendency is not currently supplied by physics, but this ! field may be later filled with IAU or other tendencies + + tend_u_phys(:,:) = 0._RKIND ! ! In case some variables are not allocated due to their associated packages,