Difference between revisions of "GEOS-Chem v8-02-01"
(→Minor error in gamap_mod.f) |
|||
Line 99: | Line 99: | ||
--[[User:Bmy|Bob Y.]] 12:28, 6 May 2009 (EDT) | --[[User:Bmy|Bob Y.]] 12:28, 6 May 2009 (EDT) | ||
+ | |||
+ | === Bug fix for Streets emissions on 0.5x0.666 nested grid === | ||
+ | |||
+ | We found the v8-01-04 of Geos-Chem doesn't read the correct Streets emission files for the 0.5x0.666 nested grid over China. Right now it would read Streets emissions at 1x1 resolution and regrid them. | ||
+ | |||
+ | To correct it, you need to change : | ||
+ | ! Read David Streets' emisisons over China / SE ASia | ||
+ | IF ( LSTREETS .and. ITS_A_NEW_MONTH() ) THEN | ||
+ | CALL EMISS_STREETS_ANTHRO | ||
+ | ENDIF | ||
+ | |||
+ | by : | ||
+ | ! Read David Streets' emisisons over China / SE ASia | ||
+ | IF ( LSTREETS .and. ITS_A_NEW_MONTH() ) THEN | ||
+ | #if defined(GRID05x0666) | ||
+ | CALL EMISS_STREETS_ANTHRO_05x0666 | ||
+ | #else | ||
+ | CALL EMISS_STREETS_ANTHRO | ||
+ | #endif | ||
+ | ENDIF | ||
+ | |||
+ | in emission_mod.f for every call to EMISS_STREETS_ANTHRO. I.e in the case of full chemistry, offline aerosol, tagged CO, CH4, CO2 and offline H2/HD. | ||
+ | |||
+ | In streets_anthro_mod.f, in the subroutine emiss_streets_anthro_05x0666, you need to change all the lines : | ||
+ | CALL READ_STREETS(... | ||
+ | by : | ||
+ | CALL READ_STREETS_05x0666(... | ||
+ | |||
+ | WARNING : make sure to change the calls to READ_STREETS only in the emiss_streets_anthro_05x0666 subroutine and NOT in emiss_streets_anthro. | ||
+ | |||
+ | --[[User:Ccarouge|Ccarouge]] 15:30, 11 March 2009 (EDT) |
Revision as of 19:43, 6 May 2009
Contents
Overview
CHEMISTRY RELEASE -- ABOUT TO BE RELEASED!!
Will contain everything in v8-01-04 plus the following:
- Glyoxal chemistry (May Fu) (NOTE: This is an option which is turned off by default)
- Updated J(O1D) cross sections (Lin Zhang)
- Updated reaction rates (Jingqiu Mao, Dylan Millet, Palmer group @ U. Edinburgh)
- HO2 uptake (Lyatt Jaegle)
- Updated dust single scattering albedo (Randall Martin)
1-year benchmarks
Run0
Three GEOS-Chem model versions were compared to each other:
Color | Quantity Plotted | Met Field Type | Anthro Emissions | Photolysis Mechanism |
Chemical Mechanism |
Annual Mean OH [105 molec/cm3] |
---|---|---|---|---|---|---|
Red | v8-01-04 Run1 | GEOS-5 4x5 version 5.0.1 w/ "quick fix" for optical depth spinup: 2005 run: 2005 |
GEIA/Piccot EDGAR emissions EMEP European emissions BRAVO Mexican emissions David Streets 2006 emissions CAC Canadian emissions EPA/NEI99 with ICARTT fix EDGAR ship emissions ARCTAS ship SO2 emissions Anthro scale year 2005 |
"old" jv_spec.dat | Mostly unchanged from v5-07-08 except for a few revisions |
10.294 |
Green | v8-01-04 Run2 | GEOS-5 4x5 version 5.1.0, "reprocessed" met fields spinup: 2004 run: 2005 |
" " | " " | " " | 11.099 |
Blue | v8-02-01 Run0 | " " | " " | "updated" jv_spec.dat | "updated" chemical mechanism (cf. J. Mao, D. Millet, T-M. Fu, Palmer Group @ U. Edinburgh) |
11.817 |
Black | Observations |
SUMMARY:
- v8-01-04 Run2 vs. v8-01-04 Run1 is a clean comparison between the GEOS-5.0.1 and GEOS-5.1.0 met products.
- v8-02-01 Run0 vs. v8-01-04 Run2 is a clean validation for the new chemistry mechanism, photolysis rates, and other updates that were added into v8-02-01.
- All runs used the same:
- A.T.E. Algorithm = RPMARES
- Advection routine: "new" tpcore
- Biogenic emissions: "new" MEGAN AEF's for isoprene
- GFED2 biomass emissions
- Aircraft NOx emissions
- Approximate random cloud overlap option in photolysis
The output plots for Run0 may be downloaded from:
ftp ftp.as.harvard.edu cd pub/geos-chem/1yr_benchmarks/v8-02-01/geos5/2005/Run0/output
Previous issues now resolved in v8-02-01
Minor error in gamap_mod.f
Noelle Selin (selin@mit.edu) wrote:
- In my testing/debugging, the compiler found another problem: in gamap_mod, line 1313 and forward (init_tracerinfo):
DO T = 1, N_TRACERS ! Store quantities for each tracer NAME (T,45) = TRACER_NAME(T) FNAME(T,45) = TRIM( NAME(T,45) ) // ' tracer' MOLC (T,45) = INT( TRACER_COEFF(T,1) ) SCALE(T,45) = 1.0e+9 INDEX(T,45) = N
- The error pulled was that N did not have a previous value. Should it be T? (I changed it to T and it seemed OK...but again it's gamap_mod so it doesn't seem like a big deal).
Philippe Le Sager (plesager@seas.harvard.edu) wrote:
- yes, it should be T. We recently found about that when debugging for Jenny, and we already put it in the next version being tested now.
--Bob Y. 12:28, 6 May 2009 (EDT)
Bug fix for Streets emissions on 0.5x0.666 nested grid
We found the v8-01-04 of Geos-Chem doesn't read the correct Streets emission files for the 0.5x0.666 nested grid over China. Right now it would read Streets emissions at 1x1 resolution and regrid them.
To correct it, you need to change :
! Read David Streets' emisisons over China / SE ASia IF ( LSTREETS .and. ITS_A_NEW_MONTH() ) THEN CALL EMISS_STREETS_ANTHRO ENDIF
by :
! Read David Streets' emisisons over China / SE ASia IF ( LSTREETS .and. ITS_A_NEW_MONTH() ) THEN #if defined(GRID05x0666) CALL EMISS_STREETS_ANTHRO_05x0666 #else CALL EMISS_STREETS_ANTHRO #endif ENDIF
in emission_mod.f for every call to EMISS_STREETS_ANTHRO. I.e in the case of full chemistry, offline aerosol, tagged CO, CH4, CO2 and offline H2/HD.
In streets_anthro_mod.f, in the subroutine emiss_streets_anthro_05x0666, you need to change all the lines :
CALL READ_STREETS(...
by :
CALL READ_STREETS_05x0666(...
WARNING : make sure to change the calls to READ_STREETS only in the emiss_streets_anthro_05x0666 subroutine and NOT in emiss_streets_anthro.
--Ccarouge 15:30, 11 March 2009 (EDT)