Hello everyone,
I gave a WRF-Chem run of the direct effect using the RADM2-MADE/SORGAM chemistry option. The thing is, the simulation of the parent domain was completed successfully without throwing any errors or warnings. But when I gave a simulation for the nested domain after including boundary conditions using ndown, the simulation was stopped at 20 minutes after printing some warnings in the rsl files (like in the following):
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 5 8 0.5000E-02 0.6806E-02
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 3 0.5000E-02 0.9308E+02
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 4 0.5000E-02 0.7421E+01
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 5 0.5000E-02 0.1127E+01
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 6 0.5000E-02 0.4981E+00
.......
And,
WARNING: Large total optical depth of 89039.00 at point i,j,ib= 102 60 11
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999
122935.4923314.1223546.7923362.27
2 4516.83 4595.96 4645.28 4605.27
3 36.71 37.17 37.76 37.78
4 32.06 32.48 32.99 32.96
5 26.29 26.62 27.05 27.05
6 4219.61 4206.24 4293.73 4398.62
726218.8626825.8926758.4227027.11
8 439.26 490.87 471.42 328.93
9 519.29 559.88 550.55 394.78
10 7258.15 7339.81 7744.52 8109.93
11 6083.74 6106.64 6368.91 6787.10
12 1380.64 1422.42 1452.39 1522.76
13 7262.47 7771.11 7491.30 7544.24
14 701.68 713.84 724.75 740.10
15 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00
25 0.62 0.64 0.66 0.67
26 0.27 0.27 0.29 0.30
27 0.04 0.04 0.05 0.05
28 0.05 0.04 0.05 0.05
29 0.01 0.01 0.01 0.01
30 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00
...
If you see the above warnings, you will notice the total optical depth is highly unrealistic. I didn't understand why these two different types of warnings were coming in the case of nested domain simulation while parent domain simulation was completed successfully without any errors. I searched online for this kind of error but didn't get any solution. Here, I want to mention that when I was giving parent domain simulation, I set seas_opt=2 instead of 1 since only by giving this option did the large total optical depth (<10, not beyond that value) warnings vanish. I noticed in the WRF-Chem user guide that for RADM2-MADE/SORGAM with aqueous chemistry, i.e., chem_opt=11 'Due to errors, dust_opt=2, seas_opt=2' has been disabled.-- is written (but not beside chem_opt=2). So, in the case of Indirect effect simulation, I used seas_opt=1, and the simulations were completed successfully for both domains (20 km and 4 km). However, only for some days, I got some large total optical depth warnings and therefore I ignored those warnings. But, in this case, I didn't understand why this kind of error is coming? Can anyone please guide me on why are the errors coming in this case and how to solve them? Could it be because of seas_opt? If so, then which seas_opt can be used? I am really looking forward to hearing from you. Thank you for your time and consideration.
Below the namelist.input files for coarsegrid run, ndown run and nestgrid run are attached.
With regards,
Ankan
I gave a WRF-Chem run of the direct effect using the RADM2-MADE/SORGAM chemistry option. The thing is, the simulation of the parent domain was completed successfully without throwing any errors or warnings. But when I gave a simulation for the nested domain after including boundary conditions using ndown, the simulation was stopped at 20 minutes after printing some warnings in the rsl files (like in the following):
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 5 8 0.5000E-02 0.6806E-02
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 3 0.5000E-02 0.9308E+02
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 4 0.5000E-02 0.7421E+01
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 5 0.5000E-02 0.1127E+01
mieaer: radius_wet set to rmax,id,i,j,k,m,rm(m,k),number 1 1 1 6 6 0.5000E-02 0.4981E+00
.......
And,
WARNING: Large total optical depth of 89039.00 at point i,j,ib= 102 60 11
Diagnostics 1: k, tauaer300, tauaer400, tauaer600, tauaer999
122935.4923314.1223546.7923362.27
2 4516.83 4595.96 4645.28 4605.27
3 36.71 37.17 37.76 37.78
4 32.06 32.48 32.99 32.96
5 26.29 26.62 27.05 27.05
6 4219.61 4206.24 4293.73 4398.62
726218.8626825.8926758.4227027.11
8 439.26 490.87 471.42 328.93
9 519.29 559.88 550.55 394.78
10 7258.15 7339.81 7744.52 8109.93
11 6083.74 6106.64 6368.91 6787.10
12 1380.64 1422.42 1452.39 1522.76
13 7262.47 7771.11 7491.30 7544.24
14 701.68 713.84 724.75 740.10
15 0.00 0.00 0.00 0.00
16 0.00 0.00 0.00 0.00
17 0.00 0.00 0.00 0.00
18 0.00 0.00 0.00 0.00
19 0.00 0.00 0.00 0.00
20 0.00 0.00 0.00 0.00
21 0.00 0.00 0.00 0.00
22 0.00 0.00 0.00 0.00
23 0.00 0.00 0.00 0.00
24 0.00 0.00 0.00 0.00
25 0.62 0.64 0.66 0.67
26 0.27 0.27 0.29 0.30
27 0.04 0.04 0.05 0.05
28 0.05 0.04 0.05 0.05
29 0.01 0.01 0.01 0.01
30 0.00 0.00 0.00 0.00
31 0.00 0.00 0.00 0.00
32 0.00 0.00 0.00 0.00
33 0.00 0.00 0.00 0.00
34 0.00 0.00 0.00 0.00
...
If you see the above warnings, you will notice the total optical depth is highly unrealistic. I didn't understand why these two different types of warnings were coming in the case of nested domain simulation while parent domain simulation was completed successfully without any errors. I searched online for this kind of error but didn't get any solution. Here, I want to mention that when I was giving parent domain simulation, I set seas_opt=2 instead of 1 since only by giving this option did the large total optical depth (<10, not beyond that value) warnings vanish. I noticed in the WRF-Chem user guide that for RADM2-MADE/SORGAM with aqueous chemistry, i.e., chem_opt=11 'Due to errors, dust_opt=2, seas_opt=2' has been disabled.-- is written (but not beside chem_opt=2). So, in the case of Indirect effect simulation, I used seas_opt=1, and the simulations were completed successfully for both domains (20 km and 4 km). However, only for some days, I got some large total optical depth warnings and therefore I ignored those warnings. But, in this case, I didn't understand why this kind of error is coming? Can anyone please guide me on why are the errors coming in this case and how to solve them? Could it be because of seas_opt? If so, then which seas_opt can be used? I am really looking forward to hearing from you. Thank you for your time and consideration.
Below the namelist.input files for coarsegrid run, ndown run and nestgrid run are attached.
With regards,
Ankan
Attachments
Last edited: