/[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.20 by mlosch, Thu Feb 28 17:25:41 2013 UTC
# Line 84  C     vectorizes even after TAFization Line 84  C     vectorizes even after TAFization
84        k = 1        k = 1
85        noSlipFac = 0. _d 0        noSlipFac = 0. _d 0
86        IF ( SEAICE_no_slip ) noSlipFac = 1. _d 0        IF ( SEAICE_no_slip ) noSlipFac = 1. _d 0
87    C     in order repoduce results before fixing a bug in r1.20 comment out
88    C     the following line
89    CML      IF ( SEAICE_no_slip ) noSlipFac = 2. _d 0
90  C  C
91        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
92         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 148  C     evaluate strain rates at Z-points Line 151  C     evaluate strain rates at Z-points
151       &         )       &         )
152       &         *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)
153       &         *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)
154       &         + 2.0 _d 0 * noSlipFac * (       &         + noSlipFac * (
155       &           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
156       &         + 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
157       &         )       &         )
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22