[Wrf-users] adaptive time step problems continue in WRF 3.3?

Zulauf, Michael Michael.Zulauf at iberdrolaren.com
Fri May 20 13:39:37 MDT 2011


Hi all. . .

While I (and others) had problems with adaptive time-stepping in WRF
3.2, I was hoping that this would be fixed in WRF 3.3.  Unfortunately,
I'm still having problems - although of a different sort.  In WRF 3.2 I
had problems with the data being output at times slightly "off" from
when they were supposed to be output (despite using the namelist options
that are supposed to fix that).  In WRF 3.3, I've gotten a crash that
appears to be related to Grid FDDA.  This is for a test case that worked
with adaptive time-stepping turned off.

The basics of the run are that there are 4 domains (27km, 9km, 3km,
1km), initialized and forced with 0.5 degree GFS, with the nested grids
being initialized at 3 hour intervals.

The crash occurs right after domain 3 is initialized (at hour 6), but
the message is that the time in the input file is not the time on the
domain for the FDDA file for domain 2.  I've included the bottom portion
of rsl.error.0000 (below), as well as my namelist.input.

Any thoughts on this?  Is this something that can be handled by changing
some of my namelist options?  Is this a problem with adaptive time
stepping and FDDA?

Thanks,
Mike



--------------------------------- rsl.error.0000
----------------------------------------------
 d02 2011-05-10_05:56:08+**/** *** Initializing nest domain # 3 from an
input file. ***
 d02 2011-05-10_05:56:08+**/** med_initialdata_input: calling
input_auxinput2
Timing for processing wrfinput file (stream 0) for domain        3:
1.82700 elapsed seconds.
 INPUT LandUse = "USGS"
 D03 3-D analysis nudging for wind is applied and Guv=   0.1000E-03
 D03 3-D analysis nudging for temperature is applied and Gt=
0.1000E-03
 D03 3-D analysis nudging for water vapor mixing ratio is applied and
Gq=   0.1000E-05
 D03 3-D analysis nudging for wind is turned off within the PBL.
 D03 3-D analysis nudging for temperature is turned off within the PBL.
 D03 3-D analysis nudging for water vapor mixing ratio is turned off
within the PBL.
 D03 3-D analysis nudging for wind is turned off below layer 19
 D03 3-D analysis nudging for temperature is turned off below layer 10
 D03 3-D analysis nudging for water vapor mixing ratio is turned off
below layer 10
 D03 analysis nudging is ramped down near the end of the nudging period,
       starting at  17.00h, ending at  18.00h.
 INPUT LandUse = "USGS"
 D03 3-D analysis nudging for wind is applied and Guv=   0.1000E-03
 D03 3-D analysis nudging for temperature is applied and Gt=
0.1000E-03
 D03 3-D analysis nudging for water vapor mixing ratio is applied and
Gq=   0.1000E-05
 D03 3-D analysis nudging for wind is turned off within the PBL.
 D03 3-D analysis nudging for temperature is turned off within the PBL.
 D03 3-D analysis nudging for water vapor mixing ratio is turned off
within the PBL.
 D03 3-D analysis nudging for wind is turned off below layer 19
 D03 3-D analysis nudging for temperature is turned off below layer 10
 D03 3-D analysis nudging for water vapor mixing ratio is turned off
below layer 10
 D03 analysis nudging is ramped down near the end of the nudging period,
       starting at  17.00h, ending at  18.00h.
Timing for main (dt= 48.00): time 2011-05-10_05:56:56 on domain   2:
11.03700 elapsed seconds.
Timing for main (dt= 48.00): time 2011-05-10_05:57:44 on domain   2:
8.30000 elapsed seconds.
Timing for main (dt= 48.00): time 2011-05-10_05:58:32 on domain   2:
4.54300 elapsed seconds.
Timing for main (dt= 48.00): time 2011-05-10_05:59:20 on domain   2:
4.84000 elapsed seconds.
Timing for Writing wrfout_d03_2011-05-10_06:00:00 for domain        3:
17.89600 elapsed seconds.
 d03 2011-05-10_06:00:00 Input data processed for aux input  10 for
domain   3
 WRF TILE   1 IS      1 IE     70 JS      1 JE     36
 WRF NUMBER OF TILES =   1
 D03 3-D analysis nudging reads new data at time =    360.000 min.
 D03 3-D analysis nudging bracketing times =   360.00  540.00 min.
Timing for main (dt= 20.00): time 2011-05-10_06:00:16 on domain   3:
27.35300 elapsed seconds.
Timing for main (dt= 48.00): time 2011-05-10_06:00:08 on domain   2:
32.13700 elapsed seconds.
Timing for main (dt=240.00): time 2011-05-10_06:00:00 on domain   1:
64.06500 elapsed seconds.
Timing for Writing wrfout_d01_2011-05-10_06:00:00 for domain        1:
13.37800 elapsed seconds.
 d01 2011-05-10_06:00:00 Input data processed for aux input  10 for
domain   1
Timing for processing lateral boundary for domain        1:    2.06900
elapsed seconds.
 D01 3-D analysis nudging reads new data at time =    360.000 min.
 D01 3-D analysis nudging bracketing times =   360.00  540.00 min.
Timing for Writing wrfout_d02_2011-05-10_06:00:00 for domain        2:
61.44100 elapsed seconds.
  Time in file: 2011-05-10_06:00:00
  Time on domain: 2011-05-10_06:00:08
 **WARNING** Time in input file not equal to time on domain **WARNING**
  **WARNING** Trying next time in file wrffdda_d02 ...
  Time in file: 2011-05-10_09:00:00
  Time on domain: 2011-05-10_06:00:08
 **WARNING** Time in input file not equal to time on domain **WARNING**
  **WARNING** Trying next time in file wrffdda_d02 ...
  Time in file: 2011-05-10_12:00:00
  Time on domain: 2011-05-10_06:00:08
 **WARNING** Time in input file not equal to time on domain **WARNING**
  **WARNING** Trying next time in file wrffdda_d02 ...
  Time in file: 2011-05-10_15:00:00
  Time on domain: 2011-05-10_06:00:08
 **WARNING** Time in input file not equal to time on domain **WARNING**
  **WARNING** Trying next time in file wrffdda_d02 ...
             3  input_wrf: wrf_get_next_time current_date:
2011-05-10_15:00:00 Status =            -4
 -------------- FATAL CALLED ---------------
 FATAL CALLED FROM FILE:  <stdin>  LINE:     751
  ... Could not find matching time in input file wrffdda_d02
------------------------------------------------------------------------
-----------------------


----------------------------------- namelist.input
--------------------------------------------
&time_control
 run_days                            = 0,
 run_hours                           = 18,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2011,2011,2011,2011,
 start_month                         = 05,05,05,05,
 start_day                           = 10,10,10,10,
 start_hour                          = 00,03,06,09,
 start_minute                        = 00,   00,   00,   00,   00,   00,
 start_second                        = 00,   00,   00,   00,   00,   00,

 end_year                            = 2011,2011,2011,2011,
 end_month                           = 05,05,05,05,
 end_day                             = 10,10,10,10,
 end_hour                            = 18,18,18,18,
 end_minute                          = 00,   00,   00,   00,   00,   00,
 end_second                          = 00,   00,   00,   00,   00,   00,
 interval_seconds                    = 10800,
 input_from_file                     =
.true.,.true.,.true.,.true.,.true.,
 fine_input_stream                   = 0, 2, 2, 2, 
 io_form_auxinput2                   = 2
 history_interval                    = 180,180,60,20,
 frames_per_outfile                  =  1,  1,  1,  1,  1,  1, 
 restart                             = .false.,
 restart_interval                    = 1440,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 debug_level                         = 0
 adjust_output_times                 = .true.
 /

 &domains
 time_step                           = 180,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 4,
 s_we                                =  1,1,1,1,
 e_we                                =  268,622,280,382,
 s_sn                                =  1,1,1,1,
 e_sn                                =  184,418,250,196,
 s_vert                              =  1,  1,  1,  1,  1, 1,
 e_vert                              = 31,  31,  31,  31,  31, 31,
 num_metgrid_levels                  =  27 ,
 eta_levels                          = 1.000, 0.993, 0.980, 0.966,
0.950, 0.933, 0.913, 0.892, 0.869, 0.844, 0.816, 0.786, 0.753, 0.718,
0.680, 0.639, 0.596, 0.550, 0.501, 0.451, 0.398, 0.345, 0.290, 0.236,
0.188, 0.145, 0.108, 0.075, 0.046, 0.021, 0.000,

 p_top_requested                     = 5000,
 dx                                  = 27000,9000,3000,1000,
 dy                                  = 27000,9000,3000,1000,
 grid_id                             = 1,  2,  3,  4,  5,  6,
 parent_id                           = 1,  1,  2,  3,  4,  5,
 i_parent_start                      = 1, 31, 91, 92,
 j_parent_start                      = 1, 23,183, 93,
 parent_grid_ratio                   = 1,  3,  3,  3,  3,  3,
 parent_time_step_ratio              = 1,  3,  3,  3,  3, 3,
 feedback                            = 0,
 smooth_option                       = 2
 use_adaptive_time_step              = .true.
 step_to_output_time                 = .true.
 target_cfl                          = 1.1,1.1,1.1,1.1,
 max_step_increase_pct               = 5, 51, 51, 51, 51, 51
 starting_time_step                  = 180, 60, 20, 6.66666667
 max_time_step                       = 240, 80, 26.66666667, 8.88888889
 min_time_step                       = 27, 9, 3, 1
 adaptation_domain                   = 4
 /

 &physics
 mp_physics                          = 5, 5, 5, 5, 
 ra_lw_physics                       = 1, 1, 1, 1, 
 ra_sw_physics                       = 1, 1, 1, 1, 
 radt                                = 30,    30,    30,    30,    30,
30,
 sf_sfclay_physics                   = 1, 1, 1, 1, 
 sf_surface_physics                  = 1, 1, 1, 1, 
 bl_pbl_physics                      = 1, 1, 1, 1, 
 bldt                                = 0,     0,     0,     0,     0,
0,
 cu_physics                          = 1,     1,     0,     0,     0,
0,
 cudt                                = 5,     5,     5,     0,     0,
0, 
 cam_abs_freq_s                      = 21600,
 levsiz                              = 59,
 paerlev                             = 29,
 cam_abs_dim1                        = 4,
 cam_abs_dim2                        = 31,
 isfflx                              = 1,
 ifsnow                              = 0,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 5,
 sf_urban_physics                    = 0,     0,     0,     0,
 mp_zero_out                         = 0,
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 slope_rad                           = 0,
 topo_shading                        = 0,
 /

 &fdda
 grid_fdda                           = 1,     1,     1,     1,
 gfdda_inname                        = "wrffdda_d<domain>",
 gfdda_interval_m                    = 180, 180,   180,   180,
 gfdda_end_h                         = 18, 18, 18, 18,
 io_form_gfdda                       = 2,
 fgdt                                = 0,     0,     0,     0,
 if_no_pbl_nudging_uv                = 1,     1,     1,     1,
 if_no_pbl_nudging_t                 = 1,     1,     1,     1,
 if_no_pbl_nudging_q                 = 1,     1,     1,     1,
 if_zfac_uv                          = 1,     1,     1,     1,
  k_zfac_uv                          = 19,   19,    19,    19,
 if_zfac_t                           = 1,     1,     1,     1,
  k_zfac_t                           = 10,   10,    10,    10,
 if_zfac_q                           = 1,     1,     1,     1,
  k_zfac_q                           = 10,   10,    10,    10,
 guv                                 = 0.0001,     0.0001,     0.0001,
0.0001,
 gt                                  = 0.0001,     0.0001,     0.0001,
0.0001,
 gq                                  = 0.000001,   0.000001,   0.000001,
0.000001,
 if_ramping                          = 1,
 dtramp_min                          = -60.0,
/

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

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /
------------------------------------------------------------------------
-----------------------

-- 

PLEASE NOTE - NEW E-MAIL ADDRESS:
michael.zulauf at iberdrolaren.com

Mike Zulauf
Meteorologist, Lead Senior
Wind Asset Management 
Iberdrola Renewables
1125 NW Couch, Suite 700
Portland, OR 97209
Office: 503-478-6304  Cell: 503-913-0403



Please be advised that email addresses for Iberdrola Renewables personnel have changed to first.last at iberdrolaREN.com effective Aug. 16, 2010.  Please make a note.  Thank you.

This message is intended for the exclusive attention of the recipient(s) indicated.  Any information contained herein is strictly confidential and privileged.  If you are not the intended recipient, please notify us by return e-mail and delete this message from your computer system. Any unauthorized use, reproduction, alteration, filing or sending of this message and/or any attached files may lead to legal action being taken against the party(ies) responsible for said unauthorized use. Any opinion expressed herein is solely that of the author(s) and does not necessarily represent the opinion of the Company. The sender does not guarantee the integrity, speed or safety of this message, and does not accept responsibility for any possible damage arising from the interception, incorporation of viruses, or any other damage as a result of manipulation.




More information about the Wrf-users mailing list