[Wrf-users] Crash at time step with new BCs

Hald, Cornelius (IMK) cornelius.hald at kit.edu
Mon Feb 20 06:26:03 MST 2017


Dear WRF-Users,

I'm running WRF3.7.1 for a single domain in LES-mode (small domain 
15x15km, 30m dx=dy). I've had similar configurations running well in the 
past. During the last two runs I ran into a problem that is new for me:

At the time step when a new set of boundary conditions (in my case 
coming from the ECMWF Interim dataset) is supposed to be handed over to 
the model, it crashes. To be more precise: The model ran from 00:00 to 
06:00 simulated time. First boundary conditions should be set on 00:00 
and the first update would be on 06:00. At time step 5:59:59 **/** 
(fraction of a second) the whole domain explodes with violations of the 
vertical cfl. This happens at the domain's borders as well as in the 
middle and at all heights. Output and restart files for 06:00 are 
written, then the model aborts. The restart file is useless, as it will 
crash immediately when using it.

I also had a model running until 12h model time which means that at 
06:00 everything worked while it crashed at the next point with new BCs, 
I don't think it has to do with faulty wrfinput-files.

Is it possible that the model solution diverges so extremely from the 
BCs that this causes a crash? Is there any way to avoid the problem?

I appreciate any help you can give. For further information I've 
attached one of my namelists to this mail.

Thanks in advance,

Cornelius

-------------- next part --------------
&time_control
 run_days                            = 0,
 run_hours                           = 42,
 run_minutes                         = 30,
 run_seconds                         = 0,
 start_year                          = 2015,
 start_month                         = 07,
 start_day                           = 15,
 start_hour                          = 05,
 start_minute                        = 30,
 start_second                        = 00,
 end_year                            = 2015,
 end_month                           = 07,
 end_day                             = 17,
 end_hour                            = 00,
 end_minute                          = 00,
 end_second                          = 00,
 interval_seconds                    = 21600,
 input_from_file                     = .true.,
 override_restart_timers             = .true.,
 history_interval_s                  = 3600,
 frames_per_outfile                  = 1,
 write_hist_at_0h_rst                = .true.,
 restart                             = .true.,
 restart_interval                    = 30,
 io_form_history                     = 2,
 io_form_restart                     = 2,
 io_form_input                       = 2,
 io_form_boundary                    = 2,
 debug_level                         = 0,
 !io_form_auxinput2                   = 2,
 !fine_input_stream                   = 2,
 !history_outname                     = 'LES_ndown_mpphys8_d<domain>_<date>'
 !auxhist23_outname                   = wrfpress_d<domain>_<date>
 !io_form_auxhist23                   = 2,
 !auxhist23_interval_s                = 10800, 3600, 3600, 3600,
 !frames_per_auxhist23                = 100, 180,180,300,
 !output_diagnostics                  = 1,
 !auxhist3_outname                    = "wrfxtrm_mpphys8_d<domain>_<date>"
 !io_form_auxhist3                    = 2,
 !auxhist3_interval                   = 10, !5,5,5,5,
 !frames_per_auxhist3                 = 1,1,1,1,
 iofields_filename                   = "excluded_output.txt","excluded_output.txt","excluded_output.txt","excluded_output.txt",
 ignore_iofields_warning             = .true.
 /

 &domains
 time_step                           = 1,
 time_step_fract_num                 = 1,
 time_step_fract_den                 = 10,
 max_dom                             = 1,
 e_we                                = 500,
 e_sn                                = 500,
 e_vert                               = 151,
 eta_levels   = 1.000, 0.9973, 0.9946, 0.9919, 0.9891, 
                0.9864, 0.9837, 0.981, 0.9783, 0.9756, 
                0.9729, 0.9703, 0.9676, 0.9649, 0.9622, 
                0.9595, 0.9568, 0.9542, 0.9515, 0.9488, 
                0.9462, 0.9435, 0.9408, 0.9382, 0.9355, 
                0.9329, 0.9302, 0.9276, 0.9249, 0.9223, 
                0.9196, 0.917, 0.9144, 0.9117, 0.9091, 
                0.9065, 0.9038, 0.9012, 0.8986, 0.896, 
                0.8934, 0.8907, 0.8881, 0.8855, 0.8829, 
                0.8803, 0.8777, 0.8751, 0.8725, 0.8699, 
                0.8673, 0.8648, 0.8622, 0.8596, 0.857, 
                0.8544, 0.8519, 0.8493, 0.8467, 0.8441, 
                0.8416, 0.839, 0.8365, 0.8339, 0.8314, 
                0.8288, 0.8263, 0.8237, 0.8212, 0.8186, 
                0.8161, 0.8135, 0.811, 0.8085, 0.8059, 
                0.8034, 0.8009, 0.7984, 0.7959, 0.7933, 
                0.7908, 0.7883, 0.7858, 0.7833, 0.7808, 
                0.7783, 0.7758, 0.7733, 0.7708, 0.7683, 
                0.7658, 0.7633, 0.7609, 0.7584, 0.7559, 
                0.7534, 0.7509, 0.7485, 0.746, 0.7435, 
                0.7411, 0.7386, 0.7361, 0.7334, 0.7306, 
                0.7276, 0.7245, 0.7211, 0.7175, 0.7136, 
                0.7094, 0.7048, 0.6998, 0.6943, 0.6883, 
                0.6816, 0.6742, 0.666, 0.6571, 0.6472, 
                0.6363, 0.6244, 0.6115, 0.5974, 0.5823, 
                0.566, 0.5487, 0.5304, 0.511, 0.4908, 
                0.4698, 0.448, 0.4256, 0.4027, 0.3794, 
                0.3557, 0.3317, 0.3075, 0.2833, 0.259, 
                0.2347, 0.2104, 0.1863, 0.1623, 0.1384, 
                0.1147, 0.0913, 0.0681, 0.0451, 0.0224, 
                0.000,
 p_top_requested                     = 50000,
 num_metgrid_levels                  = 38,
 num_metgrid_soil_levels             = 4,
 dx                                  = 30,
 dy                                  = 30,
 grid_id                             = 1,
 parent_id                           = 0,
 i_parent_start                      = 1,
 j_parent_start                      = 1,
 parent_grid_ratio                   = 1,
 parent_time_step_ratio              = 1,
 feedback                            = 0,
 smooth_option                       = 0,
 use_adaptive_time_step              = .true.,
 step_to_output_time                 = .true.,
 target_cfl                          = 0.9, !1.2,
 target_hcfl                         = 0.84,
 max_step_increase_pct               = 51,
 starting_time_step                  = 1,
 starting_time_step_den              = 10,
 max_time_step                       = 1,
 max_time_step_den                   = 4,
 min_time_step                       = 1,
 min_time_step_den                   = 16,
 adaptation_domain                   = 1,
 max_ts_level                        = 151,
 max_ts_locs                         = 3,
 !sfcp_to_sfcp                        = .true.,
 /

 &physics
 mp_physics                          = 4,
 ra_lw_physics                       = 4,
 ra_sw_physics                       = 4,
 radt                                = 4,
 sf_sfclay_physics                   = 2,
 sf_surface_physics                  = 2,
 bl_pbl_physics                      = 0,
 bldt                                = 0,
 cu_physics                          = 0,
 cudt                                = 0,
 isfflx                              = 1,
 ifsnow                              = 0,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,0,
 num_land_cat                        = 28, !28,
 topo_wind                           = 0,0,
 /

 &fdda
 grid_fdda                           = 1,
 gfdda_inname                        = "wrffdda_d<domain>",
 gfdda_interval_m                    = 360,
 gfdda_end_h                         = 2,
 io_form_gfdda                       = 2,
 fgdt                                = 0,
 if_no_pbl_nudging_uv                = 0,
 if_no_pbl_nudging_t                 = 0,
 if_no_pbl_nudging_q                 = 0,
 if_zfac_uv                          = 0,
 if_zfac_t                           = 0,
 if_zfac_q                           = 0,
 guv                                 = 0.0003,
 gt                                  = 0.0003,
 gq                                  = 0.0003,
 if_ramping                          = 0,
 /

 &dynamics
 w_damping                           = 1,
 diff_opt                            = 2,
 km_opt                              = 2,
 diff_6th_opt                        = 2,
 diff_6th_factor                     = 0.3,
 base_temp                           = 290.
 damp_opt                            = 3, !changed from 0 in the nested run
 zdamp                               = 500.,
 dampcoef                            = 0.2,
 khdif                               = 0,0,
 kvdif                               = 0,0,
 non_hydrostatic                     = .true.,
 moist_adv_opt                       = 1,
 scalar_adv_opt                      = 1,
 mix_full_fields                     = .true.,
 mix_isotropic                       = 1,
 mix_upper_bound                     = 0.1,
 time_step_sound                     = 8,
 sfs_opt                             = 2,
 epssm                               = 0.3,
 /

 &bdy_control
 spec_bdy_width                      = 15,
 spec_zone                           = 1,
 relax_zone                          = 14,
 specified                           = .true.,
 nested                              = .false.,
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 8,
 nio_groups = 1,
 /
 
 &diags
 p_lev_diags = 0,
 num_press_levels = 30,
 press_levels =100000, 99000, 98000, 97000, 96000, 95000, 94000, 93000, 92000, 91000, 90000, 89000, 88000, 87000, 86000, 85000, 84000, 83000, 82000, 81000, 80000, 79000, 78000, 77000, 76000, 75000, 74000, 73000, 72000, 71000,
 /


More information about the Wrf-users mailing list