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

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

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

revision 1.3 by adcroft, Wed Sep 26 19:05:21 2001 UTC revision 1.3.4.1 by heimbach, Wed Feb 6 15:48:08 2002 UTC
# Line 21  C \begin{equation*} Line 21  C \begin{equation*}
21  C F^x_{adv} = U \overline{ \theta  - \frac{1}{6} \delta_{ii} \theta }^i  C F^x_{adv} = U \overline{ \theta  - \frac{1}{6} \delta_{ii} \theta }^i
22  C                 + \frac{1}{12} |U| \delta_{iii} \theta  C                 + \frac{1}{12} |U| \delta_{iii} \theta
23  C \end{equation*}  C \end{equation*}
24    C Near boundaries, mask all the gradients ==> still 3rd O.
25    
26  C !USES: ===============================================================  C !USES: ===============================================================
27        IMPLICIT NONE        IMPLICIT NONE
# Line 56  CEOP Line 57  CEOP
57         uT(2-Olx,j)=0.         uT(2-Olx,j)=0.
58         uT(sNx+Olx,j)=0.         uT(sNx+Olx,j)=0.
59         DO i=1-Olx+2,sNx+Olx-1         DO i=1-Olx+2,sNx+Olx-1
60          Rjp=(tracer(i+1,j)-tracer(i,j))          Rjp=(tracer(i+1,j)-tracer(i,j))*maskW(i+1,j,k,bi,bj)
61  c    &     *maskW(i+1,j,k,bi,bj)          Rj =(tracer(i,j)-tracer(i-1,j))*maskW(i,j,k,bi,bj)
62          Rj =(tracer(i,j)-tracer(i-1,j))          Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskW(i-1,j,k,bi,bj)
63  c    &     *maskW(i,j,k,bi,bj)          Rjjp=Rjp-Rj
64          Rjm=(tracer(i-1,j)-tracer(i-2,j))          Rjjm=Rj-Rjm
 c    &     *maskW(i-1,j,k,bi,bj)  
 C- jmc: mask not needed here above if Rjjp & Rjjm & uTrans are masked  
         Rjjp=(Rjp-Rj)*maskW(i+1,j,k,bi,bj)  
         Rjjm=(Rj-Rjm)*maskW(i-1,j,k,bi,bj)  
65          uT(i,j) =          uT(i,j) =
66       &   uTrans(i,j)*(       &   uTrans(i,j)*(
67       &     Tracer(i,j)+Tracer(i-1,j)-oneSixth*( Rjjp+Rjjm )       &     Tracer(i,j)+Tracer(i-1,j)-oneSixth*( Rjjp+Rjjm )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.3.4.1
  ViewVC Help
Powered by ViewVC 1.1.22