page header
 
Case description

In this exercise you will be running verification scripts on pre-generated WRF-WRFDA cycling forecasts. This is a case from late December, 2013, when a major extratropical cyclone impacted northwestern Europe. It is the same domain as was used for the radiance practice. There are 3 series of 6-hour forecasts you will use in this practice session:

  • /home/ec2-user/test_data/verification/no_da contains a control case: 6-hour WRF forecasts initialized every 6 hours from 2013-12-23 00z to 2013-12-25 12z
  • /home/ec2-user/test_data/verification/conv_only contains a cycling experiment: 6-hour forecasts every 6 hours, with the output of the previous forecast used to initialize the next forecast after assimilating conventional obs.
  • /home/ec2-user/test_data/verification/conv_and_rad contains another cycling experiment: 6-hour forecasts every 6 hours, with the output of the previous forecast used to initialize the next forecast after assimilating conventional and AMSU-A radiance obs.
WRFDA Tools: Observation-based verification

This exercise will teach you how to run WRFDA verification tools for observation-based verification.

Reference: Download the tutorial presentation

Source code

Get the pre-compiled code, if you have not done so.

WRFDA/var/build/da_verif_obs.exe is the executable that will be used in this session. It will be called directly by scripts, so there is no need to have it in your working directory.

Choose your working directory

For this exercise you should create /home/ec2-user/workdir/verif_obs and use this as your working directory.

mkdir /home/ec2-user/workdir/verif_obs
cd /home/ec2-user/workdir/verif_obs

The observation files which will be used for this observation verification exercise have been pre-generated (using the "QC-OBS" procedure mentioned in the talk); they can be found in /home/ec2-user/test_data/arctic60/filt_ob

The easiest way to control these scripts is by creating a "wrapper" script that sets the appropriate environment variables, then calls the scripts in order. Open a new file "verif_obs_wrapper.sh" in vi (or your preferred text editor), and type the following:

> vi verif_obs_wrapper.sh
#!/bin/bash -aeu

# Wrapper script for running WRFDA obs verification package

# Settings for ./da_run_suite_verif_obs.sh

export CLEAN=false
export INITIAL_DATE=2013122318
export FINAL_DATE=2013122512
export WRFVAR_DIR=/home/ec2-user/compiled_code/WRFDA-v4.1.2
export SCRIPTS_DIR=/home/ec2-user/WRFDA_TOOLS/scripts
export OB_DIR=/home/ec2-user/test_data/arctic60/ob/
export FILTERED_OBS_DIR=/home/ec2-user/test_data/arctic60/filt_ob/
export BE_DIR=/home/ec2-user/test_data/arctic60/be/
export WINDOW_START=-3
export WINDOW_END=3
export CYCLE_PERIOD=6
export NUM_PROCS=8
export VERIFY_HOUR=6
export RUN_CMD="mpirun -np $NUM_PROCS"
export VERIFICATION_FILE_STRING=wrfout
# Here is where you set the appropriate namelist variables that the script will use to run WRFDA
export NL_ANALYSIS_TYPE=verify
export NL_E_WE=181
export NL_E_SN=181
export NL_E_VERT=41
export NL_DX=60000
export NL_DY=60000
export NL_SF_SURFACE_PHYSICS=2
export NL_NUM_LAND_CAT=24
export NL_FORCE_USE_OLD_DATA=T

export EXP=no_da
export EXP_DIR=`pwd`/${EXP}
export FC_DIR=/home/ec2-user/test_data/verification/${EXP}/fc

${SCRIPTS_DIR}/da_run_suite_verif_obs.sh

export EXP=conv_only
export EXP_DIR=`pwd`/${EXP}
export FC_DIR=/home/ec2-user/test_data/verification/${EXP}/fc

${SCRIPTS_DIR}/da_run_suite_verif_obs.sh

# Settings for da_verif_obs_plot.sh

export START_DATE=2013122318
export END_DATE=2013122512
export WDIR=/home/ec2-user/workdir/verif_obs
export NUM_EXPT=2
export EXP_NAMES='no_da conv_only'
export EXP_LEGENDS='(/"no_da", "conv_only"/)'
export EXP_DIRS="${WDIR}/no_da ${WDIR}/conv_only"
export GRAPHICS_DIR=/home/ec2-user/WRFDA_TOOLS/graphics/ncl
export WRF_FILE="/home/ec2-user/test_data/verification/conv_only/fc/2013122300/wrfout_d01_2013-12-23_00:00:00"
export Verify_Date_Range="${VERIFY_HOUR}-h FC, 12z 23 Dec - 12z 25 Dec, 2015 (${CYCLE_PERIOD}-h Cycle)"
export OBS_TYPES='synop sound'
export NUM_OBS_TYPES=2
export PLOT_WKS=pdf #"pdf" will save plots in pdf format; "x11" will display the plots and not save them

./da_verif_obs_plot.sh

Now make the script you just made executable using "chmod", then run it.

> chmod +x verif_obs_wrapper.sh

./verif_obs_wrapper.sh

Use "display" to view the plots that have been created (you can find them as .pdf files in the directory "verification"). You should see plots similar to the one at right.

You may try to add the third experiment "conv_and_rad"

WRFDA Tools: Analysis-based verification

This exercise will teach you how to run WRFDA verification tools for observation-based verification. We will be verifying the two assimilation experiments against the "control" experiment with no data assimilation.

WRFDA/var/build/da_verif_grid.exe is the executable that will be used in this session. It will be called directly by scripts, so there is no need to have it in your working directory.

Choose your working directory

For this exercise you should create /home/ec2-user/workdir/verif_grid and use this as your working directory.

mkdir /home/ec2-user/workdir/verif_grid
cd /home/ec2-user/workdir/verif_grid

The easiest way to run grid verification is by creating a "wrapper" script that sets the appropriate environment variables, then calls the script. Open a new file "verif_grid_wrapper.sh" in vi (or your preferred text editor), and type the following:

> vi verif_grid_wrapper.sh
#!/bin/bash -aeu

# Wrapper script for running WRFDA grid verification package

export WRFVAR_DIR=/home/ec2-user/compiled_code/WRFDA-v4.1.2
export SCRIPTS_DIR=/home/ec2-user/WRFDA_TOOLS/scripts
export GRAPHICS_DIR=/home/ec2-user/WRFDA_TOOLS/graphics/ncl
export NUM_EXPT=2
export CONTROL_EXP_DIR=/home/ec2-user/test_data/verification/no_da
export EXP_DIRS='/home/ec2-user/test_data/verification/conv_and_rad /home/ec2-user/test_data/verification/conv_only'
export EXP_NAMES='CONV_RAD CONV_ONLY'
export VERIFICATION_FILE_STRING='wrfout'
export EXP_LEGENDS='(/"Conv and rad obs","Conventional obs"/)'
export START_DATE=2013122312
export END_DATE=2013122512
export INTERVAL=06
export VERIFY_HOUR=06

${SCRIPTS_DIR}/da_verif_grid.sh

Now make the script you just made executable using "chmod", then run it.

> chmod +x verif_grid_wrapper.sh

./verif_grid_wrapper.sh

Use "display" to view the plots that have been created in .pdf format. To the right is one example of the plots you should see. In this example we have specified the input information for both experiments in the wrapper script, so the two experiments are plotted together.

NOTE: Some of the plots will come out rotated 90 degrees. You can click on the images within the "display" program to see a menu which has a "Transform --> Rotate Left" option.

SECOND NOTE: These experiments are very simplified for ease of use on the classroom computer. Don't take the results too seriously.


You are now done with the Tools/Verification practice case. If you still have time, you can go back and try some of the advanced exercises listed at the end of the 4dvar practice page.

 


back to top
 
practice_page_header