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

Diff of /MITgcm/pkg/timeave/timeave_cumul_dif_1t.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.3 by heimbach, Fri Aug 19 22:50:26 2005 UTC
# Line 5  C $Name$ Line 5  C $Name$
5  CStartofinterface  CStartofinterface
6        SUBROUTINE TIMEAVE_CUMUL_DIF_1T(        SUBROUTINE TIMEAVE_CUMUL_DIF_1T(
7       O   fldtave,       O   fldtave,
8       I   fld1, fld2, Ksize, dir, deltaT,       I   fld1, fld2, Ksize, dir, deltaTloc,
9       I   bi, bj, myThid )       I   bi, bj, myThid )
10  C     /==========================================================*  C     /==========================================================*
11  C     | SUBROUTINE TIMEAVE_CUMUL_DIF_1T  C     | SUBROUTINE TIMEAVE_CUMUL_DIF_1T
# Line 29  C     Ksize   - 3rd dimension of local a Line 29  C     Ksize   - 3rd dimension of local a
29        _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)        _RL fld1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize,nSx,nSy)
30        _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize)        _RL fld2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Ksize)
31        _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)
32        _RL deltaT        _RL deltaTloc
33        INTEGER bi, bj, myThid        INTEGER bi, bj, myThid
34    
35  CEndofinterface  CEndofinterface
# Line 47  C     i,j,k,bi,bj  - Loop counters Line 47  C     i,j,k,bi,bj  - Loop counters
47             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
48       &        - ( fld1(i,j,k,bi,bj) - fld1(i-1,j,k,bi,bj) )       &        - ( fld1(i,j,k,bi,bj) - fld1(i-1,j,k,bi,bj) )
49       &          * recip_dxC(i,j,bi,bj)       &          * recip_dxC(i,j,bi,bj)
50       &          * fld2(i,j,k) * deltaT       &          * fld2(i,j,k) * deltaTloc
51            ENDDO            ENDDO
52           ENDDO           ENDDO
53          ENDDO          ENDDO
# Line 60  C     i,j,k,bi,bj  - Loop counters Line 60  C     i,j,k,bi,bj  - Loop counters
60             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
61       &        - ( fld1(i,j,k,bi,bj) - fld1(i,j-1,k,bi,bj) )       &        - ( fld1(i,j,k,bi,bj) - fld1(i,j-1,k,bi,bj) )
62       &          * recip_dyC(i,j,bi,bj)       &          * recip_dyC(i,j,bi,bj)
63       &          * fld2(i,j,k) * deltaT       &          * fld2(i,j,k) * deltaTloc
64            ENDDO            ENDDO
65           ENDDO           ENDDO
66          ENDDO          ENDDO
# Line 73  C     i,j,k,bi,bj  - Loop counters Line 73  C     i,j,k,bi,bj  - Loop counters
73             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)             fldtave(i,j,k,bi,bj) = fldtave(i,j,k,bi,bj)
74       &        + ( fld1(i,j,k,bi,bj) - fld1(i,j,k-1,bi,bj) )       &        + ( fld1(i,j,k,bi,bj) - fld1(i,j,k-1,bi,bj) )
75       &          * recip_drC(k)       &          * recip_drC(k)
76       &          * fld2(i,j,k) * deltaT       &          * fld2(i,j,k) * deltaTloc
77            ENDDO            ENDDO
78           ENDDO           ENDDO
79          ENDDO          ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22