/[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.27 by jmc, Thu Sep 10 16:05:26 2009 UTC revision 1.28 by mlosch, Wed Nov 24 15:36:31 2010 UTC
# Line 63  C     average wind stress over ice and o Line 63  C     average wind stress over ice and o
63  C     stress and internal ice stresses to surface layer of ocean  C     stress and internal ice stresses to surface layer of ocean
64             areaW = 0.5 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))             areaW = 0.5 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))
65       &          * SEAICEstressFactor       &          * SEAICEstressFactor
            areaS = 0.5 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))  
      &          * SEAICEstressFactor  
66             fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)             fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)
67       &          + areaW*taux(I,J,bi,bj)       &          + areaW*taux(I,J,bi,bj)
68       &          + stressDivergenceX(I,J,bi,bj) * SEAICEstressFactor       &          + stressDivergenceX(I,J,bi,bj) * SEAICEstressFactor
69              ENDDO
70             ENDDO
71             DO J=1,sNy
72              DO I=1,sNx
73               areaS = 0.5 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))
74         &          * SEAICEstressFactor
75             fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)             fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)
76       &          + areaS*tauy(I,J,bi,bj)       &          + areaS*tauy(I,J,bi,bj)
77       &          + stressDivergenceY(I,J,bi,bj) * SEAICEstressFactor       &          + stressDivergenceY(I,J,bi,bj) * SEAICEstressFactor
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22