From 0dbd755f48e8ce2a377b34b750b231394b07e408 Mon Sep 17 00:00:00 2001 From: "L. Gwen Chen" Date: Mon, 20 May 2024 20:04:13 +0000 Subject: [PATCH] Change GRIB2 parameter names and vertical levels --- parm/post/oceanice_products.yaml | 4 ++-- ush/oceanice_nc2grib2.sh | 22 +++++++++++----------- 2 files changed, 13 insertions(+), 13 deletions(-) diff --git a/parm/post/oceanice_products.yaml b/parm/post/oceanice_products.yaml index 8d8bcd4b51..48e5a5f204 100644 --- a/parm/post/oceanice_products.yaml +++ b/parm/post/oceanice_products.yaml @@ -29,9 +29,9 @@ ocean: cosvar: "cos_rot" angvar: "" {% if model_grid == 'mx025' or model_grid == 'mx050' or model_grid == 'mx100' %} - ocean_levels: [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] + ocean_levels: [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] {% elif model_grid == 'mx500' %} - ocean_levels: [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] + ocean_levels: [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] {% endif %} subset: ['SSH', 'SST', 'SSS', 'speed', 'MLD_003', 'latent', 'sensible', 'SW', 'LW', 'LwLatSens', 'Heat_PmE', 'SSU', 'SSV', 'taux', 'tauy', 'temp', 'so', 'uo', 'vo'] data_in: diff --git a/ush/oceanice_nc2grib2.sh b/ush/oceanice_nc2grib2.sh index 2afd0e07f2..5781e06b36 100755 --- a/ush/oceanice_nc2grib2.sh +++ b/ush/oceanice_nc2grib2.sh @@ -45,23 +45,23 @@ function _ice_nc2grib2 { # tables in wgrib2 v2.0.8: # -import_netcdf "${infile}" "hs_h" "0:1:${latlon_dims}" \ -# -set_var ??? -set center 7 \ +# -set_var SNVOLSI -set center 7 \ # -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ # -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" \ # -import_netcdf "${infile}" "frzmlt_h" "0:1:${latlon_dims}" \ -# -set_var ??? -set center 7 \ +# -set_var FRZMLTPOT -set center 7 \ # -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ # -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" \ # -import_netcdf "${infile}" "albsni_h" "0:1:${latlon_dims}" \ -# -set_var ALBICE -set center 7 -rpn "100.0:/" \ +# -set_var ALBDOICE -set center 7 -rpn "100.0:/" \ # -set_date "${current_cycle}" -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:${latlon_dims}" \ -# -set_var ??? -set center 7 \ +# -set_var MLTDATE -set center 7 \ # -set_date "${current_cycle}" -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:${latlon_dims}" \ -# -set_var ??? -set center 7 \ +# -set_var FRZDATE -set center 7 \ # -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ # -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" @@ -121,11 +121,11 @@ function _ocean2D_nc2grib2 { -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" \ -import_netcdf "${infile}" "SW" "0:1:${latlon_dims}" \ - -set_var DSWRF -set center 7 \ + -set_var NSWRF -set center 7 \ -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" \ -import_netcdf "${infile}" "LW" "0:1:${latlon_dims}" \ - -set_var DLWRF -set center 7 \ + -set_var NLWRF -set center 7 \ -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ -set_scaling same same -set_grib_type c1 -grib_out "${outfile}" \ -import_netcdf "${infile}" "LwLatSens" "0:1:${latlon_dims}" \ @@ -187,22 +187,22 @@ function _ocean3D_nc2grib2 { ${WGRIB2} "${template}" \ -import_netcdf "${infile}" "temp" "0:1:${zl}:1:${latlon_dims}" \ -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 "${current_cycle}" -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:${latlon_dims}" \ -set_var SALIN -set center 7 \ - -set_lev "${depth} m below water surface" \ + -set_lev "${depth} m below sea level" \ -set_date "${current_cycle}" -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:${latlon_dims}" \ -set_var UOGRD -set center 7 \ - -set_lev "${depth} m below water surface" \ + -set_lev "${depth} m below sea level" \ -set_date "${current_cycle}" -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:${latlon_dims}" \ -set_var VOGRD -set center 7 \ - -set_lev "${depth} m below water surface" \ + -set_lev "${depth} m below sea level" \ -set_date "${current_cycle}" -set_ftime "${aperiod} hour ave fcst" \ -set_scaling same same -set_grib_type c1 -grib_out tmp.gb2