/[MITgcm]/MITgcm/pkg/seaice/seaice_ocean_stress.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_ocean_stress.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by mlosch, Thu Mar 16 14:25:40 2006 UTC revision 1.7 by mlosch, Fri Mar 17 15:53:38 2006 UTC
# Line 114  C     re-compute internal stresses with Line 114  C     re-compute internal stresses with
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)
# Line 170  C     - (d/dx) P/2 Line 170  C     - (d/dx) P/2
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  )
# Line 242  C     average wind stress over ice and o Line 242  C     average wind stress over ice and o
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
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22