Difference between revisions of "GEOS-Chem v9-02"

From Geos-chem
Jump to: navigation, search
(Bug fixes in biofuel_mod.F and emfossil.F)
(What's new in this version)
Line 56: Line 56:
 
|Science
 
|Science
 
|Yanxu Zhang (U. Washington)<br>Lyatt Jaegl&eacute; (U. Washington)
 
|Yanxu Zhang (U. Washington)<br>Lyatt Jaegl&eacute; (U. Washington)
|v9-02c<br>+ Hg benchmark (by [[Hg and POPs Working Group]])
+
|[[GEOS-Chem v9-02 benchmark history#v9-02c|v9-02c]]<br>+ Hg benchmark (by [[Hg and POPs Working Group]])
|TBD
+
|Pending approval
 
|-
 
|-
 
|[[Mercury#Nested-grid Hg simulation over North America|Nested-grid Hg simulation]]
 
|[[Mercury#Nested-grid Hg simulation over North America|Nested-grid Hg simulation]]
Line 65: Line 65:
 
|" "
 
|" "
 
|-
 
|-
|Streets future Hg emissions
+
|[[Mercury#Streets future Hg emissions|Streets future Hg emissions]]
 
|Science
 
|Science
 
|Bess Corbitt (Harvard)
 
|Bess Corbitt (Harvard)
Line 71: Line 71:
 
|" "
 
|" "
 
|-
 
|-
|Correct Hg simulation to use updated Hg(0) oxidation kinetics and Br/BrO fields
+
|[[Mercury#Updated Hg.280.29 oxidation kinetics|Correct Hg simulation to use updated Hg(0) oxidation kinetics]]
 +
|Science
 +
|Bess Corbitt(Harvard)<br>Justin Parrella (Harvard)
 +
|" "
 +
|" "
 +
|-
 +
|[[Mercury#Capability to use GEOS-Chem Br.2FBrO fields in Hg simulation|Capability to use GEOS-Chem Br/BrO fields in Hg simulation]]
 
|Science
 
|Science
 
|Bess Corbitt(Harvard)<br>Justin Parrella (Harvard)
 
|Bess Corbitt(Harvard)<br>Justin Parrella (Harvard)
Line 79: Line 85:
 
|[[POPs simulation|POPs simulation: Polycyclic aromatic hydrocarbons]]
 
|[[POPs simulation|POPs simulation: Polycyclic aromatic hydrocarbons]]
 
|Science
 
|Science
|Carey Friedman (MIT)
+
|Carey Friedman (MIT)<br>Noelle Selin (MIT)
 +
|" "
 +
|" "
 +
|-
 +
|[[#Bug fix for updated CAC emissions|Bug fix for updated CAC emissions]]
 +
|Bug fix
 +
|Aaron van Donkelaar (Dalhousie)
 
|" "
 
|" "
 
|" "
 
|" "
 
|-bgcolor="#CCFFFF"
 
|-bgcolor="#CCFFFF"
 +
|Grid-independent GEOS-Chem updates
 +
|Structural
 +
|[[GEOS-Chem Support Team]]
 +
|v9-02d
 +
|TBD
 +
|-
 
|[[Updating_standard_chemistry_with_JPL_10-6|Inorganic chemistry updates]]
 
|[[Updating_standard_chemistry_with_JPL_10-6|Inorganic chemistry updates]]
 
|Science
 
|Science
 
|Mat Evans (York)
 
|Mat Evans (York)
|v9-02d<br>+ 1 year benchmark
+
|v9-02e<br>+ 1 year benchmark
 
|TBD
 
|TBD
|-bgcolor="#CCFFFF"
+
|-
 
|Methyl peroxy nitrate (MPN) chemistry
 
|Methyl peroxy nitrate (MPN) chemistry
 
|Science
 
|Science
Line 94: Line 112:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|Fix RO2+HO2 rate constant
 
|Fix RO2+HO2 rate constant
 
|Science
 
|Science
Line 100: Line 118:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|[[ChemTelecon20111202|Adopt Paulot isoprene scheme as standard]]
 
|[[ChemTelecon20111202|Adopt Paulot isoprene scheme as standard]]
 
|Science
 
|Science
Line 106: Line 124:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|[[ChemTelecon20111202|Improved HO2 uptake by aerosol]]
 
|[[ChemTelecon20111202|Improved HO2 uptake by aerosol]]
 
|Science
 
|Science
Line 113: Line 131:
 
|" "
 
|" "
 
|-
 
|-
|-bgcolor="#CCFFFF"
 
 
|Removal of NOx-Ox partitioning
 
|Removal of NOx-Ox partitioning
 
|Science
 
|Science
Line 119: Line 136:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|[[Anthropogenic_emissions#CAC|Updated Canadian NH3 inventory]]
 
|[[Anthropogenic_emissions#CAC|Updated Canadian NH3 inventory]]
 
|Science
 
|Science
 
|Wai-Ho Lo (Dalhousie)  
 
|Wai-Ho Lo (Dalhousie)  
|v9-02e<br>+ 1-year benchmark
+
|v9-02f<br>+ 1-year benchmark
 
|TBD
 
|TBD
|-
+
|-bgcolor="#CCFFFF"
 
|RCP emission scenarios
 
|RCP emission scenarios
 
|Science
 
|Science
Line 131: Line 148:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|EDGAR v4.2 emissions for CH4 simulation
 
|EDGAR v4.2 emissions for CH4 simulation
 
|Science
 
|Science
Line 137: Line 154:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|[[Scale factors for anthropogenic emissions|Update anthropogenic scale factors through 2010]]
 
|[[Scale factors for anthropogenic emissions|Update anthropogenic scale factors through 2010]]
 
|Science
 
|Science
Line 143: Line 160:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|Update to EDGAR 4 anthropogenic emissions
 
|Update to EDGAR 4 anthropogenic emissions
 
|Science
 
|Science
Line 149: Line 166:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|Interannual lightning
 
|Interannual lightning
 
|Science
 
|Science
Line 155: Line 172:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|AEIC aircraft emissions inventory
 
|AEIC aircraft emissions inventory
 
|Science
 
|Science
Line 161: Line 178:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|[[#Bug fixes in biofuel_mod.F and emfossil.F|Bug fixes in biofuel_mod.F and emfossil.F]]
 
|[[#Bug fixes in biofuel_mod.F and emfossil.F|Bug fixes in biofuel_mod.F and emfossil.F]]
 
|Bug fix
 
|Bug fix
Line 167: Line 184:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|[[Secondary_organic_aerosols#SOA_simulation_with_semi-volatile_POA|SOA simulation with semi-volatile POA]]
 
|[[Secondary_organic_aerosols#SOA_simulation_with_semi-volatile_POA|SOA simulation with semi-volatile POA]]
 
|Science
 
|Science
 
|Havala Pye (EPA)<br>[[GEOS-Chem Support Team]]
 
|Havala Pye (EPA)<br>[[GEOS-Chem Support Team]]
|v9-02f<br>+ 1-year benchmark
+
|v9-02g<br>+ 1-year benchmark
 
|TBD
 
|TBD
|-bgcolor="#CCFFFF"
+
|-
 
|[[Mineral_dust_aerosols#Surface_chemistry_on_dust|Acid uptake on dust aerosols]]
 
|[[Mineral_dust_aerosols#Surface_chemistry_on_dust|Acid uptake on dust aerosols]]
 
|Science
 
|Science
Line 179: Line 196:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|Cloudwater pH for sulfate formation
 
|Cloudwater pH for sulfate formation
 
|Science
 
|Science
Line 185: Line 202:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|Update jv_spec.dat with better representation<br>of OC growth with RH and correction to sulfate optics
 
|Update jv_spec.dat with better representation<br>of OC growth with RH and correction to sulfate optics
 
|Science
 
|Science
Line 191: Line 208:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|Nested-grid model at 0.25&deg; x 0.3125&deg; resolution using [[GEOS-5.7.2|GEOS-5 "FP"]] met data product from GMAO
 
|Nested-grid model at 0.25&deg; x 0.3125&deg; resolution using [[GEOS-5.7.2|GEOS-5 "FP"]] met data product from GMAO
 
|Structural
 
|Structural
 
|Yuxuan Wang (Tsinghua)<br>[[GEOS-Chem Support Team]]
 
|Yuxuan Wang (Tsinghua)<br>[[GEOS-Chem Support Team]]
|v9-02g<br>+ Nested-grid benchmark (by [[Nested Model Working Group|Nested Model WG]])
+
|v9-02h<br>+ Nested-grid benchmark (by [[Nested Model Working Group|Nested Model WG]])
 
|TBD
 
|TBD
|-
+
|-bgcolor="#CCFFFF"
 
|High-resolution emissions for above nested-grid model
 
|High-resolution emissions for above nested-grid model
 
|Science
 
|Science
Line 203: Line 220:
 
|" "
 
|" "
 
|" "
 
|" "
|-
+
|-bgcolor="#CCFFFF"
 
|Global simulations with [[GEOS-5.7.2|GEOS-5 "FP"]] data product from GMAO<br>(Needed for the 0.25&deg; x 0.3125&deg; nested-grid simulations)
 
|Global simulations with [[GEOS-5.7.2|GEOS-5 "FP"]] data product from GMAO<br>(Needed for the 0.25&deg; x 0.3125&deg; nested-grid simulations)
 
|Structural
 
|Structural
Line 209: Line 226:
 
|" "
 
|" "
 
|" "
 
|" "
|-bgcolor="#CCFFFF"
+
|-
 
|Improve temporal resolution of anthropogenic CO2 sources
 
|Improve temporal resolution of anthropogenic CO2 sources
 
|Science
 
|Science
 
|Ray Nassar (Environment Canada)
 
|Ray Nassar (Environment Canada)
|v9-02h<br>+ CO2 benchmark (by [[Carbon_Gases_and_Organics_Working_Group|Carbon Gases & Organics WG]])
+
|v9-02i<br>+ CO2 benchmark (by [[Carbon_Gases_and_Organics_Working_Group|Carbon Gases & Organics WG]])
 
|TBD
 
|TBD
|-bgcolor="#CCFFFF"
+
|-
 
|Replace CASA CO2 biosphere climatology with year-specific fluxes
 
|Replace CASA CO2 biosphere climatology with year-specific fluxes
 
|Science
 
|Science
Line 223: Line 240:
 
|}
 
|}
  
--[[User:Melissa Payer|Melissa Payer]] 18:09, 16 October 2012 (EDT)<br>
+
--[[User:Melissa Payer|Melissa Payer]] 16:58, 27 November 2012 (EST)
--[[User:Bmy|Bob Y.]] 11:18, 23 October 2012 (EDT)
+
  
 
=== New data directories ===
 
=== New data directories ===

Revision as of 21:58, 27 November 2012

NOTE: In light of our new benchmarking procedure instituted at IGC5, every GEOS-Chem version released to the user community is now a public release, in the sense that we run a 1-year benchmark before making the code available. The 1-month benchmarks now replace what we used to call "beta releases". This allows us to retire the last two numbers in the version name. We shall henceforth denote GEOS-Chem versions as GEOS-Chem v9-02, etc. instead of using the previous nomenclature (e.g. GEOS-Chem v9-01-03).

Overview

What's new in this version

GEOS-Chem v9-02 shall contain the following major updates and improvements. The order has been determined by the GEOS-Chem Steering Committee. We will also continue to add source code updates for the Grid-independent GEOS-Chem project at regular intervals.

NOTE: v9-02a, v9-02b, etc. denote 1-month benchmark simulations.

Feature Type Submitted by Benchmark(s) Status
Correction for GEOS-5 PBL heights Bug Fix Sajeev Philip (Dalhousie) v9-02a Approved 16 Oct 2012
New soil NOx emission module Science Rynda Hudman Kay (EPA)
Neil Moore (formerly of Dalhousie)
Bram Maasakkers (Eindhoven)
v9-02b Approved 29 Oct 2012
Bug fix in hdf_mod.F Bug fix Sajeev Philip (Dalhousie) " " " "
Bug fixes for tagged CO simulation Bug fix Jenny Fisher (Wollongong)
GEOS-Chem Support Team
" " " "
Bug fix for nested simulation in nei_anthro_mod.F Bug fix Sajeev Philip (Dalhousie) " " " "
EPA/NEI05 North American Hg emissions Science Yanxu Zhang (U. Washington)
Lyatt Jaeglé (U. Washington)
v9-02c
+ Hg benchmark (by Hg and POPs Working Group)
Pending approval
Nested-grid Hg simulation Science Yanxu Zhang (U. Washington)
Lyatt Jaeglé (U. Washington)
" " " "
Streets future Hg emissions Science Bess Corbitt (Harvard) " " " "
Correct Hg simulation to use updated Hg(0) oxidation kinetics Science Bess Corbitt(Harvard)
Justin Parrella (Harvard)
" " " "
Capability to use GEOS-Chem Br/BrO fields in Hg simulation Science Bess Corbitt(Harvard)
Justin Parrella (Harvard)
" " " "
POPs simulation: Polycyclic aromatic hydrocarbons Science Carey Friedman (MIT)
Noelle Selin (MIT)
" " " "
Bug fix for updated CAC emissions Bug fix Aaron van Donkelaar (Dalhousie) " " " "
Grid-independent GEOS-Chem updates Structural GEOS-Chem Support Team v9-02d TBD
Inorganic chemistry updates Science Mat Evans (York) v9-02e
+ 1 year benchmark
TBD
Methyl peroxy nitrate (MPN) chemistry Science Ellie Browne (UC Berkeley) " " " "
Fix RO2+HO2 rate constant Science Fabien Paulot (Harvard) " " " "
Adopt Paulot isoprene scheme as standard Science Fabien Paulot (Harvard)
Eloise Marais (Harvard)
Jingqiu Mao (NASA/GFDL)
Mat Evans (York)
" " " "
Improved HO2 uptake by aerosol Science Jingqiu Mao (NASA/GFDL) " " " "
Removal of NOx-Ox partitioning Science GEOS-Chem Support Team " " " "
Updated Canadian NH3 inventory Science Wai-Ho Lo (Dalhousie) v9-02f
+ 1-year benchmark
TBD
RCP emission scenarios Science Chris Holmes (UC Irvine) " " " "
EDGAR v4.2 emissions for CH4 simulation Science Kevin Wecht (Harvard) " " " "
Update anthropogenic scale factors through 2010 Science Aaron van Donkelaar (Dalhousie) " " " "
Update to EDGAR 4 anthropogenic emissions Science Qiang Zhang (Tsinghua)
Sajeev Philip (Dalhousie)
" " " "
Interannual lightning Science Lee Murray (Harvard) " " " "
AEIC aircraft emissions inventory Science Steven Barret (MIT) " " " "
Bug fixes in biofuel_mod.F and emfossil.F Bug fix Christoph Keller (Harvard) " " " "
SOA simulation with semi-volatile POA Science Havala Pye (EPA)
GEOS-Chem Support Team
v9-02g
+ 1-year benchmark
TBD
Acid uptake on dust aerosols Science T. Duncan Fairlie (NASA/LARC) " " " "
Cloudwater pH for sulfate formation Science Becky Alexander (U. Washington) " " " "
Update jv_spec.dat with better representation
of OC growth with RH and correction to sulfate optics
Science David Ridley (MIT) " " " "
Nested-grid model at 0.25° x 0.3125° resolution using GEOS-5 "FP" met data product from GMAO Structural Yuxuan Wang (Tsinghua)
GEOS-Chem Support Team
v9-02h
+ Nested-grid benchmark (by Nested Model WG)
TBD
High-resolution emissions for above nested-grid model Science Yuxuan Wang (Tsinghua)
GEOS-Chem Support Team
" " " "
Global simulations with GEOS-5 "FP" data product from GMAO
(Needed for the 0.25° x 0.3125° nested-grid simulations)
Structural GEOS-Chem Support Team " " " "
Improve temporal resolution of anthropogenic CO2 sources Science Ray Nassar (Environment Canada) v9-02i
+ CO2 benchmark (by Carbon Gases & Organics WG)
TBD
Replace CASA CO2 biosphere climatology with year-specific fluxes Science Ray Nassar (Environment Canada)
Dylan Jones (U. Toronto)
" " " "

--Melissa Payer 16:58, 27 November 2012 (EST)

New data directories

The following new data directories have been added for GEOS-Chem v9-02. You will have to download the directories relevant to your simulation.

You can download these directories with anonymous FTP or the Wget utility. For instructions, please see Chapter 2.4, Downloading the GEOS-Chem shared data directories in the GEOS-Chem Online User's Guide.

GEOS_0.5x0.666_NA/RETRO_201103     # RETRO anthropogenic VOC emissions
GEOS_0.5x0.666_EU/RETRO_201103     # RETRO anthropogenic VOC emissions
GEOS_0.5x0.666_CH/RETRO_201103     # RETRO anthropogenic VOC emissions

GEOS_NATIVE/soil_NOx_201208        # Files for updated soil NOx emission scheme

--Melissa Payer 11:44, 30 October 2012 (EDT)

Previous issues now resolved in GEOS-Chem v9-02

Correction for GEOS-5 PBL heights

This update was tested in the 1-month benchmark simulation v9-02a and approved on 16 Oct 2012.

Please see this wiki post for a full description of this issue.

--Bob Y. 11:01, 5 October 2012 (EDT)

Bug fixes for tagged CO simulation

This update was tested in the 1-month benchmark simulation v9-02b and approved on 29 Oct 2012.

We have corrected the following minor issues in the Tagged CO simulation. These bugs were present in GEOS-Chem v9-01-03. The full-chemistry simulation is not affected by these issues.

(1) In emissdr.F, bracket AD46 diagnostics with IF statements in order to avoid out-of-bounds errors caused by undefined bromine tracer flags:

               ! CHBr3 emissions [kg/m2/s] -- tracer #14
               IF ( IDECHBr3 > 0 ) THEN
                  AD46(I,J,14) = AD46(I,J,14) 
     &                         + ( EMISRR(I,J,IDECHBr3)  / AREA_M2 ) 
     &                         * ( MWT_CHBr3             / AVG     )
              ENDIF
 
               ! CH2Br2 emissions [kg/m2/s] -- tracer #15
               IF ( IDECHBr3 > 0 ) THEN
                  AD46(I,J,15) = AD46(I,J,15) 
     &                         + ( EMISRR(I,J,IDECH2Br2) / AREA_M2 ) 
     &                         * ( MWT_CH2Br2            / AVG     )
               ENDIF


(2) We needed to pass the am_I_Root = .TRUE. value to routine GET_GLOBAL_CH4 from tagged_co_mod.F. We changed this line:

         CALL GET_GLOBAL_CH4( GET_YEAR(), .TRUE., 
     &                        A3090S, A0030S, A0030N, A3090N )

to:

         CALL GET_GLOBAL_CH4( GET_YEAR(), .TRUE.,             
     &                        A3090S, A0030S, A0030N, A3090N, .TRUE. )


(3) The CO_PRODS and CO_LOSSS arrays need to be made 3-dimensional for compatibility with the new GMI stratospheric chemistry data. Change these lines:

      ALLOCATE( CO_PRODS( JJPAR, LLPAR ), STAT=AS )
      . . .
      ALLOCATE( CO_LOSSS( JJPAR, LLPAR ), STAT=AS )
      . . .
      IF ( IS_PROD ) THEN
         RATE = CO_PRODS(J,L)   ! P(CO) from CH4 + OH in [v/v/s]
      ELSE
         RATE = CO_LOSSS(J,L)   ! L(CO) from CO + OH  in [s^-1]
      ENDIF

to:

      ALLOCATE( CO_PRODS( IIPAR, JJPAR, LLPAR ), STAT=AS )
      . . .
      ALLOCATE( CO_LOSSS( IIPAR, JJPAR, LLPAR ), STAT=AS )
      . . .
      IF ( IS_PROD ) THEN
         RATE = CO_PRODS(I,J,L)   ! P(CO) from CH4 + OH in [v/v/s]
      ELSE
         RATE = CO_LOSSS(I,J,L)   ! L(CO) from CO + OH  in [s^-1]
      ENDIF

Thanks to Jenny Fisher for her assistance in correcting these issues.

--Bob Y. 11:26, 23 October 2012 (EDT)

Bug fix in nei2005_anthro_mod.F

This update was tested in the 1-month benchmark simulation v9-02b and approved on 29 Oct 2012.

Sajeev Philip wrote:

Bug fix in subroutine: GET_NEI99_SEASON_05x0666

     ! Regrid to current model resolution [unitless]
     DO L = 1, 5
        INGRID = ARRAY(:,:,L)
        CALL DO_REGRID_A2A( LLFILENAME, I1x1,    J1x1,
    &                       INGRID,     OUTGRID, PERAREA=1,
    &                       netCDF=.TRUE.           

should be

    &                       INGRID,     OUTGRID, PERAREA=0,
    &                       netCDF=.TRUE.             

Bug fix in subroutine: GET_VISTAS_SEASON_05x0666

     ! Regrid to current model resolution [unitless]
     DO L = 1, 5
        INGRID = ARRAY(:,:,L)
        CALL DO_REGRID_A2A( LLFILENAME, I1x1,    J1x1,
    &                       INGRID,     OUTGRID, PERAREA=1,
    &                       netCDF=.TRUE.                   )

should be

    &                       INGRID,     OUTGRID, PERAREA=0,
    &                       netCDF=.TRUE.                   )

Otherwise, all NEI emissions will be reduced to 1/3rd of the original values.

--Bob Y. 11:07, 25 October 2012 (EDT)

Bug fix for updated CAC emissions

Aaron van Donkelaar wrote:

The code as implemented in v9-01-03 contains only few of the changes I submitted to update CAC to 2008. My guess is that somehow two versions of this file were being worked on simultaneously.
As a result of the mix-up, simulations using CAC emissions and running 2006 or 2007 will crash. Other than that it shouldn't have had any major ramifications.

--Melissa Payer 15:51, 27 November 2012 (EST)

Bug fixes in biofuel_mod.F and emfossil.F

Christoph Keller wrote:

I found two bugs in the emissions code which I think we should fix.
Firstly, the diurnal NOx scale factors should be based on UTC instead of local time. The fix is easy, we just have to use the local time for a given grid box instead of taking the UTC hour when calculating the emissions in emfossil.F
Secondly, I think that we currently double count biofuel emissions over Canada and Mexico if the respective local inventories (CAC and BRAVO) are used. Again, this is an easy fix, we just have to set the biofuel emissions for the given tracers and regions to zero if CAC and/or BRAVO are used.

--Melissa Payer 10:31, 16 November 2012 (EST)

Outstanding issues not yet resolved in GEOS-Chem v9-02

TBD