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

Diff of /MITgcm/pkg/generic_advdiff/gad_c2_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_C2_ADV_R                                  |
14    C     | o Compute vertical advective Flux of Tracer using        |
15    C     |   2nd Order Centered Scheme                              |
16    C     |==========================================================|
17        IMPLICIT NONE        IMPLICIT NONE
18    
19  C     == GLobal variables ==  C     == GLobal variables ==
# Line 29  C     == Local variables == Line 34  C     == Local variables ==
34    
35        km1=max(1,k-1)        km1=max(1,k-1)
36    
37        IF ( (rigidLid .AND. k.EQ.1) .OR. k.GT.Nr) THEN        IF ( k.EQ.1 .OR. k.GT.Nr) THEN
38         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
39          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
40           wT(i,j) = 0.           wT(i,j) = 0.
41          ENDDO          ENDDO
42         ENDDO         ENDDO
 c?    ELSEIF ( rigidLid ) THEN  
 c?     DO j=1-Oly,sNy+Oly  
 c?      DO i=1-Olx,sNx+Olx  
 c?       wT(i,j) = rTrans(i,j)*  
 c?   &       (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0  
 c?      ENDDO  
 c?     ENDDO  
43        ELSE        ELSE
44         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
45          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
46           wT(i,j) =           wT(i,j) = maskC(i,j,kM1,bi,bj)*
47       &     rTrans(i,j)       &     rTrans(i,j)*
      &      *(maskC(i,j,kM1,bi,bj)*  
48       &        (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
      &       +(maskC(i,j,k,bi,bj)-maskC(i,j,kM1,bi,bj))*  
      &        Tracer(i,j,k,bi,bj)  
      &       )  
49          ENDDO          ENDDO
50         ENDDO         ENDDO
51        ENDIF        ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22