Case 1: ARW Case with conventional data (PrepBUFR)
Case 1: ARW Case with conventional data (PrepBUFR) cindyhg Tue, 08/20/2019 - 12:41Case 1: ARW Case with conventional data (PrepBUFR)
Introduction
This case of 00z analysis on February 13th, 2014 was used here as the example case to introduce how to setup run scripts for the GSI observer and EnKF and how to check the EnKF analysis results to determine if a particular EnKF application was successful. In this case, only conventional observations from PrepBUFR are assimilated. Information on assimilating more data will be introduced in the tutorial case 02.
Case Data
The information on downloading the case data can be found from this link.
Setting Up the Case Data and Run Directory
- Create the data directory ./enkfdata and download, untar the data files. The case data is now populated in:
(PATH)/enkfdata/arw
- Create a working directory or use an existing directory as the working directory.
- In this case, we use run as working directory.
Setting Up the GSI Observer Run
Setting Up the GSI Observer Run cindyhg Tue, 08/20/2019 - 12:46Case 1: ARW Case With Conventional Data (PrepBUFR)
Setting Up the GSI Observer Run
Running GSI as observer is the same as running GSI as 3DVar analysis, except for couple of run scripts setups. Users should create a a basic GSI run script following the GSI tutorial: Prepare Base Run Script.
Based on run_gsi_regional.ksh_basic, make the following changes to run GSI as observer:
- Define run GSI observer only and ensemble members
- When set if_observer=Yes, GSI only runs observation operator part and skips the minimaztion:
if_observer=Yes no_member=20 BK_FILE_mem=${BK_ROOT}/wrfarw.mem- For this exercise, only conventional observations are assimilated. To do so, comment out the links to the other observations
# ln -s ${srcobsfile[$ii]} ${gsiobsfile[$ii]}The example script can be found : here
Run GSI observer
Run GSI observer cindyhg Tue, 08/20/2019 - 12:47Case 1: ARW Case With Conventional Data (PrepBUFR)
Run GSI observer
How to run GSI job depends on computer system. On Cheyenne (LINUX with PBS), run it by:
qsub run_gsi_regional.ksh
Check the GSI observer run results
Check the GSI observer run results cindyhg Tue, 08/20/2019 - 12:48Case 1: ARW Case With Conventional Data (PrepBUFR)
Check the GSI observer run results
cd ${WORK_ROOT}
- Check the stdout file from member 001 and compare it to the reference stdout here
- Check the files in the run directory for member 001 before the directory clean and compare them to reference file list here
- Check the files in the run directory for member 020 before the directory clean and compare them to reference file list here
- No need to check fit files and minimization
- should have diag files for the previous 19 members and mean
Setting Up the EnKF Run
Setting Up the EnKF Run cindyhg Tue, 08/20/2019 - 12:50Case 1: ARW Case With Conventional Data (PrepBUFR)
Setting Up the EnKF Run
Here, we will setup the the EnKF run script based on the basic run script from Prepare basic EnKF run script. After copy:
cp run_enkf_wrf.ksh_basic run_enkf_wrf.ksh
Make the following additional modifications to the script run_enkf_wrf.ksh:
Define ensemble parameters
NMEM_ENKF=20
BK_FILE_mem=${BK_ROOT}/wrf_en
NLONS=129
NLATS=70
NLEVS=50
Define ensemble background and observations used in analysis
IF_ARW=.true. IF_NMM=.false. list="conv"
The example script can be found : here
Run EnKF
Run EnKF cindyhg Tue, 08/20/2019 - 12:53Case 1: ARW Case With Conventional Data (PrepBUFR)
Run EnKF
It can be run with the same way as GSI runs, for example on Cheyenne (with PBS):
qsub run_enkf_wrf.ksh