/[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.15 by heimbach, Tue Apr 24 17:54:44 2007 UTC revision 1.16 by mlosch, Tue Apr 24 18:38:15 2007 UTC
# Line 81  C     surface stresses on ocean followin Line 81  C     surface stresses on ocean followin
81  C      C    
82  C     recompute strain rates, viscosities, etc. from updated ice velocities  C     recompute strain rates, viscosities, etc. from updated ice velocities
83         IF ( .NOT. SEAICEuseEVP ) THEN         IF ( .NOT. SEAICEuseEVP ) THEN
84  C     we already have the stress components and do not need to recompute them  C     only for EVP we already have the stress components otherwise we need
85    C     to recompute them here
86          CALL SEAICE_CALC_STRAINRATES(          CALL SEAICE_CALC_STRAINRATES(
87       I       uIce(1-Olx,1-Oly,1,1,1), vIce(1-Olx,1-Oly,1,1,1),       I       uIce(1-Olx,1-Oly,1,1,1), vIce(1-Olx,1-Oly,1,1,1),
88       O       e11, e22, e12,       O       e11, e22, e12,
# Line 159  C     stress and internal ice stresses t Line 160  C     stress and internal ice stresses t
160       &           + areaS*tauy(I,J,bi,bj)       &           + areaS*tauy(I,J,bi,bj)
161       &           + FY * SEAICEstressFactor       &           + FY * SEAICEstressFactor
162  C     save stress divergence for later  C     save stress divergence for later
163  #ifdef ALLOW_EVP  #ifdef SEAICE_ALLOW_EVP
164              stressDivergenceX(I,J,bi,bj) = FX              stressDivergenceX(I,J,bi,bj) = FX
165              stressDivergenceY(I,J,bi,bj) = FY              stressDivergenceY(I,J,bi,bj) = FY
166  #endif  #endif /* SEAICE_ALLOW_EVP */
167             ENDDO             ENDDO
168            ENDDO            ENDDO
169           ELSE           ELSE
170  #ifdef ALLOW_EVP  #ifdef SEAICE_ALLOW_EVP
171            DO J=1,sNy            DO J=1,sNy
172             DO I=1,sNx             DO I=1,sNx
173  C     average wind stress over ice and ocean and apply averaged wind  C     average wind stress over ice and ocean and apply averaged wind
# Line 183  C     stress and internal ice stresses t Line 184  C     stress and internal ice stresses t
184       &           + stressDivergenceY(I,J,bi,bj) * SEAICEstressFactor       &           + stressDivergenceY(I,J,bi,bj) * SEAICEstressFactor
185             ENDDO             ENDDO
186            ENDDO            ENDDO
187  #endif  #endif /* SEAICE_ALLOW_EVP */
188           ENDIF           ENDIF
189          ENDDO          ENDDO
190         ENDDO         ENDDO
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22