/[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.7 by mlosch, Fri Oct 5 10:50:47 2007 UTC
# Line 93  C     == Local variables == Line 93  C     == Local variables ==
93           MskIm   = maskLocW(i+1,j)           MskIm   = maskLocW(i+1,j)
94           MskImm  = maskLocW(i+2,j)           MskImm  = maskLocW(i+2,j)
95           MskImmm = maskLocW(i+3,j)           MskImmm = maskLocW(i+3,j)
96            ELSE
97             Qippp = 0. _d 0
98             Qipp  = 0. _d 0
99             Qip   = 0. _d 0
100             Qi    = 0. _d 0
101             Qim   = 0. _d 0
102             Qimm  = 0. _d 0
103             Qimmm = 0. _d 0
104    
105             MskIpp  = 0. _d 0
106             MskIp   = 0. _d 0
107             MskI    = 0. _d 0
108             MskIm   = 0. _d 0
109             MskImm  = 0. _d 0
110             MskImmm = 0. _d 0
111          ENDIF          ENDIF
112    
113          IF (uTrans(i,j).ne.0.) THEN          IF (uTrans(i,j).ne.0.) THEN
# Line 155  C        MP limiter Line 170  C        MP limiter
170           A = 4.*d2 - d2p1           A = 4.*d2 - d2p1
171           B = 4.*d2p1 - d2           B = 4.*d2p1 - d2
172           C = d2           C = d2
173           D = d2p1;           D = d2p1
174           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)
175           A = 4.*d2m1 - d2           A = 4.*d2m1 - d2
176           B = 4.*d2 - d2m1           B = 4.*d2 - d2m1
177           C = d2m1           C = d2m1
178           D = d2;           D = d2
179           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)
180          !qMD = 0.5*( ( Qi + Qip ) - dp1h )          !qMD = 0.5*( ( Qi + Qip ) - dp1h )
181           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.7
  ViewVC Help
Powered by ViewVC 1.1.22