56 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
57 |
DO j=1,sNy |
DO j=1,sNy |
58 |
DO i=1,sNx |
DO i=1,sNx |
59 |
distY = (yC(i,j,bi,bj)-(yC0))/lY |
c distY = (yC(i,j,bi,bj)-(yC0))/lY |
60 |
C tauX = -tauMax*cos(2. _d 0*PI*distY) |
c jmc: yC0 has been removed (not used except here, |
61 |
|
c and since it was not initialised, was generaly zero): |
62 |
|
distY = (yC(i,j,bi,bj)- 0. )/lY |
63 |
|
c tauX = -tauMax*cos(2. _d 0*PI*distY) |
64 |
tauX = tauMax*sin(PI*distY) |
tauX = tauMax*sin(PI*distY) |
65 |
fu(i,j,bi,bj) = tauX |
fu(i,j,bi,bj) = tauX |
66 |
ENDDO |
ENDDO |