/[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.3 by mlosch, Mon Mar 6 20:53:25 2006 UTC revision 1.4 by mlosch, Thu Mar 9 20:22:40 2006 UTC
# Line 39  C     i,j,bi,bj - Loop counters Line 39  C     i,j,bi,bj - Loop counters
39        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
40        _RL  SINWAT, COSWAT        _RL  SINWAT, COSWAT
41        _RL  fuIce, fvIce        _RL  fuIce, fvIce
42          _RL  areaW, areaS
43    
44  c     introduce turning angle (default is zero)  c     introduce turning angle (default is zero)
45        SINWAT=SIN(SEAICE_waterTurnAngle*deg2rad)        SINWAT=SIN(SEAICE_waterTurnAngle*deg2rad)
# Line 85  C--   Compute ice-affected wind stress ( Line 86  C--   Compute ice-affected wind stress (
86       &         +0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-GWATX(I  ,J-1,bi,bj)       &         +0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-GWATX(I  ,J-1,bi,bj)
87       &                   +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)) )
88       &         )       &         )
89            fu(I,J,bi,bj)=(ONE-AREA(I,J,1,bi,bj))*fu(I,J,bi,bj)+            areaW = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj))
90       &         AREA(I,J,1,bi,bj)*fuIce            areaS = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj))
91            fv(I,J,bi,bj)=(ONE-AREA(I,J,1,bi,bj))*fv(I,J,bi,bj)+            fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIce
92       &         AREA(I,J,1,bi,bj)*fvIce            fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIce
93           ENDDO           ENDDO
94          ENDDO          ENDDO
95         ENDDO         ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22