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

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

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

revision 1.10 by jmc, Tue Mar 29 15:47:19 2011 UTC revision 1.11 by jmc, Wed Jul 4 20:22:26 2012 UTC
# Line 63  C  wCFL              :: Courant-Friedric Line 63  C  wCFL              :: Courant-Friedric
63        _RL psiP,psiM,thetaP,thetaM        _RL psiP,psiM,thetaP,thetaM
64        _RL smallNo        _RL smallNo
65    
66        IF (inAdMode) THEN  c     IF (inAdMode) THEN
67    c      smallNo = 1.0D-20
68    c     ELSE
69         smallNo = 1.0D-20         smallNo = 1.0D-20
70        ELSE  c     ENDIF
        smallNo = 1.0D-20  
       ENDIF  
71  #endif  #endif
72    
73        km2=MAX(1,k-2)        km2=MAX(1,k-2)
74        km1=MAX(1,k-1)        km1=MAX(1,k-1)
75        kp1=MIN(Nr,k+1)        kp1=MIN(Nr,k+1)
76    
77        DO j=1-Oly,sNy+Oly        DO j=1-OLy,sNy+OLy
78         DO i=1-Olx,sNx+Olx         DO i=1-OLx,sNx+OLx
79          Rjp=(tracer(i,j,k)-tracer(i,j,kp1))          Rjp=(tracer(i,j,k)-tracer(i,j,kp1))
80       &         *maskC(i,j,kp1,bi,bj)       &         *maskC(i,j,kp1,bi,bj)
81          Rj =(tracer(i,j,km1)-tracer(i,j,k))          Rj =(tracer(i,j,km1)-tracer(i,j,k))
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22