<p><b>xinzhang</b> 2009-03-25 11:18:39 -0600 (Wed, 25 Mar 2009)</p><p>Save for Chap. 6 <br>
</p><hr noshade><pre><font color="gray">Modified: trunk/wrf/UsersGuide/users_guide_chap6.doc
===================================================================
--- trunk/wrf/UsersGuide/users_guide_chap6.doc 2009-03-23 22:47:21 UTC (rev 150)
+++ trunk/wrf/UsersGuide/users_guide_chap6.doc 2009-03-25 17:18:39 UTC (rev 151)
@@ -1,5 +1,5 @@
-ࡱ
-
T
+
./configure
serial means single processor
dmpar means Distributed Memory Parallel (MPI)
+Note: wrfplus was tested on following platforms:
+IBM AIX: xlfrte 11.1.0.5
+Linux : pgf90 6.2-5 64-bit target on x86-64 Linux
+Mac OS (Intel) : g95 0.91!
+Note: wrfplus does not support:
+Linux: Intel compiler V9.1 (not sure for higher versions, WRFPLUS can not be compiled with old version)
+Linux : gfortran (The behavior of WRFPLUS is strange)
./compile em_real
ls -ls main/*.exe
you should see wrfplus.exe
@@ -268,7 +275,7 @@
Now you have all the three input files (first guess, observation and background error statistics files in directory $DAT_DIR) required to run WRF-Var. Also, you have successfully downloaded and compiled the WRF-Var code. If this is correct, we are ready to learn how to run WRF-Var.
b. Run The Case3D-Var
The data for this case is valid at 12 UTC 5th February 2008. The first guess comes from the NCEP global final analysis system (FNL), passed through the WRF-WPS and real programs.
-To run WRF-Var, first create a working directory, for example, WRFDA/var/test, then follow the steps below:
+To run WRF-3D-Var, first create a working directory, WRFDA/var/test/tutorial is a sample, then follow the steps below:
cd WRFDA/var/test/tutorial
ln -sf WRFDA/run/LANDUSE.TBL ./LANDUSE.TBL
ln -sf $DAT_DIR/rc/2008020512/wrfinput_d01 ./fg (link first guess file as fg)
@@ -480,7 +487,108 @@
After successful completion of job, wrfvar_output (the WRF-Var analysis file, i.e. the new initial condition for WRF) should appear in the working directory along with a number of diagnostic files. Various text diagnostics output files will be explained in the next section (HYPERLINK \l "_WRF-Var_Diagnostics_1"WRF-Var Diagnostics).
In order to understand the role of various important WRF-Var options, try re-running WRF-Var by changing different namelist options. For example, making WRF-Var convergence criteria more stringent by reducing the value of EPS to e.g. 0.0001 by adding "EPS=0.0001" in namelist.input record &wrfvar6. Last section of this tutorial deals with many such WRF-Var additional exercises.
b. Run The Case4D-Var
+To run WRF-4DVar, first create a working directory, WRFDA/var/test/4dvar is a sample,
+Assume the working shell is C shell;
+Assume the directories for WRFDA , WRFNL and WRFPLUS are :
+setenv WRFDA_DIR /ptmp/$user/W
RFDA
+setenv WRFNL_DIR /ptmp/$user/WRFNL
+setenv WRFPLUS_DIR /ptmp/$user/WRFPLUS
+Assume the analysis date is 2008020512 and the test data directories are:
+setenv DATA_DIR /ptmp/$user/DATA
+> ls lr $DATA_DIR
+ob/2008020512
+ob/2008020513
+ob/2008020514
+ob/2008020515
+ob/2008020516
+ob/2008020517
+ob/2008020518
+rc/2008020512
+be
+Note: Currently, WRF-4DVar can only run with the observation data processed by OBSPROC, cannot work with PREPBUFR format data; WRF-4DVar can work with satellite radiance BUFR data, but the this capability is still under testing.
+Assume the working directory is:
+setenv WORK_DIR $WRFDA_DIR/var/test/4dvar
+Then follow the steps below:
+1) Link the executables.
+cd $WORK_DIR
+ln -fs $WRFDA_DIR/var/da/da_wrfvar.exe .
+cd $WORK_DIR/nl
+ln -fs $WRFNL_DIR/main/wrf.exe .
+cd $WORK_DIR/ad
+ln -fs $WRFPLUS_DIR/main/wrfplus.exe .
+cd $WORK_DIR/tl
+ln -fs $WRFPLUS_DIR/main/wrfplus.exe .
+2) Link the observational data, first guess and BE.
+cd $WORK_DIR
+
ln -fs $DATA_DIR/ob/2008020512/ob.ascii+ ob01.ascii
+ln -fs !
$DATA_DI
R/ob/2008020513/ob.ascii ob02.ascii
+ln -fs $DATA_DIR/ob/2008020514/ob.ascii ob03.ascii
+ln -fs $DATA_DIR/ob/2008020515/ob.ascii ob04.ascii
+ln -fs $DATA_DIR/ob/2008020516/ob.ascii ob05.ascii
+ln -fs $DATA_DIR/ob/2008020517/ob.ascii ob06.ascii
+ln -fs $DATA_DIR/ob/2008020518/ob.ascii- ob07.ascii
+ln -fs $DATA_DIR/rc/2008020512/wrfinput_d01 .
+ln -fs $DATA_DIR/rc/2008020512/wrfbdy_d01 .
+ln -fs wrfinput_d01 fg
+ln -fs wrfinput_d01 fg01
+
+ln -fs $DATA_DIR/be/be.dat .
+
+3) Establish the miscellaneous links.
+cd $WORK_DIR
+ln -fs nl/nl_d01_2008-02-05_13:00:00 fg02
+ln -fs nl/nl_d01_2008-02-05_14:00:00 fg03
+ln -fs nl/nl_d01_2008-02-05_15:00:00 fg04
+ln -fs nl/nl_d01_2008-02-05_16:00:00 fg05
+ln -fs nl/nl_d01_2008-02-05_17:00:00 fg06
+ln -fs nl/nl_d01_2008-02-05_18:00:00 fg07
+
+ln -fs ad/ad_d01_2008-02-05_12:00:00 gr01
+
+ln -fs tl/tl_d01_2008-02-05_13:00:00 tl02
+ln -fs tl/tl_d01_2008-02-05_14:00:00 tl03
+ln -fs tl/tl_d01_2008-02-05_15:00:00 tl04
+ln -fs tl/tl_d01_2008-
02-05_16:00:00 tl05
+ln -fs tl/tl_d01_2008-02-05_17:00:00 tl06
+ln -fs tl/tl_d01_2008-02-05_18:00:00 tl07
+
+cd $WORK_DIR/ad
+ln -fs ../af01 auxinput3_d01_2008-02-05_12:00:00
+ln -fs ../af02 auxinput3_d01_2008-02-05_13:00:00
+ln -fs ../af03 auxinput3_d01_2008-02-05_14:00:00
+ln -fs ../af04 auxinput3_d01_2008-02-05_15:00:00
+ln -fs ../af05 auxinput3_d01_2008-02-05_16:00:00
+ln -fs ../af06 auxinput3_d01_2008-02-05_17:00:00
+ln -fs ../af07 auxinput3_d01_2008-02-05_18:00:00
+4) Run with 1 processor (serial compilation required)
+Edit $WORK_DIR/namelist.input to match your experiment settings.
+Edit $WORK_DIR/nl/namelist.input to match your experiment settings.
+Edit $WORK_DIR/ad/namelist.input and $WORK_DIR/tl/namelist.input to match your experiment settings, but only change following variables:
+&time_control
+run_hours=06,
+start_year=2008,
+start_month=02,
+start_day=05,
+start_hour=12,
+end_year=2008,
+end_month=02,
+end_day=05,
+end_hour=18,
+......
+&domains
+time_
step=360, # NOTE:MUST BE THE SAME WITH WHICH IN $WORK_DIR/!
nl/namel
ist.input
+e_we=90,
+e_sn=60,
+e_vert=41,
+dx=60000,
+dy=60000,
+......
+
+cd $WORK_DIR
+setenv NUM_PROCS 1
+./da_wrfvar.exe >&! wrfda.log
+
WRF-Var Diagnostics
WRF-Var produces a number of diagnostic files that contain useful information on how the data assimilation has performed. This section will introduce you to some of these files, and what to look for.
Having run WRF-Var, it is important to check a number of output files to see if the assimilation appears sensible.
@@ -576,7 +684,7 @@
Under WRFDA/var/scripts/gen_be, gen_be_wrapper.ksh is used to generate the BE data, following variables need to be set to fit your case:
-WRFVAR_DIR=/karri/users/xinzhang/support/WRFDA #WRFDA path
+WRFVAR_DIR=/karri/users/noname/support/WRFDA #WRFDA path
REGION=con200 #Region ID
EXPT=expt # Experiment ID
EXP_DIR=/karri/users/xinzhang/support/$REGION/$EXPT
@@ -749,10 +857,10 @@
-WRF-ARW V3: Users Guide 6- PAGE 10
+WRF-ARW V3: Users Guide 6- PAGE 24
-WRF-ARW V3: Users Guide 6- PAGE 11
+WRF-ARW V3: Users Guide 6- PAGE 25
@@ -773,1443 +881,1451 @@
-
-h.
-
-&
+h:
+
+&
-
-h.
+
+h:
&
-&
-Px 4 #\'*.25@9
+Px 4 #\'*.25@9
& F
-& F
+& F
&
-&
-h.
-&
+&
+h:
+& F
+&
&
-& F
-& F
-h.
-&
-h.
+& F
+& F
+h:
+&
+& F
+&
+& F
+I
+h:
l
-l
+l
l
-l
+l
l
-l
+l
l
-l
hj
-t
-x
-
-
-h.
-h.
-h.
-&
-l
-
+l
J
+
+
+h:
+&
+
+
+h:
+h:
#
-l
-
-h.
+l
+&
+l
+
+h:
#
-l
+l
#
-l
+l
#
-l
-
+l
#
-l
-
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
-h.
+
#
-l
+l
+h:
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
l
+
#
-l
+l
+
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
-h.
+l
#
-l
+l
+
+h:
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
-h.
+l
#
-l
+l
#
-l
+l
+
+h:
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
+
#
-l
+l
l
+
#
-l
+l
#
-l
+l
#
-l
-h.
+l
#
-l
+l
+h:
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
-h.
+l
#
-l
+l
#
-l
+l
#
-l
+l
+h:
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
-
-h.
+l
#
-l
+l
+
+h:
#
-l
+l
#
-l
+l
+
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
#
-l
#
-l
+l
#
-l
-
+l
#
-l
+l
#
-l
+l
#
-l
-h.
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
+
+
+
+
+h:
#
l
+
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
#
l
+
+
l
-
-h.
+
+
+
+
+
+
+
+
+
#
-l
+l
+
+
+
+
+
+
+
+
+
+
+h:
+
+
+
+
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
d#
l
+
+
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
-h.
+l
d#
-l
+l
d#
-l
+l
+
+
+h:
d#
-l
+l
d#
-l
+l
+
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
+
d#
-l
-
-h.
+l
+
d#
-l
+l
d#
-l
+l
+h:
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
-
+l
l
+
d#
-l
+l
#
l
+
#
-l
+l
#
-l
-
-h.
+l
#
-l
+l
+
+h:
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
l
+
#
-l
+l
l
#
l
+
l
#
l
+
#
l
-
+
+
#
l
+
l
-
+
-
+
l
-
-
-
-
-
-
-
-
-
-
-h.
-h.
+
l
+
#
-l
+l
+h:
+h:
l
#
l
-
-
-
-
-
-
-
+
#
l
-
-
-
-
-
+
#
-l
-
-
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
-&
+l
+&
#
-l
+l
l
+
#
-l
+l
#
l
+
#
-l
+l
#
l
+
l
+
#
-l
-h.
+l
+h:
#
l
+
#
-l
+l
l
#
l
+
#
l
+
#
-l
+l
#
l
l
+
#
l
+
#
-l
+l
#
l
#
l
#
l
+
#
l
+
#
-l
+l
#
l
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
-h.
+l
#
l
#
l
+
+h:
l
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
l
#
l
+
#
l
+
#
-l
+l
m
#
@@ -2217,165 +2333,163 @@
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
-h.
+
+h:
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
-$
+
#
l
+
#
-l
-$
-h.
+l
+h:
#
l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
l
+
)].w
@@ -2586,7 +2700,7 @@
eԳw@3]#;b69ώ
=Lw(
-cl t0*V*z2 uq*L2U垧JL
w
+cl t0*V*z2 uq*L2U垧JL
w
@@ -2733,8 +2847,8 @@
^AݔV6#*
pg^%
31
-
-
+
+
bWV
g0a@ >Dž>gz_ V:\n~{Iu>f{a {y®k'\"ᠰ~hǵʰ. kk~c})l6#쨰^:n|6ʰ&8=, Wk1[aͣքa;;̵f2,õo? ,d-m؛5apͰz0 eް^kXcusV=VxG(߾{X ~a{a
;2l]ذ0o ~a&
@@ -3293,7 +3407,7 @@
q%)]? vaWh u6rpi^; m ۺpӞm - /* +, .+ ,F>B.,I?{pr9
uq};z̾c۸~|o L
F
-Yn~+GV` H`
+Yn~+GV` H`
@@ -3511,8 +3625,8 @@
l
l
-
-
+
c?7
gݾ(C+0sb!}Q1
E<f߿s <{B;wjhw9 ҥK 9hK/˗/GkQ|MKV㵠;Wx5|Chaq,1D /믑=v<)RP:tH`,1ܹ3]#iY0&ؘ _!1<mǥ>$B!B!B!$qpps?_R###C4O 320{ͺ8_>up8~|;/"|v:zOxx7 <<nAᥝG8Cq~z%qջ&ײU|);鴤XNh9G>ϥPWy%]GHIߙ;_onSs]2YZοw⳿8.ۧs>Gο+v 7չN/pݶw(ک_Ia[!\ pT ܿ? "W ,οE])ty1V#k dohp(p}HZ9}vA^!!?8 v@v]EYcG gi,sק+G οUw$yq"u}{wWv һ ɜ%MxMftz}nР߬#Q4} 8ȓQv Wѽkj8Tdʖ~Y7!H
-
-
+
+
\)QTD"tOHG/ET^"T] WyO}f7l3g flyg9sεq"E%Y#vw6 ( ɰZF;۱3qE5߶Wi#E9UH """"""r^^}UiѢmVL".\s >|X (.]RȲe9}[2k,[7Ȅ \? 6LVZe> ?>S_۷-{zwd~_ʊDtϟ
HDD(QBy]KN
.,۷os^Z=S2eȀ,Aɮ]|L
@@ -3732,8 +3845,9 @@
6<gMܻ; 4h`ˌa{I 7zY,`}}bxm]n]Sd/+R&>fr|l+:vh Y0nܸύ|_vڵkcڴi%p,3>}{mۖfK* NjSre ? _~ז:"Yw}gƶJI؎ 9bK.(Z~dsi`&f֬YXpqRYkm KL}YQog_ƾ,piuks eX%"NvC:ׯo ۗTԫW1W^f
<>M6PɯӧO}_wuRJq߃'% sޒ%K{n/66s=gْǴ#͘lwfN{\=} <xiߧ:X:}l>oʔ)祤}\ 7N: ]CE\p6֑>} ߲eo\/[/l Kl,Yg d(Qo_\f=ތh͚5~mqYccMl4-1 W͌qu2)Co?sGc3&|u<'g˗/cǎ wneI 2_5o6meAA҅][W"]NMl$EPœGA OG=xAkO B)Mư ]wx'ɼod FC QţŖ ,f dYn(C}u bEBv"8<nb*<O@*f|~4z:A~
&Zm=f kTЊMɳ8m6-ШW_uIUQr7廫%r< 7@|',+.K⛮܉Q . }QJw/ߓ_~մWځVy"MiWI
-
-
+
+
?uXMLLDFtO>EҥƘ?^oٳnkطoҏ# #Z-n߾m'OF``sPK} X`cbԨQ4+WD߾}z
@@ -4206,8 +4320,8 @@
ZoTZnc\<qR`Y CZn, %TG >VF@OA
Om=ttzfZFv=+ Mӝݺgieͷ,-(}-Ag +<ϙEd^Kd;}vV4Rn|>GLw-ޒlD&r;˱ =^܀d,`+tgY~/F >;rc Zoko,?bmW3hvU &j ~z-B嗷DJm_ R_DuDsay=kM|3J}D%-'vh&z8AT>gu-,)ҢJlÈuٲcٺqOkrΞΊVNv+y=+zW~PT2eLYdNu+7XSOtwOv88{k 5+[qr_)ЎroN8Ф jn>h )ښm+C
-
-
+
+
J&H(yA@(( >φ4@(zoRAJhUx"|o ܒ{M9Aٹ͙9sf !3'J̓
6֭GAW42D'FhC8YJ3Kͼƪc']^AA_FaaϱRrQ0VF6ꉆQiaěnF֭-PQ6.cLzfwO?Od oQT1l[TQc
Jߙ7gdct#H6\FoY!1N)#*zUS
GN]fݝ1bdy7qڬU ԝR?u
-
-
+
+
@@ -7462,11 +7576,11 @@
-
+
+
-
-
-
-
+
+
+
+
+
+
@
@
Dư
@
-
-n
-r
-}
-
)
-
-
-
+
<font color="red">
+z
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
@5
-
g
0ŀ
ŀ
0Š
0Š
0
0Š
+
+
+
+
+
+
+
+
+
iҀ
Ҁ
0Š
0
-
)
-
-
-
+
<font color="gray">
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
+
+
+
+
+
+
-
+
@
@@ -7577,112 +7825,171 @@
@
@
-@
0
0
-
-
-
-
-
+@
0
0
+
+
+
+
+
-
+
+
+
-
+
+
-&
-A
-X
-a
-j
-n
-
-
-
-
-
-
-
-
-
+&
+Y
+f
+x
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
+
+
+
-n
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
+
+
+
-
-
-
+
+
+
+
+
+
+
+
-
-
-
+
+
+
-
+
+
+(
:
l,6
^Zp;sZ+8
-AB ^bn]
+AB ^bn]
Z
@@ -7809,7 +8116,7 @@
-
+
@@ -7862,9 +8169,9 @@
-
+
^Z
-^Z
+^Z
@@ -7872,60 +8179,125 @@
-
-`GT^5bBJ>*N
-
-
-
-
+
+`GT^5bBJ>*N
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
+
+
+
+
+
+
+
-
-wwo~wo~o~ww
ZJRo{Zg9RR^cRJR^Z^kZJRo{R
+
+
+wwo~wo~o~ww
JRo{Zg9RR^cRJR^Z^kZJRo{R
RZg9kZ^ZRNswZNsVZRcwkZRNskZJR^ckZZg9s^VkZVZRs
swo{wo{o{ww
-wwwww
+wwwww
+
r
-
+
-
+
+
-
\ No newline at end of file
+
+
\ No newline at end of file
</font>
</pre>