/[MITgcm]/MITgcm/pkg/gmredi/gmredi_slope_limit.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_slope_limit.F

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

revision 1.13 by heimbach, Thu Nov 14 22:43:49 2002 UTC revision 1.14 by heimbach, Fri Nov 15 02:57:47 2002 UTC
# Line 118  cnostore CADJ STORE dSigmaDrReal(:,:) = Line 118  cnostore CADJ STORE dSigmaDrReal(:,:) =
118            IF (gradSmod(i,j) .NE. 0.) THEN            IF (gradSmod(i,j) .NE. 0.) THEN
119             dSigmaDrLtd(i,j) = -gradSmod(i,j)*GM_rMaxSlope             dSigmaDrLtd(i,j) = -gradSmod(i,j)*GM_rMaxSlope
120             IF ( dSigmaDrReal(i,j) .GE.             IF ( dSigmaDrReal(i,j) .GE.
121       &          GM_adjointRescale*dSigmaDrLtd(i,j) )       &          dSigmaDrLtd(i,j) )
122       &          dSigmaDrReal(i,j) =       &          dSigmaDrReal(i,j) =
123       &          dSigmaDrLtd(i,j)*GM_adjointRescale       &          dSigmaDrLtd(i,j)
124    cph           IF ( dSigmaDrReal(i,j) .GE.
125    cph     &          GM_adjointRescale*dSigmaDrLtd(i,j) )
126    cph     &          dSigmaDrReal(i,j) =
127    cph     &          dSigmaDrLtd(i,j)*GM_adjointRescale
128  ctest           dSigmaDrReal(i,j) = dSigmaDrLtd(i,j)  ctest           dSigmaDrReal(i,j) = dSigmaDrLtd(i,j)
129            ENDIF            ENDIF
130           ENDDO           ENDDO
# Line 236  cnostore CADJ STORE dSigmaDrReal(:,:) = Line 240  cnostore CADJ STORE dSigmaDrReal(:,:) =
240          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
241           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
242            IF ( dSigmaDrReal(i,j) .EQ. 0. ) THEN            IF ( dSigmaDrReal(i,j) .EQ. 0. ) THEN
243             IF ( dSigmaDx(i,j) .NE. 0. )             IF ( dSigmaDx(i,j) .NE. 0. ) THEN
244       &      SlopeX(i,j) = SIGN(Small_taper,dSigmaDx(i,j))              SlopeX(i,j) = SIGN(Small_taper,dSigmaDx(i,j))
245             IF ( dSigmaDy(i,j) .NE. 0. )             ELSE
246       &      SlopeY(i,j) = SIGN(Small_taper,dSigmaDy(i,j))              SlopeX(i,j) = 0. _d 0
247               ENDIF
248               IF ( dSigmaDy(i,j) .NE. 0. ) THEN
249                SlopeY(i,j) = SIGN(Small_taper,dSigmaDy(i,j))
250               ELSE
251                SlopeY(i,j) = 0. _d 0
252               ENDIF
253            ELSE            ELSE
254             dRdSigmaLtd(i,j) = 1./dSigmaDrReal(i,j)             dRdSigmaLtd(i,j) = 1./dSigmaDrReal(i,j)
255             SlopeX(i,j) = -dSigmaDx(i,j)*dRdSigmaLtd(i,j)             SlopeX(i,j) = -dSigmaDx(i,j)*dRdSigmaLtd(i,j)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22