Hello, recently I've been simulating SOA formation using the CAM-Chem module, and I have some questions regarding the mo_prod_loss.F90 file. Specifically, I'd like to ask about the origin of the coefficients used in the production code for species 165–169 (SOAG0–SOAG5). For example:
prod(k,165(SOAG0)) = (.390_r8*rxt(k,365)*y(k,194(CH3CO3)) + .310_r8*rxt(k,366)*y(k,195(CH3O2)) + &
.360_r8*rxt(k,368)*y(k,124(NO)) + .400_r8*rxt(k,369)*y(k,126(NO3)))*y(k,202(ISOPAO2)) &
+ .300_r8*rxt(k,385)*y(k,134(O3))*y(k,98(ISOP)) + .288_r8*rxt(k,49)*y(k,103(ISOPOOH))
Moreover, I noticed the production of SOAG2 mentions reactions involving S with O3 and OH:
prod(k,167(SOAG2)) = (rxt(k,464) + rxt(k,465)*y(k,134(O3)) + rxt(k,468)*y(k,211(OH)))*y(k,147(S)) &
+ rxt(k,462)*y(k,136(OCS))*y(k,133(O)) + rxt(k,112)*y(k,150(SO2))
I'm also curious about this particular reaction. Could it actually lead to SOA formation? I look forward to your response!
prod(k,165(SOAG0)) = (.390_r8*rxt(k,365)*y(k,194(CH3CO3)) + .310_r8*rxt(k,366)*y(k,195(CH3O2)) + &
.360_r8*rxt(k,368)*y(k,124(NO)) + .400_r8*rxt(k,369)*y(k,126(NO3)))*y(k,202(ISOPAO2)) &
+ .300_r8*rxt(k,385)*y(k,134(O3))*y(k,98(ISOP)) + .288_r8*rxt(k,49)*y(k,103(ISOPOOH))
Moreover, I noticed the production of SOAG2 mentions reactions involving S with O3 and OH:
prod(k,167(SOAG2)) = (rxt(k,464) + rxt(k,465)*y(k,134(O3)) + rxt(k,468)*y(k,211(OH)))*y(k,147(S)) &
+ rxt(k,462)*y(k,136(OCS))*y(k,133(O)) + rxt(k,112)*y(k,150(SO2))
I'm also curious about this particular reaction. Could it actually lead to SOA formation? I look forward to your response!