Skip to content

Commit

Permalink
Output pratemax from UPP (#332)
Browse files Browse the repository at this point in the history
* 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
  • Loading branch information
ericaligo-NOAA authored Jun 14, 2021
1 parent 7df8806 commit ce98991
Show file tree
Hide file tree
Showing 5 changed files with 94 additions and 4 deletions.
2 changes: 1 addition & 1 deletion VERSION
Original file line number Diff line number Diff line change
@@ -1 +1 @@
10.0.7
10.0.8
11 changes: 11 additions & 0 deletions parm/post_avblflds.xml
Original file line number Diff line number Diff line change
Expand Up @@ -4480,6 +4480,17 @@
<scale>-4.0</scale>
</param>

<param>
<post_avblfldidx>508</post_avblfldidx>
<shortname>MAX_PRATE_ON_SURFACE</shortname>
<longname>Maximum Precipitation Rate on surface</longname>
<pdstmpl>tmpl4_8</pdstmpl>
<pname>PRATE</pname>
<stats_proc>MAX</stats_proc>
<fixed_sfc1_type>surface</fixed_sfc1_type>
<scale>6.0</scale>
</param>

<param>
<post_avblfldidx>510</post_avblfldidx>
<shortname>MAX_QMAX_ON_SPEC_HGT_LVL_ABOVE_GRND_2m</shortname>
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-fv3lam.txt
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
2
219
259
260
PRSLEV
32769
ncep_nco
Expand Down Expand Up @@ -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
Expand Down
39 changes: 38 additions & 1 deletion parm/postxconfig-NT-hafs.txt
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
1
88
89
HURPRS
32769
ncep_nco
Expand Down Expand Up @@ -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
?
?
?
7 changes: 6 additions & 1 deletion sorc/ncep_post.fd/INITPOST_NETCDF.f
Original file line number Diff line number Diff line change
Expand Up @@ -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, &
Expand Down Expand Up @@ -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 &
Expand Down

0 comments on commit ce98991

Please sign in to comment.