Skip to content

Commit

Permalink
Per #1834, modified formatting, including adding backslashes and remo…
Browse files Browse the repository at this point in the history
…ving extra returns in command line calls, and adding indentation for config file notation, and adding a return after each 'A.'
  • Loading branch information
jprestop committed Jul 7, 2021
1 parent e705475 commit 49cbfda
Showing 1 changed file with 37 additions and 80 deletions.
117 changes: 37 additions & 80 deletions met/docs/Users_Guide/appendixA.rst
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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)
Expand All @@ -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.
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand All @@ -584,32 +562,22 @@ 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:

.. 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
Expand All @@ -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
Expand All @@ -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.
Expand Down Expand Up @@ -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:
Expand Down

0 comments on commit 49cbfda

Please sign in to comment.