/[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.19 by jmc, Tue Mar 6 16:45:20 2012 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 88  C Line 89  C
89         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
90  C     abbreviations on C-points, need to do them in separate loops  C     abbreviations on C-points, need to do them in separate loops
91  C     for vectorization  C     for vectorization
92          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
93           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
94            dudx(i,j) = _recip_dxF(i,j,bi,bj) *            dudx(i,j) = _recip_dxF(i,j,bi,bj) *
95       &         (uFld(i+1,j,bi,bj)-uFld(i,j,bi,bj))       &         (uFld(i+1,j,bi,bj)-uFld(i,j,bi,bj))
96            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))
97           ENDDO           ENDDO
98          ENDDO          ENDDO
99          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
100           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
101            dvdy(i,j) = _recip_dyF(i,j,bi,bj) *            dvdy(i,j) = _recip_dyF(i,j,bi,bj) *
102       &         (vFld(i,j+1,bi,bj)-vFld(i,j,bi,bj))       &         (vFld(i,j+1,bi,bj)-vFld(i,j,bi,bj))
103            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))
104           ENDDO           ENDDO
105          ENDDO          ENDDO
106  C     evaluate strain rates at C-points  C     evaluate strain rates at C-points
107          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
108           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
109            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)
110            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)
111           ENDDO           ENDDO
112          ENDDO          ENDDO
113  #ifndef OBCS_UVICE_OLD  #ifndef OBCS_UVICE_OLD
114  C--     for OBCS: assume no gradient beyong OB  C--     for OBCS: assume no gradient beyong OB
115          DO j=1-Oly,sNy+Oly-1          DO j=1-OLy,sNy+OLy-1
116           DO i=1-Olx,sNx+Olx-1           DO i=1-OLx,sNx+OLx-1
117            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)
118            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)
119           ENDDO           ENDDO
# Line 121  C--     for OBCS: assume no gradient bey Line 122  C--     for OBCS: assume no gradient bey
122    
123  C     abbreviations at Z-points, need to do them in separate loops  C     abbreviations at Z-points, need to do them in separate loops
124  C     for vectorization  C     for vectorization
125          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
126           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
127            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) )
128       &         * _recip_dyU(i,j,bi,bj)       &         * _recip_dyU(i,j,bi,bj)
129            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))
130           ENDDO           ENDDO
131          ENDDO          ENDDO
132          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
133           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
134            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) )
135       &         * _recip_dxV(i,j,bi,bj)       &         * _recip_dxV(i,j,bi,bj)
136            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))
137           ENDDO           ENDDO
138          ENDDO          ENDDO
139  C     evaluate strain rates at Z-points  C     evaluate strain rates at Z-points
140          DO j=1-Oly+1,sNy+Oly          DO j=1-OLy+1,sNy+OLy
141           DO i=1-Olx+1,sNx+Olx           DO i=1-OLx+1,sNx+OLx
142            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)
143            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)
144            e12Loc(i,j,bi,bj) = 0.5 _d 0 * (            e12Loc(i,j,bi,bj) = 0.5 _d 0 * (
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22