<p><b>xinzhang</b> 2008-07-02 15:12:36 -0600 (Wed, 02 Jul 2008)</p><p>Update for Chap6.<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/wrf/UsersGuide/users_guide_chap6.doc
===================================================================
--- trunk/wrf/UsersGuide/users_guide_chap6.doc 2008-07-02 01:41:37 UTC (rev 108)
+++ trunk/wrf/UsersGuide/users_guide_chap6.doc 2008-07-02 21:12:36 UTC (rev 109)
@@ -1,29 +1,28 @@
-ࡱ
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+ࡱ
+
+
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
/
&physics
mp_physics=3,
+sf_sfclay_physics=1,
sf_surface_physics=1,
num_soil_layers=5, (IMPORTANT: its essential to make sure the setting here is consistent with the number in your first guess file)
/
@@ -322,14 +322,23 @@
da_wrfvar.exe >&! wrfda.log
wrfda.log (or rsl.out.0000 if running in distributed-memory mode) contains important WRF-Var runtime log. Always check the log after a WRF-Var run:
*** VARIATIONAL ANALYSIS ***
+ Ntasks in X 1 , ntasks in Y 1
+ *************************************
+ Parent domain
+ ids,ide,jds,jde 1 45 1 45
+ ims,ime,jms,jme -4 50 -4 50
+ ips,ipe,jps,jpe 1 45 1 45
+ *************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
WRF NUMBER OF TILES = 1
Set up observations (ob)
-
+
Using ASCII format observation input
-
+
+ scan obs ascii
+ end scan obs ascii
Observation summary
- ob time 1
+ ob time 1
sound 122 global, 121 local
synop 1007 global, 995 local
pilot 112 global, 112 local
@@ -342,7 +351,7 @@
metar 2551 global, 2543 local
ship 270 global, 260 local
ssmi_rv 0 global, 0 local
- ssmi_tb 0 global, 0 local
+ ssmi_tb 0 global, 0 local
ssmt1 0 global, 0 local
ssmt2 0 global, 0 local
qscat 0 global, 0 local
@@ -351,108 +360,78 @@
bogus 0 global, 0 local
pseudo 0 global, 0 local
radar 0 global, 0 local
+ radiance 0 global, 0 local
airs retrieval 0 global, 0 local
sonde_sfc 122 global, 121 local
+ mtgirs 0 global, 0 local
Set up background errors for regional application
WRF-Var dry control variables are:psi, chi_u, t_u and psfc
Humidity control variable is q/qsg
- Using the latitude-dependent regression coefficients for unbalanced part
-
-Vertical truncation for psi = 7( 99.00%)
-Vertical truncation for chi_u = 13( 99.00%)
-Vertical truncation for t_u = 18( 99.00%)
-Vertical truncation for rh = 19( 99.00%)
+ Using the averaged regression coefficients for unbalanced part
+Vertical truncation for psi = 15( 99.00%)
+
+Vertical truncation for chi_u = 18( 99.00%)
+
+Vertical truncation for t_u = 24( 99.00%!
)
+
+Vertical truncation for rh = 21( 99.00%)
+
+
Calculate innovation vector(iv)
-
+
Minimize cost function using CG method
-For this run cost function diagnostics will not be written
-
+For this run cost function diagnostics will be written
+
Starting outer iteration : 1
-Original gradient is 1.61953446D+03
-For this outer iteration gradient target is 1.61953446D+01
-Starting cost function: 3.72932670D+04 gradient= 1.61953446D+03
+
+Original gradient is 5.11285903D+02
+For this outer iteration gradient target is 5.11285903D+00
+Starting cost function: 3.75916289D+04 gradient= 5.11285903D+02
----------------------------------------------------------
-Iter Gradient Step
- 1 1.13151571D+03 8.79101077D-04
- 2 1.46537384D+03 6.98252394D-04
- 3 1.19033459D+03 6.88681033D-04
- 4 8.78826374D+02 6.17939851D-04
- 5 7.83595468D+02 1.45421881D-03
- 6 6.9!
9681898D+02 1.29468392D-03
- 7 6.39856666D+02 !
2.002971
68D-03
- 8 5.33499758D+02 1.65057269D-03
- 9 4.69363328D+02 2.15663522D-03
- 10 4.01600476D+02 2.53583608D-03
- 11 3.73654859D+02 2.76667326D-03
- 12 3.90839239D+02 1.84643620D-03
- 13 3.23899036D+02 2.99752511D-03
- 14 2.67936599D+02 3.49609161D-03
- 15 2.64310685D+02 5.35852133D-03
- 16 2.23077349D+02 4.23941655D-03
- 17 5.68372735D+02 4.31270668D-03
- 18 2.31310832D+02 2.90459005D-04
- 19 2.09446584D+02 3.16564182D-03
- 20 1.94649895D+02 4.16554352D-03
- 21 1.86709131D+02 5.00033208D-03
- 22 3.78235215D+02 3.49485805D-03
- 23 1.69342578D+02 5.81454290D-04
- 24 1.44526289D+02 5.01503626D-03
- 25 1.53530244D+02 4.89240828D-03
- 26 1.35043290D+02 4.38122156D-03
- 27 1.43780803D+02 5.71128726D-03
- 28 1.95861290D+02 3.72046008D-03
- 29 1.10475835D+02 8.16205!
000D-04
- 30 2.87379421D+02 4.51750717D-03
- 31 1.06941836D+02 2.83923368D-04
- 32 1.04814926D+02 5.97535144D-03
- 33 1.01539664D+02 4.95983149D-03
- 34 1.06306876D+02 1.48102495D-03
- 35 7.92189989D+01 3.68566289D-03
- 36 8.36273942D+01 6.84077832D-03
- 37 8.59028495D+01 4.09124884D-03
- 38 7.84949162D+01 1.08802337D-03
- 39 1.01018622D+02 1.05878433D-03
- 40 5.37227188D+01 2.47360775D-03
- 41 7.04278915D+01 7.89824806D-03
- 42 6.62336908D+01 1.43025762D-03
- 43 9.33770584D+01 3.62044738D-03
- 44 5.55140747D+01 2.72494346D-04
- 45 3.93342642D+01 5.87168964D-03
- 46 4.94905278D+01 9.86698816D-03
- 47 4.03899127D+01 1.88442819D-03
- 48 3.81187624D+01 4.57030374D-03
- 49 2.98708908D+01 6.61640123D-03
- 50 4.86119145D+01 9.34025763D-03
- 51 3.18604235D+01 7.744!
95910D-04
- 52 2.39159452D+01 6.05519031D-03
- 53 !
2.503
98997D+01 9.18225574D-03
- 54 3.23991828D+01 8.61134331D-04
- 55 6.58894613D+01 3.48288133D-03
- 56 1.81965287D+01 2.56097798D-04
- 57 1.70228972D+01 1.21796000D-02
- 58 2.66169790D+01 6.64619203D-03
- 59 2.21045264D+01 2.20570553D-03
- 60 1.87163164D+01 1.30663254D-03
- 61 1.29292558D+01 7.10679099D-03
+Iter Cost Function Gradient Step
+ 1 3.52745747D+04 4.22893275D+02 1.77271355D-02
+ 2 3.36106840D+04 3.02695549D+02 1.86077231D-02
+ 3 3.24416329D+04 3.13509134D+02 2.55182799D-02
+ 4 3.14863535D+04 2.70049441D+02 1.94383780D-02
+ 5 3.07151345D+04 1.97752468D+02 2.11505260D-02
+ 6 3.01681604D+04 1.90679484D+02 2.79738918D-02
+ 7 2.98118845D+04 1.48417585D+02 1.95978521D-02
+ 8 2.94419853D+04 1.13477806D+02 3.35847985D-02
+ 9 2.92227775D+!
04 8.57960292D+01 3.40458452D-02
+ 10 2.90802889D+04 7.76454162D+01 3.87146721D-02
+ 11 2.89758281D+04 6.38614713D+01 3.46538671D-02
+ 12 2.89116212D+04 5.07339439D+01 3.14871956D-02
+ 13 2.88643346D+04 4.24947979D+01 3.67427178D-02
+ 14 2.88301177D+04 3.12058420D+01 3.78964627D-02
+ 15 2.88123982D+04 2.49212482D+01 3.63923716D-02
+ 16 2.87995001D+04 2.08189476D+01 4.15352477D-02
+ 17 2.87903434D+04 1.59621535D+01 4.22523996D-02
+ 18 2.87843493D+04 1.30676595D+01 4.70506154D-02
+ 19 2.87805743D+04 9.41818099D+00 4.42135584D-02
+ 20 2.87784444D+04 7.45330718D+00 4.80233115D-02
+ 21 2.87770601D+04 6.25547789D+00 4.98386142D-02
+ 22 2.87762347D+04 4.79613214D+00 4.21877119D-02
----------------------------------------------------------
-
-Inner iteration stopped after 61!
iterations
-
-Final: 61 iter, J= 2.99744265D+04, g= 1.29292!
558D+01
+
+Inner iteration stopped after 22 iterations
+
+Final: 22 iter, J= 2.87762347D+04, g= 4.79613214D+00
----------------------------------------------------------
-
+
+rsl.out.0000 (p4 of 4)
Diagnostics
- Final cost function J = 29974.43
-
- Total number of obs. = 68544
- Final value of J = 29974.42650
- Final value of Jo = 28880.81151
- Final value of Jb = 1093.61499
+ Final cost function J = 28776.23
+
+ Total number of obs. = 69249
+ Final value of J = 28776.23467
+ Final value of Jo = 26116.53564
+ Final value of Jb = 2659.69904
Final value of Jc = 0.00000
Final value of Je = 0.00000
- Final J / total num_obs = 0.43730
+ Final value of Jp = 0.00000
+ Final J / total num_obs = 0.41555
Jb factor used(1) = 1!
.00000
Jb factor used(2) = 1.00000
Jb factor used(3) = 1.00000
@@ -460,10 +439,12 @@
Jb factor used(5) = 1.00000
Jb factor used = 1.00000
Je factor used = 1.00000
-
+ VarBC factor used = 1.00000
+
SUCCESS COMPLETE WRFVAR
-
+
*** WRF-Var completed successfully ***
+
A file called namelist.output (which contains the complete namelist settings) will be generated after a successful da_wrfvar.exe run. Those settings appear in namelist.output that are not specified in your namelist.input are the default values from WRFDA/Registry/Registry.wrfvar.
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.
@@ -748,7 +729,7 @@
-WRF-ARW V3: Users Guide 6- PAGE 40
+WRF-ARW V3: Users Guide 6- PAGE 18
WRF-ARW V3: Users Guide 6- PAGE 1
@@ -772,14 +753,14 @@
-
-h
-
+h{
+
-&
+&
-h
+h{
+
&
-&
(
-Px 4 #\'*.25@9
-/
h?
-&
-&
-v
-v
-
-h
+Px 4 #\'*.25@9
+/
h?
+&
+&
+v
+v
+
+h{
l
-l
+l
l
-l
+l
l
-l
+l
l
-l
-h
-f
x
-
OJ
+h{
<font color="gray">
+
-
-
-h
-h
-&
+`tgd{
gd{
+
+h{
+h{
#
-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
#
-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
+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
+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
+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
+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
+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
+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
l
-
+
#
-l
l
-
-h
+
+h{
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
+l
#
-l
l
-
+
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
+l
d#
-l
-
-h
+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
+l
d#
-l
+l
d#
-l
-h
+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
-
+
+
#
-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
-
-
-h
-h
+
+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
-
-
-
-
-
-
-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
-
-
-
+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
+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
+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
D
)].w
@@ -2646,7 +2633,7 @@
=Lw(
cl t0*V*z2 uq*L2U垧JL
y
-
+
l
l
-
+
u(lH<CX.Ip=#H@6YYY999- _Q*<yIiRHuea{233]+5y;wѣGK"ʦV6_|^xA
-
+
@@ -3946,8 +3933,7 @@
~dCXx祗^RF{^;aY!_~Y7c \3Ƃ;N42~q ,>
? ? FcpX
4H wPs,y'Ђo]:!FN"0آbn;B1vZG܂ h O_~C{QNS (:-:?yN#
-
-
+
@@ -4047,7 +4033,8 @@
cWA t@J^f.?Pscp,R ƤyP.CYKL Cy>
{:o۶p[nՙ|{gg~vy(hs h'ӻTU@@izK!
H" " ' }" " ;@ HzĪN> " "{p0@> |@> VO
-
+
+
@@ -5126,7 +5113,7 @@
1cG]OOI-< DuvHzHi9qUV۳z<x ;5c=6 ^ >%PqZ
|z>?5Nc5ċĂd:'Q_ rm:5\?XZj%5MDI(Yczjlet~OQ͢m;ꨣbUG% b&#>qPeN0oY!gCO2LJ {= o-XiBA@A@9BHrr Dt@A@r7\ ! A@# $M BHDA@A !y A@A@9$+4 DDDDDDr;E
-
+
/
@@ -6695,8 +6682,7 @@
dT@@Q}(R@
䣀" R@
H)Q@QL@u@u@u@u@u@u:HnMu@u@u@u@u`U ]B
-
+
YU
+
+
+
+
+
ư
-
-`
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
ؒC
0
}j
x
0
r
0
㻐
㻐
p
-`
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
0
0
-
-
-
-
-
+
0
0
+
+
+
+
-
+
+
-
+
-
+
-L
-
-E
-
-
-
+
+
+)
+X
+y
+
+
+
+
+
+E
+z
+
+
+
+
+
+
+
l,6
-
+
^Z
-^Z
+^Z
@@ -10493,72 +10522,71 @@
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
+
-
-
-5O
-
+
+
+
+
wwo~wo~o~ww
ZkZVZkZ^ZcF1^ZcZZg9kZZkZZ^VNsg9kZRVc^ZVco{o{kZRkZJRc^ZR^g9VZ^Vg9R^VkZNsR
swo{wo{o{ww
-wwwww
+wwwww
-
+
+
+
+
@@ -10745,7 +10773,8 @@
-
+
+
@@ -10756,17 +10785,34 @@
-
+
+
-
-
+
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
\ No newline at end of file
+
+
+
\ No newline at end of file
</font>
</pre>