/[MITgcm]/MITgcm/pkg/matrix/matrix_init.F
ViewVC logotype

Diff of /MITgcm/pkg/matrix/matrix_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by spk, Wed Apr 20 18:05:08 2005 UTC revision 1.3 by jmc, Mon Nov 5 18:58:00 2007 UTC
# Line 7  C $Name$ Line 7  C $Name$
7        SUBROUTINE MATRIX_INIT( myThid )        SUBROUTINE MATRIX_INIT( myThid )
8    
9        IMPLICIT NONE        IMPLICIT NONE
10    
11  #include "EEPARAMS.h"  #include "EEPARAMS.h"
12  #include "SIZE.h"  #include "SIZE.h"
13  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
14  #include "PTRACERS.h"  #include "PTRACERS_PARAMS.h"
15    #include "PTRACERS_FIELDS.h"
16  #include "MATRIX.h"  #include "MATRIX.h"
17  #include "PARAMS.h"  #include "PARAMS.h"
18          
19        INTEGER myThid        INTEGER myThid
20          
21  #ifdef ALLOW_MATRIX  #ifdef ALLOW_MATRIX
22    
23  C     local variables:  C     local variables:
24        INTEGER i, j, k, bi, bj        INTEGER i, j, k, bi, bj
25        INTEGER iTracer        INTEGER iTracer
26          
27  C     Initialize counters  C     Initialize counters
28            expMatrixCounter=0        expMatrixCounter=0
29        impMatrixCounter=0        impMatrixCounter=0
30            DO iTracer=1,PTRACERS_num        DO iTracer=1,PTRACERS_num
31                  DO bj=myByLo(myThid), myByHi (myThid)          DO bj=myByLo(myThid), myByHi (myThid)
32                    DO bi=myBxLo (myThid), myBxHi (myThid)            DO bi=myBxLo (myThid), myBxHi (myThid)
33                          DO k=1,Nr              DO k=1,Nr
34                            DO j=1-Oly, sNy+Oly                DO j=1-Oly, sNy+Oly
35                          DO i=1-Olx, sNx+Olx                  DO i=1-Olx, sNx+Olx
36                                    PTRACERS_initial(i,j,k,bi,bj,iTracer)=                      PTRACERS_initial(i,j,k,bi,bj,iTracer)=
37       &                  pTracer(i,j,k,bi,bj,iTracer)       &                  pTracer(i,j,k,bi,bj,iTracer)
38                                    MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0                    MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
39                                    MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0                    MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
40                                  ENDDO                  ENDDO
41                            ENDDO                ENDDO
42                          ENDDO              ENDDO
43                    ENDDO            ENDDO
44              ENDDO          ENDDO
45            ENDDO        ENDDO
46    
47  C      DO iTracer=1,PTRACERS_numInUse  C      DO iTracer=1,PTRACERS_numInUse
48  C         DO bj=myByLo(myThid), myByHi (myThid)  C         DO bj=myByLo(myThid), myByHi (myThid)
# Line 50  C               DO k=1,Nr Line 51  C               DO k=1,Nr
51  C                  DO j=1-Oly, sNy+Oly  C                  DO j=1-Oly, sNy+Oly
52  C                     DO i=1-Olx, sNx+Olx  C                     DO i=1-Olx, sNx+Olx
53  C                        MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0  C                        MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
54  C                        MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0    C                        MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
55  C                     ENDDO  C                     ENDDO
56  C                  ENDDO  C                  ENDDO
57  C               ENDDO  C               ENDDO
58  C            ENDDO  C            ENDDO
59  C         ENDDO  C         ENDDO
60  C      ENDDO  C      ENDDO
61          
62  #endif /* ALLOW_MATRIX */  #endif /* ALLOW_MATRIX */
63        RETURN        RETURN
64        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22