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

From Geos-chem
Jump to: navigation, search
(New page: == Overview == BETA RELEASE December 2008 Contains everything in GEOS-Chem v8-02-01 plus the following updates: #Non-local Planetary Boundary Layer scheme for GEOS-5 only (Jintai Lin...)
 
 
(48 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 
== Overview ==
 
== Overview ==
BETA RELEASE December 2008
+
BETA RELEASE -- RELEASE DATE 08 Jun 2009
  
Contains everything in [[GEOS-Chem v8-02-01]] plus the following updates:
+
*Previous version: [[GEOS-Chem v8-02-01]]
 +
*Next version: [[GEOS-Chem v8-02-03]]
  
#Non-local Planetary Boundary Layer scheme for GEOS-5 only (Jintai Lin).
+
=== What's new in this version ===
#BOND inventory for BC and OC with monthly variations.
+
  
With the non-local PBL scheme, the emissions and dry deposition are not uniformly distributed in the PBL anymore. All surface emissions and dry deposition are considered as fluxes in the PBL code and not as reactions in the chemistry. Only airborne emissions (lightning and aircraft) are kept as reactions for the chemistry.
+
Contains everything in [[GEOS-Chem v8-02-01|v8-02-01]], plus the following updates:
 +
 
 +
#[http://acmg.seas.harvard.edu/publications/jlin5_2009_AE_1.pdf Non-local Planetary Boundary Layer scheme] for GEOS-5 only, developed by Jintai Lin.<br>Note: This is an option in <tt>input.geos</tt>
 +
#BOND inventory for BC and OC with monthly variations, prepared by [http://people.seas.harvard.edu/~eleibens/ Eric Leibensperger].  (Please see [[#Error in the Bond BC/OC emissons|this note]] for more information.)
 +
 
 +
With the non-local PBL scheme, the emissions and dry deposition are not uniformly distributed in the PBL anymore. All surface emissions and dry deposition are considered as fluxes in the PBL code and not as reactions in the chemistry. Only airborne emissions (lightning and aircraft) are kept as reactions for the chemistry (however, their small surface contribution is moved to the PBL scheme code as other emissions). These changes are necessary because the PBL scheme uses the emission and dry deposition fluxes to calculate the non-local terms.
 +
 
 +
=== References ===
 +
 
 +
Lin, J.-T., D. Youn, X.-Z. Liang, and D. J. Wuebbles (2008), ''Global model simulation of summertime U.S. ozone diurnal cycle and its sensitivity to PBL mixing, spatial resolution, and emissions'', <u>Atmos. Environ.</u>, doi:10.1016/j.atmosenv.2008.08.012.
  
 
== Previous issues now resolved in v8-02-02 ==
 
== Previous issues now resolved in v8-02-02 ==
 
The following user-reported bugs have now been corrected in v8-02-02:
 
The following user-reported bugs have now been corrected in v8-02-02:
 +
 +
=== Bug with ND52 diagnostic ===
 +
 +
The ND52 diagnostic for gamma(HO2) doesn't work if you ask to output a subset of levels. You need to add a check on the levels when calculating AD52.
 +
You need to change the following lines in calcrate.f:
 +
 +
              IF ( ND52 > 0 ) THEN
 +
                  ! Archive gamma HO2 in AD52
 +
                  AD52(IX,IY,IZ) =
 +
    &                AD52(IX,IY,IZ) + DUMMY3(KLOOP)
 +
              ENDIF
 +
 +
by
 +
 +
              IF ( ND52 > 0 ) THEN
 +
                  ! Archive gamma HO2 in AD52
 +
                  IF ( IZ <= LD52 ) THEN
 +
                    AD52(IX,IY,IZ) =
 +
    &                    AD52(IX,IY,IZ) + DUMMY3(KLOOP)
 +
                  ENDIF
 +
              ENDIF
 +
 +
This bug was solved in v8-02-02.
 +
 +
--[[User:Ccarouge|Ccarouge]] 16:51, 16 June 2009 (EDT)
  
 
=== EPA/NEI inventory: reset other species to zero ===
 
=== EPA/NEI inventory: reset other species to zero ===
'''''Havala O. T. Pye (havala@caltech.edu) wrote:'''''
+
 
 +
[[Image:Obsolete.jpg]]
 +
 
 +
<span style="color:red">'''''NOTE: This code was removed from [[GEOS-Chem v10-01]] and higher versions. Emissions in these versions are now handled by [[HEMCO]].'''''</span>
 +
 
 +
'''''[[User:Havala|Havala O. T. Pye]] wrote:'''''
 +
 
 
:Hi all,
 
:Hi all,
 +
 
:I've noticed a potential problem with how the EPA NEI99 inventory is used in EMFOSSIL.
 
:I've noticed a potential problem with how the EPA NEI99 inventory is used in EMFOSSIL.
:
+
 
 
:My understanding is that ANTHROEMS reads in anthropogenic emissions and puts them in the EMISR array. If a regional inventory is available, EMFOSSIL should override the EMISR values and update EMISRR.
 
:My understanding is that ANTHROEMS reads in anthropogenic emissions and puts them in the EMISR array. If a regional inventory is available, EMFOSSIL should override the EMISR values and update EMISRR.
:
+
 
 
:Emissions for ALD2 and many of the new dicarbonyl species are not in the EPA inventory so they get reset to zero over the U.S. I suggest changing the default value (lines 1269 and 1315, epa_nei_mod.f) in GET_EPA_ANTHRO to -1d0 and using a check like the other regional inventories in emfossil.f (starting line 574, excerpt below).
 
:Emissions for ALD2 and many of the new dicarbonyl species are not in the EPA inventory so they get reset to zero over the U.S. I suggest changing the default value (lines 1269 and 1315, epa_nei_mod.f) in GET_EPA_ANTHRO to -1d0 and using a check like the other regional inventories in emfossil.f (starting line 574, excerpt below).
  
Line 49: Line 90:
 
           ENDIF
 
           ENDIF
  
--[[User:Ccarouge|Ccarouge]] 11:50, 29 May 2009 (EDT)
+
--[[User:Ccarouge|Ccarouge]] 11:50, 29 May 2009 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)
  
 
=== Scale factor for oceanic acetone for GEOS5 2x2.5 ===
 
=== Scale factor for oceanic acetone for GEOS5 2x2.5 ===
 +
 +
[[Image:Obsolete.jpg]]
 +
 +
<span style="color:red">'''''NOTE: Module <tt>acetone_mod.F</tt> was removed from [[GEOS-Chem v10-01]] and higher versions. Emissions in these versions are now handled by [[HEMCO]].'''''</span>
 +
 
It seems the scale factor for the oceanic source of acetone was never defined for GEOS5 2x2.5. Old versions of the code would use a value of 1. But with version 8-02-01, the scale factor is not defined and the compilation would choke on it.
 
It seems the scale factor for the oceanic source of acetone was never defined for GEOS5 2x2.5. Old versions of the code would use a value of 1. But with version 8-02-01, the scale factor is not defined and the compilation would choke on it.
  
Line 95: Line 141:
 
This value for the SCALE_FACTOR will be added to the next version.
 
This value for the SCALE_FACTOR will be added to the next version.
  
--[[User:Ccarouge|Ccarouge]] 14:38, 28 May 2009 (EDT)
+
--[[User:Ccarouge|Ccarouge]] 14:38, 28 May 2009 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)
  
 
=== Bug with PRIVATE declaration in sulfate_mod.f ===
 
=== Bug with PRIVATE declaration in sulfate_mod.f ===
  
In SRCNH3 (in <tt>sulfate_mod.f</tt>), please replace the PRIVATE declaration line 4939 :
+
[[Image:Obsolete.jpg]]
 +
 
 +
<span style="color:red">'''''NOTE: This code was removed from [[GEOS-Chem v10-01]] and higher versions. Emissions in these versions are now handled by [[HEMCO]].'''''</span>
 +
 
 +
In SRCNH3 (in <tt>sulfate_mod.f</tt>), please replace the PRIVATE declaration line 4939 (in v8-02-01) :
 
  !$OMP+PRIVATE( I, J )
 
  !$OMP+PRIVATE( I, J )
  
Line 106: Line 156:
 
  !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF, CAC_AN )
 
  !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF, CAC_AN )
  
and in SRCSO4 please change the PRIVATE declaration in line 4690 from:
+
and in SRCSO4 please change the PRIVATE declaration in line 4690 (in v8-02-01) from:
  
 
   !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF )
 
   !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF )
Line 116: Line 166:
 
In both cases the <tt>CAC_AN</tt> variable has to be added to the PRIVATE declaration.
 
In both cases the <tt>CAC_AN</tt> variable has to be added to the PRIVATE declaration.
  
--[[User:Ccarouge|Ccarouge]] 10:22, 28 May 2009 (EDT)
+
--[[User:Ccarouge|Ccarouge]] 10:22, 28 May 2009 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)
  
 
=== Bug with online 2ndary aerosol ===
 
=== Bug with online 2ndary aerosol ===
Line 132: Line 182:
  
 
=== Bug for dust in ND48 ===
 
=== Bug for dust in ND48 ===
'''''Colette L. Heald (heald@atmos.colostate.edu) wrote:'''''
+
 
 +
<span style="color:darkorange">'''''NOTE: We are phasing out binary punch diagnostics output, in favor of netCDF output.'''''</span>
 +
 
 +
'''''[[User:Cheald|Colette L. Heald]] wrote:'''''
 +
 
 
:I was making changes in the AOD calculation in my version of the code and noticed that in diag48 the dust AOD is not properly scaled from 1000 nm to :400 nm.  I believe line 543 (in v8-01-04) should be changed from:
 
:I was making changes in the AOD calculation in my version of the code and noticed that in diag48 the dust AOD is not properly scaled from 1000 nm to :400 nm.  I believe line 543 (in v8-01-04) should be changed from:
 
                   Q(L) = ODMDUST(I,J,L,R)
 
                   Q(L) = ODMDUST(I,J,L,R)
Line 139: Line 193:
  
 
--[[User:Ccarouge|Ccarouge]] 12:07, 4 June 2009 (EDT)
 
--[[User:Ccarouge|Ccarouge]] 12:07, 4 June 2009 (EDT)
 +
 +
=== Typo in HO2 uptake ===
 +
[[User:Jmao|Jingqiu Mao]] found a typo in the function calculating HO2 uptake by aerosols. This function was added in v8-02-01 and thus the typo only affects v8-02-01 and v8-02-02.
 +
 +
In calcrate.f, line 1338 (v8-02-01) or line 1423 (v8-02-02) there is a misplaced space. Please change:
 +
            kaq = ( k1 * (A1 - 1.d0) + k 2) / (A1**2)
 +
by
 +
            kaq = ( k1 * (A1 - 1.d0) + k2) / (A1**2)
 +
 +
--[[User:Ccarouge|Ccarouge]] 17:41, 9 June 2009 (EDT)
 +
 +
== Outstanding issues not yet resolved in v8-02-02 ==
 +
 +
=== Tracers apparently omitted in ND38, ND39, ND44 diagnostics ===
 +
 +
<span style="color:darkorange">'''''NOTE: We are phasing out binary punch diagnostics output, in favor of netCDF output.'''''</span>
 +
 +
'''''[[User:Psk9|Prasad Kasibhatla]] wrote:'''''
 +
 +
:I am having an apparent problem with the ND38, ND39, and ND44 diag output in GC v8-02-02. I run fullchem simulation with thefollowing set in input.geos:
 +
 +
    ND38: Cld Conv scav loss: 30  all
 +
    ND39: Wetdep scav loss  : 30  all
 +
    ND44: Drydep flx/vel    :  1  all
 +
 +
:When I ran with these setting in v8-01-01, the <tt>ctm.bpch</tt> file contained ND38 and ND39 output for tracer numbers 7, 8, 20, 24, and 26-43. Now in v8-02-02 the output is being written out only for tracer numbers 7,8, 20, and 24.
 +
 +
:Similarly, ND44 output in v8-01-01 was being written out for tracers 1, 2, 3, 7, 8, 15, 16, 17, 20, 22, 26-43. Now in v8-02-02, output for tracers 29-43 in not being written out.
 +
 +
:Any ideas as to what might have changed?
 +
 +
'''''[[User:Bmy|Bob Yantosca]] replied:'''''
 +
 +
:The quick fix is to just list all the tracer #'s explicitly for the ND38, ND39, ND43 diagnostics:
 +
 +
    ND38: Cld Conv scav loss: 47  7 8 20 24 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
 +
    ND39: Wetdep scav loss  : 47  7 8 20 24 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
 +
    ...
 +
    ND43: Chem OH,NO,HO2,NO2:  1  1 2 3 7 8 15 16 17 20 22 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
 +
 +
:The <tt>diag3.f</tt> file was modified (by Claire Carouge) to just print out tracer #'s that are expliclitly requested for certain diagnostics.  I believe that this was in response to [[GEOS-Chem_v8-02-01#Output_selected_tracers_only_for_dry_and_wet_deposition.2C_convection_and_rain_diagnotics| another bug that was reported]] recently.
 +
 +
--[[User:Bmy|Bob Y.]] 10:32, 23 June 2009 (EDT)
 +
 +
'''''[[User:Ccarouge|Claire Carouge]] replied:'''''
 +
:For all diagnostics, when using the 'all' keyword in input.geos, GEOS-Chem "replaces" it by the N first consecutive tracers for the run, N being the maximum number of tracers for a given diagnostic. So for any diagnostic where 'all' means a subset of non-consecutive tracers, the 'all' keyword shouldn't work and the list of tracers should be given instead.
 +
 +
--[[User:Ccarouge|Ccarouge]] 10:18, 10 July 2009 (EDT)
 +
 +
=== Missing NOx data in S.E.-Asia ===
 +
 +
[[Image:Obsolete.jpg]]
 +
 +
<span style="color:green">'''''NOTE: This has now been fixed in [[GEOS-Chem v8-02-03]].'''''</span>
 +
 +
<span style="color:red">'''''NOTE: routine <tt>streets_anthro_mod.F</tt> was removed from [[GEOS-Chem v10-01]] and higher versions. Emissions in these versions are now handled by [[HEMCO]].'''''</span>
 +
 +
Data outside China from the Streets NOx inventory have been missing. The Streets_NOx_FF_2004_monthly.generic.1x1 files should be used only to provide monthly variations. Please update your code with the patch on the ftp site:
 +
  ftp ftp.as.harvard.edu
 +
  cd pub/geos-chem/patches/v8-02-02
 +
  get streets_anthro_mod.f
 +
 +
--[[User:Bmy|Bob Y.]] 16:22, 14 October 2009 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)
 +
 +
=== Mis-calculation of Courant numbers in tpcore_fvdas_mod.f90 ===
 +
 +
<span style="color:green">'''''NOTE: This has now been fixed in [[GEOS-Chem v8-02-03]].'''''</span>
 +
 +
I found that the Courant numbers where not calculated for the second band of latitude and the second last one in tpcore_fvdas_mod.f90. So crap values were used instead at these locations. Please update your code with the patch on the ftp site:
 +
  ftp ftp.as.harvard.edu
 +
  cd pub/geos-chem/patches/v8-02-02
 +
  get tpcore_fvdas_mod.f90
 +
 +
=== Error in the Bond BC/OC emissons ===
 +
 +
[[Image:Obsolete.jpg]]
 +
 +
<span style="color:red">'''''NOTE: All emissions are handled by [[HEMCO]] in [[GEOS-Chem v10-01]] and higher versions.'''''</span>
 +
 +
An [[GEOS-Chem_v8-02-03#Corrected_Bond_et_al_BC.2FOC_emissions|error in the Bond et al BC/OC emissions inventory files]] (in the <tt>carbon_200905</tt> data directories) was subsequently discovered.  The best thing to do is to upgrade to [[GEOS-Chem v8-02-03]].
 +
 +
--[[User:Bmy|Bob Y.]] 16:21, 14 October 2009 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)
 +
 +
=== Omission in Makefile.sparc ===
 +
 +
[[Image:Obsolete.jpg]]
 +
 +
<span style="color:red">'''''NOTE: This makefile has been removed from GEOS-Chem.'''''</span>
 +
 +
'''''[[User:Cheald|Colette Heald]] wrote:'''''
 +
 +
:The <tt>Makefile.sparc</tt> for v8-02-02 is incomplete.  It's missing references to the new vdiff routines:
 +
 +
    ...
 +
    drydep_mod.o                  \
 +
    vdiff_pre_mod.o              \
 +
    vdiff_mod.o                  \
 +
    ...
 +
 +
    uvalbedo_mod.o                : uvalbedo_mod.f                CMN_SIZE define.h
 +
    vdiff_pre_mod.o                : vdiff_pre_mod.f                CMN_SIZE comode.h CMN_O3 define.h
 +
    vdiff_mod.o                    : vdiff_mod.f90                define.h
 +
        $(F90) -c -xtypemap=real:64 $*.f90
 +
 +
NOTE: <tt>Makefile.sparc</tt> and the other platform-specific makefiles are now obsolete.  We suggest migrating to [[GEOS-Chem v8-02-03]], which has a totally rewritten makefile structure.
 +
 +
--[[User:Bmy|Bob Y.]] 09:48, 10 November 2009 (EST)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:42, 15 November 2017 (UTC)

Latest revision as of 17:42, 15 November 2017

Overview

BETA RELEASE -- RELEASE DATE 08 Jun 2009

What's new in this version

Contains everything in v8-02-01, plus the following updates:

  1. Non-local Planetary Boundary Layer scheme for GEOS-5 only, developed by Jintai Lin.
    Note: This is an option in input.geos
  2. BOND inventory for BC and OC with monthly variations, prepared by Eric Leibensperger. (Please see this note for more information.)

With the non-local PBL scheme, the emissions and dry deposition are not uniformly distributed in the PBL anymore. All surface emissions and dry deposition are considered as fluxes in the PBL code and not as reactions in the chemistry. Only airborne emissions (lightning and aircraft) are kept as reactions for the chemistry (however, their small surface contribution is moved to the PBL scheme code as other emissions). These changes are necessary because the PBL scheme uses the emission and dry deposition fluxes to calculate the non-local terms.

References

Lin, J.-T., D. Youn, X.-Z. Liang, and D. J. Wuebbles (2008), Global model simulation of summertime U.S. ozone diurnal cycle and its sensitivity to PBL mixing, spatial resolution, and emissions, Atmos. Environ., doi:10.1016/j.atmosenv.2008.08.012.

Previous issues now resolved in v8-02-02

The following user-reported bugs have now been corrected in v8-02-02:

Bug with ND52 diagnostic

The ND52 diagnostic for gamma(HO2) doesn't work if you ask to output a subset of levels. You need to add a check on the levels when calculating AD52. You need to change the following lines in calcrate.f:

              IF ( ND52 > 0 ) THEN
                 ! Archive gamma HO2 in AD52
                  AD52(IX,IY,IZ) =
    &                 AD52(IX,IY,IZ) + DUMMY3(KLOOP)
              ENDIF

by

              IF ( ND52 > 0 ) THEN
                 ! Archive gamma HO2 in AD52
                 IF ( IZ <= LD52 ) THEN
                    AD52(IX,IY,IZ) =
    &                    AD52(IX,IY,IZ) + DUMMY3(KLOOP)
                 ENDIF
              ENDIF

This bug was solved in v8-02-02.

--Ccarouge 16:51, 16 June 2009 (EDT)

EPA/NEI inventory: reset other species to zero

Obsolete.jpg

NOTE: This code was removed from GEOS-Chem v10-01 and higher versions. Emissions in these versions are now handled by HEMCO.

Havala O. T. Pye wrote:

Hi all,
I've noticed a potential problem with how the EPA NEI99 inventory is used in EMFOSSIL.
My understanding is that ANTHROEMS reads in anthropogenic emissions and puts them in the EMISR array. If a regional inventory is available, EMFOSSIL should override the EMISR values and update EMISRR.
Emissions for ALD2 and many of the new dicarbonyl species are not in the EPA inventory so they get reset to zero over the U.S. I suggest changing the default value (lines 1269 and 1315, epa_nei_mod.f) in GET_EPA_ANTHRO to -1d0 and using a check like the other regional inventories in emfossil.f (starting line 574, excerpt below).
         !--------------------------------------------------------------
         ! Get CO & Hydrocarbons from EPA/NEI inventory over the USA
         !--------------------------------------------------------------

         ! If we are using EPA/NEI99 emissions ...
         IF ( LNEI99 ) THEN

            ! If we are over the USA ...
            IF ( GET_USA_MASK( I, J ) > 0d0 ) THEN

               ! Get EPA/NEI emissions (and apply time-of-day factor)
               EPA_NEI = GET_EPA_ANTHRO( I, J, NN, WEEKDAY )

               ! hotp fix for species not present (hotp 5/28/09)
               IF ( .not. ( EPA_NEI < 0d0 ) ) THEN

                   ! apply time of day factor
                   EPA_NEI = EPA_NEI * TODX

                   ! Convert from molec/cm2/s to kg/box/timestep in order
                   ! to be in the proper units for EMISRR array
                   EMX(1)  = EPA_NEI * ( DTSRCE * AREA_CM2 )
     &                                / XNUMOL(NN)
               ENDIF

            ENDIF
         ENDIF

--Ccarouge 11:50, 29 May 2009 (EDT)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)

Scale factor for oceanic acetone for GEOS5 2x2.5

Obsolete.jpg

NOTE: Module acetone_mod.F was removed from GEOS-Chem v10-01 and higher versions. Emissions in these versions are now handled by HEMCO.

It seems the scale factor for the oceanic source of acetone was never defined for GEOS5 2x2.5. Old versions of the code would use a value of 1. But with version 8-02-01, the scale factor is not defined and the compilation would choke on it.

About this problem, Daniel wrote (djacob@fas.harvard.edu):

You can do the GEOS-5 scaling simply on the basis of surface area, i.e., 1 for 4x5 and 0.25 for 2x2.5. No need to worry about differences in sea surface T and wind speeds - they will have only a very small effect and this oceanic source is very uncertain anyway.

So if you want to use GEOS-5 with a resolution 2x2.5 you need to change the following lines in acetone_mod.f:

#elif defined( GEOS_5 ) 

#if defined( GRID05x0666 ) && defined ( NESTED_CH )

      ! GEOS-5 0.5 x 0.667, scaled to 4x5 (dan, 11/6/08)
      ! This scale factor produces too little acetone. (tmf, 3/05/09)
      !REAL*8, PARAMETER :: SCALE_FACTOR = 0.0008d0
      REAL*8, PARAMETER :: SCALE_FACTOR = 0.015369d0      

#elif defined( GRID4x5 )

      REAL*8, PARAMETER :: SCALE_FACTOR = 0.9551d0

#endif

by:

#elif defined( GEOS_5 ) 

#if defined( GRID05x0666 ) && defined ( NESTED_CH )

      ! GEOS-5 0.5 x 0.667, scaled to 4x5 (dan, 11/6/08)
      ! This scale factor produces too little acetone. (tmf, 3/05/09)
      !REAL*8, PARAMETER :: SCALE_FACTOR = 0.0008d0
      REAL*8, PARAMETER :: SCALE_FACTOR = 0.015369d0      

#elif defined( GRID4x5 )

      REAL*8, PARAMETER :: SCALE_FACTOR = 0.9551d0

#elif defined( GRID2x25 )

      REAL*8, PARAMETER :: SCALE_FACTOR = 0.25d0

#endif

This value for the SCALE_FACTOR will be added to the next version.

--Ccarouge 14:38, 28 May 2009 (EDT)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)

Bug with PRIVATE declaration in sulfate_mod.f

Obsolete.jpg

NOTE: This code was removed from GEOS-Chem v10-01 and higher versions. Emissions in these versions are now handled by HEMCO.

In SRCNH3 (in sulfate_mod.f), please replace the PRIVATE declaration line 4939 (in v8-02-01) :

!$OMP+PRIVATE( I, J )

with the PRIVATE declaration :

!$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF, CAC_AN )

and in SRCSO4 please change the PRIVATE declaration in line 4690 (in v8-02-01) from:

 !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF )

to

 !$OMP+PRIVATE( I, J, AREA_CM2, EPA_AN, EPA_BF, CAC_AN )

In both cases the CAC_AN variable has to be added to the PRIVATE declaration.

--Ccarouge 10:22, 28 May 2009 (EDT)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)

Bug with online 2ndary aerosol

Replace line 412 of aerosol_mod.f:

    IF ( IDTSOAM ) THEN

with:

    IF ( IDTSOAM > 0 ) THEN

This is a bug in GEOS-Chem v8-02-01 release. However, it is more of a matter of Fortran style. The expression in parentheses in the IF statement should be a logical expression and not an integer expression, as is the case above. Some compilers (e.g. PGI) will always choke on this. However, other compilers (e.g. IFORT) will evaluate the expression as TRUE if the integer IDTSOAM is nonzero. This will be fixed in the following release.

--phs 09:56, 28 May 2009 (EDT)

Bug for dust in ND48

NOTE: We are phasing out binary punch diagnostics output, in favor of netCDF output.

Colette L. Heald wrote:

I was making changes in the AOD calculation in my version of the code and noticed that in diag48 the dust AOD is not properly scaled from 1000 nm to :400 nm. I believe line 543 (in v8-01-04) should be changed from:
                 Q(L) = ODMDUST(I,J,L,R)
to:
                 Q(L) = ODMDUST(I,J,L,R) * SCALE400nm

--Ccarouge 12:07, 4 June 2009 (EDT)

Typo in HO2 uptake

Jingqiu Mao found a typo in the function calculating HO2 uptake by aerosols. This function was added in v8-02-01 and thus the typo only affects v8-02-01 and v8-02-02.

In calcrate.f, line 1338 (v8-02-01) or line 1423 (v8-02-02) there is a misplaced space. Please change:

           kaq = ( k1 * (A1 - 1.d0) + k 2) / (A1**2)

by

           kaq = ( k1 * (A1 - 1.d0) + k2) / (A1**2)

--Ccarouge 17:41, 9 June 2009 (EDT)

Outstanding issues not yet resolved in v8-02-02

Tracers apparently omitted in ND38, ND39, ND44 diagnostics

NOTE: We are phasing out binary punch diagnostics output, in favor of netCDF output.

Prasad Kasibhatla wrote:

I am having an apparent problem with the ND38, ND39, and ND44 diag output in GC v8-02-02. I run fullchem simulation with thefollowing set in input.geos:
   ND38: Cld Conv scav loss: 30   all
   ND39: Wetdep scav loss  : 30   all
   ND44: Drydep flx/vel    :  1   all
When I ran with these setting in v8-01-01, the ctm.bpch file contained ND38 and ND39 output for tracer numbers 7, 8, 20, 24, and 26-43. Now in v8-02-02 the output is being written out only for tracer numbers 7,8, 20, and 24.
Similarly, ND44 output in v8-01-01 was being written out for tracers 1, 2, 3, 7, 8, 15, 16, 17, 20, 22, 26-43. Now in v8-02-02, output for tracers 29-43 in not being written out.
Any ideas as to what might have changed?

Bob Yantosca replied:

The quick fix is to just list all the tracer #'s explicitly for the ND38, ND39, ND43 diagnostics:
   ND38: Cld Conv scav loss: 47   7 8 20 24 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
   ND39: Wetdep scav loss  : 47   7 8 20 24 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
   ...
   ND43: Chem OH,NO,HO2,NO2:  1   1 2 3 7 8 15 16 17 20 22 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43
The diag3.f file was modified (by Claire Carouge) to just print out tracer #'s that are expliclitly requested for certain diagnostics. I believe that this was in response to another bug that was reported recently.

--Bob Y. 10:32, 23 June 2009 (EDT)

Claire Carouge replied:

For all diagnostics, when using the 'all' keyword in input.geos, GEOS-Chem "replaces" it by the N first consecutive tracers for the run, N being the maximum number of tracers for a given diagnostic. So for any diagnostic where 'all' means a subset of non-consecutive tracers, the 'all' keyword shouldn't work and the list of tracers should be given instead.

--Ccarouge 10:18, 10 July 2009 (EDT)

Missing NOx data in S.E.-Asia

Obsolete.jpg

NOTE: This has now been fixed in GEOS-Chem v8-02-03.

NOTE: routine streets_anthro_mod.F was removed from GEOS-Chem v10-01 and higher versions. Emissions in these versions are now handled by HEMCO.

Data outside China from the Streets NOx inventory have been missing. The Streets_NOx_FF_2004_monthly.generic.1x1 files should be used only to provide monthly variations. Please update your code with the patch on the ftp site:

  ftp ftp.as.harvard.edu
  cd pub/geos-chem/patches/v8-02-02
  get streets_anthro_mod.f

--Bob Y. 16:22, 14 October 2009 (EDT)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)

Mis-calculation of Courant numbers in tpcore_fvdas_mod.f90

NOTE: This has now been fixed in GEOS-Chem v8-02-03.

I found that the Courant numbers where not calculated for the second band of latitude and the second last one in tpcore_fvdas_mod.f90. So crap values were used instead at these locations. Please update your code with the patch on the ftp site:

  ftp ftp.as.harvard.edu
  cd pub/geos-chem/patches/v8-02-02
  get tpcore_fvdas_mod.f90

Error in the Bond BC/OC emissons

Obsolete.jpg

NOTE: All emissions are handled by HEMCO in GEOS-Chem v10-01 and higher versions.

An error in the Bond et al BC/OC emissions inventory files (in the carbon_200905 data directories) was subsequently discovered. The best thing to do is to upgrade to GEOS-Chem v8-02-03.

--Bob Y. 16:21, 14 October 2009 (EDT)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)

Omission in Makefile.sparc

Obsolete.jpg

NOTE: This makefile has been removed from GEOS-Chem.

Colette Heald wrote:

The Makefile.sparc for v8-02-02 is incomplete. It's missing references to the new vdiff routines:
   ...
   drydep_mod.o                  \
   vdiff_pre_mod.o               \
   vdiff_mod.o                   \
   ...

   uvalbedo_mod.o                 : uvalbedo_mod.f                 CMN_SIZE define.h
   vdiff_pre_mod.o                : vdiff_pre_mod.f                 CMN_SIZE comode.h CMN_O3 define.h
   vdiff_mod.o                    : vdiff_mod.f90                define.h
        $(F90) -c -xtypemap=real:64 $*.f90 

NOTE: Makefile.sparc and the other platform-specific makefiles are now obsolete. We suggest migrating to GEOS-Chem v8-02-03, which has a totally rewritten makefile structure.

--Bob Y. 09:48, 10 November 2009 (EST)
--Bob Yantosca (talk) 17:42, 15 November 2017 (UTC)