Bob Yantosca
02 Oct 2012
GEOS-Chem Support Team
geos-chem-support@g.harvard.edu


NOTE: the Net_terrestrial_exch_5.29Pg.geos.1x1 is misnamed.
It contains data on the GENERIC 1x1 grid (360 lons x 180 lats) instead
of on the GEOS 1x1 grid (360 lons x 181 lats).

Therefore we have created a symbolic link:

   Net_terrestrial_exch_5.29Pg.generic.1x1

with the proper grid name that points to the erroneous file.

We have also modified the code in co2_mod.F (routine: READ_ANNUAL_BIONET_CO2)
to read from this linked file:

      !------------------------------------
      ! Read climatological Bionet data
      !------------------------------------       
      IF ( LBIONETCLIM ) THEN

         ! TAU value for start of "generic" year 2000
         TAU = GET_TAU0( 1, 1, 2000 )

!-----------------------------------------------------------------------------
!%%% v9-01-03 POST-RELEASE PATCH (bmy, 10/2/12)
!%%% NOTE: Someone mislabeled the Net_terrestrial_exch file, it is labeled
!%%% "geos.1x1" but is on the GENERIC 1x1 grid.  Therefore, we shall rename
!%%% this file accordingly and use the generic grid indices (I1x1,J1x1-1)
!%%% etc. below. (bmy, 10/2/12)
!%%%         ! Filename
!%%%         FILENAME = TRIM( DATA_DIR_1x1 )  //  
!%%%     &              'CO2_201003/Net_terrestrial_exch_5.29Pg.geos.1x1'
!%%%	
!%%%         ! Echo info
!%%%         WRITE( 6, 100 ) TRIM( FILENAME )
!%%%
!%%%         ! Read Net Terrestrial CO2 Exchange [molec/cm2/s]
!%%%         CALL READ_BPCH2( FILENAME, 'CO2-SRCE', 6, 
!%%%     &                    TAU,       I1x1,      J1x1      
!%%%     &                    1,         ARRAY,     QUIET=.TRUE. )
!%%%
!%%%         ! Cast to REAL*8 before regridding
!%%%         GEOS_1x1(:,:,1) = ARRAY(:,:,1)
!%%%
!%%%         ! File with lat/lon edges for regridding
!%%%         LLFILENAME = TRIM( DATA_DIR_1x1) //
!%%%     &                'MAP_A2A_Regrid_201203/MAP_A2A_latlon_geos1x1.nc'
!%%%
!%%%         ! Regrid from GEOS 1x1 --> current model resolution
!%%%         INGRID  => GEOS_1x1 (:,:,1)
!%%%         OUTGRID => GEOS_GRID(:,:,1)
!%%%         CALL DO_REGRID_A2A( LLFILENAME, I1x1,    J1x1, 
!%%%     &                       INGRID,     OUTGRID, PERAREA=0,
!%%%     &                       netCDF=.TRUE.                   )
!-----------------------------------------------------------------------------
                  
         ! Filename
         FILENAME = TRIM( DATA_DIR_1x1 )                      //  
     &              'CO2_201003/'                             //
     &              'Net_terrestrial_exch_5.29Pg.generic.1x1'
	
         ! Echo info
         WRITE( 6, 100 ) TRIM( FILENAME )

         ! Read Net Terrestrial CO2 Exchange [molec/cm2/s]
         CALL READ_BPCH2( FILENAME, 'CO2-SRCE', 6, 
     &                    TAU,       I1x1,      J1x1-1,
     &                    1,         ARRAY_1x1, QUIET=.TRUE. )

         ! Cast to REAL*8 before regridding
         GEN_1x1(:,:,1) = ARRAY_1x1(:,:,1)

         ! File with lat/lon edges for regridding
         LLFILENAME = TRIM( DATA_DIR_1x1)            //
     &                'MAP_A2A_Regrid_201203/'       //
     &                'MAP_A2A_latlon_generic1x1.nc'

         ! Regrid from GEOS 1x1 --> current model resolution
         INGRID  => GEN_1x1  (:,:,1)
         OUTGRID => GEOS_GRID(:,:,1)
         CALL DO_REGRID_A2A( LLFILENAME, I1x1,    J1x1-1, 
     &                       INGRID,     OUTGRID, PERAREA=0,
     &                       netCDF=.TRUE.                   )

         ! Free pointers
         NULLIFY( INGRID, OUTGRID )
         
         EMBIONETCO2(:,:) = GEOS_GRID(:,:,1)  !(lmw,05/16/11)     

      ENDIF