GEOS-Chem Nested Model
All users interested in the GEOS-Chem regional air quality simulations (e.g. nested-grid, GEOS-Chem/CMAQ linkages, etc.) are encouraged to subscribe to the regional air quality email list (click on the link in the contact information section below).
Contact information
| Regional Air Quality Working Group Co-Chairs | Yuxuan Wang, Jun Wang and Lin Zhang |
| Regional Air Quality Working Group email list | geos-chem-regional [at] seas.harvard.edu |
| To subscribe to email list | Send email to geos-chem-regional-join [at] seas.harvard.edu |
| To unsubscribe from email list | Send email to geos-chem-regional-leave [at] seas.harvard.edu |
Current Regional Working Group Projects (Please add yours!)
| User Group | Description | Funding | Contact Person |
|---|---|---|---|
| Tsinghua | Nested-grid GEOS-Chem model development; Regional air quality over East Asia and interactions with the globe | N/A | Yuxuan Wang |
| SNU | Regional PM simulations over East Asia | Ministry of Environment in Korea | Rokjin Park |
| US EPA | N/A | Rob Pinder | |
| U. Tennessee | N/A | Joshua Fu | |
| Nat'l Obs. Athens | N/A | Anna Protonotariou | |
| U. L'Aquila | Buondary conditions to CHIMERE model over Europe, focus on Saharan dust transport. | ASI (Italian Space Agency) | Gabriele Curci |
| Georgia Tech. | N/A | Tao Zeng | |
| U. Toronto | N/A | Jane Liu | |
| Peking University | Regional air quality in East Asia, Inversion of trace gas emissions using satellite products | N/A | May Fu |
| Dalhousie U. | Surface PM2.5 retrieval using satellite products | N/A | Aaron van Donkelaar |
| Dalhousie U. | Nested aerosol microphysics (TOMAS). NOx/ultrafine aerosol concentration relationships. Topographic influence on microphysics | N/A | Jeffrey Pierce |
| SUNY-Albany | Regional particle formation, evolution, and properties: Nested-grid GEOS-Chem (with APM) simulations and comparisons with measurements | Fangqun Yu | |
| U. Edinburgh | Quantify the impact of boreal biomass burning outflow on tropospheric chemistry over the North Atlantic: Nested grid simulations and evaluation with ground-based, aircraft, and satellite observations | Mark Parrington | |
| Emory University | Improving MISR Aerosol retrieval with GEOS-Chem simulated aerosol composition | NASA | Yang Liu |
| Peking University | Carbonaceous aerosols in China | May Fu | |
| Peking University | Two-way nesting model for GESO-Chem and impacts on regional air quality and long-range transport | Jintai Lin | |
| CSU | Investigation of inorganic aerosol and role of NH3 during CalNex | Luke Schiferl | |
| Harvard | Application to SEAC4RS aircraft campaign in Southeast Asia | Jenny Fisher, Patrick Kim, Lin Zhang | |
| Harvard | Variability of HCHO Over the Southeastern United States Observed from Space | NASA AQAST | Lei Zhu |
| UNL | GEOS-Chem Satellite Simulator for future sensors | NASA | Jun Wang |
Recent updates to nested GEOS-Chem
Fixes to Streets emissions for model years 2005 and earlier
This update was tested in the 1-month benchmark simulation v9-01-03f and approved on 17 Feb 2012.
We discovered two bugs in the implementation of Streets emissions for the nested model. These will affect all nested simulations at least through v9-01-02 for model years 2005 and earlier.
1. The nested model was using Streets 2000 emissions for years 2001-2005 instead of Streets 2006 emissions, which the global model has been using since v8-03-01. This will be updated in GEOS-Chem v9-01-03.
2. Streets 2000 anthropogenic emissions include NOx, CO, and SO2 from ships, which leads to double-counting of ship emissions as well as improper treatment of NOx emissions. This problem is described in more detail on the ship emissions page. This will be updated in GEOS-Chem v9-01-03.
--Jenny Fisher, 1 February 2012
Nested-grid capability for offline simulations
The following offline simulations have been updated for the nested grid:
| Simulation | Author | Status |
|---|---|---|
| Nested-grid CO2 simulation | Yuxuan Wang (Tsinghua U.) | Added to GEOS-Chem v9-01-03; released 14 Sep 2012 |
| CH4 simulation (North America) | Kevin Wecht (Harvard) | Added to GEOS-Chem v9-01-03; released 14 Sep 2012 |
| Hg simulation (North America) | Yanxu Zhang (U. Washington) | Slated for inclusion into GEOS-Chem v9-02 |
--Bob Y. 16:42, 14 September 2012 (EDT)
Recent linking tool updates for GEOS-Chem and CMAQ
On-going Developments
1). main.f, add following lines after DO_CHEMISTRY in line 941:
! check STT (yxw)
#if defined( GEOS_5 ) && defined( GRID05x0666 )
! Loop over grid boxes
!$OMP PARALLEL DO
!$OMP+DEFAULT( SHARED )
!$OMP+PRIVATE( I, J, L, N )
DO N = 1, N_TRACERS
DO L = 1, LLPAR
DO J = 1, JJPAR
DO I = 1, IIPAR
!---------------------------
! Check for Negatives
!---------------------------
IF ( STT(I,J,L,N) < 0d0 ) THEN
!$OMP CRITICAL
WRITE( 6, 100 ) I, J, L, N, STT(I,J,L,N)
PRINT*, 'Neg STT after chemistry ' //
& 'SET STT TO BE ZERO'
STT(I,J,L,N) = 0d0
!$OMP END CRITICAL
!---------------------------
! Check for NaN's
!---------------------------
ELSE IF ( IT_IS_NAN( STT(I,J,L,N) ) ) THEN
!$OMP CRITICAL
WRITE( 6, 100 ) I, J, L, N, STT(I,J,L,N)
PRINT*, 'NaN STT after chemistry ' //
& 'SET STT TO BE LOWER LEVEL'
STT(I,J,L,N) = STT(I,J,L-1,N)
!$OMP END CRITICAL
!----------------------------
! Check STT's for Infinities
!----------------------------
ELSE IF ( .not. IT_IS_FINITE( STT(I,J,L,N) ) ) THEN
!$OMP CRITICAL
WRITE( 6, 100 ) I, J, L, N, STT(I,J,L,N)
PRINT*, 'Inf STT after chemistry ' //
& 'SET STT TO BE LOWER LEVEL'
STT(I,J,L,N) = STT(I,J,L-1,N)
!$OMP END CRITICAL
ENDIF
ENDDO
ENDDO
ENDDO
ENDDO
!$OMP END PARALLEL DO
#endif
2). main.f, add following lines after ITS_TIME_FOR_DIAG in line 1021:
! Sometimes STT in the stratosphere can be negative at
! the nested-grid domain edges. Force them to be zero before
! CHECK_STT (yxw)
#if defined( GEOS_5 ) && defined( GRID05x0666 )
! Loop over grid boxes
!$OMP PARALLEL DO
!$OMP+DEFAULT( SHARED )
!$OMP+PRIVATE( I, J, L, N )
DO N = 1, N_TRACERS
DO L = 1, LLPAR
DO J = 1, JJPAR
DO I = 1, IIPAR
!---------------------------
! Check for Negatives
!---------------------------
IF ( STT(I,J,L,N) < 0d0 ) THEN
!$OMP CRITICAL
WRITE( 6, 100 ) I, J, L, N, STT(I,J,L,N)
PRINT*, 'Neg STT after dynamics step ' //
& 'SET STT TO BE ZERO'
STT(I,J,L,N) = 0d0
!$OMP END CRITICAL
ENDIF
ENDDO
ENDDO
ENDDO
ENDDO
!$OMP END PARALLEL DO
100 FORMAT( ' STT(',i3,',',i3,',',i3,',',i3,') = ', f13.6 )
#endif
--Wang,Chen, 05:13, 25 June 2009 (EDT) --[[User:Aaron van Donkelaar], March 4, 2010
How to run the 0.5x0.667 nested-grid for GEOS-5
We have moved this discussion to the following wiki pages:
--Bob Y. 14:07, 15 November 2010 (EST)
About linking tools
TBD