From c9cd92782b01accf26204ac32844b71bdd7b4f64 Mon Sep 17 00:00:00 2001 From: Oksana Guba Date: Thu, 14 Jan 2021 19:10:24 -0600 Subject: [PATCH 1/2] add pnh --- components/eam/src/dynamics/se/stepon.F90 | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/components/eam/src/dynamics/se/stepon.F90 b/components/eam/src/dynamics/se/stepon.F90 index 4c9c49b9366..37dff95a9fa 100644 --- a/components/eam/src/dynamics/se/stepon.F90 +++ b/components/eam/src/dynamics/se/stepon.F90 @@ -28,6 +28,7 @@ module stepon use scamMod, only: use_iop, doiopupdate, single_column, & setiopupdate, readiopdata, iop_mode use element_mod, only: element_t + use element_ops, only: get_field use shr_const_mod, only: SHR_CONST_PI use se_single_column_mod, only: scm_broadcast @@ -149,6 +150,7 @@ subroutine stepon_init(dyn_in, dyn_out ) call addfld('DYN_V' ,(/ 'lev' /), 'A', 'm/s', 'Meridional Velocity', gridname='GLL') call addfld('DYN_OMEGA',(/ 'lev' /), 'A', 'Pa/s', 'Vertical Velocity', gridname='GLL' ) call addfld('DYN_PS' ,horiz_only, 'A', 'Pa', 'Surface pressure', gridname='GLL') + call addfld('DYN_PNH' ,(/ 'lev' /), 'A', 'Pa', 'Nonhydro pressure', gridname='GLL') end subroutine stepon_init @@ -417,6 +419,15 @@ subroutine stepon_run2(phys_state, phys_tend, dyn_in, dyn_out ) call outfld('DIV',tmp_dyn(1,1,1,ie),npsq,ie) enddo endif + + if (hist_fld_active('DYN_PNH')) then + do ie=1,nelemd + ! time level ntQ is not used + call get_field(dyn_in%elem(ie),'pnh',tmp_dyn(:,:,:,ie),hvcoord,tl_f,1) + call outfld('DYN_PNH',tmp_dyn(1,1,1,ie),npsq,ie) + enddo + endif + if (hist_fld_active('FU') .or. hist_fld_active('FV') ) then do ie=1,nelemd do j=1,np From 65fe92aa56cdc5631545486b66799c6becdc5310 Mon Sep 17 00:00:00 2001 From: Oksana Guba Date: Thu, 14 Jan 2021 19:16:38 -0600 Subject: [PATCH 2/2] fix name --- components/eam/src/dynamics/se/stepon.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/components/eam/src/dynamics/se/stepon.F90 b/components/eam/src/dynamics/se/stepon.F90 index 37dff95a9fa..4dbfedce0a5 100644 --- a/components/eam/src/dynamics/se/stepon.F90 +++ b/components/eam/src/dynamics/se/stepon.F90 @@ -150,7 +150,7 @@ subroutine stepon_init(dyn_in, dyn_out ) call addfld('DYN_V' ,(/ 'lev' /), 'A', 'm/s', 'Meridional Velocity', gridname='GLL') call addfld('DYN_OMEGA',(/ 'lev' /), 'A', 'Pa/s', 'Vertical Velocity', gridname='GLL' ) call addfld('DYN_PS' ,horiz_only, 'A', 'Pa', 'Surface pressure', gridname='GLL') - call addfld('DYN_PNH' ,(/ 'lev' /), 'A', 'Pa', 'Nonhydro pressure', gridname='GLL') + call addfld('DYN_PNH' ,(/ 'lev' /), 'A', 'Pa', 'Nonhydrostatic pressure', gridname='GLL') end subroutine stepon_init