/[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.2 by jmc, Thu Jul 12 00:30:27 2001 UTC revision 1.3 by jmc, Mon Aug 20 20:45:15 2001 UTC
# Line 13  C     /================================= Line 13  C     /=================================
13  C     | SUBROUTINE GAD_C4_ADV_R                                  |  C     | SUBROUTINE GAD_C4_ADV_R                                  |
14  C     | o Compute vertical advective Flux of Tracer using        |  C     | o Compute vertical advective Flux of Tracer using        |
15  C     |   4th Order Centered Scheme                              |  C     |   4th Order Centered Scheme                              |
16    C     | o reduced to 2nd or 3rd Order near a boundary            |
17  C     |==========================================================|  C     |==========================================================|
18        IMPLICIT NONE        IMPLICIT NONE
19    
# Line 30  C     == Routine arguments == Line 31  C     == Routine arguments ==
31    
32  C     == Local variables ==  C     == Local variables ==
33        INTEGER i,j,kp1,km1,km2        INTEGER i,j,kp1,km1,km2
34        _RL Cr,Rjm,Rj,Rjp,Rjjm,Rjjp        _RL Rjm,Rj,Rjp,Rjjm,Rjjp
35  #include "GAD_FLUX_LIMITER.h"        _RL maskP1
36    
37        km2=MAX(1,k-2)        km2=MAX(1,k-2)
38        km1=MAX(1,k-1)        km1=MAX(1,k-1)
39        kp1=MIN(Nr,k+1)        kp1=MIN(Nr,k+1)
40          maskP1 = 1.
41          IF (k.GE.Nr) maskP1 = 0.
42    
43        IF ( k.EQ.1 .OR. k.GT.Nr) THEN        IF ( k.EQ.1 .OR. k.GT.Nr) THEN
44         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
# Line 43  C     == Local variables == Line 46  C     == Local variables ==
46           wT(i,j) = 0.           wT(i,j) = 0.
47          ENDDO          ENDDO
48         ENDDO         ENDDO
49          ELSEIF ( k.EQ.2 ) THEN
50             Rjjm= 0.
51           DO j=1-Oly,sNy+Oly
52            DO i=1-Olx,sNx+Olx
53             Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))
54             Rj =(tracer(i,j,k,bi,bj)-tracer(i,j,km1,bi,bj))
55             Rjjp=(Rjp-Rj)*maskC(i,j,kp1,bi,bj)*maskP1
56             wT(i,j) = maskC(i,j,km1,bi,bj)*(
57         &     rTrans(i,j)*(
58         &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,km1,bi,bj))*0.5 _d 0
59         &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )
60         &    +ABS(rTrans(i,j))*
61         &         oneSixth*(Rjjm-Rjjp)*0.5 _d 0
62         &                                  )                      
63            ENDDO
64           ENDDO
65          ELSEIF ( k.EQ.Nr ) THEN
66             Rjjp= 0.
67           DO j=1-Oly,sNy+Oly
68            DO i=1-Olx,sNx+Olx
69             Rj =(tracer(i,j,k,bi,bj)-tracer(i,j,km1,bi,bj))
70             Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,km2,bi,bj))
71             Rjjm=(Rj-Rjm)*maskC(i,j,km2,bi,bj)
72             wT(i,j) = maskC(i,j,km1,bi,bj)*(
73         &     rTrans(i,j)*(
74         &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,km1,bi,bj))*0.5 _d 0
75         &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )
76         &    +ABS(rTrans(i,j))*
77         &         oneSixth*(Rjjm-Rjjp)*0.5 _d 0
78         &                                  )                      
79            ENDDO
80           ENDDO
81        ELSE        ELSE
82         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
83          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
84           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))
      &        *maskC(i,j,kp1,bi,bj)  
85           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))
86           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,km2,bi,bj))           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,km2,bi,bj))
87       &        *maskC(i,j,km2,bi,bj)           Rjjp=(Rjp-Rj)*maskC(i,j,kp1,bi,bj)
88           Rjjp=Rjp-Rj           Rjjm=(Rj-Rjm)*maskC(i,j,km2,bi,bj)
89           Rjjm=Rj-Rjm           wT(i,j) = maskC(i,j,km1,bi,bj)*(
          wT(i,j) = maskC(i,j,kM1,bi,bj)*  
90       &     rTrans(i,j)*(       &     rTrans(i,j)*(
91       &        (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
92       &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )       &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0 )
93         &    +ABS(rTrans(i,j))*
94         &         oneSixth*(Rjjm-Rjjp)*0.5 _d 0
95         &      *( 1. _d 0 - maskC(i,j,km2,bi,bj)*maskC(i,j,kp1,bi,bj) )
96         &                                  )                      
97          ENDDO          ENDDO
98         ENDDO         ENDDO
99        ENDIF        ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22