/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_fill_state.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_fill_state.F

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

revision 1.22 by molod, Mon Jul 11 21:58:19 2005 UTC revision 1.23 by jmc, Wed Jul 13 00:46:07 2005 UTC
# Line 42  C     == Local variables == Line 42  C     == Local variables ==
42        _RL tmpMk(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL tmpMk(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43        _RL tmpMk1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1,nSx,nSy)        _RL tmpMk1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1,nSx,nSy)
44        _RL tmp1k(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL tmp1k(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45        _RL tmpFac        _RL tmpFac, uBarC, vBarC
46        INTEGER i,j,K,bi,bj        INTEGER i,j,K,bi,bj
47        INTEGER km1        INTEGER km1
48                
# Line 119  C---+----1----+----2----+----3----+----4 Line 119  C---+----1----+----2----+----3----+----4
119             DO K=1,Nr             DO K=1,Nr
120              DO j = 1,sNy              DO j = 1,sNy
121               DO i = 1,sNx               DO i = 1,sNx
122                tmpMk(i,j,K,bi,bj) = 0.25 _d 0                uBarC = 0.5 _d 0
123       &           *(uVel(i,j,K,bi,bj)+uVel(i+1,j,K,bi,bj))       &           *(uVel(i,j,K,bi,bj)+uVel(i+1,j,K,bi,bj))
124                  vBarC = 0.5 _d 0
125       &           *(vVel(i,j,K,bi,bj)+vVel(i,j+1,K,bi,bj))       &           *(vVel(i,j,K,bi,bj)+vVel(i,j+1,K,bi,bj))
126                  tmpMk(i,j,K,bi,bj) =
127         &            ( angleCosC(i,j,bi,bj)*uBarC
128         &             -angleSinC(i,j,bi,bj)*vBarC )
129         &           *( angleSinC(i,j,bi,bj)*uBarC
130         &             +angleCosC(i,j,bi,bj)*vBarC )
131               ENDDO               ENDDO
132              ENDDO              ENDDO
133             ENDDO             ENDDO
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22