/[MITgcm]/MITgcm/pkg/timeave/timeave_cumul_2vfc.F
ViewVC logotype

Diff of /MITgcm/pkg/timeave/timeave_cumul_2vfc.F

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

revision 1.2 by edhill, Thu Oct 9 04:19:20 2003 UTC revision 1.4 by heimbach, Fri Aug 19 22:50:26 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3  #include "TIMEAVE_OPTIONS.h"  #include "CPP_EEOPTIONS.h"
4    
5  CStartofinterface  CStartofinterface
6        SUBROUTINE TIMEAVE_CUMUL_2VFC(        SUBROUTINE TIMEAVE_CUMUL_2VFC(
7       O   fldtave,       O   fldtave,
8       I   fld1, fld2, fld3, Ksize, dir, deltaT,       I   fld1, fld2, fld3, Ksize, dir, deltaTloc,
9       I   bi, bj, myThid )       I   bi, bj, myThid )
10  C     /==========================================================\  C     /==========================================================\
11  C     | SUBROUTINE TIMEAVE_CUMUL_2VFC                            |  C     | SUBROUTINE TIMEAVE_CUMUL_2VFC                            |
# Line 34  C     Ksize - 3rd dimension of local arr Line 34  C     Ksize - 3rd dimension of local arr
34        _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)        _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
35        _RS fld3(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)        _RS fld3(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
36        _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)        _RL fldtave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
37        _RL deltaT        _RL deltaTloc
38        INTEGER bi, bj, myThid        INTEGER bi, bj, myThid
39    
40  CEndofinterface  CEndofinterface
# Line 50  C      DO bi = myBxLo(myThid), myBxHi(my Line 50  C      DO bi = myBxLo(myThid), myBxHi(my
50          DO k=1,Ksize          DO k=1,Ksize
51           DO j=1,sNy           DO j=1,sNy
52            DO i=1,sNx            DO i=1,sNx
53              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT*              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaTloc*
54       &         fld1(i,j,k,bi,bj)*fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)       &         fld1(i,j,k,bi,bj)*fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)
55             ENDDO             ENDDO
56            ENDDO            ENDDO
# Line 65  C      DO bi = myBxLo(myThid), myBxHi(my Line 65  C      DO bi = myBxLo(myThid), myBxHi(my
65          DO k=1,Ksize          DO k=1,Ksize
66           DO j=1,sNy           DO j=1,sNy
67            DO i=1,sNx            DO i=1,sNx
68              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT*              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaTloc*
69       &         0.5*( fld1(i-1,j,k,bi,bj) + fld1(i,j,k,bi,bj) )       &         0.5*( fld1(i-1,j,k,bi,bj) + fld1(i,j,k,bi,bj) )
70       &            *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)       &            *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)
71             ENDDO             ENDDO
# Line 81  C      DO bi = myBxLo(myThid), myBxHi(my Line 81  C      DO bi = myBxLo(myThid), myBxHi(my
81          DO k=1,Ksize          DO k=1,Ksize
82           DO j=1,sNy           DO j=1,sNy
83            DO i=1,sNx            DO i=1,sNx
84              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaT*              fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj) + deltaTloc*
85       &         0.5*( fld1(i,j-1,k,bi,bj) + fld1(i,j,k,bi,bj) )       &         0.5*( fld1(i,j-1,k,bi,bj) + fld1(i,j,k,bi,bj) )
86       &            *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)       &            *fld2(i,j,k,bi,bj)*fld3(i,j,k,bi,bj)
87             ENDDO             ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22