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

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

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

revision 1.19 by jmc, Tue Mar 6 16:45:20 2012 UTC revision 1.22 by mlosch, Fri May 26 09:08:32 2017 UTC
# Line 7  C $Name$ Line 7  C $Name$
7  #else  #else
8  # define OBCS_UVICE_OLD  # define OBCS_UVICE_OLD
9  #endif  #endif
10    #ifdef ALLOW_AUTODIFF
11    # include "AUTODIFF_OPTIONS.h"
12    #endif
13    
14  CBOP  CBOP
15  C     !ROUTINE: SEAICE_CALC_STRAINRATES  C     !ROUTINE: SEAICE_CALC_STRAINRATES
# Line 84  C     vectorizes even after TAFization Line 87  C     vectorizes even after TAFization
87        k = 1        k = 1
88        noSlipFac = 0. _d 0        noSlipFac = 0. _d 0
89        IF ( SEAICE_no_slip ) noSlipFac = 1. _d 0        IF ( SEAICE_no_slip ) noSlipFac = 1. _d 0
90    C     in order repoduce results before fixing a bug in r1.20 comment out
91    C     the following line
92    CML      IF ( SEAICE_no_slip ) noSlipFac = 2. _d 0
93  C  C
94        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
95         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 148  C     evaluate strain rates at Z-points Line 154  C     evaluate strain rates at Z-points
154       &         )       &         )
155       &         *maskC(i  ,j  ,k,bi,bj)*maskC(i-1,j  ,k,bi,bj)       &         *maskC(i  ,j  ,k,bi,bj)*maskC(i-1,j  ,k,bi,bj)
156       &         *maskC(i  ,j-1,k,bi,bj)*maskC(i-1,j-1,k,bi,bj)       &         *maskC(i  ,j-1,k,bi,bj)*maskC(i-1,j-1,k,bi,bj)
157       &         + 2.0 _d 0 * noSlipFac * (       &         + noSlipFac * (
158       &           2.0 _d 0 * uave(i,j) * _recip_dyU(i,j,bi,bj) * hFacU       &           2.0 _d 0 * uave(i,j) * _recip_dyU(i,j,bi,bj) * hFacU
159       &         + 2.0 _d 0 * vave(i,j) * _recip_dxV(i,j,bi,bj) * hFacV       &         + 2.0 _d 0 * vave(i,j) * _recip_dxV(i,j,bi,bj) * hFacV
160       &         )       &         )
# Line 159  c$$$     &         - hFacV * k1AtZ(i,j,b Line 165  c$$$     &         - hFacV * k1AtZ(i,j,b
165  c$$$     &         - hFacU * k2AtZ(i,j,bi,bj) * uave(i,j)  c$$$     &         - hFacU * k2AtZ(i,j,bi,bj) * uave(i,j)
166           ENDDO           ENDDO
167          ENDDO          ENDDO
   
168         ENDDO         ENDDO
169        ENDDO        ENDDO
170    
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22