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

Revision as of 18:46, 19 October 2009

Overview

BETA RELEASE -- Date TBA

Will contain everything in GEOS-Chem v8-02-03, plus:

• LINOZ stratospheric ozone (Dylan Jones group @ U. Toronto)
• Updated inventory of volcanic SO2 emissions from AeroCom (J. Fisher)

Previous issues now resolved in v8-02-04

Updated reactions in sulfate_mod.f

Jenny Fisher (jafisher@fas.harvard.edu) wrote:

I made two very small changes to sulfate_mod.f that should be included in the standard code. The change updates the reaction rates used for the offline aerosol simulation so that they match the rates Jingqiu implemented in the full chemistry simulation.

The changes are:

1. In CHEM_DMS, definition of RK2.
2. In CHEM_SO2, definitions of K0 and Ki. Note that I also removed the definition of Ki earlier in the subroutine (line 1497) and moved it to this section so it can be easily updated for future rate constant revisions.

--Bob Y. 09:50, 15 October 2009 (EDT)

Bug fix for EMEP ship emissions

Monika Kopacz (mkopacz@princeton.edu) wrote:

I noticed that the code crashes when EMEP is set to false but ship emissions from EMEP are true in the input file. Sure, it makes sense, but I would suggest the following addition to the input_mod.f, around say around line 1550 to prevent the code from crashing.

   ! Make sure we're not using EMEP ship emissions, if EMEP is off
   IF (.not. LEMEP ) THEN
      LEMEPSHIP = .FALSE.
   ENDIF 

Bob Yantosca (yantosca@seas.harvard.edu) replied:

Good catch. We've put this code into READ_EMISSIONS_MENU in input_mod.f:

    !%%% Bug fix!  If LEMEPSHIP is turned on but LEMEP is turned %%%
    !%%% off, this will cause an error (because arrays are not   %%%
    !%%% allocated, etc.).  For now, just turn off LEMEPSHIP     %%%
    !%%% and print a warning message.  Whoever wants to fix this %%%
    !%%% in a more robust way is welcome to do so.               %%%
    !%%% (mak, bmy, 10/19/09)                                    %%%
    IF ( LEMEPSHIP .and. ( .not. LEMEP ) ) THEN
       LEMEPSHIP = .FALSE.
       WRITE( 6, '(a)' ) REPEAT( '=', 79 )
       WRITE( 6, '(a)' ) 'WARNING! EMEP emissions are turned off,'
       WRITE( 6, '(a)' ) 'so also turn off EMEP ship emissions'
       WRITE( 6, '(a)' ) 'in order to avoid crashing GEOS-Chem!'
       WRITE( 6, '(a)' ) REPEAT( '=', 79 )
    ENDIF

As you pointed out the simple solution is to just turn off LEMEPSHIP if LEMEP=F. However, it may still be possible to keep LEMEPSHIP turned on...one would have to go thru emep_mod.f to make sure that the ship emissions arrays get allocated even if the other arrays don't. Whoever wants to can feel free to do that...

--Bob Y. 10:57, 19 October 2009 (EDT)

Minor fix for molecular weights in drydep_mod.f

Bob Yantosca (yantosca@seas.harvard.edu) replied:

I found a minor inconsistency in INIT_DRYDEP in (drydep_mod.f) which only affects the SOA tracers. In the big loop where we initialize everything, the following updates were made:

       ! ALPH (Alpha-pinene)
       ELSE IF ( N == IDTALPH ) THEN
          ...

          !----------------------------------------------------------
          ! Prior to 10/19/09:
          ! Molwt should be 136.23 not 136 even (bmy, 10/19/09)
          !XMW(NUMDEP)     = 136d-3
          !----------------------------------------------------------
          XMW(NUMDEP)     = 136.23d-3

          ...            

       ! LIMO (Limonene)
       ELSE IF ( N == IDTLIMO ) THEN
          ...

          !----------------------------------------------------------
          ! Prior to 10/19/09:
          ! Molwt should be 136.23 not 136 even (bmy, 10/19/09)
          !----------------------------------------------------------
          XMW(NUMDEP)     = 136.23d-3
  

So in other words, the mol wt of ALPH and LIMO was listed as 136d-3 but it should be 136.23d-3. It doesn't change things drastically, just in the last few decimal places of the surface resistance array RSURFC.

--Bob Y. 14:46, 19 October 2009 (EDT)

Outstanding issues not yet resolved in v8-02-04