/[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.1 by mlosch, Mon Mar 6 13:17:37 2006 UTC revision 1.2 by mlosch, Mon Mar 6 13:27:40 2006 UTC
# Line 127  C--   Compute ice-affected wind stress ( Line 127  C--   Compute ice-affected wind stress (
127       &         SINWAT *       &         SINWAT *
128       &         ( UICE(I,J,1,bi,bj)-GWATX(I,J,bi,bj) )       &         ( UICE(I,J,1,bi,bj)-GWATX(I,J,bi,bj) )
129       &         + COSWAT * 0.5 _d 0 * (       &         + COSWAT * 0.5 _d 0 * (
130       &          0.5 _d 0*(uIce(I  ,J  ,1,bi,bj)-GWATY(I  ,J  ,bi,bj)       &          0.5 _d 0*(uIce(I  ,J  ,1,bi,bj)-GWATX(I  ,J  ,bi,bj)
131       &                   +uIce(I+1,J  ,1,bi,bj)-GWATX(I+1,J  ,bi,bj))       &                   +uIce(I+1,J  ,1,bi,bj)-GWATX(I+1,J  ,bi,bj))
132       &         +0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-GWATY(I  ,J-1,bi,bj)       &         +0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-GWATX(I  ,J-1,bi,bj)
133       &                   +uIce(I+1,J-1,1,bi,bj)-GWATX(I+1,J-1,bi,bj)) )       &                   +uIce(I+1,J-1,1,bi,bj)-GWATX(I+1,J-1,bi,bj)) )
134       &         )       &         )
135            fu(I,J,bi,bj)=(ONE-AREA(I,J,1,bi,bj))*fu(I,J,bi,bj)+            fu(I,J,bi,bj)=(ONE-AREA(I,J,1,bi,bj))*fu(I,J,bi,bj)+
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22