/[MITgcm]/MITgcm/pkg/seaice/ostres.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/ostres.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.21 by mlosch, Wed Jun 24 08:23:00 2009 UTC revision 1.22 by jmc, Fri Oct 29 15:26:29 2010 UTC
# Line 5  C $Name$ Line 5  C $Name$
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE ostres( COR_ICE, myThid )        SUBROUTINE ostres( COR_ICE, myThid )
8  C     /==========================================================\  C     *==========================================================*
9  C     | SUBROUTINE ostres                                        |  C     | SUBROUTINE ostres                                        |
10  C     | o Calculate ocean surface stress                         |  C     | o Calculate ocean surface stress                         |
11  C     |==========================================================|  C     *==========================================================*
12  C     \==========================================================/  C     *==========================================================*
13        IMPLICIT NONE        IMPLICIT NONE
14    
15  C     === Global variables ===  C     === Global variables ===
# Line 43  c     introduce turning angle (default i Line 43  c     introduce turning angle (default i
43        COSWAT=COS(SEAICE_waterTurnAngle*deg2rad)        COSWAT=COS(SEAICE_waterTurnAngle*deg2rad)
44    
45  #ifdef SEAICE_ORIGINAL_BAD_ICE_STRESS  #ifdef SEAICE_ORIGINAL_BAD_ICE_STRESS
46        These lines are put here on purpose to cause the compilation        These lines are put here on purpose to cause the compilation
47        to fail because the following code is retired and can no longer        to fail because the following code is retired and can no longer
48        be used. It was probably never used anyway.        be used. It was probably never used anyway.
49  CMLC--   Following formulation is problematic and is no longer used.  CMLC--   Following formulation is problematic and is no longer used.
# Line 116  C--   Compute ice-affected wind stress Line 116  C--   Compute ice-affected wind stress
116          DO j=1,sNy          DO j=1,sNy
117           DO i=1,sNx           DO i=1,sNx
118            fuIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*(            fuIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*(
119       &         COSWAT *       &         COSWAT *
120       &         ( UICE(I,J,  bi,bj)-GWATX(I,J,  bi,bj)       &         ( UICE(I,J,  bi,bj)-GWATX(I,J,  bi,bj)
121       &         + UICE(I,J+1,bi,bj)-GWATX(I,J+1,bi,bj) )       &         + UICE(I,J+1,bi,bj)-GWATX(I,J+1,bi,bj) )
122       &         -SIGN(SINWAT,COR_ICE(I,J,bi,bj)) *       &         -SIGN(SINWAT,COR_ICE(I,J,bi,bj)) *
# Line 124  C--   Compute ice-affected wind stress Line 124  C--   Compute ice-affected wind stress
124       &         + VICE(I+1,J,bi,bj)-GWATY(I+1,J,bi,bj) )       &         + VICE(I+1,J,bi,bj)-GWATY(I+1,J,bi,bj) )
125       &         )       &         )
126            fvIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*(            fvIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*(
127       &         SIGN(SINWAT,COR_ICE(I,J,bi,bj)) *       &         SIGN(SINWAT,COR_ICE(I,J,bi,bj)) *
128       &         ( UICE(I,J,  bi,bj)-GWATX(I,J,  bi,bj)       &         ( UICE(I,J,  bi,bj)-GWATX(I,J,  bi,bj)
129       &         + UICE(I,J+1,bi,bj)-GWATX(I,J+1,bi,bj) )       &         + UICE(I,J+1,bi,bj)-GWATX(I,J+1,bi,bj) )
130       &         + COSWAT *       &         + COSWAT *
131       &         ( VICE(I,  J,i,bj)-GWATY(I,  J,bi,bj)       &         ( VICE(I,  J,bi,bj)-GWATY(I,  J,bi,bj)
132       &         + VICE(I+1,J,bi,bj)-GWATY(I+1,J,bi,bj) )       &         + VICE(I+1,J,bi,bj)-GWATY(I+1,J,bi,bj) )
133       &         )       &         )
134            fu(I,J,bi,bj)=(ONE-AREA(I,J,bi,bj))*fu(I,J,bi,bj)+            fu(I,J,bi,bj)=(ONE-AREA(I,J,bi,bj))*fu(I,J,bi,bj)+
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22