C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/ocean_inversion_project/code_ecco1x1/ptracers_read_takahashi.F,v 1.1 2004/10/06 20:26:55 dimitri Exp $
C $Name:  $

#include "PTRACERS_OPTIONS.h"

CBOP
C !ROUTINE: PTRACERS_READ_TAKAHASHI

C !INTERFACE: ==========================================================
      SUBROUTINE PTRACERS_READ_TAKAHASHI( myThid )

C !DESCRIPTION:
C     Initialize pTracerTakahashi, an array that holds the monthly-mean
C     Takahashi CO2 surface flux climatology in units of mol/m^2/yr.

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "PTRACERS.h"

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  none

#ifdef ALLOW_PTRACERS

C !LOCAL VARIABLES: ====================================================
      character*(80) infile
      INTEGER        i, j, bi, bj, icount
      INTEGER        imonth, filePrec, interp_method
      INTEGER        imt    , jmt
      PARAMETER     (imt= 72, jmt=45)
      REAL*4         lon(imt), lat(jmt)
      _RL            arrayout(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nsx,nsy)
      _RL            lon0      , lon_inc     , lat0
      PARAMETER     (lon0=2.5, lon_inc=5.0, lat0=-88.0)
      _RL            lat_inc(jmt-1)
      DATA           lat_inc / 
     &     4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,
     &     4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 /
CEOP

C==== Use exf_interp to read-in Takahashi climatology
      infile        = 'taka02_montlhy.bin'
      filePrec      = 32
      interp_method = 1
      do imonth=1,12
       call exf_interp( infile, filePrec
     &        , arrayout, imonth, xC, yC
     &        , lon0, lon_inc, lat0, lat_inc
     &        , imt, jmt, interp_method, mythid
     &        )
       _EXCH_XY_R4( arrayout,  myThid )
       do bj = mybylo(mythid), mybyhi(mythid)
        do bi = mybxlo(mythid), mybxhi(mythid)
         do j=1-OLy,sny+OLy
          do i=1-OLx,snx+OLx
           pTracerTakahashi(i,j,imonth,bi,bj) = arrayout(i,j,bi,bj)
cdb   write(10,*),imonth,i,j,pTracerTakahashi(i,j,imonth,bi,bj)
          enddo
         enddo
        enddo
       enddo
      enddo

#endif /* ALLOW_PTRACERS */

      RETURN
      END