31 |
C k :: vertical level |
C k :: vertical level |
32 |
C wFld :: vertical flow |
C wFld :: vertical flow |
33 |
C myThid :: thread number |
C myThid :: thread number |
34 |
INTEGER bi,bj,K |
INTEGER bi,bj,k |
35 |
_RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
36 |
INTEGER myThid |
INTEGER myThid |
37 |
|
|
55 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
56 |
vCoriolisTerm(i,j) = |
vCoriolisTerm(i,j) = |
57 |
& -0.5*( fCoriCos( i ,j,bi,bj)*angleSinC(i, j ,bi,bj) |
& -0.5*( fCoriCos( i ,j,bi,bj)*angleSinC(i, j ,bi,bj) |
58 |
& *0.5*( wFld( i ,j,k,bi,bj)+wFld( i ,j,kp1,bi,bj)*wMsk ) |
& *0.5*( wFld( i ,j, k ,bi,bj)*rVel2wUnit( k ) |
59 |
|
& +wFld( i ,j,kp1,bi,bj)*rVel2wUnit(kp1)*wMsk ) |
60 |
& + fCoriCos(i,j-1,bi,bj)*angleSinC(i,j-1,bi,bj) |
& + fCoriCos(i,j-1,bi,bj)*angleSinC(i,j-1,bi,bj) |
61 |
& *0.5*( wFld(i,j-1,k,bi,bj)+wFld(i,j-1,kp1,bi,bj)*wMsk ) ) |
& *0.5*( wFld(i,j-1, k ,bi,bj)*rVel2wUnit( k ) |
62 |
& *gravitySign*recip_horiVertRatio |
& +wFld(i,j-1,kp1,bi,bj)*rVel2wUnit(kp1)*wMsk ) |
63 |
|
& )*gravitySign |
64 |
ENDDO |
ENDDO |
65 |
ENDDO |
ENDDO |
66 |
|
|