/[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.9 by mlosch, Wed Jun 14 16:08:36 2006 UTC revision 1.10 by mlosch, Thu Mar 8 11:21:34 2007 UTC
# Line 225  C     - (d/dy) P/2 Line 225  C     - (d/dy) P/2
225  C      C    
226  C     recompute wind stress over ice (done already in seaice_dynsolver,  C     recompute wind stress over ice (done already in seaice_dynsolver,
227  C     but not saved)  C     but not saved)
228    #ifdef ALLOW_ATM_WIND          
229             fuIce = 0.5 _d 0 *             fuIce = 0.5 _d 0 *
230       &          ( DAIRN(I  ,J,bi,bj)*(       &          ( DAIRN(I  ,J,bi,bj)*(
231       &          COSWIN*uWind(I  ,J,bi,bj)       &          COSWIN*uWind(I  ,J,bi,bj)
# Line 241  C     but not saved) Line 242  C     but not saved)
242       &          SIGN(SINWIN, _fCori(I,J-1,bi,bj))*uWind(I,J-1,bi,bj)       &          SIGN(SINWIN, _fCori(I,J-1,bi,bj))*uWind(I,J-1,bi,bj)
243       &          +COSWIN*vWind(I,J-1,bi,bj) )       &          +COSWIN*vWind(I,J-1,bi,bj) )
244       &          )       &          )
245    #else
246               fuIce=DAIRN(I,J,bi,bj)*(
247         &         COSWIN*FU(I,J,bi,bj)
248         &         -SIGN(SINWIN, _fCori(I  ,J,bi,bj))*0.25
249         &           *(FV(I,J,  bi,bj)+FV(I-1,J,  bi,bj)
250         &            +FV(I,J+1,bi,bj)+FV(I-1,J+1,bi,bj))
251         &         )
252               fvIce=DAIRN(I,J,bi,bj)*(
253         &         SIGN(SINWIN, _fCori(I,J  ,bi,bj))*0.25
254         &         *(FU(I,J  ,bi,bj)+FU(I+1,J,  bi,bj)
255         &          +FU(I,J-1,bi,bj)+FU(I+1,J-1,bi,bj))
256         &         +COSWIN*FV(I,J,bi,bj)
257         &         )
258    #endif /* ALLOW_ATM_WIND */
259  C     average wind stress over ice and ocean and apply averaged wind  C     average wind stress over ice and ocean and apply averaged wind
260  C     stress and internal ice stresses to surface layer of ocean  C     stress and internal ice stresses to surface layer of ocean
261             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))
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22