pkg/matrix/matrix_init_varia.F

5d0c92356b Jean*0001 #include "MATRIX_OPTIONS.h"
a9544d2f65 Sama*0002
ff8fe3d29f Jean*0003       SUBROUTINE MATRIX_INIT_VARIA( myThid )
a9544d2f65 Sama*0004
                0005       IMPLICIT NONE
8f6dc1243d Jean*0006
a9544d2f65 Sama*0007 #include "EEPARAMS.h"
                0008 #include "SIZE.h"
                0009 #include "PTRACERS_SIZE.h"
8f6dc1243d Jean*0010 #include "PTRACERS_PARAMS.h"
                0011 #include "PTRACERS_FIELDS.h"
a9544d2f65 Sama*0012 #include "MATRIX.h"
                0013 #include "PARAMS.h"
8f6dc1243d Jean*0014
57d1ea6458 Sama*0015       INTEGER myThid
8f6dc1243d Jean*0016
a9544d2f65 Sama*0017 #ifdef ALLOW_MATRIX
                0018
                0019 C     local variables:
                0020       INTEGER i, j, k, bi, bj
8f6dc1243d Jean*0021       INTEGER iTracer
                0022
57d1ea6458 Sama*0023 C     Initialize counters
5d0c92356b Jean*0024       expMatrixWriteCount = 1
                0025       impMatrixWriteCount = 1
                0026       DO bj=myByLo(myThid),myByHi(myThid)
                0027        DO bi=myBxLo(myThid),myBxHi(myThid)
                0028          expMatrixCounter(bi,bj) = 0
                0029          impMatrixCounter(bi,bj) = 0
                0030        ENDDO
                0031       ENDDO
8f6dc1243d Jean*0032       DO iTracer=1,PTRACERS_num
5d0c92356b Jean*0033         DO bj=myByLo(myThid),myByHi(myThid)
                0034          DO bi=myBxLo(myThid),myBxHi(myThid)
                0035            DO k=1,Nr
                0036             DO j=1-Oly, sNy+Oly
                0037              DO i=1-Olx, sNx+Olx
                0038                PTRACERS_initial(i,j,k,bi,bj,iTracer) =
57d1ea6458 Sama*0039      &                  pTracer(i,j,k,bi,bj,iTracer)
5d0c92356b Jean*0040                MATRIX(i,j,k,bi,bj,iTracer,1) = 0. _d 0
                0041                MATRIX(i,j,k,bi,bj,iTracer,2) = 0. _d 0
                0042              ENDDO
8f6dc1243d Jean*0043             ENDDO
5d0c92356b Jean*0044            ENDDO
                0045          ENDDO
8f6dc1243d Jean*0046         ENDDO
                0047       ENDDO
a9544d2f65 Sama*0048
                0049 #endif /* ALLOW_MATRIX */
8f6dc1243d Jean*0050       RETURN
a9544d2f65 Sama*0051       END