Hi, all!
I’m doing a simulation with ‘sulfur tracking method’ (STM) in CMAQv5.3.1.
I’ve got some hints of calculating ‘aqueous-phase sulfate formation’ in aqchem.F.
But I can’t find any codes for calculating ‘sulfate from boundary and initial condition’ (ASO4ICBC) in CMAQv5.3.1 source codes.

Do you guys know where the calculation source code of ‘sulfate from boundary and initial condition’ is located ??

I wish I can get any responses from you all :slight_smile:

Hi Phy,

The source code for sulfate from boundary and initial conditions is located in two places. It is applied in CGRID_SPCS.F (CMAQ/CGRID_SPCS.F at master · USEPA/CMAQ · GitHub). Beginning on line 446 you can see the implementation of the preloaded STM species into the transported aerosols array.

Each species loaded here has every field from the AE namelist determined already in the STM module. You can find that table in STM_VARS.F , beginning on line 81 (CMAQ/STM_VARS.F at master · USEPA/CMAQ · GitHub).


Dear Dr. Murphy

Hi, I’m HyeonYeong Park, having a question of STM ICBC species.
First, Thank you for your previous response :slight_smile:

But I think that CGRID_SPCS.F and STM_VARS.F assign and define CMAQ species array and also STM species.

I also want to know the calculation routine of ICBC species concentration similar to sulfate aqueous-phase concentration routine (e.g. H2O2, O3, Fe and Mn reactions in aqchem.F).

Can you let me know another location of source code of ASO4ICBC ??

Thank you.

HyeonYeong Park

Hi HyeonYeong,

Sorry for my severely delayed response. I’m not sure I completely follow the question. The IC and BC species for STM are just assigned directly from the IC and BC file. This is done completely analogously to the other CMAQ species.

Specifically, load_cgrid.F is responsible for loading the initial conditions, and rdbcon.F is responsible for loading boundary conditions.


Hi, Dr. Murphy!
Thank you for your answering to my question.

I confirmed those source files, and I think those files what I want.

Thank you for your response again.

HyeonYeong Park