[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