/[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.18 by jmc, Fri Oct 21 17:32:01 2011 UTC revision 1.20 by mlosch, Thu Feb 28 17:25:41 2013 UTC
# Line 33  C     === Global variables === Line 33  C     === Global variables ===
33  #include "EEPARAMS.h"  #include "EEPARAMS.h"
34  #include "PARAMS.h"  #include "PARAMS.h"
35  #include "GRID.h"  #include "GRID.h"
36    #include "SEAICE_SIZE.h"
37  #include "SEAICE_PARAMS.h"  #include "SEAICE_PARAMS.h"
38  #include "SEAICE.h"  #include "SEAICE.h"
39    
# Line 51  C     iStep  :: Sub-time-step number Line 52  C     iStep  :: Sub-time-step number
52  C     myTime :: Simulation time  C     myTime :: Simulation time
53  C     myIter :: Simulation timestep number  C     myIter :: Simulation timestep number
54  C     myThid :: My Thread Id. number  C     myThid :: My Thread Id. number
55        _RL uFld   (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)        _RL uFld   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
56        _RL vFld   (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)        _RL vFld   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
57        _RL e11Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL e11Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
58        _RL e22Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL e22Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
59        _RL e12Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL e12Loc (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 83  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)
93  C     abbreviations on C-points, need to do them in separate loops  C     abbreviations on C-points, need to do them in separate loops
94  C     for vectorization  C     for vectorization
95          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
96           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
97            dudx(i,j) = _recip_dxF(i,j,bi,bj) *            dudx(i,j) = _recip_dxF(i,j,bi,bj) *
98       &         (uFld(i+1,j,bi,bj)-uFld(i,j,bi,bj))       &         (uFld(i+1,j,bi,bj)-uFld(i,j,bi,bj))
99            uave(i,j) = 0.5 _d 0 * (uFld(i,j,bi,bj)+uFld(i+1,j,bi,bj))            uave(i,j) = 0.5 _d 0 * (uFld(i,j,bi,bj)+uFld(i+1,j,bi,bj))
100           ENDDO           ENDDO
101          ENDDO          ENDDO
102          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
103           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
104            dvdy(i,j) = _recip_dyF(i,j,bi,bj) *            dvdy(i,j) = _recip_dyF(i,j,bi,bj) *
105       &         (vFld(i,j+1,bi,bj)-vFld(i,j,bi,bj))       &         (vFld(i,j+1,bi,bj)-vFld(i,j,bi,bj))
106            vave(i,j) = 0.5 _d 0 * (vFld(i,j,bi,bj)+vFld(i,j+1,bi,bj))            vave(i,j) = 0.5 _d 0 * (vFld(i,j,bi,bj)+vFld(i,j+1,bi,bj))
107           ENDDO           ENDDO
108          ENDDO          ENDDO
109  C     evaluate strain rates at C-points  C     evaluate strain rates at C-points
110          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
111           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
112            e11Loc(i,j,bi,bj) = dudx(i,j) + vave(i,j) * k2AtC(i,j,bi,bj)            e11Loc(i,j,bi,bj) = dudx(i,j) + vave(i,j) * k2AtC(i,j,bi,bj)
113            e22Loc(i,j,bi,bj) = dvdy(i,j) + uave(i,j) * k1AtC(i,j,bi,bj)            e22Loc(i,j,bi,bj) = dvdy(i,j) + uave(i,j) * k1AtC(i,j,bi,bj)
114           ENDDO           ENDDO
115          ENDDO          ENDDO
116  #ifndef OBCS_UVICE_OLD  #ifndef OBCS_UVICE_OLD
117  C--     for OBCS: assume no gradient beyong OB  C--     for OBCS: assume no gradient beyong OB
118          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
119           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
120            e11Loc(i,j,bi,bj) = e11Loc(i,j,bi,bj)*maskInC(i,j,bi,bj)            e11Loc(i,j,bi,bj) = e11Loc(i,j,bi,bj)*maskInC(i,j,bi,bj)
121            e22Loc(i,j,bi,bj) = e22Loc(i,j,bi,bj)*maskInC(i,j,bi,bj)            e22Loc(i,j,bi,bj) = e22Loc(i,j,bi,bj)*maskInC(i,j,bi,bj)
122           ENDDO           ENDDO
# Line 121  C--     for OBCS: assume no gradient bey Line 125  C--     for OBCS: assume no gradient bey
125    
126  C     abbreviations at Z-points, need to do them in separate loops  C     abbreviations at Z-points, need to do them in separate loops
127  C     for vectorization  C     for vectorization
128          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
129           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
130            dudy(i,j) = ( uFld(i,j,bi,bj) - uFld(i  ,j-1,bi,bj) )            dudy(i,j) = ( uFld(i,j,bi,bj) - uFld(i  ,j-1,bi,bj) )
131       &         * _recip_dyU(i,j,bi,bj)       &         * _recip_dyU(i,j,bi,bj)
132            uave(i,j) = 0.5 _d 0 * (uFld(i,j,bi,bj)+uFld(i  ,j-1,bi,bj))            uave(i,j) = 0.5 _d 0 * (uFld(i,j,bi,bj)+uFld(i  ,j-1,bi,bj))
133           ENDDO           ENDDO
134          ENDDO          ENDDO
135          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
136           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
137            dvdx(i,j) = ( vFld(i,j,bi,bj) - vFld(i-1,j  ,bi,bj) )            dvdx(i,j) = ( vFld(i,j,bi,bj) - vFld(i-1,j  ,bi,bj) )
138       &         * _recip_dxV(i,j,bi,bj)       &         * _recip_dxV(i,j,bi,bj)
139            vave(i,j) = 0.5 _d 0 * (vFld(i,j,bi,bj)+vFld(i-1,j  ,bi,bj))            vave(i,j) = 0.5 _d 0 * (vFld(i,j,bi,bj)+vFld(i-1,j  ,bi,bj))
140           ENDDO           ENDDO
141          ENDDO          ENDDO
142  C     evaluate strain rates at Z-points  C     evaluate strain rates at Z-points
143          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
144           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
145            hFacU = _maskW(i,j,k,bi,bj) - _maskW(i,j-1,k,bi,bj)            hFacU = _maskW(i,j,k,bi,bj) - _maskW(i,j-1,k,bi,bj)
146            hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj)            hFacV = _maskS(i,j,k,bi,bj) - _maskS(i-1,j,k,bi,bj)
147            e12Loc(i,j,bi,bj) = 0.5 _d 0 * (            e12Loc(i,j,bi,bj) = 0.5 _d 0 * (
# Line 147  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.18  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22