/[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.21 by dimitri, Thu Jan 17 23:18:39 2008 UTC revision 1.22 by mlosch, Wed Mar 18 10:26:10 2009 UTC
# Line 119  C     evaluate divergence of stress and Line 119  C     evaluate divergence of stress and
119       &           + sig12(I  ,J+1) * _dxV(I  ,J+1,bi,bj)       &           + sig12(I  ,J+1) * _dxV(I  ,J+1,bi,bj)
120       &           - sig12(I  ,J  ) * _dxV(I  ,J  ,bi,bj)       &           - sig12(I  ,J  ) * _dxV(I  ,J  ,bi,bj)
121       &           ) * recip_rAw(I,J,bi,bj)       &           ) * recip_rAw(I,J,bi,bj)
      &           -  
      &           ( sig12(I,J) + sig12(I,J+1) )  
      &           * _tanPhiAtU(I,J,bi,bj) * recip_rSphere  
      &           +  
      &           ( sig22(I,J) + sig22(I-1,J) ) * 0.5 _d 0  
      &           * _tanPhiAtU(I,J,bi,bj) * recip_rSphere  
 C     one metric term  missing for general curvilinear coordinates  
122              FY = ( sig22(I  ,J  ) * _dxF(I  ,J  ,bi,bj)              FY = ( sig22(I  ,J  ) * _dxF(I  ,J  ,bi,bj)
123       &           - sig22(I  ,J-1) * _dxF(I  ,J-1,bi,bj)       &           - sig22(I  ,J-1) * _dxF(I  ,J-1,bi,bj)
124       &           + sig12(I+1,J  ) * _dyU(I+1,J  ,bi,bj)       &           + sig12(I+1,J  ) * _dyU(I+1,J  ,bi,bj)
125       &           - sig12(I  ,J  ) * _dyU(I  ,J  ,bi,bj)       &           - sig12(I  ,J  ) * _dyU(I  ,J  ,bi,bj)
126       &           ) * recip_rAs(I,J,bi,bj)       &           ) * recip_rAs(I,J,bi,bj)
      &           -  
      &           ( sig22(I,J) + sig22(I,J-1) ) * 0.5 _d 0  
      &           * _tanPhiAtV(I,J,bi,bj) * recip_rSphere  
 C     two metric terms missing for general curvilinear coordinates  
127  C     average wind stress over ice and ocean and apply averaged wind  C     average wind stress over ice and ocean and apply averaged wind
128  C     stress and internal ice stresses to surface layer of ocean  C     stress and internal ice stresses to surface layer of ocean
129              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.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22