Problem with adding SOA variables in iofields in WRF-Chem V4.4.1

[neyestani]

Hi there,

I'm trying to add some variables to my wrfout file. So, added these two lines to my namelist:

iofields_filename = "man_output.txt","man_output.txt"
ignore_iofields_warning = .true.,

Also, I defined the variables I wanted in my wrfout in man_output.txt file by adding this line:

+:h:0:TAUAER1,TAUAER2,TAUAER3,TAUAER4,Z,Z_AT_W,soa_a01,soa_a02,soa_a03,soa_a04,ebu_pm25

After running WRF-Chem, I have all the added variables in my wrfout, except for SOA variables. Also, I received these warnings at my rsl.error files:

Domain 1 Setting history stream 0 for tauaer1
Domain 1 Setting history stream 0 for tauaer2
Domain 1 Setting history stream 0 for tauaer3
Domain 1 Setting history stream 0 for tauaer4
Domain 1 Setting history stream 0 for z
Domain 1 Setting history stream 0 for z_at_w
W A R N I N G : Unable to modify mask for soa_a01. Variable not found. File: man_output.txt at line 1
W A R N I N G : Unable to modify mask for soa_a02. Variable not found. File: man_output.txt at line 1
W A R N I N G : Unable to modify mask for soa_a03. Variable not found. File: man_output.txt at line 1
W A R N I N G : Unable to modify mask for soa_a04. Variable not found. File: man_output.txt at line 1
Domain 1 Setting history stream 0 for ebu_pm25
Ignoring problems reading man_output.txt
Continuing. To make this a fatal error, set ignore_iofields_warn to false in namelist

I'm using MOZART/MOSAIC (chem_opt=201). Please let me know if more details are required.

I appreciate any help on how to get SOA concentrations in my output.


Best,
Soroush

 

[neyestani]

Hi,

This is just a reminder to see if anyone knows how to add SOA fields in wrfout with chem_opt=201.


Best,
Soroush

 

[jordanschnell]

Hi Soroush,

It's unclear to me why the io_fields options isn't working. You can try modifying the registry instead:


state real soa_a01 ikj misc 1 - r "soa_a01" "soa_a01" "ug m^-3"

to

state real soa_a01 ikj misc 1 - rh "soa_a01" "soa_a01" "ug m^-3"

Note you will need to ./clean and recompile.

Jordan

 

[neyestani]

Thank you for the reply, Jordan.

I edited soa lines in registry.chem from:

state real soa_a01 ikj misc 1 - - "soa_a01" "soa_a01" "ug m^-3"
state real soa_a02 ikj misc 1 - - "soa_a02" "soa_a02" "ug m^-3"
state real soa_a03 ikj misc 1 - - "soa_a03" "soa_a03" "ug m^-3"
state real soa_a04 ikj misc 1 - - "soa_a04" "soa_a04" "ug m^-3"
state real soa_a05 ikj misc 1 - - "soa_a05" "soa_a05" "ug m^-3"
state real soa_a06 ikj misc 1 - - "soa_a06" "soa_a06" "ug m^-3"
state real soa_a07 ikj misc 1 - - "soa_a07" "soa_a07" "ug m^-3"
state real soa_a08 ikj misc 1 - - "soa_a08" "soa_a08" "ug m^-3"

to:

state real soa_a01 ikj misc 1 - rh "soa_a01" "soa_a01" "ug m^-3"
state real soa_a02 ikj misc 1 - rh "soa_a02" "soa_a02" "ug m^-3"
state real soa_a03 ikj misc 1 - rh "soa_a03" "soa_a03" "ug m^-3"
state real soa_a04 ikj misc 1 - rh "soa_a04" "soa_a04" "ug m^-3"
state real soa_a05 ikj misc 1 - rh "soa_a05" "soa_a05" "ug m^-3"
state real soa_a06 ikj misc 1 - rh "soa_a06" "soa_a06" "ug m^-3"
state real soa_a07 ikj misc 1 - rh "soa_a07" "soa_a07" "ug m^-3"
state real soa_a08 ikj misc 1 - rh "soa_a08" "soa_a08" "ug m^-3"

Then did "./clean -a" and recompiled. Still, I don't have soa fields in my output, though. I have same warnings in rsl.error log files.

Is there any other way to add soa into output?


Best,
Soroush

 

[jordanschnell]

Hi Soroush,

Apologies, it looks like those variables are in a package:

package diag_aero aero_diag_opt==1 - state:hoa_a01,hoa_a02,hoa_a03,hoa_a04,hoa_a05,hoa_a06,hoa_a07,hoa_a08,soa_a01,soa_a02,soa_a03,soa_a04,soa_a05,soa_a06,soa_a07,soa_a08,bboa_a01,bboa_a02,bboa_a03,bboa_a04,bboa_a05,bboa_a06,bboa_a07,bboa_a08,bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04,bbsoa_a05,bbsoa_a06,bbsoa_a07,bbsoa_a08,hsoa_a01,hsoa_a02,hsoa_a03,hsoa_a04,hsoa_a05,hsoa_a06,hsoa_a07,hsoa_a08,biog_a01,biog_a02,biog_a03,biog_a04,biog_a05,biog_a06,biog_a07,biog_a08,arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04,arosoa_a05,arosoa_a06,arosoa_a07,arosoa_a08,totoa_a01,totoa_a02,totoa_a03,totoa_a04,totoa_a05,totoa_a06,totoa_a07,totoa_a08,hsoa_c,hsoa_o,bbsoa_c,bbsoa_o,biog_v1,biog_v2,biog_v3,biog_v4,ant_v1,ant_v2,ant_v3,ant_v4,smpa_v1,smpbb_v1,asmpsoa_a01,asmpsoa_a02,asmpsoa_a03,asmpsoa_a04

Meaning you need to set `aero_diag_opt = 1`

Jordan

 

[neyestani]

Hi Jordan,

Thanks! Setting "aero_diag_opt" to 1 in the namelist worked.


Best,
Soroush