[Wrf-users] sst_update problem
Feng Liu
FLiu at azmag.gov
Mon Dec 13 16:26:27 MST 2010
Hi,
your input data have the interval of 6 hours, I think, you might use history_interval = 60 in your namelist.input for running real.exe and wrf.exe. You may have two ways to fix the problem:
1) Set the interval_seconds = 3600 in your namelist.wps, and run ungrib.exe which will create hourly intermediate files by time-interpolation, run metgrib.exe again, you will gain the hourly met_em.d0*. Then run real.exe and wrf.exe, you should fine.
2) Set end_date in namelist.wps as below, and run WPS
'2010-01-16_00:00:00', '2010-01-16_00:00:00', '2010-01-16_00:00:00',
Thanks.
Feng
From: wrf-users-bounces at ucar.edu [mailto:wrf-users-bounces at ucar.edu] On Behalf Of Gino A. Passalacqua
Sent: Monday, December 13, 2010 3:46 PM
To: wrf-users at ucar.edu; Gino A. Passalacqua
Subject: [Wrf-users] sst_update problem
Dear WRF Users,
I'm having a problem to run wrf.exe with SST_update=1 and your help and suggestions will be greatly appreciated. The SST data I'm using is the RTG_SST 0.500 degree ( http://polar.ncep.noaa.gov/sst/). For boundary conditions I'm using the NCEP FNL GFS data. I have 3 nested grids with 108, 36 and 12 km resolutions, running the model from 00:00 hrs of Jan. 01, 2010 to the 00:00 hrs of Jan. 15, 2010.
The steps followed to get to running wrf.exe are the following:
- Ran successfully geogri.exe
- Ran successfully ungrib.exe for the RTG_SST files.
- Ran successfully ungrib.exe for the FNL GFS files.
- Ran successfully metgrib.exe
- Ran successfully real.exe, creating wrfbdy_d01, wrfinput_d01, wrfinput_d02, wrfinput_d03, and wrflowinp_d01, wrflowinp_d02, wrflowinp_d03.
When I ran wrf.exe, I get the following error:
Time in file: 2010-01-15_00:00:00
Time on domain: 2010-01-01_06:07:12
**WARNING** Time in input file not equal to time on domain **WARNING**
**WARNING** Trying next time in file wrflowinp_d01 ...
d01 2010-01-01_06:07:12 module_io.F: in wrf_get_next_time
d01 2010-01-01_06:07:12 2 input_wrf: wrf_get_next_time current_date:
2010-01-15_00:00:00 Status = -4
2 input_wrf: wrf_get_next_time current_date: 2010-01-15_00:00:00 Sta
tus = -4
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 715
... Could not find matching time in input file wrflowinp_d01
-------------------------------------------
[cli_0]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
Assuming that it could be an error of the time stamp on the input files, I have checked the time stamp on the metgrid files and wrfbdy_d01, wrflowinp_d0* (1 to 3), obtaining the following:
- Metgrib files: 2010-01-01_00:00:00 to 2010-01-15_00:00:00 for nests 01, 02 & 03, every 6 hrs.
- wrfbdy_d01: 2010-01-01_00:00:00 to 2010-01-14_18:00:00 , every 6 hrs.
- wrflowinp_d0* (1 to 3): 2010-01-01_00:00:00 to 2010-01-15_00:00:00 for nests 01, 02 & 03, every 6 hrs.
I tried to change the end date in the namelist.input file to 2010-01-14_00:00:00 to run wrf.exe with the same wrfbdy_d01 and wrflowinp_d0*, but I obtained a similar error.
Am I missing something on the namelist.wps or namelist.input files (attached)?
I have search on the web, WRF user's forum, and the manual but I haven't been able to find any solutions.
I also have a couple of questions about the input of SST data:
- I have only been able to run real.exe when I start the model at 00:00:00 hrs only, if I changed the start or end time to let say 06:00:00 hrs, I get an error. Do you know why that happens? Is that a requirement to use SST data as input?
- For having a successful ungrib.exe run with SST data is it necessary 1 extra day of SST data as input, i.e. if your end date is 2010-01-15_00:00:00, you need the SST data until Jan. 16, 2010 to have a successful run. This is just out of my curiosity to try to understand better how WPS and WRF work in general.
Finally, when I remove the SST_update = 1, wrf.exe runs without any problems.
Any help would be greatly appreciated!!!
Gino
namelist.wps %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
&share
wrf_core = 'ARW',
max_dom = 3,
start_date = '2010-01-01_00:00:00','2010-01-01_00:00:00','2010-01-01_00:00:00',
end_date = '2010-01-15_00:00:00','2010-01-15_00:00:00','2010-01-15_00:00:00',
interval_seconds = 21600,
active_grid = .true.,.true.,.true.,.true.,
io_form_geogrid = 2,
opt_output_from_geogrid_path = '/LUSTRE/gpassala/WRF_DATA/Run',
debug_level = 0,
/
&geogrid
parent_id = 1,1,2,
parent_grid_ratio = 1,3,3,
i_parent_start = 1,6,10,
j_parent_start = 1,6,8,
e_we = 30,58,124,
e_sn = 30,58,130,
geog_data_res = '10m','5m','2m',
dx = 108000,
dy = 108000,
map_proj = 'lambert',
ref_lat = 24.584,
ref_lon = -91.831,
truelat1 = 24.584,
truelat2 = 24.584,
stand_lon = -91.831,
geog_data_path = '/LUSTRE/gpassala/WRF_DATA/Data/geog',
&opt_geogrid_tbl_path = '/LUSTRE/gpassala/WRF_DATA/Run/',
ref_x = 15.0,
ref_y = 15.0,
/
&ungrib
out_format = 'WPS',
prefix = '/LUSTRE/gpassala/WRF_DATA/Run/GFS',
/
&metgrid
fg_name = '/LUSTRE/gpassala/WRF_DATA/Run/GFS','/LUSTRE/gpassala/WRF_DATA/Run/RTGSST',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/LUSTRE/gpassala/WRF_DATA/Run',
&opt_metgrid_tbl_path = '/LUSTRE/gpassala/WRF_DATA/Run',
/
&mod_levs
press_pa = 201300 , 200100 , 100000 ,
95000 , 90000 ,
85000 , 80000 ,
75000 , 70000 ,
65000 , 60000 ,
55000 , 50000 ,
45000 , 40000 ,
35000 , 30000 ,
25000 , 20000 ,
15000 , 10000 ,
5000 , 1000
/
namelist.input %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
&time_control
run_days = 14,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2010, 2010, 2010,
start_month = 01, 01, 01,
start_day = 01, 01, 01,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2010, 2010, 2010,
end_month = 01, 01, 01,
end_day = 15, 15, 15,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
fine_input_stream = 0, 0, 0,
io_form_auxinput2 = 2
history_interval = 180, 180, 180,
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,
diag_print = 2,
debug_level = 500,
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360, 360, 360,
io_form_auxinput4 = 2
/
&tc
/
&fire
/
&domains
time_step = 648,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 30, 58, 124,
e_sn = 30, 58, 130,
e_vert = 35, 35, 35,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 108000, 36000, 12000,
dy = 108000, 36000, 12000,
grid_id = 1,2,3,
parent_id = 1,1,2,
i_parent_start = 1,6,10,
j_parent_start = 1,6,8,
parent_grid_ratio = 1,3,3,
parent_time_step_ratio = 1,3,3,
feedback = 1,
smooth_option = 1,
sfcp_to_sfcp = .true.,
eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
0.934, 0.909, 0.880, 0.845, 0.807,
0.765, 0.719, 0.672, 0.622, 0.571,
0.520, 0.468, 0.420, 0.376, 0.335,
0.298, 0.263, 0.231, 0.202, 0.175,
0.150, 0.127, 0.106, 0.088, 0.070,
0.055, 0.040, 0.026, 0.013, 0.000,
force_sfc_in_vinterp = 1,
interp_type = 2,
extrap_type = 2,
t_extrap_type = 2,
use_surface = .true.,
use_levels_below_ground = .true.,
lagrange_order = 1,
max_ts_locs = 20,
/
&physics
mp_physics = 16, 16, 16,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 1, 1, 1,
radt = 36, 36, 36,
sf_sfclay_physics = 5, 5, 5,
sf_surface_physics = 2, 2, 2,
sf_urban_physics = 0, 0, 0,
bl_pbl_physics = 5, 5, 5,
bldt = 0, 0, 0,
cu_physics = 5, 5, 5,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
mp_zero_out = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
sst_update = 1,
usemonalb = .false.,
tmn_update = 1,
sst_skin = 1,
/
&fdda
/
&scm
/
&dynamics
w_damping = 0,
rk_ord = 3,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 2, 2, 2,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.,
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
use_baseparam_fr_nml = .true.,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
smdiv = 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01,
epssm = 0.1, 0.1, 0.1,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 1, 1, 1,
tke_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.,
/
&dfi
&dfi_opt = 3, 3, 3,
&dfi_nfilter = 7, 7, 7,
&dfi_write_filtered_input = .true., .true., .true.,
&dfi_write_dfi_history =.false.,.false.,.false.,
&dfi_cutoff_seconds = 3600,
&dfi_time_dim = 1000,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20101213/396f36a0/attachment-0001.html
More information about the Wrf-users
mailing list