<p><b>xinzhang</b> 2009-04-09 15:23:07 -0600 (Thu, 09 Apr 2009)</p><p>Bug fix 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-04-07 23:54:30 UTC (rev 189)
+++ trunk/wrf/UsersGuide/users_guide_chap6.doc 2009-04-09 21:23:07 UTC (rev 190)
@@ -1,16 +1,7 @@
-ࡱ
+ࡱ
-
+
> mv obs_gts_2008-02-05_17:00:00.4DVAR $DAT_DIR/ob/2008020517/ob.ascii
> mv obs_gts_2008-02-05_18:00:00.4DVAR $DAT_DIR/ob/2008020518/ob.ascii-
At this pont you have three of the input files (first guess, observation and background error statistics files in directory $DAT_DIR) required to run WRF-Var, and 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
+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-3D-Var, first create and cd to a working directory, for example, WRFDA/var/test/tutorial, and then follow the steps below:
> cd WRFDA/var/test/tutorial
@@ -1032,7 +1023,7 @@
; to assimilate radiance data.use_hirs2obsfasle.true.: to read in data from hirs2.bufruse_hirs3obsfalse.true.: to read in data from hirs3.bufruse_hirs4obsfalse.true.: to read in data from hirs4.bufruse_mhsobsfalse.true.: to read in data from mhs.bufruse_msuobsfalse.true.: to read in data from msu.bufruse_amsuaobsfalse.true.: to read in data from amsua.bufruse_amsubobsfalse.true.: to read in data from amsub.bufruse_airsobsfalse.true.: to read in data from airs.bufruse_eos_amsuaobsfalse.true.: to read in data from airs.bufruse_hsbobsfalse.true.: to read in data from hsb.bufruse_ssmisobsfalse.true.: to read in data from ssmis.bufruse_obs_errfacfalse.true.: apply obs error tuning factors if errfac.dat is available for conventional data only&wrfvar5
check_max_ivtrue.true.: reject the observations whose innovations (O-B) are larger than a maximum value defined as a multiple of the observation error for each observ
ation. i.e., inv > (obs_error*factor) --> fails_error_max; the default maximum value is 5 times the observation error ; the factor of 5 can be changed through max_error_* settings.max_error_t5.0maximum check_max_iv error check factor for tmax_error_uv5.0maximum check_max_iv error check factor for u and vmax_error_pw5.0maximum check_max_iv error check factor for precipitable watermax_error_ref5.0maximum check_max_iv error check factor for gps refractivitymax_error_q5.0maximum check_max_iv error check factor for specific humiditymax_error_p5.0maximum check_max_iv error check factor for pressuremax_error_thicknessmaximum check_max_iv error check factor for thicknessmax_error_rvmaximum check_max_iv error check factor for radar radial velocitymax_error_rfmaximum check_max_iv error check factor for radar reflectivity&wrfvar6max_ext_its1number of outer loopsntmax200maximum number of iterations in an inner loopeps
0.01 (max
_ext_its)minimization convergence criterion (used dimension!
: max_ex
t_its); minimization stops when the norm of the gradient of the cost function gradient is reduced by a factor of eps. inner minimization stops either when the criterion is met or when inner iterations reach ntmax.&wrfvar7cv_options53: NCEP Background Error model
-5: NCAR Background Error model (default)as1(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 1 = stream function. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.25.as2(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 2 - unbalanced potential velocity. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.25.as3(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 3 - unbalanced temperature. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.25.as4(3) -1.0tuning factors for variance, horizontal and vertical scales for control variable 4 -
pseudo relative humidity. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.25.as5(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 5 - unbalanced surface pressure. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.25.rf_passes6number of passes of recursive filter.var_scaling11.0tuning factor of background error covariance for control variable 1 - stream function. For cv_options=5 only.var_scaling21.0tuning factor of background error covariance for control variable 2 - unbalanced velocity potential. For cv_options=5 only.var_scaling31.0tuning factor of background error covariance for control variable 3 - unbalanced temperature. For cv_options=5 only.var_scaling41.0tuning factor of background error covariance for control variable 4 - pseudo relative humidity. For cv_options=5 only.var_scaling51.0tuning factor of background error covariance for control variable 5 - unbalanced
surface pressure. For cv_options=5 only.len_scaling11.0t!
uning fa
ctor of scale-length for stream function. For cv_options=5 only.len_scaling21.0tuning factor of scale-length for unbalanced velocity potential. For cv_options=5 only.len_scaling31.0tuning factor of scale-length for unbalanced temperature. For cv_options=5 only.len_scaling41.0tuning factor of scale-length for pseudo relative humidity. For cv_options=5 only.len_scaling51.0tuning factor of scale-length for unbalanced surface pressure. For cv_options=5 only.je_factor1.0ensemble covariance weighting factor&wrfvar8 ;not used&wrfvar9for program tracing. trace_use=.true. gives additional performance diagnostics (calling tree, local routine timings, overall routine timings, memory usage) It does not change results, but does add runtime overhead.stdout6unit number for standard outputstderr0unit number for error outputtrace_unit7Unit number for tracing output note that units 10 and 9 are reserved for reading namelist.input and writing namel
ist.output respectively.trace_pe0Currently, statistics are always calculated for all processors, and output by processor 0.trace_repeat_head10the number of times any trace statement will produce output for any particular routine. This stops overwhelming trace output when a routine is called multiple times. Once this limit is reached a 'going quiet' message is written to the trace file, and no more output is produced from the routine, though statistics are still gathered.trace_repeat_body10see trace_repeat_head descriptiontrace_max_depth30define the deepest level to which tracing writes outputtrace_usetrue.true.: activate tracingtrace_use_frequentfalsetrace_use_dullfalsetrace_memorytrue.true.: calculate allocated memory using a mallinfo call. On some platforms (Cray and Mac), mallinfo is not available and no memory monitoring can be done.trace_all_pesfalse.true.: tracing is output for all pes. As stated in trace_pe, this does not change pro
cessor statistics.trace_csvtrue.true.: tracing statistic!
s are wr
itten to a xxxx.csv file in CSV formatuse_htmltrue.true.: tracing and error reporting routines will include HTML tags.warnings_are_fatalfalse.true.: warning messages that would normally allow the program to continue are treated as fatal errors.&wrfvar10 ; for code developer&wrfvar11cv_options_hum1do not changecheck_rh20 --> No supersaturation check after minimization.
+5: NCAR Background Error model (default)as1(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 1 = stream function. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.5.as2(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 2 - unbalanced potential velocity. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.5.as3(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 3 - unbalanced temperature. For cv_options=3 only. The actual defau
lt values are 0.25, 1.0, 1.5.as4(3) -1.0tuning factors for variance, horizontal and vertical scales for control variable 4 - pseudo relative humidity. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.5.as5(3)-1.0tuning factors for variance, horizontal and vertical scales for control variable 5 - unbalanced surface pressure. For cv_options=3 only. The actual default values are 0.25, 1.0, 1.5.rf_passes6number of passes of recursive filter.var_scaling11.0tuning factor of background error covariance for control variable 1 - stream function. For cv_options=5 only.var_scaling21.0tuning factor of background error covariance for control variable 2 - unbalanced velocity potential. For cv_options=5 only.var_scaling31.0tuning factor of background error covariance for control variable 3 - unbalanced temperature. For cv_options=5 only.var_scaling41.0tuning factor of background error covariance for control variable 4 - pseudo relative humidit
y. For cv_options=5 only.var_scaling51.0tuning factor of!
backgro
und error covariance for control variable 5 - unbalanced surface pressure. For cv_options=5 only.len_scaling11.0tuning factor of scale-length for stream function. For cv_options=5 only.len_scaling21.0tuning factor of scale-length for unbalanced velocity potential. For cv_options=5 only.len_scaling31.0tuning factor of scale-length for unbalanced temperature. For cv_options=5 only.len_scaling41.0tuning factor of scale-length for pseudo relative humidity. For cv_options=5 only.len_scaling51.0tuning factor of scale-length for unbalanced surface pressure. For cv_options=5 only.je_factor1.0ensemble covariance weighting factor&wrfvar8 ;not used&wrfvar9for program tracing. trace_use=.true. gives additional performance diagnostics (calling tree, local routine timings, overall routine timings, memory usage) It does not change results, but does add runtime overhead.stdout6unit number for standard outputstderr0unit number for error output
trace_unit7Unit number for tracing output note that units 10 and 9 are reserved for reading namelist.input and writing namelist.output respectively.trace_pe0Currently, statistics are always calculated for all processors, and output by processor 0.trace_repeat_head10the number of times any trace statement will produce output for any particular routine. This stops overwhelming trace output when a routine is called multiple times. Once this limit is reached a 'going quiet' message is written to the trace file, and no more output is produced from the routine, though statistics are still gathered.trace_repeat_body10see trace_repeat_head descriptiontrace_max_depth30define the deepest level to which tracing writes outputtrace_usetrue.true.: activate tracingtrace_use_frequentfalsetrace_use_dullfalsetrace_memorytrue.true.: calculate allocated memory using a mallinfo call. On some platforms (Cray and Mac), mallinfo is not available and no memory monit
oring can be done.trace_all_pesfalse.true.: tracing is o!
utput fo
r all pes. As stated in trace_pe, this does not change processor statistics.trace_csvtrue.true.: tracing statistics are written to a xxxx.csv file in CSV formatuse_htmltrue.true.: tracing and error reporting routines will include HTML tags.warnings_are_fatalfalse.true.: warning messages that would normally allow the program to continue are treated as fatal errors.&wrfvar10 ; for code developer&wrfvar11cv_options_hum1do not changecheck_rh20 --> No supersaturation check after minimization.
1 --> supersaturation (rh> 100%) and minimum rh (rh<10%) check, and make the local adjustment of q.
2 --> supersaturation (rh> 95%) and minimum rh (rh<11%) check and make the multi-level q adjustment under the constraint of conserved column integrated water vaporsfc_assi_options11 --> 1surface observations will be assimilated based on the lowest model level first guess. Observations are not used when the height difference of t
he elevation of the observing
site and the lowest model level height is larger than 100m.
@@ -1102,10 +1093,10 @@
-WRF-ARW V3: Users Guide 6- PAGE 54
+WRF-ARW V3: Users Guide 6- PAGE 52
-WRF-ARW V3: Users Guide 6- PAGE 53
+WRF-ARW V3: Users Guide 6- PAGE 43
@@ -1126,11 +1117,11 @@
-
-
-&
+
+
+&
@@ -1140,1329 +1131,1332 @@
-h^xK
+h
-
-h^xK
-h^xK
+h
+h
-Px 4 #\'*.25@9
-h^xK
-h^xK
-h^xK
-+
-h^xK
-&
-&
+Px 4 #\'*.25@9
+h
+h
+h
++
+h
+&
+&
&
+hhd
&
+
& F
+hhd
&
+hhd
&
+
& F
h(
Hd
-`Hgd^xK
+`Hgd
& F
-h^xK
-
-h^xK
+hhd
+h
+
+h
&
-
-h^xK
{{
-h^xK
-h^xK
-h^xK
+h
{{
+h
+h
+h
l
-l
+l
l
-l
-h^xK
+l
+h
l
-l
+l
l
-l
-h^xK
-h^xK
-h^xK
-
-h^xK
-h^xK
+h
+h
+h
+
+h
+h
-&
-&
-&
-h^xK
-
+&
+&
+&
+h
+
t
-l
+l
t
-l
+l
t
l
+
t
l
+
t
-l
+l
t
l
+
t
-l
-
-
-h^xK
+l
+
+
+h
t
-l
+l
t
-l
+l
t
l
t
l
t
-l
+l
t
-l
+l
t
l
t
-l
+l
t
l
t
l
t
-l
+l
t
-l
-h^xK
+l
+h
-h^xK
+
+h
-
-h^xK
-h^xK
-h^xK
-h^xK
+h
+h
+h
+h
-)
-h^xK
-h^xK
-h^xK
-:
-h^xK
-h^xK
-h^xK
-&
+)
+h
+h
+h
+:
+h
+h
+h
+&
l
+
,#
l
,#
-l
-h^xK
-h^xK
+l
+h
+h
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
-l
+l
,#
l
l
+
,#
l
-f
+f
,#
l
+
,#
-l
-h^xK
-h^xK
+l
+h
+h
l
+
,#
l
-f
-f
+f
+f
,#
l
,#
l
-f
+
+f
,#
l
+
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
-Q
-h^xK
+Q
+h
,#
l
+
,#
l
+
,#
l
+
,#
l
Q
+
,#
l
Q
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
-l
+l
,#
l
+
,#
l
-h^xK
-h^xK
+Z
+h
+h
,#
-l
+l
,#
l
l
+
,#
l
!
,#
-l
+l
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
+
,#
l
l
+
,#
-l
-a
-h^xK
+l
+h
,#
l
+
,#
l
l
+
,#
-l
+l
,#
l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
-h^xK
-h^xK
+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
-g
-h^xK
-g
+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^xK
-h^xK
-f
+l
+h
+h
+f
,#
-l
+l
,#
-l
+l
,#
l
-f
+
+f
,#
l
+
l
-
+
l
+
,#
-l
+l
,#
l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
-u
-h^xK
-h^xK
+l
+v
+h
+h
,#
-l
+l
,#
-l
+l
l
+
,#
-l
+l
l
-
+
#
-l
+l
+v
#
-l
+l
#
-l
+l
#
-l
+l
#
l
+
#
l
+
#
-l
+l
#
-l
-h^xK
-h^xK
+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
-h^xK
-h^xK
+
+
+h
+h
,#
-l
+l
,#
-l
+l
,#
-l
+l
"
,#
l
+
l
+
"
,#
-l
+l
,#
-l
+l
,#
l
+
,#
-l
+l
l
,#
-l
+l
l
+
,#
-l
+l
,#
l
-
-h^xK
+
+h
l
-
+
l
+
,#
-l
+l
,#
l
-
-
+
l
-
+
l
+
,#
-l
+l
,#
l
+
l
-
+
,#
-l
+l
,#
-l
+
,#
-l
+l
,#
-l
-h^xK
-h^xK
+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
-h^xK
-h^xK
+
+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
-h^xK
+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^xK
+
+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
-h^xK
-h^xK
+l
+h
+h
#
l
+
#
l
+
#
-l
+l
#
l
+
#
l
+
#
-l
+l
#
l
#
l
+
#
-l
+l
#
l
+
#
l
+
#
-l
+l
#
l
+
#
-l
+l
@@ -18606,13 +18600,13 @@
+
+
+
+
+
@d
@t
K
-a
+T
+p
+
+
0
0
0
- 0
0
s
0
?
?
P
p
`
+ 0
0
s
0
?
?
P
p
`
%
K
-a
-
uot;
/
+T
+p
+
+
+
uot;
/
d
"
"
quot;
4d
@
@@ -18696,7 +18694,8 @@
@
-
+
+
K
K
K
@@ -18712,32 +18711,32 @@
z
-
+
-Q
-
+Q
+
+
+
+
-
+
+
@@ -18749,7 +18748,7 @@
.
M
+c
y
B
N
-^
-g
-
g?
+T
+Y
+p
y
}
@@ -18797,12 +18797,13 @@
h[
-g
<font color="gray">
-|
-
:
+f
+p
+
+
^Zp;K`|$~
@@ -18812,31 +18813,32 @@
-^Z
+^Z
`GT^5bBJ>*N
K
-a
-6
-
+p
+
+
+
-
+
+
+
+
+
-
\ No newline at end of file
+
+
\ No newline at end of file
</font>
</pre>