[Wrf-users] Limit on Domain Size in ideal.exe?

[Anders Nottrott]

Hello,

I am trying to run ideal.exe for an ideal LES simulation. I have found 
that 'ideal.exe' crashes if I make the number of x and y grid points in 
the domain greater than 102x102. For example, with the following 
configuration 'ideal.exe' runs fine:
  max_dom                             = 1,
  s_we                                = 1,     1,     1,
  e_we                                = 102,   150,    150,
  s_sn                                = 1,     1,     1,
  e_sn                                = 102,   150,    150,
  s_vert                              = 1,     1,     1,
  e_vert                              = 134,   134,    134,
  dx                                  = 30,   50,    16.6667,
  dy                                  = 30,   50,    16.6667,
  ztop                                = 1067,  2000,  2000,

If i only change:
  e_we                                = 103,   150,    150,
  e_sn                                = 103,   150,    150,

'ideal.exe' crashes with the following error:

"
-------------- FATAL CALLED ---------------
  FATAL CALLED FROM FILE:  <stdin>  LINE:    6458
  frame/module_domain.f: Failed to allocate 
grid%rt_tendf_stoch_bye(sm31:em31,sm32:em32,spec_bdy_width).
  -------------------------------------------
wrf_abort
"

Note that any number larger than 102 causes 'ideal.exe' to crash. Has 
anyone else encountered this issue? I have attached the namelist.

Regards,

Anders
-------------- next part --------------
 &time_control
 run_days                            = 0,
 run_hours                           = 1,
 run_minutes                         = 0,
 run_seconds                         = 00,
 start_year                          = 0001, 0001, 0001,
 start_month                         = 01,   01,   01,
 start_day                           = 01,   01,   01,
 start_hour                          = 00,   01,   00,
 start_minute                        = 00,   30,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 0001, 0001, 0001,
 end_month                           = 01,   01,   01,
 end_day                             = 01,   01,   01,
 end_hour                            = 01,   02,   00,
 end_minute                          = 00,   30,   00,
 end_second                          = 00,   00,   00,
 history_interval_m                  = 10,   10,   1,
 history_interval_s                  = 00,  00,   1,
 frames_per_outfile                  = 06, 1000, 1000,
 restart                             = .false.,
 restart_interval_m                  = 60,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 debug_level                         = 0
/

 &domains
 time_step                           = 0,
 time_step_fract_num                 = 3,
 time_step_fract_den                 = 10,
 max_dom                             = 1,
 s_we                                = 1,     1,     1,
 e_we                                = 103,   150,    150,
 s_sn                                = 1,     1,     1,
 e_sn                                = 103,   150,    150,
 s_vert                              = 1,     1,     1,
 e_vert                              = 134,   134,    134,
 dx                                  = 30,   50,    16.6667, 
 dy                                  = 30,   50,    16.6667, 
 ztop                                = 1067,  2000,  2000,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 0,     10,    15,
 j_parent_start                      = 0,     10,    15,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 0,
 smooth_option                       = 0,
 ts_buf_size                         = 200,
 /

 &physics
 mp_physics                          = 0,     0,     0,
 ra_lw_physics                       = 0,     0,     0,
 ra_sw_physics                       = 0,     0,     0,
 radt                                = 0,     0,     0,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 0,     0,     0,
 bl_pbl_physics                      = 0,     0,     0,
 bldt                                = 0,     0,     0,
 cu_physics                          = 0,     0,     0,
 cudt                                = 0,     0,     0,
 isfflx                              = 1,
 ifsnow                              = 0,
 icloud                              = 0,
 num_soil_layers                     = 5,
 mp_zero_out                         = 0,
 /

 &fdda
 /

 &dynamics
 rk_ord                              = 3,
 diff_opt                            = 2,
 km_opt                              = 3,   
 damp_opt                            = 3,
 zdamp                               = 267.,  5000.,  5000.,
 dampcoef                            = 0.2,    0.2,    0.2
 khdif                               = 1.,      1.,    .05,
 kvdif                               = 1.,      1.,    .05,
 c_s                                 = 0.18
 c_k                                 = 0.10
 mix_isotropic                       = 1
 smdiv                               = 0.1,    0.1,    0.1,
 emdiv                               = 0.01,   0.01,   0.01,
 epssm                               = 0.1,    0.1,    0.1
 tke_heat_flux                       = 0.05,   0.24,   0.24,
 time_step_sound                     = 6,      6,      6,
 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,
 mix_full_fields                     = .true., .true., .true.,
 non_hydrostatic                     = .true., .true., .true.,
 pert_coriolis                       = .true., .true., .true.,
 top_lid				  = .true., .true., .true.,
 scalar_adv_opt			  = 1,      1,      1,
 tracer_adv_opt			  = 1,      1,      1,
 sfs_opt				  = 0,      0,      0,
 m_opt					  = 1,      1,      1,
 tracer_opt				  = 2,      2,      2,
 /

 &bdy_control
 periodic_x                          = .true., .false.,.false.,
 symmetric_xs                        = .false.,.false.,.false.,
 symmetric_xe                        = .false.,.false.,.false.,
 open_xs                             = .false.,.false.,.false.,
 open_xe                             = .false.,.false.,.false.,
 periodic_y                          = .true., .false.,.false.,
 symmetric_ys                        = .false.,.false.,.false.,
 symmetric_ye                        = .false.,.false.,.false.,
 open_ys                             = .false.,.false.,.false.,
 open_ye                             = .false.,.false.,.false.,
 nested                              = .false., .true., .true.,
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /