/[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.18 by jmc, Mon Jan 13 19:02:45 2003 UTC revision 1.19 by heimbach, Tue Jan 21 19:34:13 2003 UTC
# Line 195  c Line 195  c
195       &            *dRdSigmaLtd(i,j)*dSigmaDrReal(i,j)       &            *dRdSigmaLtd(i,j)*dSigmaDrReal(i,j)
196  cph             T11(i,j)=(dSigmaDrReal(i,j)**2)*dRdSigmaLtd(i,j)  cph             T11(i,j)=(dSigmaDrReal(i,j)**2)*dRdSigmaLtd(i,j)
197            ENDIF            ENDIF
198    #ifndef ALLOWW_AUTODIFF_TAMC
199  cph-- this part doesn't adjoint well  cph-- this part doesn't adjoint well
200  cph          IF ( SlopeSqr(i,j) .GT. maxSlopeSqr .AND.            IF ( SlopeSqr(i,j) .GT. maxSlopeSqr .AND.
201  cph     &         SlopeSqr(i,j) .LT. GM_slopeSqCutoff ) THEN       &         SlopeSqr(i,j) .LT. GM_slopeSqCutoff ) THEN
202  cph           taperFct(i,j) = maxSlopeSqr/SlopeSqr(i,j)             taperFct(i,j) = maxSlopeSqr/SlopeSqr(i,j)
203  cph          ELSE IF ( SlopeSqr(i,j) .GT. GM_slopeSqCutoff ) THEN            ELSE IF ( SlopeSqr(i,j) .GT. GM_slopeSqCutoff ) THEN
204  cph           taperFct(i,j) = 0. _d 0             taperFct(i,j) = 0. _d 0
205  cph          ENDIF            ENDIF
206    #endif
207           ENDDO           ENDDO
208          ENDDO          ENDDO
209    
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22