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) |
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 |