/[MITgcm]/MITgcm/verification/advect_cs/code/ini_vel.F
ViewVC logotype

Diff of /MITgcm/verification/advect_cs/code/ini_vel.F

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

revision 1.1 by adcroft, Fri Sep 28 02:30:09 2001 UTC revision 1.2 by jmc, Tue Feb 11 02:29:54 2003 UTC
# Line 22  C     myThid -  Number of this instance Line 22  C     myThid -  Number of this instance
22  C     == Local variables ==  C     == Local variables ==
23  C     bi,bj,i,j,k - Loop counters  C     bi,bj,i,j,k - Loop counters
24        INTEGER bi,bj,i,j,k        INTEGER bi,bj,i,j,k
25          INTEGER ip1,jp1
26        _RL psi,omegaprime,fac        _RL psi,omegaprime,fac
27    
28        psi(i,j,bi,bj)=fac*fCoriG(i,j,bi,bj)        psi(i,j,bi,bj)=fac*fCoriG(i,j,bi,bj)
# Line 39  C Line 40  C
40         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
41          DO k=1,Nr          DO k=1,Nr
42           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
43              jp1=MIN(j+1,sNy+Oly)
44            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
45               ip1=MIN(i+1,sNx+Olx)
46             uVel(i,j,k,bi,bj)=0.             uVel(i,j,k,bi,bj)=0.
47       &      +(psi(i,j,bi,bj)-psi(i,j+1,bi,bj))*recip_dyg(i,j,bi,bj)       &      +(psi(i,j,bi,bj)-psi(i,jp1,bi,bj))*recip_dyg(i,j,bi,bj)
48             vVel(i,j,k,bi,bj)=0.             vVel(i,j,k,bi,bj)=0.
49       &      +(psi(i+1,j,bi,bj)-psi(i,j,bi,bj))*recip_dxg(i,j,bi,bj)       &      +(psi(ip1,j,bi,bj)-psi(i,j,bi,bj))*recip_dxg(i,j,bi,bj)
50             wVel(i,j,k,bi,bj)=0.             wVel(i,j,k,bi,bj)=0.
51            ENDDO            ENDDO
52           ENDDO           ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22