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

Error in running WRF-Chem using RADM2-MADE/SORGAM after making boundary conditions for nested domain using ndown

Ankan

Member
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
 

Attachments

  • namelist.input_coarsegridrun.txt
    8.4 KB · Views: 9
  • namelist.input_ndownrun.txt
    8.4 KB · Views: 3
  • namelist.input_nestgridrun.txt
    8.3 KB · Views: 6
Last edited:
Dear all,
I gave a fresh run again but this time I set
emiss_opt = 0,
bio_emiss_opt = 0,
biomass_burn_opt = 0,
seas_opt = 0,
dust_opt = 0,
dmsemis_opt = 0,
But still, the same errors and warnings are coming, and the simulation is stopped at 00:20:20, just like before. So, I don't think this could be because of these aerosol species. There is something else wrong. But I am struggling to find the root cause of these issues.
Is there anyone who can guide me to solving this issue? That will be really helpful for me. Thank you in advance for your time and consideration.
With regards,
Ankan
 
Last edited:
Dear everyone,
After giving so many test runs, I learned that only by setting 'aer_op_opt = 0' did all the above warning messages and errors vanish. And, the simulation is also going on. I also tried by 'aer_op_opt' to 1, 2, 3, 4, and 5. But none of this worked.
I was following the tutorial exercise https://ruc.noaa.gov/wrf/wrf-chem/tutorialexercises/tutorialexercises009.html.tmp, where it was mentioned that to calculate aerosol direct effect following lines should be considered:
chem_opt=2
progn=0
aer_ra_feedback = 1
aer_op_opt = 1
So, my question is, is it okay if I consider aerosol direct effects by setting 'aer_op_opt = 0' while running WRF-Chem? Can anyone please guide me in this regard? That will be greatly appreciated. Thank you for your time and consideration.
With regards,
Ankan
 
Ankan,

Setting aer_op_opt = 0 tells the model not to calculate optical properties and thus aerosol direct effects will not be included. Something is wrong with your aerosol simulation.

Jordan
 
Hi Ankan,

Were you able to figure out a solution to successfully run the simulation with aer_op_opt > 0?
I am facing a similar issue and any input will be of great help.

Thanks
Krishna
 
Top