Difference between revisions of "Sulfate aerosols"
From Geos-chem
| Line 7: | Line 7: | ||
--[[User:Bmy|Bob Y.]] 11:05, 16 February 2010 (EST) | --[[User:Bmy|Bob Y.]] 11:05, 16 February 2010 (EST) | ||
| − | == Updated THNO3.geos4.4x5 file == | + | == References == |
| + | |||
| + | #Park, R. J., D. J. Jacob, B. D. Field, R. M. Yantosca, and M. Chin, ''Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: implications for policy'', <u>J. Geophys. Res.</u>, '''109''', D15204, 10.1029/2003JD004473, 2004. [http://acmg.seas.harvard.edu/publications/park2004.pdf PDF] | ||
| + | |||
| + | == Known issues == | ||
| + | |||
| + | === Updated THNO3.geos4.4x5 file === | ||
'''''Lyatt Jaegle (jaegle@atmos.washington.edu) wrote:''''' | '''''Lyatt Jaegle (jaegle@atmos.washington.edu) wrote:''''' | ||
| Line 19: | Line 25: | ||
--[[User:Bmy|Bob Y.]] 14:19, 24 April 2009 (EDT) | --[[User:Bmy|Bob Y.]] 14:19, 24 April 2009 (EDT) | ||
| − | == Fix for mass balance of HNO3 and NIT == | + | === Fix for mass balance of HNO3 and NIT === |
'''''Becky Alexander (beckya@atmos.washington.edu) wrote:''''' | '''''Becky Alexander (beckya@atmos.washington.edu) wrote:''''' | ||
| Line 42: | Line 48: | ||
:is appropriate as long as you also have <tt>PNIT</tt> (analogous to <tt>PNITs</tt>). <tt>PNIT</tt> is in my original code where I did all my mass balance testing. <tt>PNIT</tt> got dropped when going to the standard version. I don't recall dropping this, but my guess is that I decided it was redundant to have it when isorropia would just repartition HNO3 and NIT anyway according to thermodynamic equilibrium. But when dropping <tt>PNIT</tt>, you have to change <tt>TITR_HNO3</tt> to <tt>HNO3_SSC</tt> in the above equation in order to achieve mass balance. | :is appropriate as long as you also have <tt>PNIT</tt> (analogous to <tt>PNITs</tt>). <tt>PNIT</tt> is in my original code where I did all my mass balance testing. <tt>PNIT</tt> got dropped when going to the standard version. I don't recall dropping this, but my guess is that I decided it was redundant to have it when isorropia would just repartition HNO3 and NIT anyway according to thermodynamic equilibrium. But when dropping <tt>PNIT</tt>, you have to change <tt>TITR_HNO3</tt> to <tt>HNO3_SSC</tt> in the above equation in order to achieve mass balance. | ||
| − | '''''NOTE: This fix is | + | '''''NOTE: This fix is was standardized in [[GEOS-Chem v8-01-02]].''''' |
| − | --[[User:Bmy|Bob Y.]] | + | --[[User:Bmy|Bob Y.]] 16:36, 19 February 2010 (EST) |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
Revision as of 21:36, 19 February 2010
Contents
Emissions
Please see the following wiki pages which contain more information about sulfate aerosol emissions.
--Bob Y. 11:05, 16 February 2010 (EST)
References
- Park, R. J., D. J. Jacob, B. D. Field, R. M. Yantosca, and M. Chin, Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: implications for policy, J. Geophys. Res., 109, D15204, 10.1029/2003JD004473, 2004. PDF
Known issues
Updated THNO3.geos4.4x5 file
Lyatt Jaegle (jaegle@atmos.washington.edu) wrote:
- I was trying to run GEOS-Chem (v8-01-03) in the offline aerosol mode with GEOS-4 met fields and ran into a problem: it seems that the offline file GEOS_4x5/sulfate_sim_200508/offline/THNO3.geos4.4x5 is in units of ppbv instead of v/v as expected in routine GET_HNO3_UGM3. This leads to issues in RPMARES which thinks that HNO3 is very large. I checked all the other files, and they are in v/v, as expected.
Bob Yantosca (yantosca@seas.harvard.edu) wrote:
- Thanks Lyatt. I've downloaded the file and updated the README in the GEOS_4x5/sulfate_sim_200508/offline directory.
--Bob Y. 14:19, 24 April 2009 (EDT)
Fix for mass balance of HNO3 and NIT
Becky Alexander (beckya@atmos.washington.edu) wrote:
- We need to make a change in sulfate_mod in order to have mass balance for HNO3 and NIT. Duncan Fairlie noticed the bug. There is a simple change:
- In routine SEASALT_CHEM in sulfate_mod.f: In order to have mass balance, you need to change:
!HNO3 lost [eq/timestep] converted back to [v/v/timestep] HNO3_ss = TITR_HNO3 * 0.063 * TCVV(IDTHNO3)/AD(I,J,L)
- to:
!HNO3 lost [eq/timestep] converted back to [v/v/timestep] HNO3_ss = HNO3_SSC * 0.063 * TCVV(IDTHNO3)/AD(I,J,L)
- In my original code where I added isorropia and the new tracers, NITs and SO4s, the line above:
!HNO3 lost [eq/timestep] converted back to [v/v/timestep] HNO3_ss = TITR_HNO3 * 0.063 * TCVV(IDTHNO3)/AD(I,J,L)
- is appropriate as long as you also have PNIT (analogous to PNITs). PNIT is in my original code where I did all my mass balance testing. PNIT got dropped when going to the standard version. I don't recall dropping this, but my guess is that I decided it was redundant to have it when isorropia would just repartition HNO3 and NIT anyway according to thermodynamic equilibrium. But when dropping PNIT, you have to change TITR_HNO3 to HNO3_SSC in the above equation in order to achieve mass balance.
NOTE: This fix is was standardized in GEOS-Chem v8-01-02.
--Bob Y. 16:36, 19 February 2010 (EST)
