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

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

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

revision 1.1 by adcroft, Wed May 30 19:34:48 2001 UTC revision 1.2 by jmc, Thu Jul 12 00:30:27 2001 UTC
# Line 9  C $Name$ Line 9  C $Name$
9       I           tracer,       I           tracer,
10       O           wT,       O           wT,
11       I           myThid )       I           myThid )
12    C     /==========================================================\
13    C     | SUBROUTINE GAD_C4_ADV_R                                  |
14    C     | o Compute vertical advective Flux of Tracer using        |
15    C     |   4th Order Centered Scheme                              |
16    C     |==========================================================|
17        IMPLICIT NONE        IMPLICIT NONE
18    
19  C     == GLobal variables ==  C     == GLobal variables ==
# Line 32  C     == Local variables == Line 37  C     == Local variables ==
37        km1=MAX(1,k-1)        km1=MAX(1,k-1)
38        kp1=MIN(Nr,k+1)        kp1=MIN(Nr,k+1)
39    
40        IF ( k.GT.Nr) THEN        IF ( k.EQ.1 .OR. k.GT.Nr) THEN
41         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
42          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
43           wT(i,j) = 0.           wT(i,j) = 0.
# Line 43  C     == Local variables == Line 48  C     == Local variables ==
48          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
49           Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))           Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))
50       &        *maskC(i,j,kp1,bi,bj)       &        *maskC(i,j,kp1,bi,bj)
51           Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,kM1,bi,bj))           Rj =(tracer(i,j,k,bi,bj)-tracer(i,j,km1,bi,bj))
52       &        *maskC(i,j,k,bi,bj)           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,km2,bi,bj))
53           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,kM2,bi,bj))       &        *maskC(i,j,km2,bi,bj)
      &        *maskC(i,j,km1,bi,bj)  
54           Rjjp=Rjp-Rj           Rjjp=Rjp-Rj
55           Rjjm=Rj-Rjm           Rjjm=Rj-Rjm
56           wT(i,j) =           wT(i,j) = maskC(i,j,kM1,bi,bj)*
57       &     rTrans(i,j)*(       &     rTrans(i,j)*(
      &       (maskC(i,j,kM1,bi,bj)*  
58       &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0       &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0
59       &       +(maskC(i,j,k,bi,bj)-maskC(i,j,kM1,bi,bj))*       &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )
      &        Tracer(i,j,k,bi,bj)  
      &       )  
      &     -maskC(i,j,km1,bi,bj)*oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )  
60          ENDDO          ENDDO
61         ENDDO         ENDDO
62        ENDIF        ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22