Vert_cfl exceed after wrf 3dvar

Issues with running 3DVAR.

Vert_cfl exceed after wrf 3dvar

Postby alexiceview » Tue Jan 19, 2010 12:15 pm

Hi:

I have run a fog case with two domains (parent domain and nested domain). Everything is successful to run ./real.exe and ./wrf.exe, and the namelist.input is specified as below:

&time_control
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2007, 2007, 2007,
start_month = 10, 10, 10,
start_day = 24, 24, 24,
start_hour = 12, 12, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2007, 2007, 2007,
end_month = 10, 10, 10,
end_day = 25, 25, 25,
end_hour = 12, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 180, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/

&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 95, 142, 139,
e_sn = 95, 136, 145,
e_vert = 51, 51, 51,
eta_levels = 1.000, 0.997, 0.994, 0.991, 0.988, 0.985, 0.982, 0.976, 0.970, 0.964,
0.956, 0.948, 0.940, 0.932, 0.924, 0.916, 0.908, 0.900, 0.890, 0.880,
0.865, 0.850, 0.835, 0.820, 0.800, 0.780, 0.760, 0.740, 0.720, 0.700,
0.675, 0.650, 0.625, 0.600, 0.575, 0.550, 0.525, 0.500, 0.450, 0.400,
0.350, 0.300, 0.250, 0.200, 0.150, 0.100, 0.050, 0.030, 0.020, 0.005, 0,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000,
dy = 27000, 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 26, 47,
j_parent_start = 1, 27, 47,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0,
use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.2, 1.2, 1.2,
max_step_increase_pct = 5, 51, 51,
starting_time_step = -1, -1, -1,
max_time_step = -1, -1, -1,
min_time_step = -1, -1, -1,
/

&physics
mp_physics = 6, 6, 6,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 27, 27, 27,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 1,
num_land_cat = 20,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/

&fdda
/

&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/



As it seems to be OK,then, I assimilate wrfinput_d01 and update wrfbdy_d01 with surface observations and sounding data and rerun the ./wrf.exe. However, vert_cfl exceeding happen at the very beginning of running ./wrf.exe in many points as showed below:

Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND 20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICAION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
Timing for Writing wrfout_d01_2007-10-24_12:00:00 for domain 1: 0.57100 elapsed seconds.
Timing for processing lateral boundary for domain 1: 0.05200 elapsed seconds.
WRF NUMBER OF TILES = 1
Timing for Writing wrfout_d02_2007-10-24_12:00:00 for domain 2: 0.87700 elapsed seconds.
WRF NUMBER OF TILES = 1
Timing for main: time 2007-10-24_12:00:40 on domain 2: 5.40100 elapsed seconds.
d02 2007-10-24_12:00:40 52 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:00:40 hours
d02 2007-10-24_12:00:40 MAX AT i,j,k: 2 2 17 vert_cfl,w,d(eta)= 3.199152 7.351084 8.0000162E-03
d02 2007-10-24_12:00:40 5355 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:00:40 hours
d02 2007-10-24_12:00:40 MAX AT i,j,k: 2 2 17 vert_cfl,w,d(eta)= 4.306213 7.319943 8.0000162E-03
d02 2007-10-24_12:00:40 5558 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:00:40 hours
d02 2007-10-24_12:00:40 MAX AT i,j,k: 2 2 37 vert_cfl,w,d(eta)= 4.203724 33.97660 2.4999976E-02
Timing for main: time 2007-10-24_12:01:20 on domain 2: 2.88900 elapsed seconds.
d02 2007-10-24_12:01:20 5785 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:01:20 hours
d02 2007-10-24_12:01:20 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 4.329267 -2.794484 2.4999976E-02
d02 2007-10-24_12:01:20 5722 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:01:20 hours
d02 2007-10-24_12:01:20 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 4.678005 -3.241263 2.4999976E-02
d02 2007-10-24_12:01:20 5663 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:01:20 hours
d02 2007-10-24_12:01:20 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 4.836264 -3.462708 2.4999976E-02
Timing for main: time 2007-10-24_12:02:00 on domain 2: 2.91700 elapsed seconds.
Timing for main: time 2007-10-24_12:02:00 on domain 1: 18.58100 elapsed seconds.
d01 2007-10-24_12:02:00 419 points exceeded cfl=2 in domain d01 at time 2007-10-24_12:02:00 hours
d01 2007-10-24_12:02:00 MAX AT i,j,k: 27 43 37 vert_cfl,w,d(eta)= 2.418804 0.5768565 2.4999976E-02
d01 2007-10-24_12:02:00 958 points exceeded cfl=2 in domain d01 at time 2007-10-24_12:02:00 hours
d01 2007-10-24_12:02:00 MAX AT i,j,k: 27 28 37 vert_cfl,w,d(eta)= 3.330005 -14.06733 2.4999976E-02
d01 2007-10-24_12:02:00 1003 points exceeded cfl=2 in domain d01 at time 2007-10-24_12:02:00 hours
d01 2007-10-24_12:02:00 MAX AT i,j,k: 27 28 37 vert_cfl,w,d(eta)= 3.946744 -2.834586 2.4999976E-02
d02 2007-10-24_12:02:00 6069 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:02:00 hours
d02 2007-10-24_12:02:00 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 5.348660 -4.119401 2.4999976E-02
d02 2007-10-24_12:02:00 5816 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:02:00 hours
d02 2007-10-24_12:02:00 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 5.557653 -3.853061 2.4999976E-02
d02 2007-10-24_12:02:00 5729 points exceeded cfl=2 in domain d02 at time 2007-10-24_12:02:00 hours
d02 2007-10-24_12:02:00 MAX AT i,j,k: 1 5 37 vert_cfl,w,d(eta)= 5.628192 -3.719646 2.4999976E-02
....................................................................................................................................................
....................................................................................................................................................
....................................................................................................................................................
....................................................................................................................................................


This problem is still even I decrease the time_step and turn on the w_damping. It may be something wrong about WRF 3DVAR process. However, everything shows successfully complete in the process. The version of wrf and wrfvar I use are V3.1.1. And the set of namelist.obsproc is showed as:

&record1
obs_gts_filename = '../obs_gts/obs_gts.07102412',
fg_format = 'WRF',
obs_err_filename = '../obserr.txt',
/

&record2
time_window_min = '2007-10-24_11:00:00',
time_analysis = '2007-10-24_12:00:00',
time_window_max = '2007-10-24_13:00:00',
/

&record3
max_number_of_obs = 400000,
fatal_if_exceed_max_obs = .TRUE.,
/

&record4
qc_test_vert_consistency = .TRUE.,
qc_test_convective_adj = .TRUE.,
qc_test_above_lid = .TRUE.,
remove_above_lid = .TRUE.,
domain_check_h = .TRUE.,
Thining_SATOB = .FALSE.,
Thining_SSMI = .FALSE.,
Thining_QSCAT = .FALSE.,
/

&record5
print_gts_read = .FALSE.,
print_gpspw_read = .FALSE.,
print_recoverp = .FALSE.,
print_duplicate_loc = .FALSE.,
print_duplicate_time = .FALSE.,
print_recoverh = .FALSE.,
print_qc_vert = .FALSE.,
print_qc_conv = .FALSE.,
print_qc_lid = .FALSE.,
print_uncomplete = .FALSE.,
/

&record6
ptop = 1000.0,
base_pres = 100000.0,
base_temp = 290.0,
base_lapse = 50.0,
base_strat_temp = 215.0,
base_tropo_pres = 20000.0
/

&record7
IPROJ = 1,
PHIC = 39.0,
XLONC = 117.0,
TRUELAT1= 30.0,
TRUELAT2= 60.0,
MOAD_CEN_LAT = 39.0,
STANDARD_LON = 117.0,
/

&record8
IDD = 1,
MAXNES = 3,
NESTIX = 95, 136, 145, 181, 211,
NESTJX = 95, 142, 139, 196, 211,
DIS = 27, 9., 3., 1.1, 1.1,
NUMC = 1, 1, 2, 3, 4,
NESTI = 1, 26, 47, 35, 45,
NESTJ = 1, 27, 47, 65, 55,
/

&record9
PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output_filename',
PREPBUFR_TABLE_FILENAME = '../prepbufr_table_filename',
OUTPUT_OB_FORMAT = 2
use_for = '3DVAR',
num_slots_past = 3,
num_slots_ahead = 3,
write_synop = .true.,
write_ship = .true.,
write_metar = .true.,
write_buoy = .true.,
write_pilot = .true.,
write_sound = .true.,
write_amdar = .true.,
write_satem = .true.,
write_satob = .true.,
write_airep = .true.,
write_gpspw = .true.,
write_gpsztd= .true.,
write_gpsref= .true.,
write_gpseph= .true.,
write_ssmt1 = .true.,
write_ssmt2 = .true.,
write_ssmi = .true.,
write_tovs = .true.,
write_qscat = .true.,
write_profl = .true.,
write_bogus = .true.,
write_airs = .true.,
/

And the namelist.input for da_wrfvar.exe are specified as:

&wrfvar1
print_detail_grad=false,
/
&wrfvar2
/
&wrfvar3
/
&wrfvar4
/
&wrfvar5
/
&wrfvar6
/
&wrfvar7
cv_options=3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2007-10-24_12:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2007-10-24_11:00:00.0000",
/
&wrfvar22
time_window_max="2007-10-24_13:00:00.0000",
/
&wrfvar23
/
&time_control
start_year=2007,
start_month=10,
start_day=24,
start_hour=12,
end_year=2007,
end_month=10,
end_day=24,
end_hour=12,
/
&fdda
/
&domains
e_we=95,
e_sn=95,
e_vert=51,
dx=27000,
dy=27000,
/
&dfi_control
/
&tc
/
&physics
mp_physics=6,
ra_lw_physics=1,
ra_sw_physics=1,
radt=27,
sf_sfclay_physics=1,
sf_surface_physics=2,
bl_pbl_physics=1,
cu_physics=1,
cudt=5,
num_soil_layers=4,
num_land_cat=20,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&fire
/
&namelist_quilt
/



The set of parame.in of da_update_bc.exe is:

&control_param
wrfvar_output_file = '/home/alex/WRFVAR/WRFDA_V3.1.1/var/test/tutorial/fog071024-25fg01/wrfvar_output'
wrf_bdy_file = './wrfbdy_d01'
wrf_input = '/home/alex/wrf/WRFV3/test/em_real/wrfinput_d01'
domain_id = 1
cycling = .false.
debug = .true.
low_bdy_only = .false.
update_lsm = .false.
/

Some one can help me on this problem? Any suggestion or solution about this problem is welcomed!
Thanks in advance!

Yours Sincerely
Alex
alexiceview
 
Posts: 1
Joined: Sun Jan 03, 2010 11:24 pm

Re: Vert_cfl exceed after wrf 3dvar

Postby piyush » Thu May 06, 2010 2:56 am

Reference to your query help we are also facing the same problem in 3dvar.
Hopefully you might have got the solution for that.
Can you suggest waht could be done to resolve this problem?????

Regards
sravan (sravan_sase@yahoo.in)
Piyush
piyush
 
Posts: 2
Joined: Tue May 04, 2010 12:52 am

Re: Vert_cfl exceed after wrf 3dvar

Postby enjoywill » Sun May 16, 2010 8:19 am

Hi alexiceview
Is it necessary to update your wrfinput_d02 and wrfinput_d03?
enjoywill
 
Posts: 3
Joined: Sun May 16, 2010 8:13 am


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 2 guests