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

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

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

revision 1.5 by adcroft, Fri May 11 18:24:31 2007 UTC revision 1.8 by mlosch, Fri Feb 29 01:30:59 2008 UTC
# Line 57  C     == Local variables == Line 57  C     == Local variables ==
57         uT(sNx+Olx-2,j)=0. _d 0         uT(sNx+Olx-2,j)=0. _d 0
58         uT(sNx+Olx-1,j)=0. _d 0         uT(sNx+Olx-1,j)=0. _d 0
59         uT(sNx+Olx,j)=0. _d 0         uT(sNx+Olx,j)=0. _d 0
60          ENDDO
61          DO j=1-Oly,sNy+Oly
62         DO i=1-Olx+4,sNx+Olx-3         DO i=1-Olx+4,sNx+Olx-3
63    
64          uLoc = uFld(i,j)          uLoc = uFld(i,j)
# Line 93  C     == Local variables == Line 95  C     == Local variables ==
95           MskIm   = maskLocW(i+1,j)           MskIm   = maskLocW(i+1,j)
96           MskImm  = maskLocW(i+2,j)           MskImm  = maskLocW(i+2,j)
97           MskImmm = maskLocW(i+3,j)           MskImmm = maskLocW(i+3,j)
98            ELSE
99             Qippp = 0. _d 0
100             Qipp  = 0. _d 0
101             Qip   = 0. _d 0
102             Qi    = 0. _d 0
103             Qim   = 0. _d 0
104             Qimm  = 0. _d 0
105             Qimmm = 0. _d 0
106    
107             MskIpp  = 0. _d 0
108             MskIp   = 0. _d 0
109             MskI    = 0. _d 0
110             MskIm   = 0. _d 0
111             MskImm  = 0. _d 0
112             MskImmm = 0. _d 0
113          ENDIF          ENDIF
114    
115          IF (uTrans(i,j).ne.0.) THEN          IF (uTrans(i,j).ne.0.) THEN
# Line 155  C        MP limiter Line 172  C        MP limiter
172           A = 4.*d2 - d2p1           A = 4.*d2 - d2p1
173           B = 4.*d2p1 - d2           B = 4.*d2p1 - d2
174           C = d2           C = d2
175           D = d2p1;           D = d2p1
176           dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)           dp1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)
177           A = 4.*d2m1 - d2           A = 4.*d2m1 - d2
178           B = 4.*d2 - d2m1           B = 4.*d2 - d2m1
179           C = d2m1           C = d2m1
180           D = d2;           D = d2
181           dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)           dm1h = max(min(A,B,C,D),0. _d 0)+min(max(A,B,C,D),0. _d 0)
182          !qMD = 0.5*( ( Qi + Qip ) - dp1h )          !qMD = 0.5*( ( Qi + Qip ) - dp1h )
183           qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h )           qMD = 0.5*( ( 2.*Qi + DelIp ) - dp1h )
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22