/[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.1 by spk, Mon Apr 18 19:40:31 2005 UTC revision 1.2 by edhill, Tue Apr 19 18:44:50 2005 UTC
# Line 31  C to file. Line 31  C to file.
31        DATA expMatrixWriteCount /0/        DATA expMatrixWriteCount /0/
32        DATA impMatrixWriteCount /0/        DATA impMatrixWriteCount /0/
33    
34        IF ((mod(myTime-startTime,expMatrixWriteTime).EQ.(0.0 _d 0))) THEN                IF ( (mod(myTime-startTime,expMatrixWriteTime)
35          recipExpDeltaTtracer = (1. _d 0)/(expMatrixCounter*dTtracerLev(1))       &     .EQ. (0.0 _d 0))) THEN
36            recipExpDeltaTtracer =
37         &       (1. _d 0)/(expMatrixCounter*dTtracerLev(1))
38          IF (expMatrixWriteCount.EQ.0) expMatrixWriteCount=1          IF (expMatrixWriteCount.EQ.0) expMatrixWriteCount=1
39          iRec=expMatrixWriteCount          iRec=expMatrixWriteCount
40          DO iTracer=1,PTRACERS_numInUse          DO iTracer=1,PTRACERS_numInUse
# Line 42  C to file. Line 44  C to file.
44                  DO j=1-Oly, sNy+Oly                  DO j=1-Oly, sNy+Oly
45                    DO i=1-Olx, sNx+Olx                    DO i=1-Olx, sNx+Olx
46                      MATRIX(i,j,k,bi,bj,iTracer,1) =                      MATRIX(i,j,k,bi,bj,iTracer,1) =
47       &                   MATRIX(i,j,k,bi,bj,iTracer,1)*recipExpDeltaTtracer       &                   MATRIX(i,j,k,bi,bj,iTracer,1)
48         &                   *recipExpDeltaTtracer
49                    ENDDO                    ENDDO
50                  ENDDO                  ENDDO
51                ENDDO                ENDDO
52              ENDDO              ENDDO
53            ENDDO            ENDDO
54            WRITE(suff,'(A9,I2.2)') 'MATRIXEXP',iTracer            WRITE(suff,'(A9,I2.2)') 'MATRIXEXP',iTracer
55            write(*,*)'Writing explicit matrix :',iTracer, expMatrixWriteCount, expMatrixCounter            write(*,*)'Writing explicit matrix :',iTracer,
56         &         expMatrixWriteCount, expMatrixCounter
57            CALL WRITE_REC_XYZ_RL(suff,            CALL WRITE_REC_XYZ_RL(suff,
58       &         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)
59          ENDDO          ENDDO
# Line 71  C       Reset explicit matrix Line 75  C       Reset explicit matrix
75          ENDDO          ENDDO
76        ENDIF        ENDIF
77                    
78        IF ((mod(myTime-startTime,impMatrixWriteTime).EQ.(0.0 _d 0))) THEN                IF ( (mod(myTime-startTime,impMatrixWriteTime)
79         &     .EQ.(0.0 _d 0)) ) THEN
80          recipImpMatrixCounter = (1. _d 0)/dble(impMatrixCounter)          recipImpMatrixCounter = (1. _d 0)/dble(impMatrixCounter)
81          IF (impMatrixWriteCount.EQ.0) impMatrixWriteCount=1          IF (impMatrixWriteCount.EQ.0) impMatrixWriteCount=1
82          iRec=impMatrixWriteCount          iRec=impMatrixWriteCount
# Line 82  C       Reset explicit matrix Line 87  C       Reset explicit matrix
87                  DO j=1-Oly, sNy+Oly                  DO j=1-Oly, sNy+Oly
88                    DO i=1-Olx, sNx+Olx                    DO i=1-Olx, sNx+Olx
89                      MATRIX(i,j,k,bi,bj,iTracer,2) =                      MATRIX(i,j,k,bi,bj,iTracer,2) =
90       &                   MATRIX(i,j,k,bi,bj,iTracer,2)*recipImpMatrixCounter       &                   MATRIX(i,j,k,bi,bj,iTracer,2)
91         &                   *recipImpMatrixCounter
92                    ENDDO                    ENDDO
93                  ENDDO                  ENDDO
94                ENDDO                ENDDO
95              ENDDO              ENDDO
96            ENDDO            ENDDO
97            WRITE(suff,'(A9,I2.2)') 'MATRIXIMP',iTracer            WRITE(suff,'(A9,I2.2)') 'MATRIXIMP',iTracer
98            write(*,*)'Writing implicit matrix :',iTracer, impMatrixWriteCount, impMatrixCounter            write(*,*)'Writing implicit matrix :',iTracer,
99         &         impMatrixWriteCount, impMatrixCounter
100            CALL WRITE_REC_XYZ_RL(suff,            CALL WRITE_REC_XYZ_RL(suff,
101       &         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)
102          ENDDO            ENDDO  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22