WRF-Chem predicted aerosol fields are all zero

[Aman Khan]

Hello Everyone,

I'm using the WRF-Chem v4.6 with MOZART-MOSIAC option (chem_opt=202). I have updated the IC/BC conditions from CAM-CHEM model using MOZBC utility. The model is initialized correctly and runs successfully for the complete duration of the simulation period, but the aerosol fields are all zero in the prediction. The rsl.out.0000 show the following warning at regular intervals:

leaving mosaic_aerchem_driver - ktau =    915
 wetscav_mozcart: percent columns with precip =    2.551020     %
Timing for main: time 2019-04-13_15:15:00 on domain   1:    4.39074 elapsed seconds
entering mosaic_cloudchem_driver - ktau =    916
 dvode--  maxord (=i1) .lt. 0
      in above message,  i1 =-546816000
 *** module_cmuaq_bulk, aubr aqintegr1 -- svode failed, istate =          -3
 *** mosaic_cloudchem_driver, subr interface_to_aqoperator1
    id,it,jt,kt, istat_fatal, warn =    1   14   14    8       -10         0
leaving  mosaic_cloudchem_driver - ktau =    916      1

Can anyone please help me identify the issue?

I have attached my namelist and rsl.out.0000 files. 

Thanks,

Aman

[Mary Barth]

Hi Aman,

It would be good to have your time step be a factor of the hourly time step of reading emissions. I suggest changing your time step to 360 s for domain 1. That's a 6-minute time step. Therefore, set chemist = 6  (or simply set it to 0 so that it uses the meteorology time step).

Check the other time steps in the physics and chem sections. With a 6 minute time step, radiation dt should be 12 or 18 minutes, photdt should also be 12 or 18 (match the radiation dt). 

Also, your namelist shows chem_opt = 112 (not 202). Option 112 is good for getting a simulation working for the first time -- and cheaper. 

- Mary

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

Mary Barth 

Senior Scientist   

phone: 303-497-8186    email: bar...@ucar.edu

NSF - National Center for Atmospheric Research

P.O. Box 3000

Boulder, CO 80307

Staff web page

SIMA Web Site

MUSICA Web Site

DC3 Web Site 

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

[Mary Barth]

Hi Aman,

Please disregard my last message. I was looking at the wrong namelist.input file. I'm obviously not removing namelists that others have sent in the past. 

I'll take a look at yours later today.

Mary

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

Mary Barth 

Senior Scientist   

phone: 303-497-8186    email: bar...@ucar.edu

NSF - National Center for Atmospheric Research

P.O. Box 3000

Boulder, CO 80307

Staff web page

SIMA Web Site

MUSICA Web Site

DC3 Web Site 

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

[Mary Barth]

Hi Aman,

Here are some things I see in your namelist:

 auxinput6_inname                    = 'wrfbiochemi_d01',

 auxinput7_inname                    = 'wrffirechemi_d<domain>',

!auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',

!auxinput12_inname                   = 'wrf_chem_input',

!auxinput13_inname                   = 'wrfchemv_d<domain>',

 auxinput5_interval_m                = 86400, 86400, 60,

 auxinput7_interval_m                = 86400, 86400, 60,

!auxinput8_interval_m                = 86400, 86400, 60,

 auxinput13_interval_m               = 86400, 86400, 60,

 io_form_auxinput2                   = 2,

 io_form_auxinput5                   = 2,

 io_form_auxinput6                   = 2,

 io_form_auxinput7                   = 0,

!io_form_auxinput8                   = 2,

 io_form_auxinput12                  = 0,

 io_form_auxinput13                  = 0,

auxinput5_inname is not set, so there are no anthropogenic emissions. Although I see the file being opened in rsl.out.0000 for 2019-04-13_00:00:00 (only). 

auxinput5_interval_m and auxinput7_interval_m have values = 86400 which is high. The units are minutes, not seconds. 

io_form_auxinput7 = 0, Should this be = 2 because the fire emissions are using netcdf format for input?

Also, be sure wrfbiochemi information is in the wrfinput file.

The rsl.out.0000 says there is an error opening wrffirechemi_d01. 

Do the initial conditions in wrfinput_d01 look correct? 

Consider using a shorter time step for biogenic emissions (currently bioemdt  =60,). I'd suggest 10 or 15 minutes. 

In summary, I'm wondering if the emissions files are being read properly. I'd suggest checking and running for a few hours and review results before running the entire simulation. 

- Mary

And sorry I confused anyone by looking at the wrong file earlier. Good luck!

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

Mary Barth 

Senior Scientist   

phone: 303-497-8186    email: bar...@ucar.edu

NSF - National Center for Atmospheric Research

P.O. Box 3000

Boulder, CO 80307

Staff web page

SIMA Web Site

MUSICA Web Site

DC3 Web Site 

^--^--^--^--^--^--^--^--^--^--^--^--^--^--^-^--^--^--^-^--^

On Mon, Nov 4, 2024 at 1:26 AM Aman Khan <amanwah...@gmail.com> wrote:

[MarwaClimate2024]

Hello, Professor Aman. I hope you are well. I am a researcher and I am doing research on the simulation model you are talking about. Can you explain to me what this model is and how it works?

[Georgy Nerobelov]

Hi Aman,

Have the same issue with chem_opt=32

Could you resolve it?

Respectfully,

Georgii.