Reference: Download the tutorial presentation
Source code
Get the pre-compiled code, if you have not done so.
WRFDA/var/build/da_wrfvar_3dvar.exe is the executable that will be used in this session.
Important note: If you downloaded the pre-compiled code the first day, it was not compiled with RTTOV. If you wish to do this exercise with RTTOV you should get the most recent version of the code at /kumquat/wrfhelp/SOURCE_CODE/WRFDA/WRFDA_compiled.tar.gz
Choice of your working directory
You should run each lesson in a separate directory, so that it will be easier to check for the necessary input files and look for what output files are created after a successful run.
We recommend you create /kumquat/users/${USER}/DA/radiance to be your working directory for this session.
mkdir /kumquat/users/${USER}/DA/radiance
cd /kumquat/users/${USER}/DA/radiance
Input files
The radiance case is the same CONUS60 case you have used in obsproc, running WRFDA-3DVAR and gen_be sessions. Keep in mind that when compiling WRFDA yourself, additional settings (specifically, environment variables for BUFR and an RTM) must be specified at compile time to utilize radiance data. This has been done for you in the pre-compiled code.
Radiance assimilation for 3DVAR is much the same procedure as Running WRFDA-3DVAR with conventional observations. To begin this exercise, link/copy the basic input files above that you should have been familiar with by now:
ln -sf /kumquat/users/${USER}/DA/WRFDA/run/LANDUSE.TBL ./LANDUSE.TBL
ln -sf /kumquat/users/${USER}/DA/obsproc/obs_gts_2008-02-05_12:00:00.3DVAR ./ob.ascii
ln -sf /kumquat/users/${USER}/DA/gen_be/be.dat ./be.dat
(use /kumquat/wrfhelp/DATA/WRFDA/CONUS60/be/be.dat if you have not done the GEN_BE exercise)
cp -p /kumquat/wrfhelp/DATA/WRFDA/CONUS60/rc/2008020512/wrfinput_d01 ./fg
In addition, the following extra files are required for radiance assimilation:
- Radiance data in NCEP BUFR format
ln -fs /kumquat/wrfhelp/DATA/WRFDA/CONUS60/ob/2008020512/gdas1.t12z.1bamua.tm00.bufr_d.linux ./amsua.bufr
ln -fs /kumquat/wrfhelp/DATA/WRFDA/CONUS60/ob/2008020512/gdas1.t12z.1bamub.tm00.bufr_d.linux ./amsub.bufr
- radiance_info files (you must link the whole directory)
ln -fs /kumquat/users/${USER}/DA/WRFDA/var/run/radiance_info ./radiance_info
- File for variational bias correction
ln -fs /kumquat/users/${USER}/DA/WRFDA/var/run/VARBC.in ./VARBC.in
- Radiative Transfer Model (CRTM or RTTOV) coefficient files (you must link the whole directory)
ln -fs /kumquat/users/${USER}/DA/WRFDA/var/run/crtm_coeffs ./crtm_coeffs
ln -fs /kumquat/wrfhelp/WRFDA/external/rttov11/rtcoef_rttov11/rttov7pred54L ./rttov_coeffs
Edit namelist.input
A sample namelist.input is provided in the WRFDA package as WRFDA/var/test/radiance/namelist.input
cp /kumquat/users/${USER}/DA/WRFDA/var/test/radiance/namelist.input ./namelist.input
vi namelist.input (pay special attention to the following radiance-related settings)
&wrfvar4
use_amsuaobs = .true.
use_amsubobs = .true.
/
&wrfvar14
rtminit_nsensor = 6,
rtminit_platform = 1,1,1,1,1,1
rtminit_satid = 15,16,18,15,16,17
rtminit_sensor = 3,3,3,4,4,4
thinning_mesh = 120.0,120.0,120.0,120.0,120.0,120.0,
thinning = .true.
qc_rad = .true.
write_iv_rad_ascii = .false.
write_oa_rad_ascii = .true.
rtm_option = 2 (set to "1" for RTTOV)
only_sea_rad = .false.
use_varbc = .true.
use_crtm_kmatrix = .true. (set "use_rttov_kmatrix" instead for RTTOV)
varbc_nbgerr = 5000
/
&wrfvar21
time_window_min = "2008-02-05_11:00:00"
/
&wrfvar22
time_window_max = "2008-02-05_13:00:00"
/
Run WRFDA
mpirun -np 4 /kumquat/users/${USER}/DA/WRFDA/var/build/da_wrfvar_3dvar.exe
Check output
Check analysis increments.
Check statistics.
There are extra diagnostic files when assimilating radiance data:
- 01_oma_*amsu*: Observation minus analysis statistics for each radiance instrument.
- 01_qcstat_*amsu*: Quality control statistics for each radiance instrument.
- VARBC.out: Output file for Variational Bias Correction. This process is described in further detail in the User's Guide
Graphics
WRFDA/var/scripts/da_rad_diags.ksh
export START_DATE=2008020512
export END_DATE=2008020512
export CYCLE_PERIOD=12
#
export NUM_PROCS=4 # number of processors used in running WRF-Var
export WRFVAR_DIR=WRFVAR_DIR=/kumquat/users/($USER)/DA/WRFDA
export BUILD_DIR=${WRFVAR_DIR}/var/build
export EXP_DIR=/kumquat/users/$USER/DA/radiance
export VAR_RUN_DIR1=$EXP_DIR
export VAR_RUN_DIR2=none
...
export FGFILE=$VAR_RUN_DIR1/fg
...
export FILE_PREFIX=01_oma
Also change the following line (around line 150):
ln -sf $VAR_RUN_DIR1/$DATE/$FILE_PREFIX* . --> ln -sf $VAR_RUN_DIR1/$FILE_PREFIX* .
The .ncgm plots will be created in your working directory. Use the command "idt" to view these .ncgm files.
Other practice
- Try linking /kumquat/wrfhelp/DATA/WRFDA/CONUS60/ob/2008020512/gdas1.t12z.1bamua.tm00.bufr_d as amsua.bufr in your working directory and run WRFDA to see what log messages you will get.
Advanced
