/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_u3c4_impl_r.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_u3c4_impl_r.F

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

revision 1.10 by jmc, Thu Dec 1 14:14:44 2011 UTC revision 1.11 by jmc, Fri Apr 4 20:29:08 2014 UTC
# Line 70  C     rC4km, rC4kp :: high order contrib Line 70  C     rC4km, rC4kp :: high order contrib
70  C     rHigh        :: high order term factor  C     rHigh        :: high order term factor
71        LOGICAL flagC4        LOGICAL flagC4
72        INTEGER i,j,kp1,km2        INTEGER i,j,kp1,km2
73  #if (defined ALLOW_AUTODIFF_TAMC && defined TARGET_NEC_SX)  #if (defined ALLOW_AUTODIFF && defined TARGET_NEC_SX)
74        _RL rC4km2D  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rC4km2D  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
75        _RL rC4kp2D  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rC4kp2D  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
76        _RL rCenter2D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rCenter2D(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 96  C--   process interior interface only: Line 96  C--   process interior interface only:
96    
97  C--    Add centered, upwind and high-order contributions  C--    Add centered, upwind and high-order contributions
98         deltaTcfl = deltaTarg(k)         deltaTcfl = deltaTarg(k)
99  #if (defined ALLOW_AUTODIFF_TAMC && defined TARGET_NEC_SX)  #if (defined ALLOW_AUTODIFF && defined TARGET_NEC_SX)
100         DO j=jMin,jMax         DO j=jMin,jMax
101          DO i=iMin,iMax          DO i=iMin,iMax
102           rCenter2D(i,j) =           rCenter2D(i,j) =
# Line 122  c         rUpwind2D(i,j) = (2. _d 0*rHig Line 122  c         rUpwind2D(i,j) = (2. _d 0*rHig
122           ENDIF           ENDIF
123          ENDDO          ENDDO
124         ENDDO         ENDDO
125  #endif /* ALLOW_AUTODIFF_TAMC and TARGET_NEC_SX */  #endif /* ALLOW_AUTODIFF and TARGET_NEC_SX */
126         DO j=jMin,jMax         DO j=jMin,jMax
127           DO i=iMin,iMax           DO i=iMin,iMax
128  #if (defined ALLOW_AUTODIFF_TAMC && defined TARGET_NEC_SX)  #if (defined ALLOW_AUTODIFF && defined TARGET_NEC_SX)
129             rC4km   = rC4km2D  (i,j)             rC4km   = rC4km2D  (i,j)
130             rC4kp   = rC4kp2D  (i,j)             rC4kp   = rC4kp2D  (i,j)
131             rCenter = rCenter2D(i,j)             rCenter = rCenter2D(i,j)
# Line 151  c           rUpwind= (2. _d 0*rHigh - wC Line 151  c           rUpwind= (2. _d 0*rHigh - wC
151              rC4km  = oneSixth*(rCenter+ABS(rCenter))*mskM              rC4km  = oneSixth*(rCenter+ABS(rCenter))*mskM
152              rC4kp  = oneSixth*(rCenter-ABS(rCenter))*mskP              rC4kp  = oneSixth*(rCenter-ABS(rCenter))*mskP
153             ENDIF             ENDIF
154  #endif /* ALLOW_AUTODIFF_TAMC and TARGET_NEC_SX */  #endif /* ALLOW_AUTODIFF and TARGET_NEC_SX */
155             a5d(i,j,k)   = a5d(i,j,k)             a5d(i,j,k)   = a5d(i,j,k)
156       &                  + rC4km       &                  + rC4km
157       &                   *deltaTarg(k)       &                   *deltaTarg(k)
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22