/[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.2 by adcroft, Thu Sep 20 20:07:06 2001 UTC revision 1.3 by jmc, Fri Sep 24 16:55:44 2004 UTC
# 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
45    c     _RL dfx(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46    c     _RL dfy(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47  CEOP  CEOP
48    
49  C     Difference of zonal fluxes ...  C--   Initialise to zero.
50        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
51         DO i=1-Olx,sNx+Olx-1         DO i=1-Olx,sNx+Olx
52          del2(i,j)=dTdx(i+1,j)-dTdx(i,j)          del2(i,j)= 0. _d 0
53         ENDDO         ENDDO
        del2(sNx+Olx,j)=0.  
54        ENDDO        ENDDO
55    
56  C     ... add difference of meridional fluxes and divide by volume  C     Difference of zonal fluxes & meridional fluxes
57    C      and divide by grid-cell volume
58        DO j=1-Oly,sNy+Oly-1        DO j=1-Oly,sNy+Oly-1
59         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx-1
60          del2(i,j)=recip_rA(i,j,bi,bj)          del2(i,j)=recip_rA(i,j,bi,bj)
61       &           *recip_drF(k)       &           *recip_drF(k)
62       &           *_recip_hFacC(i,j,k,bi,bj)       &           *_recip_hFacC(i,j,k,bi,bj)
63       &           *(       &           *( ( dTdx(i+1,j)-dTdx(i,j) )
64       &              del2(i,j)       &             +( dTdy(i,j+1)-dTdy(i,j) )
      &            +( dTdy(i,j+1)-dTdy(i,j) )  
65       &            )       &            )
66         ENDDO         ENDDO
67        ENDDO        ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22