<p><b>xinzhang</b> 2009-04-01 14:58:42 -0600 (Wed, 01 Apr 2009)</p><p>Save for Chap6 changes.<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-01 17:57:26 UTC (rev 161)
+++ trunk/wrf/UsersGuide/users_guide_chap6.doc        2009-04-01 20:58:42 UTC (rev 162)
@@ -1,31 +1,5 @@
-ࡱ
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
Description of Namelist Variables
WRF-Var namelist variables.
-Variable NamesDefault ValueDescription&wrfvar1write_incrementsfalse.true.: write out the analysis increment file.wrfvar_mem_model0not usedvar4dfalse.false.: 3D-Var mode; .true.: 4D-Var modemulti_inc0> 0: multi-incremental runvar4d_coupling2not usedprint_detail_radarfalseprint_detail_xxx: output extra (sometimes can be too many) diagnostics for debugging; not recommended to turn them on for production runs.print_detail_xafalseprint_detail_xbfalseprint_detail_obsfalseprint_detail_gradfalsethe purpose of print_detail_grad is changed in V3.1
-.true.: to print out detailed gradient of each observation type at each iteration and write out detailed cost function and gradient into files called cost_fn and grad_fn.check_max_iv_printtrueobsolete (only used by Radar)&wrfvar2analysis_accu900Specify the accuracy of the time difference in seconds between the
analysis time and the first guess time. calc_w_incrementfalse.TRUE.: the increment of the vertical velocity W will be diagnosed based on the increments of other fields. If there is information of the W from observations assimilated, such as the Radar radial velocity, the W increments ar always computed, no matter CALC_W_INCREMENTS=.TRUE. or .FALSE.
-.FALSE.: the increment of the vertical velocity W is zero if no W information assimilated.
-dt_cloud_modelfalsewrite_qcwfalsewrite_qrnfalsewrite_qcifalsewrite_qsnfalsewrite_qgrfalsewrite_filtered_obsfalse&wrfvar3fg_format1Indicates the first guess format: 1 = WRF; ( 2 = MM5; 3 = KMA, not tested)ob_format2Indicates the observation file format: 1 = BURF; 2 = ASCIInum_fgat_time11 = 3DVar, >1 apply to FGAT and 4DVar&wrfvar4use_synopobstrue.TRUE.: assimilate the SYNOP data if any available.
-.FALSE.: not assimilate the SYNOP data even if available from OBS file.use_s
hipsobstrue.TRUE.: assimilate the SHIPS data if any avail!
able.
-.
FALSE.: not assimilate the SHIPS data even if available from OBS file.use_metarobstrue.TRUE.: assimilate the METAR data if any available.
-.FALSE.: not assimilate the METAR data even if available from OBS file.use_soundobstrue.TRUE.: assimilate the SOUND data if any available.
-.FALSE.: not assimilate the SOUND data even if available from OBS file.use_mtgirsobstrue.TRUE.: assimilate the MTGIRS data if any available.
-.FALSE.: not assimilate the MTGIRS data even if available from OBS file.use_pilotobstrue.TRUE.: assimilate the PILOT data if any available.
-.FALSE.: not assimilate the PILOT data even if available from OBS file.use_airepobstrue.TRUE.: assimilate the AIREP data if any available.
-.FALSE.: not assimilate the AIREP data even if available from OBS file.use_geoamvobstrue.TRUE.: assimilate the GEOAMV data if any available.
-.FALSE.: not assimilate the GEOAMV data even if available from OBS file.use_polaramvobstrue.TRUE.: assimilate t
he POLARAMV data if any available.
-.FALSE.: not assimilate the POLARAMV data even if available from OBS file.use_bogusobstrue.TRUE.: assimilate the BOGUS data if any available.
-.FALSE.: not assimilate the BOGUS data even if available from OBS file.use_buoyobstrue.TRUE.: assimilate the BUOY data if any available.
-.FALSE.: not assimilate the BUOY data even if available from OBS file.use_profilerobstrue.TRUE.: assimilate the PROFILER data if any available.
-.FALSE.: not assimilate the PROFILER data even if available from OBS file.use_satemobstrue.TRUE.: assimilate the SATEM data if any available.
-.FALSE.: not assimilate the SATEM data even if available from OBS file.
-** Currently only the retrieved SATEM thickness data can be assimilateduse_gpsztdobsfalse.TRUE.: assimilate the GPSZTD data if any available.
-.FALSE.: not assimilate the GPSZTD data even if available from OBS file.use_gpspwobstrue.TRUE.: assimilate the GPSPW data if any available
.
-.FALSE.: not assimilate the GPSPW data even if available !
from OBS
file.use_gpsrefobstrue.TRUE.: assimilate the GPSREF data if any available.
-.FALSE.: not assimilate the GPSREF data even if available from OBS file.use_ssmiretrievalobsfalse.TRUE.: assimilate the SSMI retrieval data (PW and surface wind speed) if any available.
-.FALSE.: not assimilate the SSMI retrieval data (PW and surface wind
-speed) even if available from OBS file.use_ssmitbobsfalse.TRUE.: assimilate the SSMI Brightness temperature data if any available.
-.FALSE.: not assimilate the SSMI Brightness temperature data even if available from OBS file.
-use_ssmt1obsfalse.TRUE.: assimilate the SSMT1 (temperature) data if any available.
-.FALSE.: not assimilate the SSMT1 (temperature) data even if available
-from OBS file.use_ssmt2obsfalse.TRUE.: assimilate the SSMT2 (relative humidity) data if any available.
-.FALSE.: not assimilate the SSMT2 (relative humidity) data even if available from OBS file.
-use_qscatobstrue.TRUE.: assimilate the Quik SC
AT (SeaWinds) data if any available.
-.FALSE.: not assimilate the Quik SCAT (SeaWinds) data even if available
-** Note at this moment, assimilation of only the SeaWinds Level 2B data from JPL.
-use_radarobsfalse.TRUE.: assimilate the RADAR data if any available.
-.FALSE.: not assimilate the RADAR data even if available from OBS file.
-use_radar_rvfalseuse_hirs2obsfasleuse_hirs3obsfalseuse_hirs4obsfalseuse_mhsobsfalseuse_msuobsfalseuse_amsuaobsfalseuse_amsubobsfalseuse_airsobsfalseuse_airsretobstrueuse_eos_amsuaobsfalseuse_eos_radobsfalseuse_hsbobsfalseuse_ssmisobsfalseuse_kma1dvarfalseuse_filtered_radfalseuse_obs_errfacfalse.TRUE.: the OBS error tuning factors need to be applied.
-.FALSE.: no the OBS error tuning factors need to be applied. At this moment, always set use_obs_errfac =.FALSE., not change.
-&wrfvar5
-check_max_ivtrue.TRUE.: Go through the Maximum OBS error check based on the
innovations (O-B). For safety, in general, it should set ch!
eck_max_
iv = .TRUE.
-.FALSE.: no Maximum OBS error check. Only when you sure that all of the observation data have the good quality, such as in the SINGLE-OBS tests, you can set check_max_iv = .FALSE.
-put_rand_seedfalseset put_rand_seed = .FALSE., not changeomb_set_randfalseset omb_set_rand = .FALSE., not change omb_add_noisefalseset omb_add_noise = .FALSE., not change.position_lev_dependantfalseobs_qc_pointer0Indicates the quality flag (qc) criterion, integer number, default is 0:
-the data with qc >= obs_qc_pointer are assimilated.
-the data with qc < obs_qc_pointer are discarded.max_sound_input5000max_mtgirs_input50000max_synop_input50000max_geoamv_input300000max_polaramv_input300000max_airep_input200000max_satem_input50000max_pilot_input30000max_radar_input10000max_metar_input30000max_gpspw_input4000max_ships_input20000max_profiler_input30000max_bogus_input5000max_buoy_input30000max_ssmi_rv_
input30000max_ssmi_tb_input30000max_ssmt1_input30000max_ssmt2_input30000max_qscat_input50000max_gpsref_input800000max_airsr_input100000max_tovs_input10000000max_ssmis_input10000000report_start1report_end10000000tovs_start1tovs_end10000000&wrfvar6max_ext_its1Used to assign the number of the outer loops. For example, 1 = a single outer loop, 3 = 3 outer loops, etc., now the maximum
-number of outer loops is allowed to be 10.ntmax200The maximum number of iterations of the minimization in internal loopnsave4Not usedwrite_interval5Not usedeps
-0.01 (max_ext_its)Used to assign the minimization convergence criterion for the internal loop
-within each of the outer loops.precondition_cgfalseprecondition_factor1.0&wrfvar7rf_passes6Used to assign the number of passes of recursive filter.var_scaling11.0The background error variance tuning factor for control variable 1 --- stream function PS
I.var_scaling21.0The background error variance tuning fa!
ctor for
control variable 2 --- potential velocity CHI.var_scaling31.0The background error variance tuning factor for control variable 3 --- unbalanced pressure Pu.var_scaling41.0The background error variance tuning factor for control variable 4 --- specific humidity q or relative humidity RH.var_scaling51.0Not used currentlylen_scaling11.0The scale-length (used in the recursive filter) tuning factor for control variable 1 --- Stream function PSI.len_scaling21.0The scale-length (used in the recursive filter) tuning factor for control variable 2 --- potential velocity CHI.len_scaling31.0The scale-length (used in the recursive filter) tuning factor for control variable 3 --- unbalanced pressure Pu.len_scaling41.0The scale-length (used in the recursive filter) tuning factor for control variable 4 --- specific humidity q or relative humidity RH.len_scaling51.0Not used currentlyjb_factor1.0je_factor1.0power_truncation1.0&wrfvar8def_sub_d
omainfalsex_start_sub_domain1.0y_start_sub_domain1.0x_end_sub_domain1.0e10y_end_sub_domain1.0e10&wrfvar9stdout6stderr0trace_unit7trace_pe0trace_repeat_head10trace_repeat_body10trace_max_depth30trace_usetruetrace_use_frequentfalsetrace_use_dullfalsetrace_memorytruetrace_all_pesfalsetrace_csvtrueuse_htmltruewarnings_are_fatalfalse&wrfvar10test_transformsfalseTRUE.: Do the VToX transform correctness check. Here V represents the control variables, X represents the analyzed variables. VToX means the transform formulas and their adjoint. When you change the V or X, you must do the VToX check first to guarantee the correctness. Program stopped after check finished.
-.FALSE. Not do adjoint check, but do the data assimilation.test_statisticsfalseinterpolate_statstrueNot usedtest_dm_exactfalse&wrfvar11cv_options_hum11 --> Use the specific humidity q as the contro
l variable.
-2 --> Use the relative humidity RH as the co!
ntrol va
riable.check_rh20 --> No supersaturation check after minimization.
-1 --> With the supersaturation (rh>100%) and minimum rh (rh<10%) check, and make the local adjustment of thespecific humidity.
-2 --> With the maximum (rh>95%) and minimum (rh<11%) rh check, and make the multi-level q adjustment under the constraints of integrated water vapor in column conserved.set_omb_rand_fac1.0Scaling factor for random O-B (*sigma_o). Not change.seed_array10The first element of the seed for subroutine random_seed. Not change.seed_array2 0The second element of the seed for subroutine random_seed. Not change.sfc_assi_options11 --> Before the minimization, the surface observed parameters (wind, temperature, moisture and pressure) are reduced to the lowest model level. This means that the original observations are modified based on the information of the first guess fields. The observations are discarded when the elevations of the observing sites ar
e too far from the lowest model level height (> 100m).
-2 --> Based on the surface similarity theory in PBL, some of the observation operators (Hx) are developed to get the analysis equivalent of the surface observed parameters, then the innovations are computed.
-The original values of the surface observation are never changed. Regardless of the elevations of the observing sites, the analysis 10-m wind and 2-m temperature, moisture are computed. Only the analysis pressure computation need to consider the station heights.calculate_cg_cost_fnfalse.TRUE.: Compute the cost functions every iteration and output.
-.FALSE.: Only the initial and final cost functions are computed and output.lat_stats_optionfalseTRUE. The I-dependent regression coefficients for balance part of pressure are used. If the input BES file does not haveI-dependent regression coefficients, program stopped.
-.FALSE.: The domain averaged regression coefficients for balance part of pressure ar
e used.&wrfvar12balance_type1Momentum-mass balance!
type:
-
1 = Geostrophic;
-2 = Cyclostrophic;
-3 = Geo/cyclostrophic.&wrfvar13vert_corr2Currently you must set vert_corr = 2, i.e. the eigenvector decomposition
-method is applied in vertical in modeling the vertical correlationship.vertical_ip0The type of the vertical inner product:
-0=Euclidean,
-1=SQRT(pressure),
-2=pressurevert_evalue11 --> use the global eigenvector/eigenvalue. Now program hardwired vert_evalue = 1, other values will stop the job. Not change it.
-2 --> use the I-dependent eigenvector/eigenvalue.
-Note: Although the I-dependent (approximately latitude-dependent) eigenvector and eigenvalue are available in BE file, the GLOBAL eigenvector/eigenvalue gave smoother analysis results.max_vert_var199.0Specify the maximum truncation value (percentage) to explain the variance of the control variable 1 -- Stream function, in eigenvector decompositionmax_vert_var299.0Specify the maximum truncation value (percentage) to explain the variance o
f the control variable 2 -- Potential velocity, in eigenvector decompositionmax_vert_var399.0Specify the maximum truncation value (percentage) to explain the variance of the control variable 3 -- Unbalanced part of pressure, in eigenvector decompositionmax_vert_var499.0Specify the maximum truncation value (percentage) to explain the variance of the control variable 4 -- Specify humidity or relative humidity, in eigenvector decompositionmax_vert_var599.0Not used&wrfvar14rtminit_print1rtminit_nsensor1rtminit_platform-1 (max_instruments)rtminit_satid-1.0 (max_instruments)rtminit_sensor-1.0 (max_instruments)rad_monitoring0 (max_instruments)thinning_mesh60.0 (max_instruments)thinning_meshfalseread_biascoeftruebiascorrfalsebiasprepfalsewrite_profilefalsewrite_jacobianfalseqc_radtruewrite_iv_rad_asciifalsewrite_oa_rad_asciifalsewrite_filtered_radfalseuse_error_factor_radfalseuse_landemfals
emw_emis_sea1tovs_min_transfer10tovs_batchfalse!
rtm_op
tion1use_crtm_kmatrixfalseuse_crtm_kmatrix_fastfalsecrtm_cloudfalseonly_sea_radfalseuse_simulated_radfalseuse_pseudo_radfalsepseudo_rad_platid1pseudo_rad_satid16pseudo_rad_senid3pseudo_rad_ichan6pseudo_rad_lat35.0pseudo_rad_lon130.0pseudo_rad_tb258.0pseudo_rad_err1.0spccoeff_fileCRTM_SpcCoeff.bintaucoeff_fileCRTM_TauCoeff.binaerosolcoeff_fileCRTM_AerosolCoeff.bincloudcoeff_fileCRTM_CldCoeff.binemiscoeff_fileCRTM_EmisCoeff.binuse_varbcfalsefreeze_varbcfalsevarbc_factor1.0use_airs_mmrfalseairs_warmest_fovfalsecrtm_atmosphere0&wrfvar15num_pseudo0Set the number of pseudo observations. Currently only the value of 1 is allowed, i.e. num_pseudo = 1 --- the SINGLE-OBS testpseudo_x1.0Set the x-position (I) of the OBS in unit of grid-point.pseudo_y1.0Set the y-position (J) of the OBS in unit of grid-point.pseudo_z1.0Set the z-position (ZK) of OBS with the vertical
level index.pseudo_val1.0Set the innovation of the OBS variable: u/v at m/s, p at Pa, t at degree, q at kg/kg.pseudo_err1.0Set the OBS error of the OBS variable.&wrfvar16alphacv_method1ensdim_alpha0alpha_truncation1alpha_corr_type3alpha_corr_scale3alpha_std_dev16 &wrfvar17analysis_type3D-VARIndicates the job type, now only '3D-VAR', not change.&wrfvar18analysis_date2002-08-03_00:00:00.0000Specify the analysis time.&wrfvar19pseudo_vartSet the name of the OBS variable:
+Variable NamesDefault ValueDescription&wrfvar1write_incrementsfalse.true.: write out a binary analysis increment file.var4dfalse.true.: 4D-Var modemulti_inc0> 0: multi-incremental runvar4d_coupling21: var4d_coupling_disk_linear, 2: var4d_coupling_disk_simulprint_detail_radarfalseprint_detail_xxx: output extra (sometimes can be too many) diagnostics for debugging; not recommended to turn them on for product
ion runs.print_detail_xafalseprint_detail_xbfalsep!
rint_det
ail_obsfalseprint_detail_gradfalsethe purpose of print_detail_grad is changed in V3.1
+.true.: to print out detailed gradient of each observation type at each iteration and write out detailed cost function and gradient into files called cost_fn and grad_fn.check_max_iv_printtrueobsolete (only used by Radar)&wrfvar2analysis_accu900seconds, if time difference between namelist setting
+ (analysis_date) and date info read in from first guess is larger than analysis_accu, WRF-Var will issue a warning message ("=======> Wrong xb time found???"), but won't abort.
+calc_w_incrementfalsedt_cloud_modelfalse&wrfvar3fg_format1 1: fg_format_wrf_arw_regional (default)
+ 2: fg_format_wrf_nmm_regional
+ 3: fg_format_wrf_arw_global
+ 4: fg_format_kma_global
+ob_format21: ob_format_bufr (NCEP PREPBUFR), read in data from ob.bufr (not fully tested)
+2: ob_format_ascii (output from obsproc), read in data from ob.ascii (default)
+3: ob_format
_madis (not tested)
+num_fgat_time11: 3DVar
+> 1: number of time slots for FGAT and 4DVAR (for ob_format=2 and radiance only)&wrfvar4thin_convtruefor ob_format=1 (NCEP PREPBUFR) only. thinning is mandatory for ob_format=1 as time-duplicate data are "thinned" within thinning routine, however, thin_conv can be set to .false. for debugging purpose.thin_mesh_conv 20. (max_instruments)for ob_format=1 (NCEP PREPBUFR) only.
+km, each observation type can set its thinning mesh and the index/order follows the definition in
+WRFDA/var/da/da_control/da_control.f90use_synopobstrueuse_xxxobs - .true.: assimilate xxx obs if available
+.false.: not assimilate xxx obs even available
+.
+.use_shipsobstrueuse_metarobstrueuse_soundobstrueuse_pilotobstrueuse_airepobstrueuse_geoamvobstrueuse_polaramvobstrueuse_bogusobstrueuse_buoyobstrueuse_profilerobstrueuse_satemobstrueuse_gpspwobstrueuse_gpsrefobstrueuse_
qscatobstrueuse_radarobsfalseuse_radar_rvfalseus!
e_radar_
rffalseuse_airsretobstrue ; use_hirs2obs, use_hirs3obs, use_hirs4obs, use_mhsobs
+ ; use_msuobs, use_amsuaobs, use_amsubobs, use_airsobs,
+ ; use_eos_amsuaobs, use_hsbobs, use_ssmisobs are
+ ; radiance-related variables that only control if reading
+ ; in corresponding BUFR files into WRF-Var or not, but
+ ; do not control if assimilate the data or not.
+ ; Some more variables have to be set in &wrfvar14 in order
+ ; 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.t
rue.: 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 observation. 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 f
actor for specific humiditymax_error_p5.0maximum check_m!
ax_iv er
ror 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_ext_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.as2(3)1.0tuning factors for variance, horizontal and vertical scales for cont
rol variable 2 - unbalanced potential velocity. For cv_options=3 only.as3(3)1.0tuning factors for variance, horizontal and vertical scales for control variable 3 - unbalanced temperature. For cv_options=3 only.as4(3) 1.0tuning factors for variance, horizontal and vertical scales for control variable 4 - pseudo relative humidity. For cv_options=3 only.as5(3)1.0tuning factors for variance, horizontal and vertical scales for control variable 5 - unbalanced surface pressure. For cv_options=3 only.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.0tun
ing factor of background error covariance for control varia!
ble 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.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 o
verhead.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 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 mal
linfo call. On some platforms (Cray and Mac), mallinfo is no!
t availa
ble 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 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 the elevation of the observing
+site and the lowest model level height is larger than 100m.
+2 --> surface observations will be assimilated based on surface similarity theory in PBL. Innovations are computed based on 10-m wind, 2-m temperature and 2-m moisture.calculate_cg_cost_fnfalsethe purpose of calculate_cg_cost_fn is changed.
+use print_detail_grad=.true. to dump cost function and gradient of each iteration to cost_fn and grad_fn. conjugate gradient algorithm does not require the computation of cost function at every iteration during minimization..true.: cost function is printed out and is directly derived from the gradient using the fully linear properties inside the inner-loop..false.: Only the initial and final cost functions are computedlat_stats_optionfalsedo not change&wrfvar12balance_type1obsolete&wrfvar13vert_corr2do not changevertical_ip
0obsoletevert_evalue1do not changemax_vert_var199.0!
specify
the maximum truncation value (in percentage) to explain the variance of stream function in eigenvector decompositionmax_vert_var299.0specify the maximum truncation value (in percentage) to explain the variance of unbalanced potential velocity in eigenvector decompositionmax_vert_var399.0specify the maximum truncation value (in percentage) to explain the variance of the unbalanced temperature in eigenvector decompositionmax_vert_var499.0specify the maximum truncation value (percentage) to explain the variance of pseudo relative humidity in eigenvector decompositionmax_vert_var599.0for unbalanced surface pressure no effect except for max_vert_var5=0.0&wrfvar14the following 4 variables (rtminit_nsensor, rtminit_platform, rtminit_satid, rtminit_sensor) together control what sensors to be assimilated.rtminit_nsensor1total number of sensors to be assimilatedrtminit_platform-1
+(max_instruments)platforms IDs array (used dimension: rtminit_nsensor);
e.g., 1 for NOAA, 9 for EOS, 10 for METOP and 2 for DMSPrtminit_satid-1.0
+(max_instruments)satellite IDs array (used dimension: rtminit_nsensor)rtminit_sensor-1.0
+(max_instruments)sensor IDs array (used dimension: rtminit_nsensor); e.g., 0 for HIRS, 3 for AMSU-A, 4 for AMSU-B, 15 for MHS, 10 for SSMIS, 11 for AIRSrad_monitoring0
+(max_instruments)integer array (used dimension: rtminit_nsensor); 0: assimilating mode
+; 1: monitoring mode (only calculate innovations)thinning_mesh60.0
+(max_instruments)real array (used dimension: rtminit_nsensor); specify thinning mesh size (in KM) for different sensors.thinningfalse.true.: perform thinning on radiance dataqc_radtrue.true.: perform quality control. always .true.write_iv_rad_asciifalse.true.: output radiance Observation minus Background files, which are in ASCII format and separated by sensors and processors.write_oa_rad_asciifalse.true.: output radiance Observation minus Background files,
which are in ASCII format and separated by sensors and proc!
essors.
use_error_factor_radfalse.true.: use a radiance error tuning factor file "radiance_error.factor", which can be created with empirical values or generated using variational tuning method (Desroziers and Ivanov, 2001)use_antcorrfalse
+(max_instruments).true.: perform Antenna Correction in CRTMrtm_option1what RTM (Radiative Transfer Model) to use 1: RTTOV (WRF-Var needs to compile with RTTOV) 2: CRTM (WRF-Var needs to compile with CRTM)only_sea_radfalse.true.: assimilate radiance over water onlyuse_varbcfalse.true.: perform Variational Bias Correction. A parameter file in ASCII format called VARBC.in (a template is provided with the source code tar ball)
+is required.freeze_varbcfalse.true: together with use_varbc=.false., keep the VarBC bias parameters constant in time. In this case, the bias correction is read and applied to the
+innovations, but it is not updated during the minimization.varbc_factor1.0for scaling the VarBC preconditionin
gvarbc_nobsmin10defines the minimum number of observations required for the computation of the predictor statistics during the first assimilation cycle. If there are not enough data (according to "VARBC_NOBSMIN") on the first cycle, the next cycle will perform a coldstart again.airs_warmest_fovfalse.true.: uses the observation brightness temperature forAIRS Window channel #914 as criterion for GSI thinning (with a higher amplitude than the distance
+from the observation location to the nearest grid point).crtm_atmosphere0climatology reference profile used above model top for CRTM Radiative Transfer Model (up to 0.01hPa
+ 0: Invalid (default, use U.S. Standard Atmosphere)
+1: Tropical
+2: Midlatitude summer
+3: Midlatitude winter
+4: Subarctic summer
+5: Subarctic winter
+6: U.S. Standard Atmosphere&wrfvar15 (needs to be set together with &wrfvar19)num_pseudo0Set the number of pseudo observations, either 0 or 1 (single ob)pseudo_x1.0Se
t the x-position (I) of the OBS in unit of grid-point.pseu!
do_y1.0
Set the y-position (J) of the OBS in unit of grid-point.pseudo_z1.0Set the z-position (K) of OBS with the vertical level index, in bottom-up order.pseudo_val1.0Set the innovation of the ob; wind in m/s, pressure in Pa, temperature in K,; specific humidity in kg/kg
+pseudo_err1.0set the error of the pseudo ob. Unit the same as pseudo_val.; if pseudo_var="q", pseudo_err=0.001 is more reasonable.&wrfvar16 (for hybrid WRF-Var/ensemble)alphacv_method21: ensemble perturbations in control variable space
+2: ensemble perturbations in model variable spaceensdim_alpha0ensemble sizealpha_corr_type31: alpha_corr_type_exp
+2: alpha_corr_type_soar
+3: alpha_corr_type_gaussian (default)alpha_corr_scale1500.0km&wrfvar17analysis_type3D-VAR"3D-VAR": 3D-VAR mode (default);
+ "QC_OBS": 3D-VAR mode plus extra filtered_obs output;
+"VERIFY": verification mode. WRF-Var resets check_max_iv=.false. and ntmax=0;
"RANDOMCV": for creating ensemble perturbations&wrfvar18 (needs to set &wrfvar21 and &wrfvar22 as well if ob_format=1 and/or radiances are used)analysis_date2002-08-03_00:00:00.0000specify the analysis time. It should be consistent with the first guess time. However, if time difference between analysis_date and date info read in from first guess is larger than analysis_accu, WRF-Var will issue a warning message ("=======> Wrong xb time found???"), but won't abort.&wrfvar19 (needs to be set together with &wrfvar15)pseudo_vartSet the name of the OBS variable:
'u' = X-direction component of wind,
'v' = Y-direction component of wind,
't' = Temperature,
'p' = Prerssure,
-'q' = Specific humidity&wrfvar20documentation_urlhttp://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk&wrfvar21time_window_min"2002-08-02_21:00:00.0000"&wrfvar22time_window_max"2002-08-03_03:00:00.0000&
quot;&wrfvar23jcdfi_usefalse.True.: Include JcDF !
term in
cost function.
-.False.: Ignore JcDF term in cost function.jcdfi_iofalse.True.: Read JcDF output from WRF+. Even jcdfi_use= false. Used for diagnosis.
-.False.: Ignore the JcDF output from WRF+jcdfi_tauc10800Filter time window second.jcdfi_gama1.0Scaling number used to tune the weight of JcDF termjcdfi_error_wind3.0Error for Wind (m/s) used in JcDFjcdfi_error_t1.0Error for Temperature (k) used in JcDFjcdfi_error_q0.001Error for Specific Humidity (kg/kg) used in JcDFjcdfi_error_mu1000.Error for MU (Pa) used in JcDF
+'q' = Specific humidity
+"pw": total precipitable water
+"ref": refractivity
+"ztd": zenith total delay&wrfvar20documentation_urlhttp://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk&wrfvar21time_window_min"2002-08-02_21:00:00.0000"start time of assimilation time window used for ob_format=1 and radiances to select observations inside the defined time_window. Note: new in
V3.1, this variable is also used for ob_format=2 to double-check if the obs are within the specified time window.
+&wrfvar22time_window_max"2002-08-03_03:00:00.0000"end time of assimilation time window used for ob_format=1 and radiances to select observations inside the defined time_window. Note: new in V3.1, this variable is also used for ob_format=2 to double-check if the obs are within the specified time window.&wrfvar23 (settings related to the 4D-Var penalty term option, which controls the high-frequency gravity waves using a digital filter)jcdfi_usefalse.true.: Include JcDF term in cost function.
+.False.: Ignore JcDF term in cost function.jcdfi_iofalse.true.: Read JcDF output from WRF+. Even jcdfi_use= false. Used for diagnosis.
+.False.: Ignore the JcDF output from WRF+jcdfi_tauc10800seconds, filter time window second.jcdfi_gama1.0Scaling number used to tune the weight of JcDF termjcdfi_error_wind3.0m/s, wind error used
in JcDFjcdfi_error_t1.0K, temperature error used in JcDF!
jcdfi_
error_q0.001kg/kg, specific humidity error used in JcDFjcdfi_error_mu1000.Pa, perturbation pressure (mu) error used in JcDF
OBSPROC namelist variables.
Variable NamesDescription&record1obs_gts_filenamename and path of decoded observation filefg_format'MM5' for MM5 application, 'WRF' for WRF applicationobserr.txtname and path of observational error filefirst_guess_filename and path of the first guess file&record2time_window_minBeginning of time window (included) as ccyy-mm-dd_hh:mn:sstime_analysisAnanlysis time as ccyy-mm-dd_hh:mn:sstime_window_maxEnd of time window (included) as ccyy-mm-dd_hh:mn:ss
** Note : Only observations between [time_window_min, time_window_max] will kept.&record3max_number_of_obsMaximum number of observations to be loaded, ie in domain and time window, this is independent of the number of obs actually read.fatal_if_exceed_max_obs.TRUE.: will stop when more than max_number_of_obs are loaded
@@ -
1094,10 +1066,10 @@
-WRF-ARW V3: Users Guide        6- PAGE 38
+WRF-ARW V3: Users Guide        6- PAGE 50
-WRF-ARW V3: Users Guide        6- PAGE 1
+WRF-ARW V3: Users Guide        6- PAGE 51
@@ -1118,11 +1090,11 @@
-
-
-&
+
+
+&
@@ -1134,1900 +1106,1489 @@
-hk
-
-&
-&
+h
+&
+
-Px4 #\'*.25@9
-&F
-&
-&
+Px4 #\'*.25@9
+&F
+&
&F
-&
-&
-&F
-&F
-hk
-&
-&F
-&
-&F
-&F
-hk
`gdk
-hk
+
+&
+&F
+h
+h
+&
+&F
+&
d
+h
+e
l
-l
+l
l
l
+i
l
-l
+l
l
-l
J
-hk
-&
+l
J
+h
+&
t
-l
+l
t
l
+
t
-l
+l
t
l
+
t
-l
+l
t
l
+
t
l
$hz
+
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
t
-l
+l
+
+
"
-hk
-&
-hk
-hk
+
ot;
+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
+
+f
,#
-l
+l
+h
+
+l
+
,#
-l
+l
+f
+f
,#
-l
-hk
+l
,#
-l
+l
+f
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-l
-
-
-,#
-
-l
-
-,#
-
l
+
,#
-l
+l
l
-
+
,#
-l
-
-,#
-
l
-
+
+
,#
-
-l
-
+
+l
+\
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-6
-8
-hk
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+h
+
,#
-
-l
-
+
+l
+
,#
-
-l
-8
-
+
+l
+
,#
-
-l
-
+
+l
+3
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
-
+
+l
+
,#
-
-l
+
+l
,#
-l
+l
+
+l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
-hk
+l
,#
-l
+l
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
+l
,#
-
-l
+
+l
,#
-l
+l
+h
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
+l
l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
-hk
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
+
+l
+
,#
-l
+l
,#
-l
+l
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-E
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-hk
+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
-hk
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
+l
l
+
,#
-l
+l
,#
l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
-N
+l
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
+l
+h
l
+
,#
-l
-hk
+l
,#
l
-
-,#
-
-l
-
-,#
-
-l
+
l
-
+
,#
-l
-
-,#
-
l
+
,#
-l
+l
+f
+f
,#
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
-
-,#
-
-l
+l
l
+
,#
-l
+l
+
+l
+
,#
l
-
-l
-
+
,#
-l
-Z
+l
,#
-l
-Z
-_
-hk
+l
,#
-l
+l
,#
-l
+l
+U
+W
+h
,#
-l
+l
,#
-l
+l
,#
-l
+l
l
+
,#
-l
+l
-d#
-
-
+#
+
l
+
+
+l
+
+W
-d#
-
-
-l
+#
+
+l
-d#
-
-
-l
-_
+#
+
+l
-d#
-
-
-l
-_
+#
+
+l
-d#
-
-
-l
+#
+
+l
-d#
-
-
-l
+#
+
+l
-d#
-
-
-l
+#
+
+l
-d#
-
-
-l
+#
+
+l
-d#
-
-
-l
-hk
+#
+
+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
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-hk
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
-
-d#
-
-
-l
+#
+
+l
l
+
-d#
-
-
-l
+#
+
+l
,#
l
+
,#
-l
+l
,#
-l
+l
,#
-l
+l
,#
-l
+l
l
+
,#
-l
-f
-hk
+l
+h
,#
l
+
,#
-l
+l
,#
-l
-
-,#
-
-l
-
-,#
-
-l
+l
l
+
,#
-l
-
-l
-
+l
,#
l
+
l
-
+
+
+l
+
,#
-l
+l
,#
l
-g
+
,#
l
+
l
-
-
+
+
l
+
l
+
,#
-l
+l
l
-
+
,#
l
-hk
-hk
+
+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
+h
#
l
+
#
-l
+l
l
-
+
#
l
-hk
+
#
l
+
#
-l
+l
#
-l
+l
l
+
#
l
+
#
-l
+l
#
-l
+l
#
l
-p
+
#
l
+
#
l
+
#
-l
+l
#
-l
+l
l
+
#
l
+
#
-l
+l
#
-l
+l
#
l
+
#
l
+
#
l
+
#
-l
+l
#
-l
+l
#
l
+
+h
l
-
+
#
l
-v
-hk
+
#
l
+
#
-l
+l
#
-l
+l
#
l
-v
+
#
l
+
#
l
+
#
-l
+l
#
-l
+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
-}
-hk
+l
#
-l
+l
#
l
+
#
l
-}
+
+
#
l
+
#
-l
+l
#
-l
+l
#
l
+
+
#
l
+
+
#
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
-hk
-hk
-
+
+
#
l
-
+
+
#
-l
+l
+
#
-l
+l
+
#
l
+
+
#
l
+
+
-
+
+#
+#?
Zcd
@@ -3522,10 +3083,11 @@
aoq[;kw}SO>ؖ*''f:S~0`]^yUZ73nEc~Wн>
'=)|L9/f}%io}]>h|O)~.K闯_>w~տ~:[w;yoH        s|]H90;X'5/4HOvwH2쬬D
-
-
-
+n/I
+epConف(,@E@|@ZU!439;$37uBb^Bb%c5
+)EE)ELb`yĢTԔ|TbҢT=NL
+
+
aߗ!,3@E
hA[^Ql^*Q,Z@*(KAKā0@        3d 94ɽI$|2mړON}r\2bKթ+-b8WZ,o˯[ZܮMGK8KX3X?秶,"cϚ>bv
UW]U|̇
@@ -3668,10 +3230,10 @@
>.$vܜR@̓';;$qU~PvVVW
pg^%
-31
-
-
-
+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Ͱz0eް^kXcusV=VxG(߾{X        ~a{a
;2l]ذ0o~a&
@@ -4230,9 +3792,19 @@
q%)]?vaWhu6rpi^;mۺpӞm-/*+,.+,F>B.,I?{pr9
uq};z̾c۸~|oL
F
-Yn~+GV`H`
-
-
+Yn~+GV`H`
+)EE)zb`]E
+y))
+%
+ũ
+%
+%Ez@63r/
+
+
+
+
)wsCO>c8㚚j-ղQ9{>3f}
&lM
@@ -4446,40 +4018,44 @@
l
l
-
-t
-t
-t
-t
-t
-t
-t
+
+t
+t
+t
+t
+t
+t
+t
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
+6`6
t
-
-
+
+
+
+
+{{
+{{
qxX_{N?f1_q3
c?7
gݾ(C+0sb!}Q1
-E<f߿s<{B;wjhw9ҥK9hK/˗/GkQ|MKV㵠;Wx5|Chaq,1D/믑=v<)RP:tH`,1ܹ3]#iY0&ؘ_!1<mǥ>$B!B!B!$qpps?_R###C4O320{ͺ8_>up8~|;/"|v:zOxx7 <<nAᥝG8Cq~z%qջ&ײU|);鴤XNh9G>ϥPWy%]GHIߙ;_onSs]2YZοw⳿8.ۧs>Gο+v7չN/pݶw(ک_Ia[!\pTܿ?        "W,οE])ty1V#kdohp(p}HZ9}vA^!!?8v@v]EYcGgi,sק+GοUw$yq"u}{wWvһ ɜ%MxMftz}nР߬#Q4}8ȓQvWѽkj8Tdʖ~Y7!H
-
-
+E<f߿s<{B;wjhw9ҥK9hK/˗/GkQ|MKV㵠;Wx5|Chaq,1D/믑=v<)RP:tH`,1ܹ3]#iY0&ؘ_!1<mǥ>$B!B!B!$qpps?_R###C4O320{ͺ8_>up8~|;/"|v:zOxx7 <<nAᥝG8Cq~z%qջ&ײU|);鴤XNh9G>ϥPWy%]GHIߙ;_onSs]2YZοw⳿8.ۧs>Gο+v7չN/pݶw(ک_Ia[!\pTܿ?        "W,οE])ty1V#kdohp(p}HZ9}vA^!!?8v@v]EYcGgi,sק+GοUw$yq"u}{wWvһ ɜ%MxMftz}nР߬#Q4}8ȓQvWѽkj8Tdʖ~Y7!I
+
+
+
+
+)EE){0eRRSSRSJSSJ22
+2KJRN20lRg^
\)QTD"tOHG/ET^"T]WyO}f7l3gflyg9sε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
@@ -4698,9 +4278,13 @@
*wȜ9Y+\r{.Ç;!,uDs1S̾3gΘHU"oxb5ҥKsu;{iᇾc$VGX9sfxoҤo<3xB6l@ddYɞ={ӧO7WDժU3        wǺŋm/1w\E>}\[>ϟ={o_r~0]vMX\Ovq^9揙żI>ID뒓9J>H"3XWs5'NUЪUTe+b~#kך
6<gMܻ; 4h`ˌa{I        7zY,`}}bxm]n]Sd/+R&>fr|l+:vhY0nܸύ|_vڵkcڴi%p,3>}{mۖfK*NjSre?_~ז:"Yw}gƶJI؎9bK.(Z~dsi`&f֬YXpqRYkmKL}YQog_ƾ,piukseX%"NvC:ׯoۗTԫW1W^f
<>M6PɯӧO}_wuRJq߃'%sޒ%K{n/66s=gْǴ#͘lwfN{\=}<xiߧ:X:}l>oʔ)祤}\7N:        ]CE\p6֑>}߲eo\/[/lKl,Ygd(Qo_\f=ތh͚5~mqYccMl4-1W͌qu2)Co?sGc3&|u<'g˗/cǎwneI2_5o6meAA҅][W"]NMl$EPœGAOG=xAkO        B)Mư ]wx'ɼod FCQţŖ,fdYn(C}ubEBv"8<nb*<O@*f|~4z:A~
-&Zm=fkTЊMɳ8m6-ШW_uIUQr7廫%r< 7@|',+.K⛮܉Q.}QJw/ߓ_~մWځVy"MiWI
-
-
+
+
+
+
+y)
+v2\ˠ[PZ\_r XF!(U!/55%5E$_85U$#X 3(U:DM
?uXMLLDFtO>EҥƘ?^oٳnkطoҏ# #Z-n߾m'OF``sPK}X`cbԨQ4+WD߾}z
@@ -5172,9 +4756,13 @@
@        jߒznMŔρ}q$:_uS-/Y$Ěm6UySiX<ĄQ߄L!I0^CR$=VSK
ZoTZnc\<qR`YCZn, %TG>VF@OA
-Om=ttzfZFv=+Mӝݺgieͷ,-(}-Ag+<ϙEd^Kd;}vV4Rn|>GLw-ޒlD&r;˱=^܀d,`+tgY~/F>;rcZoko,?bmW3hvU&j~z-B嗷DJm_R_DuDsay=kM|3J}D%-'vh&z8AT>gu-,)ҢJlÈuٲcٺqOkrΞΊVNv+y=+zW~PT2eLYdNu+7XSOtwOv88{k5+[qr_)ЎroN8Фjn>h)ښm+C
-
-
+Om=ttzfZFv=+Mӝݺgieͷ,-(}-Ag+<ϙEd^Kd;}vV4Rn|>GLw-ޒlD&r;˱=^܀d,`+tgY~/F>;rcZoko,?bmW3hvU&j~z-B嗷DJm_R_DuDsay=kM|3J}D%-'vh&z8AT>gu-,)ҢJlÈuٲcٺqOkrΞΊVNv+y=+zW~PT2eLYdNu+7XSOtwOv88{k5+[qr_)ЎroN8Фjn>h)ښm+E
+
+
+
+
+y)
+1\ˠ[PZ\_r4XF!(U!/55%5E$_85U$#X 3(UlM
J&H(yA@(( >φ4@(zoRAJhUx"|oܒ{M9Aٹ͙9sf !3'J̓
ƾEw{ly&cQL0窲
6֭GAW42D'FhC8YJ3Kͼƪc']^AA_FaaϱRrQ0VF6ꉆQiaěnF֭-PQ6.cLzfwO?OdoQT1l[TQc
Jߙ7gdct#H6\FoY!1N)#*zUS
-GN]fݝ1bdy7qڬUԝR?ux
-
-
+GN]fݝ1bdy7qڬUԝR?uy
+
+
+
+
+0f&gdTHKQH|Ƹ
+EV^܂ORRSSRSJSSJ22
+2KJR13lRg^
,(K,1(DDD9xL0Avk֬u+WdDDD9rX,={E!""2̬Ytz͚5>|8KFDDDݻwɓٞ "u߾}%<<\bbbdDDDDDDD)5|DEEɭ[dȐ!{|7ސ+WJ:u$oC[jo^oߞ!N:%ӦM#F3QFɤIuqٳgܾ}EGG{vڑ#GkDDDd,ժU5kBѢEH"Æ
sپ[nz}̙rRL*˖-slDD4H2d*r_~l۶˶ҥ7~x
<^t~
@@ -7168,27 +6761,24 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
2
55#v#v#v
:V
2
55#v#v#v
:V
-2
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
+:V
2
55#v#v#v
:V
@@ -7205,24 +6795,6 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -7240,64 +6812,65 @@
2
55#v#v#v
:V
-2
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
55#v#v#v
-:V
-2
+:V
+2
+
55#v#v#v
:V
2
@@ -7321,36 +6894,18 @@
2
55#v#v#v
:V
-2
+2
55#v#v#v
-:V
+:V
2
55#v#v#v
:V
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -7382,7 +6937,8 @@
2
55#v#v#v
:V
-2
+2
+
55#v#v#v
:V
2
@@ -7391,69 +6947,6 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -7476,10 +6969,6 @@
2
55#v#v#v
:V
-2
-
-55#v#v#v
-:V
2
55#v#v#v
:V
@@ -7513,12 +7002,6 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -7568,22 +7051,8 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
+
55#v#v#v
:V
@@ -7596,31 +7065,9 @@
2
55#v#v#v
:V
-2
-
-55#v#v#v
-:V
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -7651,294 +7098,67 @@
55#v#v#v
:V
2
+
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
-2
-5.
-5U#vU#v.
-#v
-:V
-5
+5^5%#v%#v^#v
+:V
2
55#v#v#v
@@ -7971,12 +7191,6 @@
2
55#v#v#v
:V
-2
-55#v#v#v
-:V
-2
-55#v#v#v
-:V
2
55#v#v#v
@@ -8308,9 +7522,9 @@
-
-&F^
-l
+&F
+l
Px4 #\'*.25@9
-&
-&
-
+&
+
dh
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-@
-@
-Dư
-@
-
+@d
+@t
-
-2
-?
-H
-T
-Y
-j
-p
-w
-
-
|
|
0
U#
U#
+-
+Q
+|
+}
+
K
0
0
0
U#
-
-2
-?
-H
-T
-Y
-j
-p
-w
-
-
-
-
-
-@
-
-
-
-
-@
-
-
-
-
-
-@
-@
-@
-
-
-
0
0
-_
-
-
+
+-
+Q
+|
+}
+
$d
$d
$d
$d
$d
$d
x@t
+
+
+@
+@
+@
}@d
+K
+K
+K
+K
+K
+
+
+
+
+
x
+
+
-
+
-
+
-
-
Y
-r
-{
-
-
-
-
-
-
-2
d
-k
-{
+'
+2
+H
+}
+
-2
+\
o
-
+
+
-2
-8
-J
-`
-d
-j
-p
:
-
-l,6
-)
"G)'
-
-
-^Zp;sZ+8
-AB^bn]
-Z
-n4
-b
-*o
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
P
OJ
^`OJ
`
-
-
-
-
-
-p
-^
-`ph
-
-
-
-
-
-
-
-
-
-
-
`
-
-
-
-
-
-
-
-
p^p`OJ
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
^p`
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-^
-`
-
-^2
-`OJ
OJ
-
-^
-`
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
^`
OJ
^`
.
^`
-
-
-
-
-
-
-
-
-
-
-^
-`
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+(
P
:
+^Zp;K`|$~
-
-
-
-
-
-
-
-
-
-
-^Z
-^Z
-
-n
-
-
-`GT^5bBJ>*N
-2
-?
-H
-T
-Y
-p
-w
N
-
-
-
-
-{{w{ww{{
-RZg9kZ^ZRNswZNsVZRcwkZRNskZJR^ckZZg9s^VkZVZRs
-swo{wo{o{ww
-wwwww
-
+
+^Z
+`GT^5bBJ>*N
+-
+Q
+|
+}
+
+
+
+
r
-
+
+
+
-        
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
\ No newline at end of file
+
+
+
\ No newline at end of file
</font>
</pre>