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

Diff of /MITgcm/pkg/generic_advdiff/gad_u3_adv_y.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.4 by jmc, Tue Jan 8 21:43:24 2002 UTC
# Line 21  C \begin{equation*} Line 21  C \begin{equation*}
21  C F^y_{adv} = V \overline{ \theta  - \frac{1}{6} \delta_{jj} \theta }^j  C F^y_{adv} = V \overline{ \theta  - \frac{1}{6} \delta_{jj} \theta }^j
22  C                 + \frac{1}{12} |V| \delta_{jjj} \theta  C                 + \frac{1}{12} |V| \delta_{jjj} \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 58  CEOP Line 59  CEOP
59        ENDDO        ENDDO
60        DO j=1-Oly+2,sNy+Oly-1        DO j=1-Oly+2,sNy+Oly-1
61         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
62          Rjp=(tracer(i,j+1)-tracer(i,j))          Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj)
63  c    &     *maskS(i,j+1,k,bi,bj)          Rj =(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj)
64          Rj =(tracer(i,j)-tracer(i,j-1))          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskS(i,j-1,k,bi,bj)
65  c    &     *maskS(i,j,k,bi,bj)          Rjjp=Rjp-Rj
66          Rjm=(tracer(i,j-1)-tracer(i,j-2))          Rjjm=Rj-Rjm
 c    &     *maskS(i,j-1,k,bi,bj)  
 C- jmc: mask not needed here above if Rjjp & Rjjm & vTrans are masked  
         Rjjp=(Rjp-Rj)*maskS(i,j+1,k,bi,bj)  
         Rjjm=(Rj-Rjm)*maskS(i,j-1,k,bi,bj)  
67          vT(i,j) =          vT(i,j) =
68       &   vTrans(i,j)*(       &   vTrans(i,j)*(
69       &     Tracer(i,j)+Tracer(i,j-1)-oneSixth*( Rjjp+Rjjm )       &     Tracer(i,j)+Tracer(i,j-1)-oneSixth*( Rjjp+Rjjm )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22