""" Grid-Stat: Verification of TC forecasts against merged TDR data =============================================================== model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF.conf """ ############################################################################## # .. contents:: # :depth: 1 # :local: # :backlinks: none ############################################################################## # Scientific Objective # -------------------- # To provide useful statistical information on the relationship between merged # Tail Doppler Radar (TDR) data in NetCDF format to a gridded forecast. These # values can be used to assess the skill of the prediction. The TDR data is # available every 0.5 km AGL. So, the TC forecasts need to be in height coordinates # to compare with the TDR data. ############################################################################## # Version Added # ------------- # # METplus version 4.0 ############################################################################## # Datasets # -------- # # **Forecast:** NOAA Hurricane Analysis and Forecast System (HAFS) zonal wind # # **Observation:** NOAA Hurricane Research Division (HRD) Tail Doppler Radar (TDR) merged_zonal_wind, # a subset of the Merged Analysis (v2d_combined_xy_rel_merged_ships.nc). # # **Climatology:** None # # **Location:** All of the input data required for this use case can be # found in a sample data tarball. Each use case category will have # one or more sample data tarballs. It is only necessary to download # the tarball with the use case’s dataset and not the entire collection # of sample data. Click here to access the METplus releases page and download sample data # for the appropriate release: https://github.com/dtcenter/METplus/releases # This tarball should be unpacked into the directory that you will # set the value of INPUT_BASE. See :ref:`running-metplus` section for more information. ############################################################################## # METplus Components # ------------------ # # The observations in the use case contains data mapped into Cartesian Grids with a horizontal grid spacing of 2 km and vertical grid spacing of 0.5 km. Hence the model output needs to be in height (km) (vertical coordinates) instead of pressure levels. Both observation and model output are available with the release. The instructions below tells how the input to the use case was prepared. # The Hurricane Analysis and Forecast System (HAFS) (pressure levels in GRIB2 format) outputs are converted to height level (in NetCDF4 format) using METcalcpy vertical interpolation routine. Under METcalcpy/examples directory user can modify the vertical_interp_hwrf.sh or create a similar file for their own output. The $DATA_DIR is the top level output directory where the pressure level data resides. The --input and --output should point to the input and output file names resp. The --config points to a yaml file. Users should edit the yaml file, if needed. For this use case only zonal wind (u) at 4 (200m, 2000m, 4000m and 6000m) vertical levels are provided. The use case will compare the HAFS 2 km zonal wind (u) data against TDR's merged_zonal_wind at 2km. The user need to run the shell script to get the height level output in NetCDF4 format. # This use case utilizes the METplus python embedding to read the TDR data and compare them to gridded forecast data using GridStat. ############################################################################## # METplus Workflow # ---------------- # # **Beginning time (VALID_BEG):** 2019082912 # # **End time (VALID_END):** 2019082912 # # **Increment between beginning and end times (VALID_INCREMENT):** 21600 # # **Sequence of forecast leads to process (LEAD_SEQ):** 0,6,12,18 # # The use case runs the python embedding scripts (GridStat_fcstHAFS_obsTDR_NetCDF/read_tdr.py: to read the TDR data) and run Grid-Stat (compute statistics against HAFS model output, in height coordinates), called in this example. # # It processes the following run times: Valid at 2019-08-29 12Z # # Forecast lead times: 0,6,12 and 18 UTC # # The mission number (e.g CUSTOM_LOOP_LIST = 190829H1) # # Height level (for TDR: OBS_VERT_LEVEL_KM = 2, HAFS: FCST_VAR1_LEVELS = "(0,1,*,*)") ############################################################################## # METplus Configuration # --------------------- # # METplus first loads all of the configuration files found in parm/metplus_config, # then it loads any configuration files passed to METplus via the command line, # i.e. parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF.conf # # .. highlight:: bash # .. literalinclude:: ../../../../parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF.conf # ############################################################################## # MET Configuration # ----------------- # # METplus sets environment variables based on user settings in the METplus # configuration file. See :ref:`How METplus controls MET config file settings` for more details. # # **YOU SHOULD NOT SET ANY OF THESE ENVIRONMENT VARIABLES YOURSELF! THEY WILL BE OVERWRITTEN BY METPLUS WHEN IT CALLS THE MET TOOLS!** # # If there is a setting in the MET configuration file that is currently # not supported by METplus you’d like to control, please refer to: # :ref:`Overriding Unsupported MET config file settings` # # .. dropdown:: GridStatConfig_wrapped # # .. literalinclude:: ../../../../parm/met_config/GridStatConfig_wrapped ############################################################################## # Python Embedding # ---------------- # # This use case uses a Python embedding script to read input data. # # .. dropdown:: parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF/read_tdr.py # # .. highlight:: python # .. literalinclude:: ../../../../parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF/read_tdr.py # # The read_tdr.py script above imports another script called tdr_utils.py in the same directory: # # .. dropdown:: parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF/tdr_utils.py # # .. highlight:: python # .. literalinclude:: ../../../../parm/use_cases/model_applications/tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF/tdr_utils.py # # For more information on the basic requirements to utilize Python Embedding in METplus, # please refer to the MET User’s Guide section on # `Python embedding `_. ############################################################################## # Python Embedding # ---------------- # # This use case does not use Python embedding. ############################################################################## # User Scripting # -------------- # # User Scripting is not used in this use case. ############################################################################## # Running METplus # --------------- # # Pass the use case configuration file to the run_metplus.py script along # with any user-specific system configuration files if desired:: # # run_metplus.py /path/to/METplus/parm/use_cases/model_applications//tc_and_extra_tc/GridStat_fcstHAFS_obsTDR_NetCDF.conf /path/to/user_system.conf # # See :ref:`running-metplus` for more information. ############################################################################## # Expected Output # --------------- # # A successful run will output the following both to the screen and to the logfile:: # # INFO: METplus has successfully finished running. # # Refer to the value set for **OUTPUT_BASE** to find where the output data was generated. # Output for this use case will be found in nam (relative to **OUTPUT_BASE**) # and will contain the following files: # # * grid_stat_HAFS_vs_TDR_000000L_20190829_120000V_fho.txt # * grid_stat_HAFS_vs_TDR_000000L_20190829_120000V_pairs.nc # * grid_stat_HAFS_vs_TDR_000000L_20190829_120000V.stat # * The use case is run for 4 lead times valid at 2019081912, so four directories will be generated which contains similar files as above. ############################################################################## # Keywords # -------- # # .. note:: # # * TCandExtraTCAppUseCase # * GridStatToolUseCase # * TropicalCycloneUseCase # # Navigate to the :ref:`quick-search` page to discover other similar use cases. # # # # sphinx_gallery_thumbnail_path = '_static/tc_and_extra_tc-GridStat_fcstHAFS_obsTDR_NetCDF.png' #