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

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

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

revision 1.13 by jmc, Tue Mar 29 15:47:19 2011 UTC revision 1.14 by jmc, Wed Jul 4 20:22:26 2012 UTC
# Line 63  C  vCFL              :: Courant-Friedric Line 63  C  vCFL              :: Courant-Friedric
63        _RL smallNo        _RL smallNo
64  c     _RL Rjjm,Rjjp  c     _RL Rjjm,Rjjp
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        DO i=1-Olx,sNx+Olx        DO i=1-OLx,sNx+OLx
74         vT(i,1-Oly)=0.         vT(i,1-OLy)=0.
75         vT(i,2-Oly)=0.         vT(i,2-OLy)=0.
76         vT(i,sNy+Oly)=0.         vT(i,sNy+OLy)=0.
77        ENDDO        ENDDO
78        DO j=1-Oly+2,sNy+Oly-1        DO j=1-OLy+2,sNy+OLy-1
79         DO i=1-Olx,sNx+Olx         DO i=1-OLx,sNx+OLx
80          Rjp=(tracer(i,j+1)-tracer(i, j ))*maskLocS(i,j+1)          Rjp=(tracer(i,j+1)-tracer(i, j ))*maskLocS(i,j+1)
81          Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )          Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )
82          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22