diff --git a/met/docs/Users_Guide/appendixA.rst b/met/docs/Users_Guide/appendixA.rst index e84b703211..f5c379ef3d 100644 --- a/met/docs/Users_Guide/appendixA.rst +++ b/met/docs/Users_Guide/appendixA.rst @@ -95,13 +95,10 @@ Re-run but pointing only to the admin_China_data: ${MET_BUILD_BASE}/bin/plot_data_plane china_tmp_2m.nc china_tmp_2m_admin.ps \ - 'name="TMP_Z2"; level="(*,*)"; - map_data = { source = [ - { file_name = - "${MET_BUILD_BASE}/data/map/admin_by_country/admin_China_data"; - } - ]; - }' + 'name="TMP_Z2"; level="(*,*)"; \ + map_data = { source = [ { file_name = \ + "${MET_BUILD_BASE}/data/map/admin_by_country/admin_China_data"; } \ + ]; }' An arbitrary number of map_data "file_name" entries can be listed. However, using "country_data" doesn't @@ -255,9 +252,9 @@ Gen_Vx_Mask **Q. Gen_Vx_Mask - How do I Mask Region Intersection between Stations and Polyline?** -I have a list of stations to use for verification. I also have a poly +**I have a list of stations to use for verification. I also have a poly region defined. If I specify both of these should the result -be a union of them? +be a union of them?** A. These settings are defined in the "mask" section of the Point-Stat @@ -449,7 +446,6 @@ be used to aggregate them together into a single FSS value, like this: .. code-block:: ini stat_analysis -job aggregate -line_type NBRCNT \ - -lookin out/grid_stat Be sure to pick thresholds (e.g. for the thunderstorms and monsoons) @@ -463,9 +459,11 @@ generated by a MET tool: .. code-block:: ini - fcst = { field = [ - - { name = "HGT_P500"; level = [ "(*,*)" ]; } ]; } + fcst = { field = + [ + { name = "HGT_P500"; level = [ "(*,*)" ]; } + ]; + } Do not use numbers, such as "(181,360)", please use "(\*,\*)" instead. NetCDF variables can have an arbitrary number of dimensions. @@ -493,46 +491,27 @@ The config file should look something like this... .. code-block:: ini fcst = { - - wind_thresh = [ NA ]; - - field = [ - - { - - name = "LCDC"; - - level = [ "L0" ]; - - prob = TRUE; - - cat_thresh = [ >=0.0, >=0.1, >=0.2, >=0.3, >=0.4, >=0.5, >=0.6, >=0.7, >=0.8, >=0.9]; - - } - - ]; - - }; + wind_thresh = [ NA ]; + field = [ + { + name = "LCDC"; + level = [ "L0" ]; + prob = TRUE; + cat_thresh = [ >=0.0, >=0.1, >=0.2, >=0.3, >=0.4, >=0.5, >=0.6, >=0.7, >=0.8, >=0.9]; + } + ]; + }; obs = { - - wind_thresh = [ NA ]; - - field = [ - - { - - name = "WIND"; - - level = [ "Z2" ]; - - cat_thresh = [ >=34 ]; - - } - - ]; - - }; + wind_thresh = [ NA ]; + field = [ + { + name = "WIND"; + level = [ "Z2" ]; + cat_thresh = [ >=34 ]; + } + ]; + }; Without seeing how it's encoded in the GRIB file, it is unclear how to handle “name” in the forecast section. The PROB flag is set to TRUE @@ -569,7 +548,6 @@ Note in there the following... .. code-block:: ini "prob = TRUE;" # tells MET to interpret this data a probability field. - "cat_thresh = [ >=0.0, >=0.1, >=0.2, >=0.3, >=0.4, >=0.5, >=0.6, >=0.7, >=0.8, >=0.9]; " Here the thresholds are used to fully partition the probability space @@ -584,15 +562,10 @@ Run Grid-Stat using the following commands and the attached config file .. code-block:: ini mkdir out - ${MET_BUILD_BASE}/bin/grid_stat \ - gfs_4_20160220_0000_012.grb2 \ - ST4.2016022012.06h \ - - GridStatConfig \ - + GridStatConfig \ -outdir out Note the following two sections of the Grid-Stat config file: @@ -600,16 +573,11 @@ Note the following two sections of the Grid-Stat config file: .. code-block:: ini regrid = { - - to_grid = OBS; - - vld_thresh = 0.5; - - method = BUDGET; - - width = 2; - - } + to_grid = OBS; + vld_thresh = 0.5; + method = BUDGET; + width = 2; + } This tells Grid-Stat to do verification on the "observation" grid. Grid-Stat reads the GFS and Stage4 data and then automatically regrids @@ -621,8 +589,7 @@ regrid both to NCEP Grid 212 before comparing them. .. code-block:: ini - mask = { grid = [ "FULL" ]; - + mask = { grid = [ "FULL" ]; poly = [ "MET_BASE/poly/CONUS.poly" ]; } This will compute statistics over the FULL model domain as well @@ -634,11 +601,8 @@ commands to plot its NetCDF matched pairs output file: .. code-block:: ini ${MET_BUILD_BASE}/bin/plot_data_plane \ - out/grid_stat_120000L_20160220_120000V_pairs.nc \ - out/DIFF_APCP_06_A06_APCP_06_A06_CONUS.ps \ - 'name="DIFF_APCP_06_A06_APCP_06_A06_CONUS"; level="(*,*)";' Examine the resulting plot of that difference field. @@ -671,19 +635,12 @@ point where the temperature is less than 290 K to a value of 0: .. code-block:: ini {MET_BUILD_BASE}/bin/gen_vx_mask \ - data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \ - data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212 \ - APCP_03_where_2m_TMPge290.nc \ - -type data \ - -input_field 'name="APCP"; level="A3";' \ - -mask_field 'name="TMP"; level="Z2";' \ - -thresh 'lt290&&ne-9999' -v 4 -value 0 So this is a bit confusing. Here's what is happening: