/[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.1 by spk, Mon Apr 18 19:42:56 2005 UTC revision 1.2 by spk, Wed Apr 20 18:05:08 2005 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7        SUBROUTINE MATRIX_INIT( iFirstTime, myThid )        SUBROUTINE MATRIX_INIT( myThid )
8    
9        IMPLICIT NONE        IMPLICIT NONE
10    
# Line 15  C $Name$ Line 15  C $Name$
15  #include "MATRIX.h"  #include "MATRIX.h"
16  #include "PARAMS.h"  #include "PARAMS.h"
17                
18        INTEGER iFirstTime, myThid        INTEGER myThid
19                
20  #ifdef ALLOW_MATRIX  #ifdef ALLOW_MATRIX
21    
# Line 23  C     local variables: Line 23  C     local variables:
23        INTEGER i, j, k, bi, bj        INTEGER i, j, k, bi, bj
24        INTEGER iTracer        INTEGER iTracer
25                
26        IF (iFirstTime.EQ.1) THEN  C     Initialize counters
27           expMatrixCounter=0            expMatrixCounter=0
28           impMatrixCounter=0        impMatrixCounter=0
29           DO iTracer=1,PTRACERS_num            DO iTracer=1,PTRACERS_num
30              DO bj=myByLo(myThid), myByHi (myThid)                  DO bj=myByLo(myThid), myByHi (myThid)
31                 DO bi=myBxLo (myThid), myBxHi (myThid)                    DO bi=myBxLo (myThid), myBxHi (myThid)
32                    DO k=1,Nr                          DO k=1,Nr
33                       DO j=1-Oly, sNy+Oly                            DO j=1-Oly, sNy+Oly
34                          DO i=1-Olx, sNx+Olx                          DO i=1-Olx, sNx+Olx
35                             PTRACERS_initial(i,j,k,bi,bj,iTracer)=                                    PTRACERS_initial(i,j,k,bi,bj,iTracer)=
36       &                          pTracer(i,j,k,bi,bj,iTracer)       &                  pTracer(i,j,k,bi,bj,iTracer)
37                             MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0                                    MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
38                             MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0                                    MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
39                          ENDDO                                  ENDDO
40                       ENDDO                            ENDDO
41                    ENDDO                          ENDDO
42                 ENDDO                    ENDDO
43              ENDDO              ENDDO
44           ENDDO            ENDDO
       ENDIF  
45    
46  C      DO iTracer=1,PTRACERS_numInUse  C      DO iTracer=1,PTRACERS_numInUse
47  C         DO bj=myByLo(myThid), myByHi (myThid)  C         DO bj=myByLo(myThid), myByHi (myThid)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22