Skip to content

Home

Jump to bottom
YaliMao-NOAA edited this page Feb 22, 2024 · 11 revisions

Welcome to the EVS WAFS wiki!

The plotting web site is: https://www.emc.ncep.noaa.gov/users/verification/aviation/wafs/prod/

Plot script changes not merged in EVS yet: https://github.com/NOAA-EMC/EVS/compare/develop...YaliMao-NOAA:EVS:develop

Basic scripts:

Stat:

dev/drivers/scripts/stats/wafs/jevs_wafs_atmos_stats.sh
jobs/JEVS_WAFS_ATMOS_STATS
scripts/stats/wafs/exevs_wafs_atmos_stats.sh

Plot:

dev/drivers/scripts/plots/wafs/jevs_wafs_atmos_plots.sh
jobs/JEVS_WAFS_ATMOS_PLOTS
scripts/plots/wafs/exevs_wafs_atmos_plots.sh

Formulas:

ush/wafs/evs_wafs_atmos_plot_util.py

To turn on WDIR plotting:

export WDIR=yes; qsub dev/drivers/scripts/plots/wafs/jevs_wafs_atmos_plots.sh

WDIR control how to assign a value to 'VAR_NAME' in the parm file:

parm/evs_config/wafs/config.evs.stats.wafs.atmos.standalone

        if [[ $WDIR = 'yes' ]] ; then
            export VAR_NAME="TMP, WIND, WIND80, WDIR" # For plotting only                                                                                                                     
        else
            export VAR_NAME="TMP, WIND, WIND80" # For plotting only                                                                                                                           
        fi

Plot over more than 90 days of data

Specify DAYS_LIST and COMIN

export WDIR=yes; export DAYS_LIST=400; export COMIN=/lfs/h2/emc/vpppg/noscrub/yali.mao/stats_from_vsdb; qsub dev/drivers/scripts/plots/wafs/jevs_wafs_atmos_plots.sh

To change the pressure levels of observation data:

Stat (for icing severity):

  • Modify VAR1_LEVELS_OBS:
parm/evs_config/wafs/config.evs.wafs.standalone

 export VAR1_LEVELS_OBS="P812, P696.8, P595.2, P506, P392.7"

Stat (for u/v/t):

  • Modify VAR1_LEVELS :
parm/evs_config/wafs/config.evs.wafs.standalone

 export VAR_LEVELS="P843.1, P696.8, P595.2, P506, P392.7, P300.9, P250, P196.8, P147.5, P100.4"
  • Modify OBS_GRID_STAT_INPUT_TEMPLATE from pgrb2.0p25.anl to wafs.0p25.anl:
parm/metplus_config/stats/wafs/atmos/GridStat_fcstWAFS_obsGFS.conf

 OBS_GRID_STAT_INPUT_TEMPLATE = gfs.{valid?fmt=%Y%m%d}/{valid?fmt=%H}/atmos/gfs.t{valid?fmt=%H}z.wafs.0p25.anl
  • Modify sourcefile from pgrb2.0p25.anl to wafs.0p25.anl:
ush/wafs/evs_wafs_atmos_stats_preparedata.sh

 sourcefile=$sourcedir/gfs.t${cc}z.wafs.0p25.anl

Plot:

ush/wafs/evs_wafs_atmos_plot_settings.py

Under 'var_dict' 'ICESEV', modify obs_var_levels:

'obs_var_levels': ['P812','P696.8','P595.2','P506','P392.7']

Under 'var_dict' 'TMP/WIND/WIND80/WDIR', modify fcst_var_levels and obs_var_levels:

'fcst_var_levels': ['P843.1','P696.8','P595.2','P506','P392.7','P300.9','P250','P196.8','P147.5','P100.4']
'obs_var_levels': ['P843.1','P696.8','P595.2','P506','P392.7','P300.9','P250','P196.8','P147.5','P100.4']
Clone this wiki locally