59 |
DO j=1,sNy |
DO j=1,sNy |
60 |
DO i=1,sNx |
DO i=1,sNx |
61 |
distY = yC(i,j,bi,bj)/lY |
distY = yC(i,j,bi,bj)/lY |
62 |
C tauX = -tauMax*cos(2. _d 0*PI*distY) |
tauX = -tauMax*cos(2. _d 0*PI*distY) |
63 |
tauX = tauMax*sin(PI*distY) |
c tauX = tauMax*sin(PI*distY) |
64 |
fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil) |
fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil) |
65 |
C fu(i,j,bi,bj) = -tauMax |
C fu(i,j,bi,bj) = -tauMax |
66 |
ENDDO |
ENDDO |
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |
69 |
ENDDO |
ENDDO |
70 |
|
fu(Nx/2,Ny/2,1,1) = fu(Nx/2,Ny/2,1,1)*0.917d0 |
71 |
C |
C |
72 |
_EXCH_XY_R4(fu , myThid ) |
_EXCH_XY_R4(fu , myThid ) |
73 |
_EXCH_XY_R4(fv , myThid ) |
_EXCH_XY_R4(fv , myThid ) |