Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

ICE Thickness/Accumulation - CANICE variable in the WRF

utkarsh-7

New member
Hi,
I was actually facing difficulty in figuring out whether the CANICE(mm) output in WRF is incremental or cumulative(total ice accumulated from start till that particular time instant)?
The reason why I am asking is because I ran a simulations for an ice storms and then after observing the output this variable was showing very small values when actually it should be showing very large values of ice accumulation.

Also if you can suggest some changes in the namelist files to get accurate output for ice accumulation.

Output -
Timesouth_northwest_eastXLATXLONGXTIMECANICE
0​
1​
0​
32.68152​
-96.9915​
03-02-2022 09:20​
0​
0​
1​
1​
32.68532​
-96.9635​
03-02-2022 09:20​
0.034570932​
0​
1​
2​
32.68915​
-96.9354​
03-02-2022 09:20​
0.035885155​
0​
1​
3​
32.69304​
-96.9074​
03-02-2022 09:20​
0.047002126​
0​
1​
4​
32.69687​
-96.8793​
03-02-2022 09:20​
0.054122966​
0​
1​
5​
32.7007​
-96.8513​
03-02-2022 09:20​
0.08609927​
0​
1​
6​
32.7045​
-96.8232​
03-02-2022 09:20​
0.1374387​
0​
1​
7​
32.70839​
-96.7952​
03-02-2022 09:20​
0.3059743​

Namelist files -
wps -
&share
nocolons = .true.
max_dom = 2
start_date = '2022-02-02_00:00:00', '2022-02-02_00:00:00'
end_date = '2022-02-03_21:00:00', '2022-02-03_21:00:00'
interval_seconds = 10800
/

&geogrid
parent_id = 1, 1
parent_grid_ratio = 1, 3
i_parent_start = 1, 3
j_parent_start = 1, 4
e_we = 10, 10
e_sn = 10, 10
map_proj = 'lambert'
dx = 9000
dy = 9000
ref_lon = -96.808891
ref_lat = 32.779167
geog_data_res = 'lowres', 'lowres'
geog_data_path = '/wrf_14/DATA/geog/WPS_GEOG_LOW_RES'
truelat1 = 3.5
truelat2 = 7.0
stand_lon = 4.0
/

&metgrid
fg_name = 'FILE'
/
input -
&time_control
start_year = 2022, 2022
start_month = 2, 2
start_day = 02, 02
start_hour = 00, 00
end_year = 2022, 2022
end_month = 2, 2
end_day = 03, 03
end_hour = 21, 21
interval_seconds = 10800
input_from_file = .true., .true.
history_interval = 10, 10
frames_per_outfile = 100, 100
restart = .false.
restart_interval = 7200
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
start_minute = 0, 0
start_second = 0, 0
end_minute = 0, 0
end_second = 0, 0
nocolons = .true.
/

&domains
time_step = 40
time_step_fract_num = 0
time_step_fract_den = 1
max_dom = 2
e_we = 10, 10
e_sn = 10, 10
e_vert = 44, 44
p_top_requested = 1000
num_metgrid_levels = 34
num_metgrid_soil_levels = 4
dx = 9000.0, 3000.0
dy = 9000.0, 3000.0
grid_id = 1, 2
parent_id = 1, 1
i_parent_start = 1, 3
j_parent_start = 1, 4
parent_grid_ratio = 1, 3
parent_time_step_ratio = 1, 3
feedback = 1
smooth_option = 0
/

&physics
physics_suite = 'CONUS'
sf_surface_physics = 4,4
mp_physics = 4, 4
sf_sfclay_physics = 2, 2
mp_zero_out = 2
cu_physics = 1, 1
ra_lw_physics = 1,1
ra_sw_physics = 1,1
bl_pbl_physics = 2,2
isfflx = 1
radt = 9, 3
bldt = 0, 0
cudt = 0, 0
icloud = 2
num_land_cat = 21
sf_urban_physics = 1, 1
/

&noah_mp
dveg = 4,
opt_crs = 1,
opt_btr = 1,
opt_run = 3,
opt_sfc = 1,
opt_frz = 1,
opt_inf = 1,
opt_rad = 3,
opt_alb = 2,
opt_snf = 1,
opt_tbot = 2,
opt_stc = 1,
opt_gla = 1,
opt_rsf = 1,
opt_soil = 1,
opt_pedo = 1,
opt_crop = 0,
/

&dynamics
hybrid_opt = 2
w_damping = 0
diff_opt = 1, 1
km_opt = 4, 4
diff_6th_opt = 0, 0, 0
diff_6th_factor = 0.12, 0.12
base_temp = 290.0
damp_opt = 3
zdamp = 5000.0, 5000.0
dampcoef = 0.2, 0.2
khdif = 0, 0
kvdif = 0, 0
/

&bdy_control
spec_bdy_width = 5
specified = .true.
/

&fdda
/

&namelist_quilt
/


Thank you for your help,
Utkarsh
 
Top