From ce989911a7a09a2e2a0e61b3acc87588b5b9fc26 Mon Sep 17 00:00:00 2001 From: ericaligo-NOAA <48365233+ericaligo-NOAA@users.noreply.github.com> Date: Mon, 14 Jun 2021 11:06:45 -0400 Subject: [PATCH] Output pratemax from UPP (#332) * add pratemax to xml and create new flat file for fv3lam * pratemax in mm/hr added to initpost_netcdf.f * switched pratemax position in file and created new flat file * include hafs flat file for pratemax change * change version from 10.0.7 to 10.0.8 --- VERSION | 2 +- parm/post_avblflds.xml | 11 ++++++++ parm/postxconfig-NT-fv3lam.txt | 39 ++++++++++++++++++++++++++++- parm/postxconfig-NT-hafs.txt | 39 ++++++++++++++++++++++++++++- sorc/ncep_post.fd/INITPOST_NETCDF.f | 7 +++++- 5 files changed, 94 insertions(+), 4 deletions(-) diff --git a/VERSION b/VERSION index 9380cfccb8..5219a0df74 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -10.0.7 +10.0.8 diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml index 2ca5b9087a..2bdb0f090f 100755 --- a/parm/post_avblflds.xml +++ b/parm/post_avblflds.xml @@ -4480,6 +4480,17 @@ -4.0 + + 508 + MAX_PRATE_ON_SURFACE + Maximum Precipitation Rate on surface + tmpl4_8 + PRATE + MAX + surface + 6.0 + + 510 MAX_QMAX_ON_SPEC_HGT_LVL_ABOVE_GRND_2m diff --git a/parm/postxconfig-NT-fv3lam.txt b/parm/postxconfig-NT-fv3lam.txt index bfed37d887..a78cfffe27 100644 --- a/parm/postxconfig-NT-fv3lam.txt +++ b/parm/postxconfig-NT-fv3lam.txt @@ -1,6 +1,6 @@ 2 219 -259 +260 PRSLEV 32769 ncep_nco @@ -9156,6 +9156,43 @@ isobaric_sfc ? ? ? +508 +MAX_PRATE_ON_SURFACE +Maximum Precipitation Rate on surface +1 +tmpl4_8 +PRATE +? +MAX +surface +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? 444 LTNG_ON_SURFACE lightning diff --git a/parm/postxconfig-NT-hafs.txt b/parm/postxconfig-NT-hafs.txt index c3f237ce32..d60fabd68e 100644 --- a/parm/postxconfig-NT-hafs.txt +++ b/parm/postxconfig-NT-hafs.txt @@ -1,5 +1,5 @@ 1 -88 +89 HURPRS 32769 ncep_nco @@ -3272,3 +3272,40 @@ entire_atmos_single_lyr ? ? ? +508 +MAX_PRATE_ON_SURFACE +Maximum Precipitation Rate on surface +1 +tmpl4_8 +PRATE +? +MAX +surface +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? diff --git a/sorc/ncep_post.fd/INITPOST_NETCDF.f b/sorc/ncep_post.fd/INITPOST_NETCDF.f index 70622e21da..891893c7b5 100644 --- a/sorc/ncep_post.fd/INITPOST_NETCDF.f +++ b/sorc/ncep_post.fd/INITPOST_NETCDF.f @@ -64,7 +64,7 @@ SUBROUTINE INITPOST_NETCDF(ncid3d) avgedir,avgecan,avgetrans,avgesnow,avgprec_cont,avgcprate_cont,rel_vort_max, & avisbeamswin,avisdiffswin,airbeamswin,airdiffswin,refdm10c_max,wspd10max, & alwoutc,alwtoac,aswoutc,aswtoac,alwinc,aswinc,avgpotevp,snoavg, & - ti,aod550,du_aod550,ss_aod550,su_aod550,oc_aod550,bc_aod550 + ti,aod550,du_aod550,ss_aod550,su_aod550,oc_aod550,bc_aod550,prate_max use soil, only: sldpth, sh2o, smc, stc use masks, only: lmv, lmh, htm, vtm, gdlat, gdlon, dx, dy, hbm2, sm, sice use physcons_post, only: grav => con_g, fv => con_fvirt, rgas => con_rd, & @@ -1308,6 +1308,11 @@ SUBROUTINE INITPOST_NETCDF(ncid3d) ! GFS does not have accumulated total, gridscale, and convective precip, will use inst precip to derive in SURFCE.f +! max hourly surface precipitation rate + VarName='pratemax' + call read_netcdf_2d_scatter(me,ncid2d,1,im,jm,jsta,jsta_2l & + ,jend_2u,MPI_COMM_COMP,icnt,idsp,spval,VarName,prate_max) + if(debugprint)print*,'sample ',VarName,' = ',prate_max(isa,jsa) ! max hourly 1-km agl reflectivity VarName='refdmax' call read_netcdf_2d_scatter(me,ncid2d,1,im,jm,jsta,jsta_2l &