Skip to content

Add ifi products to post #418

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
hu5970 merged 8 commits into from
Sep 26, 2022
Merged

Add ifi products to post #418

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
f54f8ea
add ifi to post files
SamuelTrahanNOAA Aug 19, 2022
ef67c01
correct name of IFI grib file
SamuelTrahanNOAA Aug 19, 2022
c91ec5b
Updating exregional_run_post.sh for the new IFI GRIB2 files.
EricJames-NOAA Aug 24, 2022
a105018
Updating exregional_run_wgrib2 to handle new icing fields. Using nea…
EricJames-NOAA Aug 24, 2022
8846f32
update IFI variable numbers
SamuelTrahanNOAA Aug 25, 2022
5fe6ea2
Merge branch 'ifi' of ssh://github.com/SamuelTrahanNOAA/regional_work…
SamuelTrahanNOAA Aug 25, 2022
eff5143
Updating fix file to be consistent with latest UPP. It has been test…
EricJames-NOAA Sep 14, 2022
0573bef
Merge branch 'feature/RRFS_dev1' into ifi
hu5970 Sep 26, 2022
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
45 changes: 45 additions & 0 deletions fix/upp/fv3lam_rrfs.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3539,4 +3539,49 @@
</param>

</paramset>

<paramset>
<datset>IFIFIP</datset>
<grid_num>4</grid_num>
<sub_center>ncep_emc</sub_center>
<version_no>v2003</version_no>
<local_table_vers_no>local_tab_yes1</local_table_vers_no>
<sigreftime>fcst</sigreftime>
<prod_status>oper</prod_status>
<data_type>fcst</data_type>
<gen_proc_type>fcst</gen_proc_type>
<time_range_unit>hour</time_range_unit>
<orig_center>nws_ncep</orig_center>
<gen_proc>hrrr</gen_proc>
<packing_method>complex_packing_spatial_diff</packing_method>
<order_of_sptdiff>2nd_ord_sptdiff</order_of_sptdiff>
<field_datatype>fltng_pnt</field_datatype>
<comprs_type>lossless</comprs_type>

<!-- New IFI products, only available if libIFI is linked to a post with USE_IFI enabled -->

<param>
<shortname>ICE_PROB_IFI_FLIGHT_LEVEL</shortname>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>

<param>
<shortname>SLD_IFI_FLIGHT_LEVEL</shortname>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>

<param>
<shortname>ICE_SEV_CAT_IFI_FLIGHT_LEVEL</shortname>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>
</paramset>
</postxml>
43 changes: 40 additions & 3 deletions fix/upp/post_avblflds.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6299,16 +6299,17 @@

<param>
<post_avblfldidx>685</post_avblfldidx>
<shortname>DU_CR_AER_SFC_MASS_CON</shortname>
<shortname>DUST10_SFC_MASS_CON</shortname>
<stats_proc>AVE</stats_proc>
<pdstmpl>tmpl4_48</pdstmpl>
<pname>MASSDEN</pname>
<pname>PMTC</pname>
<fixed_sfc1_type>surface</fixed_sfc1_type>
<aerosol_type>dust_dry</aerosol_type>
<typ_intvl_size>smaller_than_first_limit</typ_intvl_size>
<scale_fact_1st_size>6</scale_fact_1st_size>
<scale_val_1st_size>10</scale_val_1st_size>
<fixed_sfc1_type>surface</fixed_sfc1_type>
<scale> 9.0</scale>
<scale>9.0</scale>
</param>

<param>
Expand Down Expand Up @@ -8023,5 +8024,41 @@
<fixed_sfc1_type>surface</fixed_sfc1_type>
<scale>3.0</scale>
</param>

<param>
<post_avblfldidx>1003</post_avblfldidx>
<shortname>ICE_PROB_IFI_FLIGHT_LEVEL</shortname>
<pname>ICPRB</pname>
<table_info>NCEP</table_info>
<fixed_sfc1_type>spec_alt_above_mean_sea_lvl</fixed_sfc1_type>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>

<param>
<post_avblfldidx>1004</post_avblfldidx>
<shortname>SLD_IFI_FLIGHT_LEVEL</shortname>
<pname>SIPD</pname>
<table_info>NCEP</table_info>
<fixed_sfc1_type>spec_alt_above_mean_sea_lvl</fixed_sfc1_type>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>

<param>
<post_avblfldidx>1005</post_avblfldidx>
<shortname>ICE_SEV_CAT_IFI_FLIGHT_LEVEL</shortname>
<pname>ICSEV</pname>
<table_info>NCEP</table_info>
<fixed_sfc1_type>spec_alt_above_mean_sea_lvl</fixed_sfc1_type>
<scale>4.0</scale>
<!-- Level units are in 0.1 feet, so final output will be in units of 0.1 meters -->
<scale_fact_fixed_sfc1>1</scale_fact_fixed_sfc1>
<level>5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.</level>
</param>
</post_avblflds>
</postxml>
130 changes: 129 additions & 1 deletion fix/upp/postxconfig-NT-fv3lam_rrfs.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,5 @@
2
3
3
239
276
PRSLEV
Expand Down Expand Up @@ -19088,3 +19089,130 @@ top_of_atmos
?
?
?
IFIFIP
4
ncep_emc
v2003
local_tab_yes1
fcst
oper
fcst
fcst
hour
nws_ncep
hrrr
complex_packing_spatial_diff
2nd_ord_sptdiff
fltng_pnt
lossless
1003
ICE_PROB_IFI_FLIGHT_LEVEL
?
1
tmpl4_0
ICPRB
NCEP
?
spec_alt_above_mean_sea_lvl
1
1
60
5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
1004
SLD_IFI_FLIGHT_LEVEL
?
1
tmpl4_0
SIPD
NCEP
?
spec_alt_above_mean_sea_lvl
1
1
60
5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
1005
ICE_SEV_CAT_IFI_FLIGHT_LEVEL
?
1
tmpl4_0
ICSEV
NCEP
?
spec_alt_above_mean_sea_lvl
1
1
60
5000. 10000. 15000. 20000. 25000. 30000. 35000. 40000. 45000. 50000. 55000. 60000. 65000. 70000. 75000. 80000. 85000. 90000. 95000. 100000. 105000. 110000. 115000. 120000. 125000. 130000. 135000. 140000. 145000. 150000. 155000. 160000. 165000. 170000. 175000. 180000. 185000. 190000. 195000. 200000. 205000. 210000. 215000. 220000. 225000. 230000. 235000. 240000. 245000. 250000. 255000. 260000. 265000. 270000. 275000. 280000. 285000. 290000. 295000. 300000.
?
0
?
0
?
?
?
0
0.0
0
0.0
?
0
0.0
0
0.0
1
4.0
0
0
0
?
?
?
2 changes: 2 additions & 0 deletions scripts/exregional_run_post.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -342,9 +342,11 @@ fi
bgdawp=${postprd_dir}/${NET}.t${cyc}z.bgdawpf${subh_fhr}.${tmmark}.grib2
bgrd3d=${postprd_dir}/${NET}.t${cyc}z.bgrd3df${subh_fhr}.${tmmark}.grib2
bgsfc=${postprd_dir}/${NET}.t${cyc}z.bgsfcf${subh_fhr}.${tmmark}.grib2
bgifi=${postprd_dir}/${NET}.t${cyc}z.bgifif${subh_fhr}.${tmmark}.grib2

wgrib2 PRSLEV.GrbF${post_fhr} -set center 7 -grib ${bgdawp}
wgrib2 NATLEV.GrbF${post_fhr} -set center 7 -grib ${bgrd3d}
wgrib2 IFIFIP.GrbF${post_fhr} -set center 7 -grib ${bgifi}

#
#-----------------------------------------------------------------------
Expand Down
12 changes: 10 additions & 2 deletions scripts/exregional_run_wgrib2.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -216,6 +216,7 @@ fhr=${subh_fhr}

bgdawp=${postprd_dir}/${NET}.t${cyc}z.bgdawpf${fhr}.${tmmark}.grib2
bgrd3d=${postprd_dir}/${NET}.t${cyc}z.bgrd3df${fhr}.${tmmark}.grib2
bgifi=${postprd_dir}/${NET}.t${cyc}z.bgifif${fhr}.${tmmark}.grib2
bgsfc=${postprd_dir}/${NET}.t${cyc}z.bgsfcf${fhr}.${tmmark}.grib2

# extract the output fields for the testbed
Expand Down Expand Up @@ -245,14 +246,17 @@ fi
basetime=$( date +%y%j%H%M -d "${yyyymmdd} ${hh}" )
cp_vrfy ${bgdawp} ${comout}/${NET}.t${cyc}z.bgdawpf${fhr}.${tmmark}.grib2
cp_vrfy ${bgrd3d} ${comout}/${NET}.t${cyc}z.bgrd3df${fhr}.${tmmark}.grib2
cp_vrfy ${bgifi} ${comout}/${NET}.t${cyc}z.bgifif${fhr}.${tmmark}.grib2
cp_vrfy ${bgsfc} ${comout}/${NET}.t${cyc}z.bgsfcf${fhr}.${tmmark}.grib2
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgdawpf${fhr}.${tmmark}.grib2 ${comout}/BGDAWP_${basetime}${post_fhr}
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgrd3df${fhr}.${tmmark}.grib2 ${comout}/BGRD3D_${basetime}${post_fhr}
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgifif${fhr}.${tmmark}.grib2 ${comout}/BGIFI_${basetime}${post_fhr}
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgsfcf${fhr}.${tmmark}.grib2 ${comout}/BGSFC_${basetime}${post_fhr}

net4=$(echo ${NET:0:4} | tr '[:upper:]' '[:lower:]')
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgdawpf${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.prslev.f${fhr}.conus_3km.grib2
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgrd3df${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.natlev.f${fhr}.conus_3km.grib2
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgifif${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.ififip.f${fhr}.conus_3km.grib2
ln_vrfy -sf --relative ${comout}/${NET}.t${cyc}z.bgsfcf${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.testbed.f${fhr}.conus_3km.grib2
# Remap to additional output grids if requested
if [ ${#ADDNL_OUTPUT_GRIDS[@]} -gt 0 ]; then
Expand All @@ -272,7 +276,7 @@ if [ ${#ADDNL_OUTPUT_GRIDS[@]} -gt 0 ]; then

for grid in ${ADDNL_OUTPUT_GRIDS[@]}
do
for leveltype in dawp rd3d sfc
for leveltype in dawp rd3d ifi sfc
do

eval grid_specs=\$grid_specs_${grid}
Expand All @@ -287,7 +291,7 @@ if [ ${#ADDNL_OUTPUT_GRIDS[@]} -gt 0 ]; then
-new_grid_vectors "UGRD:VGRD:USTM:VSTM:VUCSH:VVCSH" \
-new_grid_interpolation bilinear \
-if ":(WEASD|APCP|NCPCP|ACPCP|SNOD):" -new_grid_interpolation budget -fi \
-if ":(NCONCD|NCCICE|SPNCR|CLWMR|CICE|RWMR|SNMR|GRLE|PMTF|PMTC|REFC|CSNOW|CICEP|CFRZR|CRAIN|LAND|ICEC|TMP:surface|VEG|CCOND|SFEXC|MSLMA|PRES:tropopause|LAI|HPBL|HGT:planetary boundary layer):" -new_grid_interpolation neighbor -fi \
-if ":(NCONCD|NCCICE|SPNCR|CLWMR|CICE|RWMR|SNMR|GRLE|PMTF|PMTC|REFC|CSNOW|CICEP|CFRZR|CRAIN|LAND|ICEC|TMP:surface|VEG|CCOND|SFEXC|MSLMA|PRES:tropopause|LAI|HPBL|HGT:planetary boundary layer):|ICPRB|SIPD|ICSEV" -new_grid_interpolation neighbor -fi \
-new_grid ${grid_specs} ${subdir}/${fhr}/tmp_${grid}.grib2 &
else
wgrib2 ${infile} -set_bitmap 1 -set_grib_type c3 -new_grid_winds grid \
Expand Down Expand Up @@ -316,6 +320,10 @@ if [ ${#ADDNL_OUTPUT_GRIDS[@]} -gt 0 ]; then
ln_vrfy -fs --relative ${comout}/${grid}_grid/${NET}.t${cyc}z.bg${leveltype}f${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.natlev.f${fhr}.conus_3km.grib2
fi

if [ $leveltype = 'ifi' ]; then
ln_vrfy -fs --relative ${comout}/${grid}_grid/${NET}.t${cyc}z.bg${leveltype}f${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.ififip.f${fhr}.conus_3km.grib2
fi

if [ $leveltype = 'sfc' ]; then
ln_vrfy -fs --relative ${comout}/${grid}_grid/${NET}.t${cyc}z.bg${leveltype}f${fhr}.${tmmark}.grib2 ${comout}/${net4}.t${cyc}z.testbed.f${fhr}.conus_3km.grib2
fi
Expand Down