""" UserScript: Make ERA OMI plot from calculated MJO indices ========================================================= model_applications/s2s_mjo/UserScript_obsERA_obsOnly_OMI.py """ ############################################################################## # .. contents:: # :depth: 1 # :local: # :backlinks: none ############################################################################## # Scientific Objective # -------------------- # # The Madden-Julian Oscillation (MJO) is the largest element of intraseasonal variability in the # tropics and is characterized by eastward moving regions of enhanced and suppressed rainfall. These # phases are typically grouped into numbers 1 - 8 based on the geographic location of the enhanced and # suppressed rainfall. The MJO affects global weather including summer monsoons, tropical cyclone # development, and sudden stratospheric warming events, and has teleconnections to mid latitude weather # systems. # # This use case uses outgoing longwave radiation (OLR) to compute the OLR based MJO Index (OMI), which is # a convective index of the MJO. OMI is computed for the ERA observations and then displayed on a phase # diagram to evaluate the model reprentation of this important oscillation. The code for computing OMI # came from Maria Gehne at PSL. ############################################################################## # Version Added # ------------- # # METplus version 4.1.0 ############################################################################## # Datasets # -------- # # **Forecast:** None # # **Observation:** ERA Reanlaysis Outgoing Longwave Radiation, 1979 - 2012 # # **Climatology:** None # # **EOFs:** Observed OMI EOF1 and EOF2 patterns from the PSL Website (https://psl.noaa.gov/mjo/mjoindex/) # # **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 # ------------------ # # This use case calls UserScript twice, first to create a list of EOF files and then to run # the OMI calculation. In addition, there is one optional pre-processing steps, using # Regrid-Data-Plane. # # This use case requires METcalcpy, METplotpy, and METdataio to run. The METcalcpy scripts accessed include the following: # # * metcalcpy/contributed/rmm_omi/compute_mjo_indices.py # # The METplotpy scripts accessed include the following: # # * metplotpy/contributed/mjo_rmm_omi/plot_mjo_indioces.py # # The METdataio scripts accessed include the following: # # * METreadnc/util/read_netcdf.py ############################################################################## # METplus Workflow # ---------------- # # **Beginning time calculation (VALID_BEG):** 01-01-1979 # # **End time calculation (VALID_END):** 12-30-2012 # # **Increment between beginning and end times (VALID_INCREMENT):** 1 day # # **Sequence of forecast leads to process (LEAD_SEQ):** 0 hour # # This use case does not loop, but the UserScript to create and EOF filelist is run once and the OMI driver script is # run once. The EOF filelist is done separately since the EOF files are needed for each day of the year while the OMI # calculation is on a separate time frame. The optional pre-processing step loops by valid time. The optional # pre-processing step uses Regrid-Data-Plane to cut the observation grid to only include -20 to 20 latitude. This # omitted step can be turned back on by using the PROCESS_LIST that is commented out in the file: # # PROCESS_LIST = RegridDataPlane(regrid_obs_olr), UserScript(create_eof_filelist), UserScript(script_omi) # # Settings for the optional pre-processing step can be found in the regrid_obs_olr section of the configuration. # Data is not provided in the tarball to run this steps, but the configuration is provided for reference on # how to set up this step. # # The Phase diagram plot is created over a different time frame than the calculation, 10-01-2012 to 03-30-2012. ############################################################################## # 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/s2s_mjo/UserScript_obsERA_obsOnly_OMI.conf. # # .. highlight:: bash # .. literalinclude:: ../../../../parm/use_cases/model_applications/s2s_mjo/UserScript_obsERA_obsOnly_OMI.conf ############################################################################## # MET Configuration # ----------------- # # There are no MET configuration files used in this use case. ############################################################################## # Python Embedding # ---------------- # # This use case does not use python embedding ############################################################################## # Python Scripting # ---------------- # # This use case runs the OMI driver which computes OMI and creates a phase diagram. Inputs to the # OMI driver include netCDF files formatted in MET's netCDF version. In addition, a text file containing # the listing of these input netCDF files is required, as well as text file listings of the EOF1 and # EOF2 files. These text files can be generated using the USER_SCRIPT_INPUT_TEMPLATES in the # [create_eof_filelist] and [script_omi] sections. # # For the OMI calculation, the OLR data are then projected onto Empirical Orthogonal Function (EOF) # data that is computed for each day of the year, latitude, and longitude. The OLR is then filtered # for 20 - 96 days, and regressed onto the daily EOFs. Finally, it's normalized and these normalized # components are plotted on a phase diagram. The OMI driver script orchestrates the calculation of the # MJO indices and the generation of a phase diagram OMI plot. Variables for the OMI calculation are set # in the [user_env_vars] section of the .conf file. # # .. dropdown:: parm/use_cases/model_applications/s2s_mjo/common/OMI_driver.py # # .. highlight:: python # .. literalinclude:: ../../../../parm/use_cases/model_applications/s2s_mjo/common/OMI_driver.py ############################################################################## # 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/s2s_stratosphere/UserScript_obsERA_obsOnly_OMI.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 {OUTPUT_BASE}/model_applications/s2s_mjo/UserScript_fcstGFS_obsERA_OMI/plots. # A Phase diagram plots will be generated:: # # * obs_OMI_comp_phase.png # # If the pre-processing steps are turned on, the output will include the regridded data and daily averaged files. ############################################################################## # Keywords # -------- # # .. note:: # # * S2SAppUseCase # * S2SMJOAppUseCase # * RegridDataPlaneUseCase # * PCPCombineUseCase # * UserScriptUseCase # * METdataioUseCase # * METcalcpyUseCase # * METplotpyUseCase # # Navigate to :ref:`quick-search` to discover other similar use cases. # # sphinx_gallery_thumbnail_path = '_static/s2s_mjo-UserScript_obsERA_obsOnly_OMI.png'