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

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

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

revision 1.3 by jmc, Fri Sep 24 16:55:44 2004 UTC revision 1.4 by jmc, Tue Dec 5 05:26:46 2006 UTC
# Line 7  CBOP Line 7  CBOP
7  C !ROUTINE: GAD_DEL2  C !ROUTINE: GAD_DEL2
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DEL2(        SUBROUTINE GAD_DEL2(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           dTdx,dTdy,       I           dTdx,dTdy,
13       O           del2,       O           del2,
# Line 42  C  del2                 :: Laplacian of Line 42  C  del2                 :: Laplacian of
42  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
43  C  i,j                  :: loop indices  C  i,j                  :: loop indices
44        INTEGER i,j        INTEGER i,j
 c     _RL dfx(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
 c     _RL dfy(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
45  CEOP  CEOP
46    
47  C--   Initialise to zero.  C--   Initialise to zero.
# Line 57  C     Difference of zonal fluxes & merid Line 55  C     Difference of zonal fluxes & merid
55  C      and divide by grid-cell volume  C      and divide by grid-cell volume
56        DO j=1-Oly,sNy+Oly-1        DO j=1-Oly,sNy+Oly-1
57         DO i=1-Olx,sNx+Olx-1         DO i=1-Olx,sNx+Olx-1
58          del2(i,j)=recip_rA(i,j,bi,bj)          del2(i,j)=recip_rA(i,j,bi,bj)*recip_deepFacC(k)
59       &           *recip_drF(k)       &           *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj)
      &           *_recip_hFacC(i,j,k,bi,bj)  
60       &           *( ( dTdx(i+1,j)-dTdx(i,j) )       &           *( ( dTdx(i+1,j)-dTdx(i,j) )
61       &             +( dTdy(i,j+1)-dTdy(i,j) )       &             +( dTdy(i,j+1)-dTdy(i,j) )
62       &            )       &            )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22