GEOS-Chem v8-03-02
Contents
Overview
PUBLIC RELEASE -- Summer 2010
Will contain everything in GEOS-Chem v8-03-01, plus:
- Terrestrial and deep ocean mercury (N. Smith-Downey, E. Sunderland)
- Updated CO2 simulation (R. Nassar)
- Liquid water content taken from GEOS-5 met fields (J. Fisher)
- Wet scavenging by snow
--Bob Y. 14:36, 3 August 2010 (EDT)
Previous issues now resolved in v8-03-02
Patches from v8-03-01
The following bug fix patches are now standardized in v8-03-02:
- Patch to enable call to SOA_PARA_INIT
- Patch to fix declaration of IONIC in ISORROPIA II
- Patch to fix bug in Aromatic SOA
- Patch to scale AOD diagnostic output to other wavelengths
--Bob Y. 12:58, 30 June 2010 (EDT)
Fix for EPA/NEI 2005 emissions
Aaron van Donkelaar submitted updated EPA/NEI05 emissions files in order to fix abnormally low VOC emissions. This fix will be standardized in v8-03-02. Please see this wiki page for more information.
--Bob Y. 15:38, 12 July 2010 (EDT)
Minor fixes in gamap_mod.f
Prasad Kasibhatla reported a bug in gamap_mod.f that was caused by the definition of the parameters:
INTEGER, PARAMETER :: MAXDIAG = 70 INTEGER, PARAMETER :: MAXTRACER = 120
The problem is that when the actual number of tracers is increased beyond 120, this results in an array-out-of-bounds error.
The fix is to use the MAX_DIAG and MAX_TRACER parameters from the CMN_DIAG header file to define the MAXDIAG and MAXTRACER parameters in gamap_mod.f, as follows:
INTEGER, PARAMETER :: MAXDIAG = MAX_DIAG INTEGER, PARAMETER :: MAXTRACER = MAX_TRACER
and also to define MAX_TRACER (in the CMN_DIAG header file) as follows:
INTEGER, PARAMETER :: MAX_TRACER = NNPAR+6 ! For non-TOMAS simulations
The NNPAR+6 is required for compatibility with the ND09 diagnostic. Therefore, now whenever NNPAR is updated, the values of MAXTRACER and MAXDIAG in gamap_mod.f will also be updated.
--Bob Y. 14:35, 3 August 2010 (EDT)
Bug fix for ND42 diagnostic and ISORROPIA II
Eric Sofen reported a minor bug fix for the ND42 diagnostic in isoropiaII_mod.f. The following line of code:
PH_SAV = -999d0
should have been:
PH_SAV(I,J,L) = -999d0
This only affects ND42 diagnostic output and not the computation of the ATE diagnostic.
--Bob Y. 14:42, 3 August 2010 (EDT)
Minor fix for AOD calculation
Aaron van Donkelaar wrote:
- I've come across a minor bug in the calculation of AOD within aerosol_mod.f. At around line 950, the lines:
! Loop over relative humidity bins DO R = 1, NRH ! Wet radius in "jv_spec.dat" RW(R) = RAA(4,IND(N)+R-1) ! Wet frac of aerosol FWET = (RW(R)**3 - RW(1)**3) / RW(R)**3 ! Extinction efficiency Q for each RH bin QW(R) = QAA(4,IND(N)+R-1)*FWET + QAA(4,IND(N))*(1.d0-FWET) ENDDO
- should be replaced with:
! Loop over relative humidity bins DO R = 1, NRH ! Wet radius in "jv_spec.dat" RW(R) = RAA(4,IND(N)+R-1) ! Extinction efficiency Q for each RH bin QW(R) = QAA(4,IND(N)+R-1) ENDDO
- As a result, the variable FWET is no longer used and no longer needs to be declared in routine RDAER.
- Basically the extinction efficiency for each RH bin is currently being calculated as the water volume-weighted average of each bin's RH and dry value, when the value corresponding to each bin's RH should simply be directly used. The implications should not be large (Colette estimated between 4-8%), but this ought to be corrected.
--Bob Y. 14:48, 3 August 2010 (EDT)