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

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

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

revision 1.3.2.1 by adcroft, Tue Feb 26 16:04:48 2002 UTC revision 1.5 by jmc, Mon Aug 18 12:22:46 2014 UTC
# Line 7  CBOP Line 7  CBOP
7  C !ROUTINE: GAD_U3_ADV_R  C !ROUTINE: GAD_U3_ADV_R
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_U3_ADV_R(        SUBROUTINE GAD_U3_ADV_R(
11       I           bi,bj,k,       I           bi, bj, k,
12       I           rTrans,       I           rTrans,
13       I           tracer,       I           tracer,
14       O           wT,       O           wT,
# Line 37  C  tracer               :: tracer field Line 37  C  tracer               :: tracer field
37  C  myThid               :: thread number  C  myThid               :: thread number
38        INTEGER bi,bj,k        INTEGER bi,bj,k
39        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
41        INTEGER myThid        INTEGER myThid
42    
43  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
# Line 60  CEOP Line 60  CEOP
60        kp1=MIN(Nr,k+1)        kp1=MIN(Nr,k+1)
61    
62        IF ( k.EQ.1 .OR. k.GT.Nr ) THEN        IF ( k.EQ.1 .OR. k.GT.Nr ) THEN
63         DO j=1-Oly,sNy+Oly         DO j=1-OLy,sNy+OLy
64          DO i=1-Olx,sNx+Olx          DO i=1-OLx,sNx+OLx
65           wT(i,j) = 0.           wT(i,j) = 0.
66          ENDDO          ENDDO
67         ENDDO         ENDDO
68        ELSE        ELSE
69         DO j=1-Oly,sNy+Oly         DO j=1-OLy,sNy+OLy
70          DO i=1-Olx,sNx+Olx          DO i=1-OLx,sNx+OLx
71           Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))           Rjp = (tracer(i,j,kp1)-tracer(i,j, k ))*maskC(i,j,kp1,bi,bj)
72       &        *maskC(i,j,kp1,bi,bj)           Rj  = (tracer(i,j, k )-tracer(i,j,km1))
73           Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,km1,bi,bj))           Rjm = (tracer(i,j,km1)-tracer(i,j,km2))*maskC(i,j,km2,bi,bj)
74           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,km2,bi,bj))           Rjjp = Rjp-Rj
75       &        *maskC(i,j,km2,bi,bj)           Rjjm = Rj-Rjm
          Rjjp=Rjp-Rj  
          Rjjm=Rj-Rjm  
76           wT(i,j) = maskC(i,j,km1,bi,bj)*(           wT(i,j) = maskC(i,j,km1,bi,bj)*(
77       &     rTrans(i,j)*(       &     rTrans(i,j)*( (tracer(i,j,k)+tracer(i,j,km1))*0.5 _d 0
78       &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,km1,bi,bj))*0.5 _d 0       &                  -oneSixth*(Rjjm+Rjjp)*0.5 _d 0  )
      &        -oneSixth*(Rjjm+Rjjp)*0.5 _d 0  )  
79       &    +ABS(rTrans(i,j))*       &    +ABS(rTrans(i,j))*
80       &         oneSixth*(Rjjm-Rjjp)*0.5 _d 0       &                   oneSixth*(Rjjm-Rjjp)*0.5 _d 0
81       &                                  )       &                                  )
82          ENDDO          ENDDO
83         ENDDO         ENDDO
Legend:
Removed from v.1.3.2.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22