Hi Ming!
I didnot find "R"THCUTEN in my output files. However, there are "A"THCUTEN, "A"THMPTEN and "A"THSHTEN .
Is there a spell error?
I check the value of the 3 variables. Surprisingly, the ATHCUTEN and ATHSHTEN are all 0 during the entire simulation. Only ATHMPTEN is not zero. Can you tell me the reason?
Thank you very much! I upload my namelist file below:
&time_control
run_days = 6,
run_hours = 1,
run_minutes = 0,
run_seconds = 0,
start_year = 2024,
start_month = 01,
start_day =21,
start_hour = 0,
start_minute = 00,
start_second = 00,
history_interval = 60,
frames_per_outfile = 1,
restart = .true.,
restart_interval = 2880,
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 = 385,
s_sn = 1,
e_sn = 385,
s_vert = 1,
e_vert = 50,
dx = 3000,
dy = 3000,
ztop = 28000.,
eta_levels = 1.000, 0.995, 0.990, 0.983, 0.975,
0.965, 0.953, 0.938, 0.920, 0.898,
0.876, 0.854, 0.832, 0.810, 0.788,
0.766, 0.744, 0.722, 0.700, 0.678,
0.656, 0.634, 0.612, 0.590, 0.570,
0.550, 0.530, 0.510, 0.490, 0.470,
0.450, 0.430, 0.410, 0.390, 0.370,
0.350, 0.330, 0.310, 0.290, 0.270,
0.250, 0.225, 0.195, 0.160, 0.120,
0.070, 0.035, 0.015, 0.005, 0.000,
/
&physics
sf_ocean_physics=1,
oml_hml0=20,
oml_relaxation_time=1296000,
oml_gamma=0.14,
acc_phy_tend = 1,
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 = 13000.,
w_damping=1,
khdif = 0,
kvdif = 0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.4,
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 = .true.,
/
&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
/