/[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.5 by jmc, Wed Jun 22 00:27:47 2005 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 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)*recip_drC(k)       &            *recip_drC(k)
66       &     *(Tracer(i,j,k,bi,bj)-Tracer(i,j,km1,bi,bj))*rkSign       &            *(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.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22