/[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.24 by mlosch, Fri May 29 14:51:21 2009 UTC revision 1.25 by mlosch, Wed Jun 24 08:01:43 2009 UTC
# Line 59  C Line 59  C
59            DO I=1,sNx            DO I=1,sNx
60  C     average wind stress over ice and ocean and apply averaged wind  C     average wind stress over ice and ocean and apply averaged wind
61  C     stress and internal ice stresses to surface layer of ocean  C     stress and internal ice stresses to surface layer of ocean
62             areaW = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj))             areaW = 0.5 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))
63       &          * SEAICEstressFactor       &          * SEAICEstressFactor
64             areaS = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj))             areaS = 0.5 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))
65       &          * SEAICEstressFactor       &          * 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)
# Line 86  C     ice cover Line 86  C     ice cover
86           DO i=1,sNx           DO i=1,sNx
87            fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) )*            fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) )*
88       &         COSWAT *       &         COSWAT *
89       &         ( UICE(I,J,1,bi,bj)-uVel(I,J,1,bi,bj) )       &         ( UICE(I,J,bi,bj)-uVel(I,J,1,bi,bj) )
90       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *
91       &         ( DWATN(I  ,J,bi,bj) *       &         ( DWATN(I  ,J,bi,bj) *
92       &         0.5 _d 0*(vIce(I  ,J  ,1,bi,bj)-vVel(I  ,J  ,1,bi,bj)       &         0.5 _d 0*(vIce(I  ,J  ,bi,bj)-vVel(I  ,J  ,1,bi,bj)
93       &                  +vIce(I  ,J+1,1,bi,bj)-vVel(I  ,J+1,1,bi,bj))       &                  +vIce(I  ,J+1,bi,bj)-vVel(I  ,J+1,1,bi,bj))
94       &         + DWATN(I-1,J,bi,bj) *       &         + DWATN(I-1,J,bi,bj) *
95       &         0.5 _d 0*(vIce(I-1,J  ,1,bi,bj)-vVel(I-1,J  ,1,bi,bj)       &         0.5 _d 0*(vIce(I-1,J  ,bi,bj)-vVel(I-1,J  ,1,bi,bj)
96       &                  +vIce(I-1,J+1,1,bi,bj)-vVel(I-1,J+1,1,bi,bj))       &                  +vIce(I-1,J+1,bi,bj)-vVel(I-1,J+1,1,bi,bj))
97       &         )       &         )
98            fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) )*            fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) )*
99       &         COSWAT *       &         COSWAT *
100       &         ( VICE(I,J,1,bi,bj)-vVel(I,J,1,bi,bj) )       &         ( VICE(I,J,bi,bj)-vVel(I,J,1,bi,bj) )
101       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *
102       &         ( DWATN(I,J  ,bi,bj) *       &         ( DWATN(I,J  ,bi,bj) *
103       &         0.5 _d 0*(uIce(I  ,J  ,1,bi,bj)-uVel(I  ,J  ,1,bi,bj)       &         0.5 _d 0*(uIce(I  ,J  ,bi,bj)-uVel(I  ,J  ,1,bi,bj)
104       &                  +uIce(I+1,J  ,1,bi,bj)-uVel(I+1,J  ,1,bi,bj))       &                  +uIce(I+1,J  ,bi,bj)-uVel(I+1,J  ,1,bi,bj))
105       &         + DWATN(I,J-1,bi,bj) *       &         + DWATN(I,J-1,bi,bj) *
106       &         0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-uVel(I  ,J-1,1,bi,bj)       &         0.5 _d 0*(uIce(I  ,J-1,bi,bj)-uVel(I  ,J-1,1,bi,bj)
107       &                  +uIce(I+1,J-1,1,bi,bj)-uVel(I+1,J-1,1,bi,bj))       &                  +uIce(I+1,J-1,bi,bj)-uVel(I+1,J-1,1,bi,bj))
108       &         )       &         )
109            areaW = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj))            areaW = 0.5 _d 0 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))
110       &         * SEAICEstressFactor       &         * SEAICEstressFactor
111            areaS = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj))            areaS = 0.5 _d 0 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))
112       &         * SEAICEstressFactor       &         * SEAICEstressFactor
113            fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIceLoc            fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIceLoc
114            fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIceLoc            fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIceLoc
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22