How to suppress false convection in the upper layer of the model and avoid w-cfl problem

[husiyang22]

hello everyone! I run an idealized case for several times and the w-cfl problems always happen... I don't know if it is related to my eta_levels sets?
below are my namelist:
&time_control
run_days = 100,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2024,
start_month = 01,
start_day = 01,
start_hour = 0,
start_minute = 00,
start_second = 00,
history_interval = 60,
frames_per_outfile = 1,
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
/

&domains
time_step = 10,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
s_we = 1,
e_we = 513,
s_sn = 1,
e_sn = 513,
s_vert = 1,
e_vert = 65,
dx = 3000,
dy = 3000,
ztop = 28000.,
eta_levels = 1.000000, 0.996000, 0.992074, 0.987000, 0.982000, 0.976000, 0.969000, 0.962800, 0.957000, 0.950000, 0.942500, 0.934000, 0.925000,0.917000,
0.909300, 0.887500, 0.875600, 0.863600, 0.851500, 0.839300, 0.817000, 0.804600,
0.792100, 0.779500, 0.766800, 0.753000, 0.741100, 0.728000, 0.715000,0.701800, 0.686500, 0.671000,
0.658000, 0.634400, 0.615600, 0.596700, 0.574600, 0.554400, 0.534000, 0.511500,
0.491800, 0.470000, 0.45000, 0.429800, 0.409400, 0.388800, 0.365000, 0.340000,
0.325000, 0.308000, 0.280000, 0.261000, 0.241000, 0.220000, 0.200000,
0.180000, 0.160000, 0.140000, 0.120000, 0.100000, 0.080000,
0.060000, 0.040000, 0.020000, 0.000000, !

/

&physics
mp_physics = 6,
ra_lw_physics = 4,
ra_sw_physics = 4,
radt = 3,
sf_sfclay_physics = 1,
sf_surface_physics = 1,
bl_pbl_physics = 1,
bldt = 0,
cu_physics = 0,
cudt = 0,
bucket_mm =10.,
bucket_J =1.e8,
/

&fdda
/

&dynamics
hybrid_opt = 0,
rk_ord = 3,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
damp_opt = 3,
dampcoef = .2,
zdamp = 12000.,
w_damping=1,
khdif = 0,
kvdif = 0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.2,
time_step_sound = 6,
h_mom_adv_order = 5,
v_mom_adv_order = 3,
h_sca_adv_order = 5,
v_sca_adv_order = 3,
mix_full_fields = .true.,
non_hydrostatic = .true.,
pert_coriolis = .false.,
/

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

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&ideal
ideal_case = 3
/
I would appreciate a lot if you can give some suggestions!

 

[Ming Chen]

Which ideal case are you running? How did you create initial condition for this case? How long did the model run before it crashed? In addiiton to CFL violation, did you see other error messages?

 

[husiyang22]

Dear Chen Ming:

Which ideal case are you running?

I'm running the em_convrad

How did you create initial condition for this case?

I use a temperature and humidity profile got by a 100-days small region simulation of radiation convection equilibrium.
Besides, I set the solar radiation fixed at 413.9W/m2, and change the delt to 0.1K instead of 1K, which is the amplitude of the initial random temperature pertubation of lowest 10 levels.

How long did the model run before it crashed?

The mode stopped running after 11 hours. I obtained approximately 17days simulation results

In addiiton to CFL violation, did you see other error messages?

so far, I did not see other error