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

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

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

revision 1.4 by jmc, Fri Aug 26 19:47:00 2011 UTC revision 1.5 by jmc, Mon Apr 2 21:26:12 2012 UTC
# Line 15  C to file. Line 15  C to file.
15  #include "PARAMS.h"  #include "PARAMS.h"
16  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
17  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
 c#include "PTRACERS_FIELDS.h"  
18  #include "MATRIX.h"  #include "MATRIX.h"
19    
20        _RL myTime        _RL myTime
# Line 45  C---+----1----+----2----+----3----+----4 Line 44  C---+----1----+----2----+----3----+----4
44            DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
45             DO bi=myBxLo (myThid),myBxHi(myThid)             DO bi=myBxLo (myThid),myBxHi(myThid)
46               DO k=1,Nr               DO k=1,Nr
47                DO j=1-Oly, sNy+Oly                DO j=1-OLy, sNy+OLy
48                 DO i=1-Olx, sNx+Olx                 DO i=1-OLx, sNx+OLx
49                   MATRIX(i,j,k,bi,bj,iTracer,1) =                   MATRIX(i,j,k,bi,bj,iTracer,1) =
50       &                   MATRIX(i,j,k,bi,bj,iTracer,1)       &                   MATRIX(i,j,k,bi,bj,iTracer,1)
51       &                   *recipExpDeltaTtracer       &                   *recipExpDeltaTtracer
# Line 55  C---+----1----+----2----+----3----+----4 Line 54  C---+----1----+----2----+----3----+----4
54               ENDDO               ENDDO
55             ENDDO             ENDDO
56            ENDDO            ENDDO
57              _BEGIN_MASTER(myThid)
58            bj=myByLo(myThid)            bj=myByLo(myThid)
59            bi=myBxLo(myThid)            bi=myBxLo(myThid)
60            WRITE(msgBuf,'(A,I4,2I6)') 'Writing explicit matrix :',            WRITE(msgBuf,'(A,I4,2I6)') 'Writing explicit matrix :',
61       &       iTracer, expMatrixWriteCount, expMatrixCounter(bi,bj)       &       iTracer, expMatrixWriteCount, expMatrixCounter(bi,bj)
62            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
63       &                        SQUEEZE_RIGHT, myThid )       &                        SQUEEZE_RIGHT, myThid )
64              _END_MASTER(myThid)
65            WRITE(suff,'(A9,I2.2)') 'MATRIXEXP',iTracer            WRITE(suff,'(A9,I2.2)') 'MATRIXEXP',iTracer
66            CALL WRITE_REC_XYZ_RL(suff,            CALL WRITE_REC_XYZ_RL(suff,
67       &         MATRIX(1-Olx,1-Oly,1,1,1,iTracer,1),iRec,myIter,myThid)       &         MATRIX(1-OLx,1-OLy,1,1,1,iTracer,1),iRec,myIter,myThid)
68          ENDDO          ENDDO
69          _BARRIER          _BARRIER
70          _BEGIN_MASTER(myThid)          _BEGIN_MASTER(myThid)
# Line 80  C       Reset explicit matrix Line 81  C       Reset explicit matrix
81            DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
82             DO bi=myBxLo (myThid),myBxHi(myThid)             DO bi=myBxLo (myThid),myBxHi(myThid)
83               DO k=1,Nr               DO k=1,Nr
84                DO j=1-Oly, sNy+Oly                DO j=1-OLy, sNy+OLy
85                 DO i=1-Olx, sNx+Olx                 DO i=1-OLx, sNx+OLx
86                   MATRIX(i,j,k,bi,bj,iTracer,1) = 0. _d 0                   MATRIX(i,j,k,bi,bj,iTracer,1) = 0. _d 0
87                 ENDDO                 ENDDO
88                ENDDO                ENDDO
# Line 102  C       Reset explicit matrix Line 103  C       Reset explicit matrix
103            DO bj=myByLo(myThid), myByHi (myThid)            DO bj=myByLo(myThid), myByHi (myThid)
104             DO bi=myBxLo (myThid), myBxHi (myThid)             DO bi=myBxLo (myThid), myBxHi (myThid)
105               DO k=1,Nr               DO k=1,Nr
106                DO j=1-Oly, sNy+Oly                DO j=1-OLy, sNy+OLy
107                 DO i=1-Olx, sNx+Olx                 DO i=1-OLx, sNx+OLx
108                   MATRIX(i,j,k,bi,bj,iTracer,2) =                   MATRIX(i,j,k,bi,bj,iTracer,2) =
109       &                   MATRIX(i,j,k,bi,bj,iTracer,2)       &                   MATRIX(i,j,k,bi,bj,iTracer,2)
110       &                   *recipImpMatrixCounter       &                   *recipImpMatrixCounter
# Line 112  C       Reset explicit matrix Line 113  C       Reset explicit matrix
113               ENDDO               ENDDO
114             ENDDO             ENDDO
115            ENDDO            ENDDO
116              _BEGIN_MASTER(myThid)
117            bj=myByLo(myThid)            bj=myByLo(myThid)
118            bi=myBxLo(myThid)            bi=myBxLo(myThid)
119            WRITE(msgBuf,'(A,I4,2I6)') 'Writing implicit matrix :',            WRITE(msgBuf,'(A,I4,2I6)') 'Writing implicit matrix :',
120       &       iTracer, impMatrixWriteCount, impMatrixCounter(bi,bj)       &       iTracer, impMatrixWriteCount, impMatrixCounter(bi,bj)
121            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
122       &                        SQUEEZE_RIGHT, myThid )       &                        SQUEEZE_RIGHT, myThid )
123              _END_MASTER(myThid)
124            WRITE(suff,'(A9,I2.2)') 'MATRIXIMP',iTracer            WRITE(suff,'(A9,I2.2)') 'MATRIXIMP',iTracer
125            CALL WRITE_REC_XYZ_RL(suff,            CALL WRITE_REC_XYZ_RL(suff,
126       &         MATRIX(1-Olx,1-Oly,1,1,1,iTracer,2),iRec,myIter,myThid)       &         MATRIX(1-OLx,1-OLy,1,1,1,iTracer,2),iRec,myIter,myThid)
127          ENDDO          ENDDO
128          _BARRIER          _BARRIER
129          _BEGIN_MASTER(myThid)          _BEGIN_MASTER(myThid)
# Line 137  C       Reset implicit matrix Line 140  C       Reset implicit matrix
140            DO bj=myByLo(myThid), myByHi (myThid)            DO bj=myByLo(myThid), myByHi (myThid)
141             DO bi=myBxLo (myThid), myBxHi (myThid)             DO bi=myBxLo (myThid), myBxHi (myThid)
142               DO k=1,Nr               DO k=1,Nr
143                DO j=1-Oly, sNy+Oly                DO j=1-OLy, sNy+OLy
144                 DO i=1-Olx, sNx+Olx                 DO i=1-OLx, sNx+OLx
145                   MATRIX(i,j,k,bi,bj,iTracer,2) = 0. _d 0                   MATRIX(i,j,k,bi,bj,iTracer,2) = 0. _d 0
146                 ENDDO                 ENDDO
147                ENDDO                ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22