Skip to content

Commit

Permalink
Change GRIB2 parameter names and vertical levels
Browse files Browse the repository at this point in the history
  • Loading branch information
GwenChen-NOAA committed May 20, 2024
1 parent 0cdc279 commit 803a772
Show file tree
Hide file tree
Showing 3 changed files with 17 additions and 14 deletions.
11 changes: 6 additions & 5 deletions ush/nc2grb/nc2grb.ice.sh
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,7 @@ set -x
# input netCDF file has a different resolution, use wgrib2 to regrid
# the template file.
# Author: L. Gwen Chen ([email protected]), 11/02/2023
# updated on 05/20/2024

module load intel/19.1.3.304
module load wgrib2/2.0.8
Expand Down Expand Up @@ -47,23 +48,23 @@ wgrib2 template.global.0p25.gb2 \
# tables in wgrib2 v2.0.8:
#
# -import_netcdf $infile "hs_h" "0:1:0:721:0:1440" \
# -set_var ??? -set center 7 \
# -set_var SNVOLSI -set center 7 \
# -set_date $ddate -set_ftime "$aperiod hour ave fcst" \
# -set_scaling same same -set_grib_type c1 -grib_out $outfile \
# -import_netcdf $infile "frzmlt_h" "0:1:0:721:0:1440" \
# -set_var ??? -set center 7 \
# -set_var FRZMLTPOT -set center 7 \
# -set_date $ddate -set_ftime "$aperiod hour ave fcst" \
# -set_scaling same same -set_grib_type c1 -grib_out $outfile \
# -import_netcdf $infile "albsni_h" "0:1:0:721:0:1440" \
# -set_var ALBICE -set center 7 -rpn "100.0:/" \
# -set_var ALBDOICE -set center 7 -rpn "100.0:/" \
# -set_date $ddate -set_ftime "$aperiod hour ave fcst" \
# -set_scaling same same -set_grib_type c1 -grib_out $outfile \
# -import_netcdf $infile "mlt_onset_h" "0:1:0:721:0:1440" \
# -set_var ??? -set center 7 \
# -set_var MLTDATE -set center 7 \
# -set_date $ddate -set_ftime "$aperiod hour ave fcst" \
# -set_scaling same same -set_grib_type c1 -grib_out $outfile \
# -import_netcdf $infile "frz_onset_h" "0:1:0:721:0:1440" \
# -set_var ??? -set center 7 \
# -set_var FRZDATE -set center 7 \
# -set_date $ddate -set_ftime "$aperiod hour ave fcst" \
# -set_scaling same same -set_grib_type c1 -grib_out $outfile \

Expand Down
7 changes: 4 additions & 3 deletions ush/nc2grb/nc2grb.ocean_2D.sh
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,7 @@ set -x
# input netCDF file has a different resolution, use wgrib2 to regrid
# the template file.
# Author: L. Gwen Chen ([email protected]), 11/02/2023
# updated on 05/20/2024

module load intel/19.1.3.304
module load wgrib2/2.0.8
Expand All @@ -30,7 +31,7 @@ wgrib2 template.global.0p25.gb2 \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out $outfile \
-import_netcdf $infile "SSS" "0:1:0:721:0:1440" \
-set_var SALTY -set center 7 -rpn "1000.0:/" \
-set_var SALIN -set center 7 \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out $outfile \
-import_netcdf $infile "speed" "0:1:0:721:0:1440" \
Expand All @@ -54,11 +55,11 @@ wgrib2 template.global.0p25.gb2 \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out $outfile \
-import_netcdf $infile "SW" "0:1:0:721:0:1440" \
-set_var DSWRF -set center 7 \
-set_var NSWRF -set center 7 \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out $outfile \
-import_netcdf $infile "LW" "0:1:0:721:0:1440" \
-set_var DLWRF -set center 7 \
-set_var NLWRF -set center 7 \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out $outfile \
-import_netcdf $infile "LwLatSens" "0:1:0:721:0:1440" \
Expand Down
13 changes: 7 additions & 6 deletions ush/nc2grb/nc2grb.ocean_3D.sh
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,7 @@ set -x
# input netCDF file has a different resolution, use wgrib2 to regrid
# the template file.
# Author: L. Gwen Chen ([email protected]), 11/02/2023
# updated on 05/20/2024

module load intel/19.1.3.304
module load wgrib2/2.0.8
Expand All @@ -19,30 +20,30 @@ infile=ocean.0p25.nc
outfile=ocean_3D.0p25.gb2
ddate=2013100100
aperiod="0-6"
levlst="5 15 25 35 45 55 65 75 85 95 105 115 125 135 145 155 165 175 185 195 205 215 225.86945 241.06255 266.5239 308.7874 373.9288 467.3998 593.87915 757.1453 959.97325 1204.059 1489.9735 1817.1455 2183.879 2587.3995 3023.9285 3488.7875 3976.524 4481.0625"
levlst="5 15 25 35 45 55 65 75 85 95 105 115 125 135 145 155 165 175 185 195 205 215 226 241 267 309 374 467 594 757 960 1204 1490 1817 2184 2587 3024 3489 3977 4481"

zl=0
for depth in $levlst; do

wgrib2 template.global.0p25.gb2 \
-import_netcdf $infile "temp" "0:1:$zl:1:0:721:0:1440" \
-set_var WTMP -set center 7 -rpn "273.15:+" \
-set_lev "$depth m below water surface" \
-set_lev "$depth m below sea level" \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out tmp.gb2 \
-import_netcdf $infile "so" "0:1:$zl:1:0:721:0:1440" \
-set_var SALTY -set center 7 -rpn "1000.0:/" \
-set_lev "$depth m below water surface" \
-set_var SALIN -set center 7 \
-set_lev "$depth m below sea level" \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out tmp.gb2 \
-import_netcdf $infile "uo" "0:1:$zl:1:0:721:0:1440" \
-set_var UOGRD -set center 7 \
-set_lev "$depth m below water surface" \
-set_lev "$depth m below sea level" \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out tmp.gb2 \
-import_netcdf $infile "vo" "0:1:$zl:1:0:721:0:1440" \
-set_var VOGRD -set center 7 \
-set_lev "$depth m below water surface" \
-set_lev "$depth m below sea level" \
-set_date $ddate -set_ftime "$aperiod hour ave fcst" \
-set_scaling same same -set_grib_type c1 -grib_out tmp.gb2

Expand Down

0 comments on commit 803a772

Please sign in to comment.