From 218242f736ff15e5d5722f3767b5713c2b8e3284 Mon Sep 17 00:00:00 2001 From: jatkinson1000 Date: Mon, 30 Sep 2024 10:05:18 -0600 Subject: [PATCH] Add code to produce output required by the testcase. --- src/physics/cam/gw_drag.F90 | 13 +++++++++++++ 1 file changed, 13 insertions(+) diff --git a/src/physics/cam/gw_drag.F90 b/src/physics/cam/gw_drag.F90 index ae7bbfca6f..cc39d3bad6 100644 --- a/src/physics/cam/gw_drag.F90 +++ b/src/physics/cam/gw_drag.F90 @@ -1063,6 +1063,15 @@ subroutine gw_init() call add_default('EKGW', 1, ' ') end if + call addfld ('RHOI', (/ 'ilev' /), 'A', 'kg/m3', & + 'density at interfaces') + + call addfld ('DSE', (/ 'lev' /), 'I', 'J/kg', & + 'dry static energy') + + call addfld ('NMBV', (/ 'lev' /), 'I', 'J/kg', & + 'Brunt Vaisala Frequency') + call addfld ('UTGW_TOTAL', (/ 'lev' /), 'A','m/s2', & 'Total U tendency due to gravity wave drag') call addfld ('VTGW_TOTAL', (/ 'lev' /), 'A','m/s2', & @@ -2096,6 +2105,10 @@ subroutine gw_tend(state, pbuf, dt, ptend, cam_in, flx_heat) call outfld('CLDLIQTGW', ptend%q(:,:,ixcldliq), pcols, lchnk) call outfld('CLDICETGW', ptend%q(:,:,ixcldice), pcols, lchnk) + call outfld('RHOI', rhoi, pcols, lchnk) + call outfld('DSE', dse, ncol, lchnk) + call outfld('NMBV', nm, ncol, lchnk) + ! Destroy objects. call p%finalize()