Skip to content

Enabling GSL p-type, winter QPF and buckets, and lightning threat indices for RRFS #695

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
WenMeng-NOAA merged 6 commits into from
Apr 25, 2023
Merged

Enabling GSL p-type, winter QPF and buckets, and lightning threat indices for RRFS #695

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
d723ed1
Enabling GSL precip type
EricJames-NOAA Apr 20, 2023
e74e60d
Merge remote-tracking branch 'origin/develop' into winter_qpf
EricJames-NOAA Apr 20, 2023
98e3c68
Trying to output lightning threats.
EricJames-NOAA Apr 21, 2023
51dae2c
Getting lightning threat indices output.
EricJames-NOAA Apr 21, 2023
5207ffb
Updating control files
EricJames-NOAA Apr 24, 2023
1fffd20
Adding "tmpl4_8" back into lightning variables.
EricJames-NOAA Apr 24, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 17 additions & 13 deletions parm/fv3lam_rrfs.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -726,12 +726,6 @@
<scale>-4.0</scale>
</param>

<param>
<shortname>ACM_WEASD_ON_SURFACE</shortname>
<pname>WEASD</pname>
<scale>4.0</scale>
</param>

<param>
<shortname>ACM_SNOM_ON_SURFACE</shortname>
<pname>SNOM</pname>
Expand Down Expand Up @@ -759,28 +753,28 @@
</param>

<param>
<shortname>INST_CRAIN_ON_SURFACE</shortname>
<shortname>GSD_INST_CRAIN_ON_SURFACE</shortname>
<pname>CRAIN</pname>
<table_info>NCEP</table_info>
<scale>1.0</scale>
</param>

<param>
<shortname>CSNOW_ON_SURFACE</shortname>
<shortname>GSD_CSNOW_ON_SURFACE</shortname>
<pname>CSNOW</pname>
<table_info>NCEP</table_info>
<scale>1.0</scale>
</param>

<param>
<shortname>CICEP_ON_SURFACE</shortname>
<shortname>GSD_CICEP_ON_SURFACE</shortname>
<pname>CICEP</pname>
<table_info>NCEP</table_info>
<scale>1.0</scale>
</param>

<param>
<shortname>CFRZR_ON_SURFACE</shortname>
<shortname>GSD_CFRZR_ON_SURFACE</shortname>
<pname>CFRZR</pname>
<table_info>NCEP</table_info>
<scale>1.0</scale>
Expand Down Expand Up @@ -2396,8 +2390,8 @@
</param>

<param>
<shortname>ACM_WEASD_ON_SURFACE</shortname>
<pname>WEASD</pname>
<shortname>ACM_SNOWFALL_ON_SURFACE</shortname>
<pname>TSNOWP</pname>
<scale>4.0</scale>
</param>

Expand Down Expand Up @@ -3550,10 +3544,20 @@
<scale>4.0</scale>
</param>

<param>
<shortname>GSD_MAX_LTG_THREAT1_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>GSD_MAX_LTG_THREAT2_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
<scale>4.0</scale>
</param>

<param>
<shortname>GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS</shortname>
<pname>LTNG</pname>
<scale>3.0</scale>
<scale>4.0</scale>
</param>

<param>
Expand Down
16 changes: 8 additions & 8 deletions parm/post_avblflds.xml
View file
Open in desktop
Comment thread
WenMeng-NOAA marked this conversation as resolved.
Original file line number Diff line number Diff line change
Expand Up @@ -6561,23 +6561,23 @@
<!-- no grib2 definition for 702-710 -->
<param>
<post_avblfldidx>702</post_avblfldidx>
<shortname>GSD_MAX_LTG_THREAT1_ON_ENTIRE_ATMOS</shortname>
<shortname>GSD_MAX_LTG_THREAT1_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
<pdstmpl>tmpl4_8</pdstmpl>
Comment thread
WenMeng-NOAA marked this conversation as resolved.
<pname>LTGTHREAT1</pname>
<table_info>NCEP</table_info>
<pname>LTNGSD</pname>
<stats_proc>MAX</stats_proc>
<fixed_sfc1_type>entire_atmos</fixed_sfc1_type>
<fixed_sfc1_type>spec_hgt_lvl_above_grnd</fixed_sfc1_type>
<level>1.</level>
<scale>3.0</scale>
</param>

<param>
<post_avblfldidx>703</post_avblfldidx>
<shortname>GSD_MAX_LTG_THREAT2_ON_ENTIRE_ATMOS</shortname>
<shortname>GSD_MAX_LTG_THREAT2_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
<pdstmpl>tmpl4_8</pdstmpl>
Comment thread
WenMeng-NOAA marked this conversation as resolved.
<pname>LTGTHREAT2</pname>
<table_info>NCEP</table_info>
<pname>LTNGSD</pname>
<stats_proc>MAX</stats_proc>
<fixed_sfc1_type>entire_atmos</fixed_sfc1_type>
<fixed_sfc1_type>spec_hgt_lvl_above_grnd</fixed_sfc1_type>
<level>2.</level>
<scale>3.0</scale>
</param>

Expand Down
147 changes: 92 additions & 55 deletions parm/postxconfig-NT-fv3lam_rrfs.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
3
4
235
288
237
287
PRSLEV
32769
ncep_nco
Expand Down Expand Up @@ -3718,43 +3718,6 @@ surface
?
?
?
35
ACM_WEASD_ON_SURFACE
?
1
tmpl4_8
WEASD
?
ACM
surface
0
?
0
?
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
121
ACM_SNOM_ON_SURFACE
?
Expand Down Expand Up @@ -3903,9 +3866,9 @@ surface
?
?
?
160
INST_CRAIN_ON_SURFACE
?
407
GSD_INST_CRAIN_ON_SURFACE
GSD_instant precipitation type on surface
1
tmpl4_0
CRAIN
Expand Down Expand Up @@ -3940,9 +3903,9 @@ surface
?
?
?
551
CSNOW_ON_SURFACE
Categorical snow on surface
559
GSD_CSNOW_ON_SURFACE
GSD_Categorical snow on surface
1
tmpl4_0
CSNOW
Expand Down Expand Up @@ -3977,9 +3940,9 @@ surface
?
?
?
552
CICEP_ON_SURFACE
Categorical ice pellets on surface
560
GSD_CICEP_ON_SURFACE
GSD_Categorical ice pellets on surface
1
tmpl4_0
CICEP
Expand Down Expand Up @@ -4014,9 +3977,9 @@ surface
?
?
?
553
CFRZR_ON_SURFACE
Categorical freezing rain on surface
561
GSD_CFRZR_ON_SURFACE
GSD_Categorical freezing rain on surface
1
tmpl4_0
CFRZR
Expand Down Expand Up @@ -13021,12 +12984,12 @@ surface
?
?
?
35
ACM_WEASD_ON_SURFACE
1004
ACM_SNOWFALL_ON_SURFACE
?
1
tmpl4_8
WEASD
TSNOWP
?
ACM
surface
Expand Down Expand Up @@ -19311,6 +19274,80 @@ top_of_atmos
?
?
?
702
GSD_MAX_LTG_THREAT1_ON_SPEC_HGT_LVL_ABOVE_GRND
?
1
tmpl4_8
LTNGSD
?
MAX
spec_hgt_lvl_above_grnd
0
?
1
1.
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
703
GSD_MAX_LTG_THREAT2_ON_SPEC_HGT_LVL_ABOVE_GRND
?
1
tmpl4_8
LTNGSD
?
MAX
spec_hgt_lvl_above_grnd
0
?
1
2.
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
704
GSD_MAX_LTG_THREAT3_ON_ENTIRE_ATMOS
?
Expand Down Expand Up @@ -19341,7 +19378,7 @@ entire_atmos
0
0.0
1
3.0
4.0
0
0
0
Expand Down
37 changes: 20 additions & 17 deletions sorc/ncep_post.fd/INITPOST_NETCDF.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@
!> 2023-04-04 |Li(Kate Zhang) |Add namelist optoin for CCPP-Chem(UFS-Chem)
! and 2D diag. output (d2d_chem) for GEFS-Aerosols and CCPP-Chem model.
!> 2023-04-17 | Eric James | Read in unified ext550 extinction (and remove aodtot) for RRFS
!> 2023-04-21 | Eric James | Read in / calculate some fields needed for GSL p-type diagnosis for RRFS
!>
!> @author Hui-Ya Chuang @date 2016-03-04
SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
Expand Down Expand Up @@ -64,7 +65,7 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
uz0, vz0, ptop, htop, pbot, hbot, ptopl, pbotl, ttopl, ptopm, pbotm, ttopm, &
ptoph, pboth, pblcfr, ttoph, runoff, tecan, tetran, tedir, twa, maxtshltr, &
mintshltr, maxrhshltr, fdnsst, acgraup, graup_bucket, acfrain, frzrn_bucket, &
snow_acm, snow_bkt, &
snow_acm, snow_bkt, snownc, graupelnc, qrmax, &
minrhshltr, dzice, smcwlt, suntime, fieldcapa, htopd, hbotd, htops, hbots, &
cuppt, dusmass, ducmass, dusmass25, ducmass25, aswintoa,rel_vort_maxhy1, &
maxqshltr, minqshltr, acond, sr, u10h, v10h,refd_max, w_up_max, w_dn_max, &
Expand Down Expand Up @@ -877,10 +878,18 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
spval,recname(18),ext550(ista_2l,jsta_2l,1),lm)
endif

! Compute max QRAIN in the column to be used later in precip type computation
do j = jsta_2l, jend_2u
do i = 1, im
qrmax(i,j)=0.
end do
end do

! calculate CWM from FV3 output
do l=1,lm
do j=jsta,jend
do i=ista,iend
qrmax(i,j)=max(qrmax(i,j),qqr(i,j,l))
cwm(i,j,l)=qqg(i,j,l)+qqs(i,j,l)+qqr(i,j,l)+qqi(i,j,l)+qqw(i,j,l)
enddo
enddo
Expand Down Expand Up @@ -1709,6 +1718,16 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,frzrn_bucket)

! time step snow (in m)
VarName='snow'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,snownc)

! time step graupel (in m)
VarName='graupel'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,graupelnc)

Comment thread
WenMeng-NOAA marked this conversation as resolved.
! aerodynamic conductance
VarName='acond'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
Expand Down Expand Up @@ -2872,22 +2891,6 @@ SUBROUTINE INITPOST_NETCDF(ncid2d,ncid3d)
smstav(i,j) = buf(i,j)
enddo
enddo
VarName='accswe_land'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,buf)
VarName='accswe_ice'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,buf2)
!$omp parallel do private(i,j)
do j = jsta_2l, jend_2u
do i=ista,iend
if(buf(i,j)<spval .and. buf2(i,j)<spval) then
acsnow(i,j) = buf(i,j) + buf2(i,j)
else
acsnow(i,j) = spval
endif
enddo
enddo
VarName='snacc_land'
call read_netcdf_2d_para(ncid2d,ista,ista_2l,iend,iend_2u,jsta,jsta_2l,jend,jend_2u, &
spval,VarName,buf)
Expand Down
Loading