WRF-Chem : WARNING: Large total sw optical depth of
27 views
Skip to first unread message
kamo karapetian
Dec 2, 2025, 5:24:38 AM12/2/25
to wrf-chem-run
Hello everyone,
I am currently experiencing some difficulties with WRF-Chem(chem_opt=202). I am using WRF v4.5.2 and WPS v4.5.
I first ran a 1 hour simulation without any problem, then extended it to 3 hours, again without any issues. However, when I extended it to a 24 hour run, I encountered a problem. I would like to emphasize that for each case I adapted and downloaded the appropriate input data.
The message is going to be long, but it will be very detailed and complete in order to clearly explain the situation and the origin of the problem.
We are using a CESM file as GFS chemical boundary data for ungrib.
The global file was obtained here:
https://www.acom.ucar.edu/cesm/subset.shtml
We filled in the following information for our project:
Full name: Kamo
Affiliation: LATMOS
Address:
Save contact: on
lat2: 52.3
lon1: 356
lon2: 10
lat1: 45
Start year: 2021
Start month: 07
Start day: 01
End year: 2021
End month: 07
End day: 02
simulation: cesm
I am also using EDGAR+anthro, MEGAN, and FINN fires. I am attaching all these files for clarity.
I have already seen this problem mentioned in different threads, but no clear answer or fix has been provided, so I now find myself in the same situation as the others.
My rsl.error and rsl.out files are very large because of the extensive debugging, so I am only including the relevant part related to the problem.
Here is the context of the error message:
WARNING: Large total sw optical depth of 48.79 at point i,j,nb= 113 46 1
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1 8.88 13.52 15.18 8.33 12.05
2 22.99 34.20 39.99 23.81 36.41
3 0.16 0.21 0.24 0.14 0.18
4 0.09 0.14 0.17 0.10 0.15
5 0.00 0.00 0.00 0.00 0.00
6 0.00 0.00 0.00 0.00 0.00
7 0.00 0.00 0.00 0.00 0.00
8 0.00 0.00 0.00 0.00 0.00
9 0.00 0.00 0.00 0.00 0.00
10 0.00 0.00 0.00 0.00 0.00
11 0.00 0.00 0.00 0.00 0.00
12 0.00 0.00 0.00 0.00 0.00
13 0.00 0.00 0.00 0.00 0.00
14 0.00 0.00 0.00 0.00 0.00
15 0.00 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00 0.00
25 0.00 0.00 0.00 0.00 0.00
26 0.00 0.00 0.00 0.00 0.00
27 0.00 0.00 0.00 0.00 0.00
28 0.00 0.00 0.00 0.00 0.00
29 0.00 0.00 0.00 0.00 0.00
30 0.00 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00 0.00
35 0.00 0.00 0.00 0.00 0.00
36 0.00 0.00 0.00 0.00 0.00
37 0.00 0.00 0.00 0.00 0.00
38 0.00 0.00 0.00 0.00 0.00
39 0.00 0.00 0.00 0.00 0.00
40 0.00 0.00 0.00 0.00 0.00
41 0.00 0.00 0.00 0.00 0.00
42 0.00 0.00 0.00 0.00 0.00
43 0.00 0.00 0.00 0.00 0.00
44 0.00 0.00 0.00 0.00 0.00
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.62 0.76 0.79 0.68
2 0.66 0.73 0.80 0.68
3 0.69 0.75 0.79 0.69
4 0.65 0.73 0.79 0.67
5 0.73 0.73 0.72 0.58
6 0.68 0.71 0.74 0.62
7 0.69 0.70 0.72 0.62
8 0.69 0.70 0.72 0.63
9 0.67 0.71 0.73 0.64
10 0.64 0.71 0.77 0.65
11 0.65 0.71 0.77 0.64
12 0.64 0.71 0.76 0.63
13 0.64 0.72 0.76 0.64
14 0.64 0.74 0.78 0.67
15 0.66 0.74 0.78 0.71
16 0.66 0.74 0.77 0.70
17 0.66 0.72 0.77 0.71
18 0.67 0.74 0.78 0.74
19 0.65 0.72 0.78 0.71
20 0.55 0.61 0.69 0.60
21 0.66 0.73 0.77 0.66
22 0.66 0.72 0.77 0.65
23 0.66 0.74 0.78 0.66
24 0.66 0.74 0.77 0.66
25 0.64 0.74 0.77 0.68
26 0.66 0.73 0.77 0.66
27 0.63 0.74 0.77 0.67
28 0.63 0.72 0.79 0.72
29 0.56 0.62 0.69 0.60
30 0.56 0.62 0.69 0.60
31 0.59 0.58 0.67 0.59
32 0.59 0.50 0.61 0.57
33 0.61 0.40 0.39 0.51
34 0.61 0.40 0.39 0.51
35 0.61 0.40 0.39 0.51
36 0.61 0.40 0.39 0.51
37 0.61 0.40 0.39 0.50
38 0.61 0.40 0.39 0.50
39 0.61 0.40 0.39 0.50
40 0.61 0.41 0.40 0.51
41 0.61 0.40 0.39 0.51
42 0.61 0.41 0.39 0.51
43 0.61 0.41 0.40 0.51
44 0.61 0.41 0.40 0.51
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 0.91 0.94 0.96 0.95
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 0.99 0.99 1.00 1.00
29 0.91 0.93 0.96 0.95
30 0.91 0.94 0.96 0.95
31 0.92 0.93 0.95 0.93
32 0.94 0.94 0.94 0.93
33 0.95 0.93 0.88 0.85
34 0.95 0.93 0.88 0.85
35 0.95 0.93 0.88 0.85
36 0.95 0.93 0.88 0.84
37 0.95 0.93 0.88 0.85
38 0.95 0.93 0.88 0.84
39 0.95 0.93 0.88 0.85
40 0.95 0.93 0.88 0.85
41 0.95 0.93 0.88 0.85
42 0.95 0.93 0.88 0.85
43 0.95 0.93 0.88 0.85
44 0.95 0.93 0.88 0.85
Diagnostics 4: k, ssaal, asyal, taual
0 1.00 1.00 0.00
1 1.00 1.00 0.00
2 1.00 1.00 0.00
3 1.00 1.00 0.00
4 1.00 1.00 0.00
5 1.00 1.00 0.00
6 1.00 0.00 0.00
7 1.00 0.00 0.00
8 1.00 0.00 0.00
9 1.00 0.00 0.00
10 1.00 0.00 0.00
11 1.00 0.00 0.00
12 1.00 0.00 0.00
13 1.00 0.00 0.00
14 1.00 0.00 0.00
15 1.00 0.00 0.00
16 1.00 0.00 0.00
17 1.00 0.00 0.00
18 1.00 0.00 0.00
19 1.00 0.00 0.00
20 1.00 0.00 0.00
21 1.00 0.00 0.00
22 1.00 0.00 0.00
23 1.00 0.00 0.00
24 1.00 1.00 12.39
25 1.00 1.00 35.86
26 1.00 1.00 0.19
27 1.00 1.00 0.15
28 1.00 0.67 0.00
29 1.00 0.95 0.00
30 1.00 0.85 0.00
31 1.00 0.84 0.00
32 1.00 0.87 0.00
33 1.00 1.00 0.00
34 1.00 1.00 0.00
35 1.00 1.00 0.00
36 1.00 1.00 0.00
37 1.00 1.00 0.00
38 1.00 1.00 0.00
39 1.00 1.00 0.00
40 1.00 1.00 0.00
41 1.00 1.00 0.00
42 1.00 1.00 0.00
43 1.00 0.00 0.00
44 1.00 1.00 0.00
I also see this at different points. Then, at the end, I get:
WARNING: Large total sw optical depth of 25.91 at point i,j,nb= 117 77 14
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1 5.83 6.02 6.51 6.37 7.53
2 7.51 7.56 7.97 8.13 9.22
3 6.99 6.99 6.61 7.15 7.40
4 1.94 2.05 1.93 1.79 1.68
5 0.10 0.09 0.09 0.09 0.08
6 0.01 0.01 0.00 0.00 0.00
7 0.00 0.00 0.00 0.00 0.00
8 0.00 0.00 0.00 0.00 0.00
9 0.00 0.00 0.00 0.00 0.00
10 0.00 0.00 0.00 0.00 0.00
11 0.00 0.00 0.00 0.00 0.00
12 0.00 0.00 0.00 0.00 0.00
13 0.00 0.00 0.00 0.00 0.00
14 0.00 0.00 0.00 0.00 0.00
15 0.00 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00 0.00
25 0.00 0.00 0.00 0.00 0.00
26 0.00 0.00 0.00 0.00 0.00
27 0.00 0.00 0.00 0.00 0.00
28 0.00 0.00 0.00 0.00 0.00
29 0.00 0.00 0.00 0.00 0.00
30 0.00 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00 0.00
35 0.00 0.00 0.00 0.00 0.00
36 0.00 0.00 0.00 0.00 0.00
37 0.00 0.00 0.00 0.00 0.00
38 0.00 0.00 0.00 0.00 0.00
39 0.00 0.00 0.00 0.00 0.00
40 0.00 0.00 0.00 0.00 0.00
41 0.00 0.00 0.00 0.00 0.00
42 0.00 0.00 0.00 0.00 0.00
43 0.00 0.00 0.00 0.00 0.00
44 0.00 0.00 0.00 0.00 0.00
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.84 0.82 0.82 0.83
2 0.85 0.83 0.81 0.81
3 0.83 0.84 0.83 0.80
4 0.83 0.82 0.82 0.80
5 0.84 0.83 0.84 0.82
6 0.81 0.79 0.78 0.78
7 0.74 0.72 0.73 0.62
8 0.65 0.72 0.77 0.71
9 0.65 0.72 0.78 0.71
10 0.66 0.73 0.77 0.69
11 0.65 0.74 0.78 0.69
12 0.66 0.73 0.78 0.70
13 0.66 0.73 0.78 0.71
14 0.66 0.73 0.78 0.71
15 0.64 0.74 0.78 0.73
16 0.68 0.70 0.79 0.77
17 0.72 0.64 0.62 0.73
18 0.70 0.73 0.79 0.76
19 0.66 0.73 0.79 0.76
20 0.66 0.74 0.77 0.70
21 0.66 0.74 0.78 0.69
22 0.66 0.74 0.77 0.66
23 0.67 0.73 0.77 0.64
24 0.64 0.71 0.77 0.64
25 0.65 0.74 0.78 0.69
26 0.65 0.73 0.78 0.71
27 0.65 0.73 0.78 0.74
28 0.65 0.73 0.78 0.74
29 0.56 0.62 0.69 0.60
30 0.56 0.62 0.69 0.60
31 0.57 0.62 0.69 0.60
32 0.58 0.57 0.66 0.59
33 0.60 0.39 0.37 0.49
34 0.61 0.40 0.38 0.50
35 0.61 0.40 0.39 0.51
36 0.61 0.40 0.39 0.50
37 0.61 0.40 0.39 0.51
38 0.61 0.40 0.39 0.51
39 0.61 0.40 0.39 0.51
40 0.61 0.41 0.39 0.51
41 0.61 0.41 0.40 0.51
42 0.61 0.41 0.40 0.51
43 0.61 0.41 0.40 0.51
44 0.61 0.41 0.40 0.51
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 1.00 1.00 1.00 1.00
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 1.00 1.00 1.00 1.00
29 0.91 0.93 0.96 0.95
30 0.91 0.93 0.96 0.95
31 0.91 0.93 0.96 0.95
32 0.93 0.94 0.95 0.94
33 0.95 0.94 0.90 0.85
34 0.95 0.93 0.88 0.84
35 0.95 0.93 0.88 0.85
36 0.95 0.93 0.88 0.84
37 0.95 0.93 0.88 0.85
38 0.95 0.93 0.88 0.85
39 0.95 0.93 0.88 0.85
40 0.95 0.93 0.88 0.85
41 0.95 0.93 0.88 0.85
42 0.95 0.93 0.88 0.85
43 0.95 0.93 0.88 0.85
44 0.95 0.93 0.88 0.85
Diagnostics 4: k, ssaal, asyal, taual
0 0.00 NaN 0.00
1 0.00 NaN 0.00
2 0.00 0.00 0.00
3 0.00 0.00 0.00
4 0.00 0.00 0.00
5 0.00 0.00 0.00
6 0.00 0.00 0.00
7 0.00 0.00 0.00
8 0.00 0.05 0.00
9 0.00 0.00 0.00
10 0.00 0.02 0.00
11 0.00 0.00 0.00
12 0.00 0.00 0.00
13 0.00 0.00 0.00
14******** 0.00 0.00
15 0.00 0.00 0.00
16 13.78******** 0.00
17 14.13 0.00 0.00
18 15.54 0.00 0.00
19 19.00 0.00 0.00
20 19.34 0.00 0.00
21 19.33 0.00 0.00
22 18.70 0.00 0.01
23 20.38 0.00 0.00
24 39.02 0.00 0.00
25 47.78 0.00 0.00
26 51.14 0.00 0.00
27 54.01 0.00 0.00
28 55.83 0.00 0.00
29 99.42 0.00 0.00
30 114.35 10.12 0.00
31 130.16 0.00 0.00
32 145.19 0.00 0.00
33 158.71 0.00 0.00
34 171.34 0.00 0.00
35 182.15 0.00 0.00
36 190.74 0.00 0.00
37 197.49 0.00 0.00
38 202.89 0.00 0.00
39 206.31 0.00 0.00
40 209.10******** 0.00
41 214.38 0.00 0.00
42 219.37 0.00 0.00
43 220.22 0.00 0.00
44 221.02 0.00 0.00
-------------------------
d01 2021-07-01_03:45:00 calling inc/HALO_PWP_inline.inc
d01 2021-07-01_03:45:00 call surface_driver
d01 2021-07-01_03:45:00 in SFCLAY
d01 2021-07-01_03:45:00 in NOAH DRV
d01 2021-07-01_03:45:00 call pbl_driver
d01 2021-07-01_03:45:00 in MYNNPBL
d01 2021-07-01_03:45:00 call HALO CHEM AFTER PBL
d01 2021-07-01_03:45:00 calling inc/HALO_EM_CHEM_E_5_inline.inc
d01 2021-07-01_03:45:00 call cumulus_driver
d01 2021-07-01_03:45:00 calling inc/HALO_CUP_G3_IN_inline.inc
d01 2021-07-01_03:45:00 in grelldrv
d01 2021-07-01_03:45:00 returning from cumulus_driver
d01 2021-07-01_03:45:00 call shallow_cumulus_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_FDDA_SFC_inline.inc
d01 2021-07-01_03:45:00 call fddagd_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SCALAR_E_5_inline.inc
d01 2021-07-01_03:45:00 call calculate_phy_tend
d01 2021-07-01_03:45:00 call compute_diff_metrics
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_A1_inline.inc
d01 2021-07-01_03:45:00 call bc for diffusion_metrics
d01 2021-07-01_03:45:00 call cal_deform_and_div
d01 2021-07-01_03:45:00 calling inc/HALO_EM_HELICITY_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_D_inline.inc
d01 2021-07-01_03:45:00 call calculate_km_kh
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_E_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_PHY_BC_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_CHEM_inline.inc
d01 2021-07-01_03:45:00 call phy_bc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_PBL_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_DIFFUSION_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_5_inline.inc
d01 2021-07-01_03:45:00 call update_phy_ten
d01 2021-07-01_03:45:00 call horizontal_diffusion_2
d01 2021-07-01_03:45:00 call rk_tendency
d01 2021-07-01_03:45:00 call small_step_prep
d01 2021-07-01_03:45:00 calling inc/HALO_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C2_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B3_inline.inc
d01 2021-07-01_03:45:00 call rk_small_finish
d01 2021-07-01_03:45:00 calling inc/HALO_EM_D_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_EM_DA_inline.inc
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend in chem_tile_loop_1
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call advance_ppt
d01 2021-07-01_03:45:00 call phy_prep_part2
d01 2021-07-01_03:45:00 call moist_physics_prep
d01 2021-07-01_03:45:00 call microphysics_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SBM_inline.inc
d01 2021-07-01_03:45:00 microphysics_driver: calling morrison two moment
d01 2021-07-01_03:45:00 microphysics_driver: returning from
d01 2021-07-01_03:45:00 call moist_physics_finish
d01 2021-07-01_03:45:00 calling inc/HALO_EM_THETAM_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_EM_THETAM_inline.inc
d01 2021-07-01_03:45:00 call chem polar filter
d01 2021-07-01_03:45:00 calling inc/HALO_EM_D3_5_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_D3_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_MOIST_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_CHEM_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_TRACER_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_SCALAR_inline.inc
d01 2021-07-01_03:45:00 call set_phys_bc_dry_2
d01 2021-07-01_03:45:00 call spec_bdy_final
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_E_inline.inc
d01 2021-07-01_03:45:00 call set_w_surface
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_W_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_HCW_inline.inc
d01 2021-07-01_03:45:00 --> TOP OF AFTER ALL RK STEPS
d01 2021-07-01_03:45:00 --> CALLING DIAGNOSTICS DRIVER
d01 2021-07-01_03:45:00 --> TOP OF DIAGNOSTICS PACKAGE
d01 2021-07-01_03:45:00 --> CALL DIAGNOSTICS PACKAGE: ACCUMULATED AND BUCKET DIAGNOSTICS
d01 2021-07-01_03:45:00 call HALO_RK_E
d01 2021-07-01_03:45:00 calling inc/HALO_EM_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_MOIST
d01 2021-07-01_03:45:00 calling inc/HALO_EM_MOIST_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_CHEM
d01 2021-07-01_03:45:00 calling inc/HALO_EM_CHEM_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_SCALAR
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SCALAR_E_5_inline.inc
d01 2021-07-01_03:45:00 call end of solve_em
d01 2021-07-01_03:45:00 DEBUG wrf_timetoa(): returning with str = [2021-07-01_03:45:00]
d01 2021-07-01_03:45:00 calling mozart_mosaic_4bin_aq_kpp from chem_driver
d01 2021-07-01_03:45:00 calling chem_opt=? from chem_driver
d01 2021-07-01_03:45:00 call chem_prep
d01 2021-07-01_03:45:00 calling emissions driver
d01 2021-07-01_03:45:00 No tracer option selected
d01 2021-07-01_03:45:00 MOSAIC or SORGAM sea salt emissions
d01 2021-07-01_03:45:00 Gocart dust emissions
d01 2021-07-01_03:45:00 no dms emissions
d01 2021-07-01_03:45:00 fire emissions: adding biomassb emissions
d01 2021-07-01_03:45:00 nothing done with burn emissions for tracers here
d01 2021-07-01_03:45:00 biogenic emissions: calling megan v2.04
calculate MEGAN emissions at ktau, gmtp, tmidh = 901 3.00000000 3.75208330
par out of range at i=149 ,j= 6 par = -2.214
d01 2021-07-01_03:45:00 emissions_driver calling add_anthropogenics
d01 2021-07-01_03:45:00 add_anthropogenics
d01 2021-07-01_03:45:00 emissions_driver calling add_biogenics
d01 2021-07-01_03:45:00 emissions_driver not calling add_emiss_cptec
d01 2021-07-01_03:45:00 emissions_driver calling mosaic_addemiss
d01 2021-07-01_03:45:00 mosaic_addemiss: emiss_opt = eptec
d01 2021-07-01_03:45:00 mosaic_addemiss: gocart speciation being mapped to mosaic
d01 2021-07-01_03:45:00 mosaic_addemiss: adding fire emissions to MOSAIC
d01 2021-07-01_03:45:00 calling optical driver
d01 2021-07-01_03:45:00 optical driver: process aerosols true
d01 2021-07-01_03:45:00 optical averaging
d01 2021-07-01_03:45:00 optical driver: call optical averaging
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 16 3 29 4 0.7404E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 17 3 29 4 0.8284E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 18 3 29 4 0.8648E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 19 3 29 4 0.8779E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 3 5 4 0.7080E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 4 29 4 0.1978E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 4 5 4 0.7746E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 7 28 4 0.1094E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 8 2 4 0.7275E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 27 10 3 4 0.7014E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 20 12 7 4 0.8006E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 24 14 5 4 0.8569E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 4 16 24 4 0.1029E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 18 28 4 0.1344E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 19 28 4 0.1506E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 20 28 4 0.1459E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 21 28 4 0.8573E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 22 15 4 0.4402E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 23 29 4 0.1691E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 24 29 4 0.1435E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 25 29 4 0.1943E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 26 29 3 0.7767E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 26 29 4 0.2388E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 109 26 1 3 0.7794E+00
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 27 29 3 0.1015E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 27 29 4 0.3145E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 14 27 24 4 0.1712E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 18 4 0.2673E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 29 3 0.9201E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 29 4 0.2595E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 31 18 4 0.2107E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 31 29 4 0.1959E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 32 29 4 0.1464E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 33 29 4 0.1458E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 34 29 4 0.1309E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 35 34 28 4 0.1809E-02
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0 0xfe45244251f in ???
#1 0x57b90986cdb3 in ???
#2 0x57b909871fb3 in ???
#3 0x57b909872b64 in ???
#4 0x57b909872dc6 in ???
#5 0x57b9098737e3 in ???
#6 0x57b9098aa680 in ???
#7 0x57b908aefb1a in ???
#8 0x57b9089ee84d in ???
#9 0x57b9087e9ae9 in ???
#10 0x57b906d35fe8 in ???
#11 0x57b906c98be7 in ???
#12 0x57b906c9809e in ???
#13 0xfe452429d8f in ???
#14 0xfe452429e3f in ???
#15 0x57b906c980d4 in ???
#16 0xffffffffffffffff in ???
I have checked my input files and there does not seem to be any problem; the values look consistent. The issue appears during the run itself, and more specifically at this model time. I only noticed that in mozbc I mapped to 8 bins instead of 4 bins, but this is not the cause of the crash.
I also examined my output file at 03:00:00, and here is what I found:
Sulfates (so4_a01 to so4_a04):
Maximum values remain between ~0.25 and ~2.35, with low mean values (~0.004–0.08).
No abnormal behavior is observed.
Sulfates do not seem to be responsible for the large optical depths.
Aerosol number (num_a01 to num_a04):
num_a01 to num_a03 already show high values (up to ~10¹³).
num_a04 is particularly problematic:
max ≈ 1.1 × 10¹⁷
mean ≈ 5.4 × 10¹¹
These orders of magnitude are unrealistic for aerosol number and indicate strong instability or inconsistency in this mode.
SOA (glysoa_*):
For the fine mode (a01), values remain low and consistent.
From a02 onwards, and especially a03 and a04, the maxima become extremely large. For example:
glysoa_sfc_a03 max ≈ 2.0 × 10⁶
glysoa_nh4_a03 max ≈ 2.9 × 10⁷
glysoa_oh_a03 max ≈ 5.1 × 10⁷
glysoa_r1_a04 max ≈ 7.4 × 10⁹
glysoa_nh4_a04 max ≈ 3.1 × 10¹³
glysoa_oh_a04 max ≈ 2.7 × 10¹³
Several of these fields also have very high mean values (up to ~10⁸).
The aqueous terms (glysoa_*_cw0*) also show abnormally large values (max up to ~10⁶–10⁷).
I then examined the outputs at 00:00:00, 01:00:00 and 02:00:00, and here is what I observe:
wrfout_d01_2021-07-01_00:00:00
Sulfates (so4_a0*):
Reasonable values (max ~2.4 µg/m³, mean ~0.08–0.12).
A few negative values (e.g. so4_a03 min = -0.1559), which are nonphysical, but with moderate amplitude. This is already a small numerical issue but not catastrophic.
Aerosol number (num_a0*):
Again, some negative values (e.g. num_a01 min = -1.7×10⁷).
Maxima remain of the usual order of magnitude (10⁹–10¹⁰ for num_a0*).
This is still “manageable”: sign errors, but not an explosion yet.
SOA (glysoa*):
Maxima remain reasonable (often << 10³, sometimes up to 4×10⁴ for glysoa_oh_a03).
No huge values and no absurdly large means.
File 01:00:00 :
Here the real drift begins:
num_a0*:
num_a02 max ≈ 5.5×10¹²
num_a04 max ≈ 1.3×10¹⁶, mean ≈ 1.1×10¹⁰
We gain several orders of magnitude compared to 00:00. This is clearly suspicious.
SOA (glysoa*):
The aerosol and surface modes (a02, a03, a04, sfc) blow up:
glysoa_oh_a04 max ≈ 1.0×10¹², mean ≈ 8.7×10⁵
glysoa_sfc_a04 max ≈ 4.4×10⁸
glysoa_nh4_a04 max ≈ 2.2×10⁹
These values are completely unrealistic for aerosol/trace gas concentrations.
File 02:00:00 :
Here, it is a total blow-up:
num_a03 max ≈ 1.7×10¹⁴
num_a04 max ≈ 2.3×10¹⁷
glysoa_r1_a03 max ≈ 1.4×10⁹
glysoa_r1_a04 max ≈ 8.0×10¹⁰
glysoa_sfc_a04 max ≈ 7.0×10⁹
glysoa_nh4_a04 max ≈ 6.4×10¹³
glysoa_oh_a04 max ≈ 5.6×10¹³
The mean values themselves also become huge (for example, glysoa_r1_a04 mean ≈ 2.9×10⁴).
However, when I plot the time series for my given point, the behavior still looks consistent.
I am also attaching plots for mode 4, where the problem is mainly located.
Thank you in advance to anyone who can provide further insight to help me resolve these issues.
PS :
In the wrfbiochemi_d01 file there are special characters, so instead of attaching it, here is how I generated it in essence:
We prepare the input for MEGAN:
cat > ~/WRFChem/megan_bio_emiss/megan_bio.inp << 'EOF'
&control
domains = 1,
start_lai_mnth = 7,
end_lai_mnth = 7,
wrf_dir = '/home/kamo/WRFChem/WRF/run',
megan_dir = '/home/kamo/WRFChem/WRF-Chem-data/MEGAN'
/
EOF
We run MEGAN:
cd ~/WRFChem/megan_bio_emiss
./megan_bio_emiss < megan_bio.inp 2>&1 | tee run_megan.log
We move the file into the WRF/run directory:
mv ~/WRFChem/megan_bio_emiss/wrfbiochemi_d01 ~/WRFChem/WRF/run/
Similarly, here is how I generate the EDGAR emissions:
We link all the objects into an executable:
$FC -o anthro_emis *.o $NETCDF_LIB
ls -lh anthro_emis
We copy the example input file:
cd ~/WRFChem/ANTHRO/src
cp ~/WRFChem/EDGAR-HTAP/INP-Examples/MOZART-MOSAIC.inp anthro_edgar.inp
Then we open the file and edit it:
nano anthro_edgar.inp
We go in and manually modify:
anthro_dir = '/home/kamo/WRFChem/EDGAR-HTAP/MOZART_MOSAIC'
wrf_dir = '/home/kamo/WRFChem/WRF/run'
We run:
export OMP_NUM_THREADS=1
export HDF5_USE_FILE_LOCKING=FALSE
ulimit -s unlimited
./anthro_emis < anthro_edgar.inp | tee run_edgar.log
We move our files into the run directory:
cd ~/WRFChem/ANTHRO/src
cp wrfchemi_00z_d01 wrfchemi_12z_d01 /home/kamo/WRFChem/WRF/run
We create dated links that will be used during the simulation:
ln -sf wrfchemi_00z_d01 wrfchemi_d01_2021-07-01_00:00:00
ln -sf wrfchemi_12z_d01 wrfchemi_d01_2021-07-01_12:00:00
For mozbc, I adapt the mapping to my specific case.
I am also attaching the namelist.output file, although I am not sure if it will be useful.
If you need any further details or even specific checks, please let me know.
Best regards,
Kamo
LATMOS
I am currently experiencing some difficulties with WRF-Chem(chem_opt=202). I am using WRF v4.5.2 and WPS v4.5.
I first ran a 1 hour simulation without any problem, then extended it to 3 hours, again without any issues. However, when I extended it to a 24 hour run, I encountered a problem. I would like to emphasize that for each case I adapted and downloaded the appropriate input data.
The message is going to be long, but it will be very detailed and complete in order to clearly explain the situation and the origin of the problem.
We are using a CESM file as GFS chemical boundary data for ungrib.
The global file was obtained here:
https://www.acom.ucar.edu/cesm/subset.shtml
We filled in the following information for our project:
Full name: Kamo
Affiliation: LATMOS
Address:
Save contact: on
lat2: 52.3
lon1: 356
lon2: 10
lat1: 45
Start year: 2021
Start month: 07
Start day: 01
End year: 2021
End month: 07
End day: 02
simulation: cesm
I am also using EDGAR+anthro, MEGAN, and FINN fires. I am attaching all these files for clarity.
I have already seen this problem mentioned in different threads, but no clear answer or fix has been provided, so I now find myself in the same situation as the others.
My rsl.error and rsl.out files are very large because of the extensive debugging, so I am only including the relevant part related to the problem.
Here is the context of the error message:
WARNING: Large total sw optical depth of 48.79 at point i,j,nb= 113 46 1
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1 8.88 13.52 15.18 8.33 12.05
2 22.99 34.20 39.99 23.81 36.41
3 0.16 0.21 0.24 0.14 0.18
4 0.09 0.14 0.17 0.10 0.15
5 0.00 0.00 0.00 0.00 0.00
6 0.00 0.00 0.00 0.00 0.00
7 0.00 0.00 0.00 0.00 0.00
8 0.00 0.00 0.00 0.00 0.00
9 0.00 0.00 0.00 0.00 0.00
10 0.00 0.00 0.00 0.00 0.00
11 0.00 0.00 0.00 0.00 0.00
12 0.00 0.00 0.00 0.00 0.00
13 0.00 0.00 0.00 0.00 0.00
14 0.00 0.00 0.00 0.00 0.00
15 0.00 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00 0.00
25 0.00 0.00 0.00 0.00 0.00
26 0.00 0.00 0.00 0.00 0.00
27 0.00 0.00 0.00 0.00 0.00
28 0.00 0.00 0.00 0.00 0.00
29 0.00 0.00 0.00 0.00 0.00
30 0.00 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00 0.00
35 0.00 0.00 0.00 0.00 0.00
36 0.00 0.00 0.00 0.00 0.00
37 0.00 0.00 0.00 0.00 0.00
38 0.00 0.00 0.00 0.00 0.00
39 0.00 0.00 0.00 0.00 0.00
40 0.00 0.00 0.00 0.00 0.00
41 0.00 0.00 0.00 0.00 0.00
42 0.00 0.00 0.00 0.00 0.00
43 0.00 0.00 0.00 0.00 0.00
44 0.00 0.00 0.00 0.00 0.00
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.62 0.76 0.79 0.68
2 0.66 0.73 0.80 0.68
3 0.69 0.75 0.79 0.69
4 0.65 0.73 0.79 0.67
5 0.73 0.73 0.72 0.58
6 0.68 0.71 0.74 0.62
7 0.69 0.70 0.72 0.62
8 0.69 0.70 0.72 0.63
9 0.67 0.71 0.73 0.64
10 0.64 0.71 0.77 0.65
11 0.65 0.71 0.77 0.64
12 0.64 0.71 0.76 0.63
13 0.64 0.72 0.76 0.64
14 0.64 0.74 0.78 0.67
15 0.66 0.74 0.78 0.71
16 0.66 0.74 0.77 0.70
17 0.66 0.72 0.77 0.71
18 0.67 0.74 0.78 0.74
19 0.65 0.72 0.78 0.71
20 0.55 0.61 0.69 0.60
21 0.66 0.73 0.77 0.66
22 0.66 0.72 0.77 0.65
23 0.66 0.74 0.78 0.66
24 0.66 0.74 0.77 0.66
25 0.64 0.74 0.77 0.68
26 0.66 0.73 0.77 0.66
27 0.63 0.74 0.77 0.67
28 0.63 0.72 0.79 0.72
29 0.56 0.62 0.69 0.60
30 0.56 0.62 0.69 0.60
31 0.59 0.58 0.67 0.59
32 0.59 0.50 0.61 0.57
33 0.61 0.40 0.39 0.51
34 0.61 0.40 0.39 0.51
35 0.61 0.40 0.39 0.51
36 0.61 0.40 0.39 0.51
37 0.61 0.40 0.39 0.50
38 0.61 0.40 0.39 0.50
39 0.61 0.40 0.39 0.50
40 0.61 0.41 0.40 0.51
41 0.61 0.40 0.39 0.51
42 0.61 0.41 0.39 0.51
43 0.61 0.41 0.40 0.51
44 0.61 0.41 0.40 0.51
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 0.91 0.94 0.96 0.95
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 0.99 0.99 1.00 1.00
29 0.91 0.93 0.96 0.95
30 0.91 0.94 0.96 0.95
31 0.92 0.93 0.95 0.93
32 0.94 0.94 0.94 0.93
33 0.95 0.93 0.88 0.85
34 0.95 0.93 0.88 0.85
35 0.95 0.93 0.88 0.85
36 0.95 0.93 0.88 0.84
37 0.95 0.93 0.88 0.85
38 0.95 0.93 0.88 0.84
39 0.95 0.93 0.88 0.85
40 0.95 0.93 0.88 0.85
41 0.95 0.93 0.88 0.85
42 0.95 0.93 0.88 0.85
43 0.95 0.93 0.88 0.85
44 0.95 0.93 0.88 0.85
Diagnostics 4: k, ssaal, asyal, taual
0 1.00 1.00 0.00
1 1.00 1.00 0.00
2 1.00 1.00 0.00
3 1.00 1.00 0.00
4 1.00 1.00 0.00
5 1.00 1.00 0.00
6 1.00 0.00 0.00
7 1.00 0.00 0.00
8 1.00 0.00 0.00
9 1.00 0.00 0.00
10 1.00 0.00 0.00
11 1.00 0.00 0.00
12 1.00 0.00 0.00
13 1.00 0.00 0.00
14 1.00 0.00 0.00
15 1.00 0.00 0.00
16 1.00 0.00 0.00
17 1.00 0.00 0.00
18 1.00 0.00 0.00
19 1.00 0.00 0.00
20 1.00 0.00 0.00
21 1.00 0.00 0.00
22 1.00 0.00 0.00
23 1.00 0.00 0.00
24 1.00 1.00 12.39
25 1.00 1.00 35.86
26 1.00 1.00 0.19
27 1.00 1.00 0.15
28 1.00 0.67 0.00
29 1.00 0.95 0.00
30 1.00 0.85 0.00
31 1.00 0.84 0.00
32 1.00 0.87 0.00
33 1.00 1.00 0.00
34 1.00 1.00 0.00
35 1.00 1.00 0.00
36 1.00 1.00 0.00
37 1.00 1.00 0.00
38 1.00 1.00 0.00
39 1.00 1.00 0.00
40 1.00 1.00 0.00
41 1.00 1.00 0.00
42 1.00 1.00 0.00
43 1.00 0.00 0.00
44 1.00 1.00 0.00
I also see this at different points. Then, at the end, I get:
WARNING: Large total sw optical depth of 25.91 at point i,j,nb= 117 77 14
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999, tauaer
1 5.83 6.02 6.51 6.37 7.53
2 7.51 7.56 7.97 8.13 9.22
3 6.99 6.99 6.61 7.15 7.40
4 1.94 2.05 1.93 1.79 1.68
5 0.10 0.09 0.09 0.09 0.08
6 0.01 0.01 0.00 0.00 0.00
7 0.00 0.00 0.00 0.00 0.00
8 0.00 0.00 0.00 0.00 0.00
9 0.00 0.00 0.00 0.00 0.00
10 0.00 0.00 0.00 0.00 0.00
11 0.00 0.00 0.00 0.00 0.00
12 0.00 0.00 0.00 0.00 0.00
13 0.00 0.00 0.00 0.00 0.00
14 0.00 0.00 0.00 0.00 0.00
15 0.00 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00 0.00
25 0.00 0.00 0.00 0.00 0.00
26 0.00 0.00 0.00 0.00 0.00
27 0.00 0.00 0.00 0.00 0.00
28 0.00 0.00 0.00 0.00 0.00
29 0.00 0.00 0.00 0.00 0.00
30 0.00 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00 0.00
35 0.00 0.00 0.00 0.00 0.00
36 0.00 0.00 0.00 0.00 0.00
37 0.00 0.00 0.00 0.00 0.00
38 0.00 0.00 0.00 0.00 0.00
39 0.00 0.00 0.00 0.00 0.00
40 0.00 0.00 0.00 0.00 0.00
41 0.00 0.00 0.00 0.00 0.00
42 0.00 0.00 0.00 0.00 0.00
43 0.00 0.00 0.00 0.00 0.00
44 0.00 0.00 0.00 0.00 0.00
Diagnostics 2: k, gaer300, gaer400, gaer600, gaer999
1 0.84 0.82 0.82 0.83
2 0.85 0.83 0.81 0.81
3 0.83 0.84 0.83 0.80
4 0.83 0.82 0.82 0.80
5 0.84 0.83 0.84 0.82
6 0.81 0.79 0.78 0.78
7 0.74 0.72 0.73 0.62
8 0.65 0.72 0.77 0.71
9 0.65 0.72 0.78 0.71
10 0.66 0.73 0.77 0.69
11 0.65 0.74 0.78 0.69
12 0.66 0.73 0.78 0.70
13 0.66 0.73 0.78 0.71
14 0.66 0.73 0.78 0.71
15 0.64 0.74 0.78 0.73
16 0.68 0.70 0.79 0.77
17 0.72 0.64 0.62 0.73
18 0.70 0.73 0.79 0.76
19 0.66 0.73 0.79 0.76
20 0.66 0.74 0.77 0.70
21 0.66 0.74 0.78 0.69
22 0.66 0.74 0.77 0.66
23 0.67 0.73 0.77 0.64
24 0.64 0.71 0.77 0.64
25 0.65 0.74 0.78 0.69
26 0.65 0.73 0.78 0.71
27 0.65 0.73 0.78 0.74
28 0.65 0.73 0.78 0.74
29 0.56 0.62 0.69 0.60
30 0.56 0.62 0.69 0.60
31 0.57 0.62 0.69 0.60
32 0.58 0.57 0.66 0.59
33 0.60 0.39 0.37 0.49
34 0.61 0.40 0.38 0.50
35 0.61 0.40 0.39 0.51
36 0.61 0.40 0.39 0.50
37 0.61 0.40 0.39 0.51
38 0.61 0.40 0.39 0.51
39 0.61 0.40 0.39 0.51
40 0.61 0.41 0.39 0.51
41 0.61 0.41 0.40 0.51
42 0.61 0.41 0.40 0.51
43 0.61 0.41 0.40 0.51
44 0.61 0.41 0.40 0.51
Diagnostics 3: k, waer300, waer400, waer600, waer999
1 1.00 1.00 1.00 1.00
2 1.00 1.00 1.00 1.00
3 1.00 1.00 1.00 1.00
4 1.00 1.00 1.00 1.00
5 1.00 1.00 1.00 1.00
6 1.00 1.00 1.00 1.00
7 1.00 1.00 1.00 1.00
8 1.00 1.00 1.00 1.00
9 1.00 1.00 1.00 1.00
10 1.00 1.00 1.00 1.00
11 1.00 1.00 1.00 1.00
12 1.00 1.00 1.00 1.00
13 1.00 1.00 1.00 1.00
14 1.00 1.00 1.00 1.00
15 1.00 1.00 1.00 1.00
16 1.00 1.00 1.00 1.00
17 1.00 1.00 1.00 1.00
18 1.00 1.00 1.00 1.00
19 1.00 1.00 1.00 1.00
20 1.00 1.00 1.00 1.00
21 1.00 1.00 1.00 1.00
22 1.00 1.00 1.00 1.00
23 1.00 1.00 1.00 1.00
24 1.00 1.00 1.00 1.00
25 1.00 1.00 1.00 1.00
26 1.00 1.00 1.00 1.00
27 1.00 1.00 1.00 1.00
28 1.00 1.00 1.00 1.00
29 0.91 0.93 0.96 0.95
30 0.91 0.93 0.96 0.95
31 0.91 0.93 0.96 0.95
32 0.93 0.94 0.95 0.94
33 0.95 0.94 0.90 0.85
34 0.95 0.93 0.88 0.84
35 0.95 0.93 0.88 0.85
36 0.95 0.93 0.88 0.84
37 0.95 0.93 0.88 0.85
38 0.95 0.93 0.88 0.85
39 0.95 0.93 0.88 0.85
40 0.95 0.93 0.88 0.85
41 0.95 0.93 0.88 0.85
42 0.95 0.93 0.88 0.85
43 0.95 0.93 0.88 0.85
44 0.95 0.93 0.88 0.85
Diagnostics 4: k, ssaal, asyal, taual
0 0.00 NaN 0.00
1 0.00 NaN 0.00
2 0.00 0.00 0.00
3 0.00 0.00 0.00
4 0.00 0.00 0.00
5 0.00 0.00 0.00
6 0.00 0.00 0.00
7 0.00 0.00 0.00
8 0.00 0.05 0.00
9 0.00 0.00 0.00
10 0.00 0.02 0.00
11 0.00 0.00 0.00
12 0.00 0.00 0.00
13 0.00 0.00 0.00
14******** 0.00 0.00
15 0.00 0.00 0.00
16 13.78******** 0.00
17 14.13 0.00 0.00
18 15.54 0.00 0.00
19 19.00 0.00 0.00
20 19.34 0.00 0.00
21 19.33 0.00 0.00
22 18.70 0.00 0.01
23 20.38 0.00 0.00
24 39.02 0.00 0.00
25 47.78 0.00 0.00
26 51.14 0.00 0.00
27 54.01 0.00 0.00
28 55.83 0.00 0.00
29 99.42 0.00 0.00
30 114.35 10.12 0.00
31 130.16 0.00 0.00
32 145.19 0.00 0.00
33 158.71 0.00 0.00
34 171.34 0.00 0.00
35 182.15 0.00 0.00
36 190.74 0.00 0.00
37 197.49 0.00 0.00
38 202.89 0.00 0.00
39 206.31 0.00 0.00
40 209.10******** 0.00
41 214.38 0.00 0.00
42 219.37 0.00 0.00
43 220.22 0.00 0.00
44 221.02 0.00 0.00
-------------------------
d01 2021-07-01_03:45:00 calling inc/HALO_PWP_inline.inc
d01 2021-07-01_03:45:00 call surface_driver
d01 2021-07-01_03:45:00 in SFCLAY
d01 2021-07-01_03:45:00 in NOAH DRV
d01 2021-07-01_03:45:00 call pbl_driver
d01 2021-07-01_03:45:00 in MYNNPBL
d01 2021-07-01_03:45:00 call HALO CHEM AFTER PBL
d01 2021-07-01_03:45:00 calling inc/HALO_EM_CHEM_E_5_inline.inc
d01 2021-07-01_03:45:00 call cumulus_driver
d01 2021-07-01_03:45:00 calling inc/HALO_CUP_G3_IN_inline.inc
d01 2021-07-01_03:45:00 in grelldrv
d01 2021-07-01_03:45:00 returning from cumulus_driver
d01 2021-07-01_03:45:00 call shallow_cumulus_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_FDDA_SFC_inline.inc
d01 2021-07-01_03:45:00 call fddagd_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SCALAR_E_5_inline.inc
d01 2021-07-01_03:45:00 call calculate_phy_tend
d01 2021-07-01_03:45:00 call compute_diff_metrics
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_A1_inline.inc
d01 2021-07-01_03:45:00 call bc for diffusion_metrics
d01 2021-07-01_03:45:00 call cal_deform_and_div
d01 2021-07-01_03:45:00 calling inc/HALO_EM_HELICITY_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_D_inline.inc
d01 2021-07-01_03:45:00 call calculate_km_kh
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_E_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_PHY_BC_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_CHEM_inline.inc
d01 2021-07-01_03:45:00 call phy_bc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_PBL_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_DIFFUSION_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_TKE_5_inline.inc
d01 2021-07-01_03:45:00 call update_phy_ten
d01 2021-07-01_03:45:00 call horizontal_diffusion_2
d01 2021-07-01_03:45:00 call rk_tendency
d01 2021-07-01_03:45:00 call small_step_prep
d01 2021-07-01_03:45:00 calling inc/HALO_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C2_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_B3_inline.inc
d01 2021-07-01_03:45:00 call rk_small_finish
d01 2021-07-01_03:45:00 calling inc/HALO_EM_D_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_EM_DA_inline.inc
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call rk_scalar_tend in chem_tile_loop_1
d01 2021-07-01_03:45:00 call rk_update_scalar
d01 2021-07-01_03:45:00 call advance_ppt
d01 2021-07-01_03:45:00 call phy_prep_part2
d01 2021-07-01_03:45:00 call moist_physics_prep
d01 2021-07-01_03:45:00 call microphysics_driver
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SBM_inline.inc
d01 2021-07-01_03:45:00 microphysics_driver: calling morrison two moment
d01 2021-07-01_03:45:00 microphysics_driver: returning from
d01 2021-07-01_03:45:00 call moist_physics_finish
d01 2021-07-01_03:45:00 calling inc/HALO_EM_THETAM_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_EM_THETAM_inline.inc
d01 2021-07-01_03:45:00 call chem polar filter
d01 2021-07-01_03:45:00 calling inc/HALO_EM_D3_5_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_D3_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_MOIST_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_CHEM_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_TRACER_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_SCALAR_inline.inc
d01 2021-07-01_03:45:00 call set_phys_bc_dry_2
d01 2021-07-01_03:45:00 call spec_bdy_final
d01 2021-07-01_03:45:00 calling inc/HALO_EM_C_inline.inc
d01 2021-07-01_03:45:00 calling inc/PERIOD_BDY_EM_E_inline.inc
d01 2021-07-01_03:45:00 call set_w_surface
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_W_inline.inc
d01 2021-07-01_03:45:00 calling inc/HALO_EM_PHYS_HCW_inline.inc
d01 2021-07-01_03:45:00 --> TOP OF AFTER ALL RK STEPS
d01 2021-07-01_03:45:00 --> CALLING DIAGNOSTICS DRIVER
d01 2021-07-01_03:45:00 --> TOP OF DIAGNOSTICS PACKAGE
d01 2021-07-01_03:45:00 --> CALL DIAGNOSTICS PACKAGE: ACCUMULATED AND BUCKET DIAGNOSTICS
d01 2021-07-01_03:45:00 call HALO_RK_E
d01 2021-07-01_03:45:00 calling inc/HALO_EM_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_MOIST
d01 2021-07-01_03:45:00 calling inc/HALO_EM_MOIST_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_CHEM
d01 2021-07-01_03:45:00 calling inc/HALO_EM_CHEM_E_5_inline.inc
d01 2021-07-01_03:45:00 call HALO_RK_SCALAR
d01 2021-07-01_03:45:00 calling inc/HALO_EM_SCALAR_E_5_inline.inc
d01 2021-07-01_03:45:00 call end of solve_em
d01 2021-07-01_03:45:00 DEBUG wrf_timetoa(): returning with str = [2021-07-01_03:45:00]
d01 2021-07-01_03:45:00 calling mozart_mosaic_4bin_aq_kpp from chem_driver
d01 2021-07-01_03:45:00 calling chem_opt=? from chem_driver
d01 2021-07-01_03:45:00 call chem_prep
d01 2021-07-01_03:45:00 calling emissions driver
d01 2021-07-01_03:45:00 No tracer option selected
d01 2021-07-01_03:45:00 MOSAIC or SORGAM sea salt emissions
d01 2021-07-01_03:45:00 Gocart dust emissions
d01 2021-07-01_03:45:00 no dms emissions
d01 2021-07-01_03:45:00 fire emissions: adding biomassb emissions
d01 2021-07-01_03:45:00 nothing done with burn emissions for tracers here
d01 2021-07-01_03:45:00 biogenic emissions: calling megan v2.04
calculate MEGAN emissions at ktau, gmtp, tmidh = 901 3.00000000 3.75208330
par out of range at i=149 ,j= 6 par = -2.214
d01 2021-07-01_03:45:00 emissions_driver calling add_anthropogenics
d01 2021-07-01_03:45:00 add_anthropogenics
d01 2021-07-01_03:45:00 emissions_driver calling add_biogenics
d01 2021-07-01_03:45:00 emissions_driver not calling add_emiss_cptec
d01 2021-07-01_03:45:00 emissions_driver calling mosaic_addemiss
d01 2021-07-01_03:45:00 mosaic_addemiss: emiss_opt = eptec
d01 2021-07-01_03:45:00 mosaic_addemiss: gocart speciation being mapped to mosaic
d01 2021-07-01_03:45:00 mosaic_addemiss: adding fire emissions to MOSAIC
d01 2021-07-01_03:45:00 calling optical driver
d01 2021-07-01_03:45:00 optical driver: process aerosols true
d01 2021-07-01_03:45:00 optical averaging
d01 2021-07-01_03:45:00 optical driver: call optical averaging
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 16 3 29 4 0.7404E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 17 3 29 4 0.8284E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 18 3 29 4 0.8648E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 19 3 29 4 0.8779E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 3 5 4 0.7080E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 4 29 4 0.1978E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 4 5 4 0.7746E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 7 28 4 0.1094E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 26 8 2 4 0.7275E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 27 10 3 4 0.7014E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 20 12 7 4 0.8006E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 24 14 5 4 0.8569E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 4 16 24 4 0.1029E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 18 28 4 0.1344E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 19 28 4 0.1506E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 20 28 4 0.1459E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 21 28 4 0.8573E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 22 15 4 0.4402E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 23 29 4 0.1691E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 24 29 4 0.1435E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 25 29 4 0.1943E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 26 29 3 0.7767E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 26 29 4 0.2388E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 109 26 1 3 0.7794E+00
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 27 29 3 0.1015E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 27 29 4 0.3145E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 14 27 24 4 0.1712E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 18 4 0.2673E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 29 3 0.9201E-03
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 30 29 4 0.2595E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 31 18 4 0.2107E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 31 29 4 0.1959E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 32 29 4 0.1464E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 33 29 4 0.1458E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 3 34 29 4 0.1309E-02
mieaer_sc: radius_wet set to rmax,id,i,j,k,m,rm(m,k) 1 35 34 28 4 0.1809E-02
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0 0xfe45244251f in ???
#1 0x57b90986cdb3 in ???
#2 0x57b909871fb3 in ???
#3 0x57b909872b64 in ???
#4 0x57b909872dc6 in ???
#5 0x57b9098737e3 in ???
#6 0x57b9098aa680 in ???
#7 0x57b908aefb1a in ???
#8 0x57b9089ee84d in ???
#9 0x57b9087e9ae9 in ???
#10 0x57b906d35fe8 in ???
#11 0x57b906c98be7 in ???
#12 0x57b906c9809e in ???
#13 0xfe452429d8f in ???
#14 0xfe452429e3f in ???
#15 0x57b906c980d4 in ???
#16 0xffffffffffffffff in ???
I have checked my input files and there does not seem to be any problem; the values look consistent. The issue appears during the run itself, and more specifically at this model time. I only noticed that in mozbc I mapped to 8 bins instead of 4 bins, but this is not the cause of the crash.
I also examined my output file at 03:00:00, and here is what I found:
Sulfates (so4_a01 to so4_a04):
Maximum values remain between ~0.25 and ~2.35, with low mean values (~0.004–0.08).
No abnormal behavior is observed.
Sulfates do not seem to be responsible for the large optical depths.
Aerosol number (num_a01 to num_a04):
num_a01 to num_a03 already show high values (up to ~10¹³).
num_a04 is particularly problematic:
max ≈ 1.1 × 10¹⁷
mean ≈ 5.4 × 10¹¹
These orders of magnitude are unrealistic for aerosol number and indicate strong instability or inconsistency in this mode.
SOA (glysoa_*):
For the fine mode (a01), values remain low and consistent.
From a02 onwards, and especially a03 and a04, the maxima become extremely large. For example:
glysoa_sfc_a03 max ≈ 2.0 × 10⁶
glysoa_nh4_a03 max ≈ 2.9 × 10⁷
glysoa_oh_a03 max ≈ 5.1 × 10⁷
glysoa_r1_a04 max ≈ 7.4 × 10⁹
glysoa_nh4_a04 max ≈ 3.1 × 10¹³
glysoa_oh_a04 max ≈ 2.7 × 10¹³
Several of these fields also have very high mean values (up to ~10⁸).
The aqueous terms (glysoa_*_cw0*) also show abnormally large values (max up to ~10⁶–10⁷).
I then examined the outputs at 00:00:00, 01:00:00 and 02:00:00, and here is what I observe:
wrfout_d01_2021-07-01_00:00:00
Sulfates (so4_a0*):
Reasonable values (max ~2.4 µg/m³, mean ~0.08–0.12).
A few negative values (e.g. so4_a03 min = -0.1559), which are nonphysical, but with moderate amplitude. This is already a small numerical issue but not catastrophic.
Aerosol number (num_a0*):
Again, some negative values (e.g. num_a01 min = -1.7×10⁷).
Maxima remain of the usual order of magnitude (10⁹–10¹⁰ for num_a0*).
This is still “manageable”: sign errors, but not an explosion yet.
SOA (glysoa*):
Maxima remain reasonable (often << 10³, sometimes up to 4×10⁴ for glysoa_oh_a03).
No huge values and no absurdly large means.
File 01:00:00 :
Here the real drift begins:
num_a0*:
num_a02 max ≈ 5.5×10¹²
num_a04 max ≈ 1.3×10¹⁶, mean ≈ 1.1×10¹⁰
We gain several orders of magnitude compared to 00:00. This is clearly suspicious.
SOA (glysoa*):
The aerosol and surface modes (a02, a03, a04, sfc) blow up:
glysoa_oh_a04 max ≈ 1.0×10¹², mean ≈ 8.7×10⁵
glysoa_sfc_a04 max ≈ 4.4×10⁸
glysoa_nh4_a04 max ≈ 2.2×10⁹
These values are completely unrealistic for aerosol/trace gas concentrations.
File 02:00:00 :
Here, it is a total blow-up:
num_a03 max ≈ 1.7×10¹⁴
num_a04 max ≈ 2.3×10¹⁷
glysoa_r1_a03 max ≈ 1.4×10⁹
glysoa_r1_a04 max ≈ 8.0×10¹⁰
glysoa_sfc_a04 max ≈ 7.0×10⁹
glysoa_nh4_a04 max ≈ 6.4×10¹³
glysoa_oh_a04 max ≈ 5.6×10¹³
The mean values themselves also become huge (for example, glysoa_r1_a04 mean ≈ 2.9×10⁴).
However, when I plot the time series for my given point, the behavior still looks consistent.
I am also attaching plots for mode 4, where the problem is mainly located.
Thank you in advance to anyone who can provide further insight to help me resolve these issues.
PS :
In the wrfbiochemi_d01 file there are special characters, so instead of attaching it, here is how I generated it in essence:
We prepare the input for MEGAN:
cat > ~/WRFChem/megan_bio_emiss/megan_bio.inp << 'EOF'
&control
domains = 1,
start_lai_mnth = 7,
end_lai_mnth = 7,
wrf_dir = '/home/kamo/WRFChem/WRF/run',
megan_dir = '/home/kamo/WRFChem/WRF-Chem-data/MEGAN'
/
EOF
We run MEGAN:
cd ~/WRFChem/megan_bio_emiss
./megan_bio_emiss < megan_bio.inp 2>&1 | tee run_megan.log
We move the file into the WRF/run directory:
mv ~/WRFChem/megan_bio_emiss/wrfbiochemi_d01 ~/WRFChem/WRF/run/
Similarly, here is how I generate the EDGAR emissions:
We link all the objects into an executable:
$FC -o anthro_emis *.o $NETCDF_LIB
ls -lh anthro_emis
We copy the example input file:
cd ~/WRFChem/ANTHRO/src
cp ~/WRFChem/EDGAR-HTAP/INP-Examples/MOZART-MOSAIC.inp anthro_edgar.inp
Then we open the file and edit it:
nano anthro_edgar.inp
We go in and manually modify:
anthro_dir = '/home/kamo/WRFChem/EDGAR-HTAP/MOZART_MOSAIC'
wrf_dir = '/home/kamo/WRFChem/WRF/run'
We run:
export OMP_NUM_THREADS=1
export HDF5_USE_FILE_LOCKING=FALSE
ulimit -s unlimited
./anthro_emis < anthro_edgar.inp | tee run_edgar.log
We move our files into the run directory:
cd ~/WRFChem/ANTHRO/src
cp wrfchemi_00z_d01 wrfchemi_12z_d01 /home/kamo/WRFChem/WRF/run
We create dated links that will be used during the simulation:
ln -sf wrfchemi_00z_d01 wrfchemi_d01_2021-07-01_00:00:00
ln -sf wrfchemi_12z_d01 wrfchemi_d01_2021-07-01_12:00:00
For mozbc, I adapt the mapping to my specific case.
I am also attaching the namelist.output file, although I am not sure if it will be useful.
If you need any further details or even specific checks, please let me know.
Best regards,
Kamo
LATMOS
Reply all
Reply to author
Forward
0 new messages