C $Header: /u/gcmpack/MITgcm/pkg/matrix/matrix_init.F,v 1.2 2005/04/20 18:05:08 spk Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

      SUBROUTINE MATRIX_INIT( myThid )

      IMPLICIT NONE
 
#include "EEPARAMS.h"
#include "SIZE.h"
#include "PTRACERS_SIZE.h"
#include "PTRACERS.h"
#include "MATRIX.h"
#include "PARAMS.h"
      
      INTEGER myThid
      
#ifdef ALLOW_MATRIX

C     local variables:
      INTEGER i, j, k, bi, bj
      INTEGER iTracer 
      
C     Initialize counters
	  expMatrixCounter=0
      impMatrixCounter=0
 	  DO iTracer=1,PTRACERS_num
	 	DO bj=myByLo(myThid), myByHi (myThid)
		  DO bi=myBxLo (myThid), myBxHi (myThid)
		  	DO k=1,Nr
		  	  DO j=1-Oly, sNy+Oly
		        DO i=1-Olx, sNx+Olx
		  		  PTRACERS_initial(i,j,k,bi,bj,iTracer)= 
     &                  pTracer(i,j,k,bi,bj,iTracer)
				  MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
				  MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
				ENDDO
			  ENDDO
			ENDDO
		  ENDDO
	    ENDDO
	  ENDDO

C      DO iTracer=1,PTRACERS_numInUse
C         DO bj=myByLo(myThid), myByHi (myThid)
C            DO bi=myBxLo (myThid), myBxHi (myThid)
C               DO k=1,Nr
C                  DO j=1-Oly, sNy+Oly
C                     DO i=1-Olx, sNx+Olx
C                        MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
C                        MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0  
C                     ENDDO
C                  ENDDO
C               ENDDO
C            ENDDO
C         ENDDO
C      ENDDO
      
#endif /* ALLOW_MATRIX */
      RETURN 
      END