C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/ocean_inversion_project/code_bombC14/TR1.h,v 1.2 2006/05/25 17:16:56 dimitri Exp $
C $Name:  $
C
CBOP
C    !ROUTINE: TR1.h
C    !INTERFACE:
C    include TR1.h
C    !DESCRIPTION: \bv
C     *==========================================================*
C     | TR1.h                                                     
C     | o Tracer 1 prognostic state.                              
C     *==========================================================*
C     | The value and two levels of time tendency are held for    
C     | each prognostic variable.                                 
C     | This file is hardcoded for bombC14 experiment.
C     *==========================================================*
C     \ev
CEOP
C
C     TR1     :: tracer concentration (tr per unit volume).
C     GTR1    :: tracer concentration (TR1 per second).
C     GTR1NM1 :: tracer concentration (TR1 per second).

      COMMON /TR1_R/ tr1, gTr1, gTr1NM1, surfaceTendencyTr1
      _RL  tr1     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL  gTr1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL  gTr1NM1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RS  surfaceTendencyTr1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

      COMMON /TR1_S/ fice, fice0, fice1,
     &                xkw,  xkw0,  xkw1,
     &               patm, patm0, patm1,
     &                WOAsss,  WOAsss0,  WOAsss1,
     &                WOAsst,  WOAsst0,  WOAsst1,
     &      SurfDicPreind, SurfDicAnthro, SurfNatD14C
      _RL  fice    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  fice0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  fice1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  xkw     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  xkw0    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  xkw1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  patm    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  patm0   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  patm1   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsss     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsss0    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsss1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsst     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsst0    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  WOAsst1    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL SurfDicPreind (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL SurfDicAnthro (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL SurfNatD14C (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

      COMMON /TR1_C/ FiceFile,XkwFile,PatmFile,WOAsssFile,WOAsstFile,
     &       SurfDicPreindFile, SurfDicAnthroFile, SurfNatD14CFile
      CHARACTER*(MAX_LEN_FNAM)
     &        FiceFile, XkwFile, PatmFile, WOAsssFile, WOAsstFile,
     &       SurfDicPreindFile, SurfDicAnthroFile, SurfNatD14CFile

C--   Parameters that specify size of input files for
C     FiceFile, XkwFile, PatmFile, WOAsssFile, and WOAsstFile
      INTEGER    nlon     , nlat
      PARAMETER( nlon=360 , nlat=180 )
      real*8     lon0     , lon_inc     , lat0
      PARAMETER( lon0=0.5 , lon_inc=1.0 , lat0=-89.5 )
      real*8 lat_inc(nlat)
      DATA   lat_inc /1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
     &                1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1/

      COMMON /TR1_DC14/ 
     &     DC14_atm_year     , DC14_atm_n     ,
     &     DC14_atm_eq       , DC14_atm_s
      _RL  DC14_atm_year(240), DC14_atm_n(240)
      _RL  DC14_atm_eq(240)  , DC14_atm_s(240)

      COMMON /TR1_CO2/ CO2_atm_year     , CO2_atm
      _RL              CO2_atm_year(478), CO2_atm(478)

      COMMON /TR1_INIT/ monthly_atm_co2_init
      _RL               monthly_atm_co2_init