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

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

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

revision 1.3 by adcroft, Fri Sep 21 13:11:43 2001 UTC revision 1.6 by jmc, Tue Dec 5 05:26:46 2006 UTC
# Line 7  CBOP Line 7  CBOP
7  C !ROUTINE: GAD_DIFF_R  C !ROUTINE: GAD_DIFF_R
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DIFF_R(        SUBROUTINE GAD_DIFF_R(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           KappaR,       I           KappaR,
13       I           tracer,       I           tracer,
# Line 35  C  KappaR               :: vertical diff Line 35  C  KappaR               :: vertical diff
35  C  tracer               :: tracer field  C  tracer               :: tracer field
36  C  myThid               :: thread number  C  myThid               :: thread number
37        INTEGER bi,bj,k        INTEGER bi,bj,k
38        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40        INTEGER myThid        INTEGER myThid
41    
# Line 49  C  km1                  :: =k+1 for k<Nr Line 49  C  km1                  :: =k+1 for k<Nr
49        INTEGER i,j,km1        INTEGER i,j,km1
50  CEOP  CEOP
51    
52        km1=max(1,k-1)        km1=MAX(1,k-1)
53    
54        IF (k.eq.1 .OR. k.gt.Nr) THEN        IF (k.EQ.1 .OR. k.GT.Nr) THEN
55         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
56          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
57           dfr(i,j) = 0.           dfr(i,j) = 0.
# Line 60  CEOP Line 60  CEOP
60        ELSE        ELSE
61         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
62          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
63           dfr(i,j) =           dfr(i,j) = -KappaR(i,j)
64       &    -_rA(i,j,bi,bj)       &            *_rA(i,j,bi,bj)*deepFac2F(k)*rhoFacF(k)
65       &     *KappaR(i,j,k)*recip_drC(k)       &            *recip_drC(k)
66       &     *(Tracer(i,j,km1,bi,bj)-Tracer(i,j,k,bi,bj))*rkFac       &            *(tracer(i,j,k,bi,bj)-tracer(i,j,km1,bi,bj))*rkSign
67          ENDDO          ENDDO
68         ENDDO         ENDDO
69        ENDIF        ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22