114 |
ENDDO |
ENDDO |
115 |
DO J = 1,sNy |
DO J = 1,sNy |
116 |
DO I = 1,sNx |
DO I = 1,sNx |
117 |
C First FX = (d/dx)*sigam |
C First FX = (d/dx)*sigma |
118 |
C + d/dx[ eta+zeta d/dx ] U |
C + d/dx[ eta+zeta d/dx ] U |
119 |
FX = _recip_dxC(I,J,bi,bj) * |
FX = _recip_dxC(I,J,bi,bj) * |
120 |
& ( etaPlusZeta(I ,J) * dUdx(I ,J) |
& ( etaPlusZeta(I ,J) * dUdx(I ,J) |
170 |
FX = _maskW(I,J,1,bi,bj) * ( FX - _recip_dxC(I,J,bi,bj) |
FX = _maskW(I,J,1,bi,bj) * ( FX - _recip_dxC(I,J,bi,bj) |
171 |
& * ( press(I,J,bi,bj) - press(I-1,J,bi,bj) ) ) |
& * ( press(I,J,bi,bj) - press(I-1,J,bi,bj) ) ) |
172 |
C |
C |
173 |
C then FY = (d/dy)*sigam |
C then FY = (d/dy)*sigma |
174 |
C + d/dy [(eta+zeta) d/dy] V |
C + d/dy [(eta+zeta) d/dy] V |
175 |
FY = _recip_dyC(I,J,bi,bj) * |
FY = _recip_dyC(I,J,bi,bj) * |
176 |
& ( dVdy(I,J ) * etaPlusZeta(I,J ) |
& ( dVdy(I,J ) * etaPlusZeta(I,J ) |
242 |
C stress and internal ice stresses to surface layer of ocean |
C stress and internal ice stresses to surface layer of ocean |
243 |
areaW = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj)) |
areaW = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj)) |
244 |
areaS = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj)) |
areaS = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj)) |
245 |
fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*(fuIce + FX) |
fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIce + FX |
246 |
fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*(fvIce + FY) |
fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIce + FY |
247 |
END DO |
END DO |
248 |
END DO |
END DO |
249 |
ENDDO |
ENDDO |