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 and Jun Wang
Regional Air Quality Working Group email list	geos-chem-regional@seas.harvard.edu
To subscribe to email list	Send email to geos-chem-regional-join@seas.harvard.edu
To unsubscribe from email list	Send email to geos-chem-regional-leave@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

On-going Developments

1. Modified for Code.v8-02-01

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:

1. GEOS-Chem nested grid simulations
2. Setting up GEOS-Chem nested grid simulations

--Bob Y. 14:07, 15 November 2010 (EST)

About linking tools

TBD