/[MITgcm]/MITgcm/model/src/calc_div_ghat.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_div_ghat.F

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

revision 1.7 by adcroft, Tue Dec 15 00:20:34 1998 UTC revision 1.8 by adcroft, Mon Mar 22 15:54:03 1999 UTC
# Line 22  C     == Global variables == Line 22  C     == Global variables ==
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "CG2D.h"  #include "CG2D.h"
25    #ifdef ALLOW_NONHYDROSTATIC
26    #include "CG3D.h"
27    #endif
28    
29  C     == Routine arguments ==  C     == Routine arguments ==
30  C     pH - Hydrostatic pressure  C     pH - Hydrostatic pressure
# Line 69  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM Line 72  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM
72        ENDDO        ENDDO
73        DO j=1,sNy        DO j=1,sNy
74         DO i=1,sNx         DO i=1,sNx
75          cg2d_b(i,j,bi,bj) =  cg2d_b(i,j,bi,bj) +          cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
76       &   pf(i+1,j)-pf(i,j)       &   pf(i+1,j)-pf(i,j)
77         ENDDO         ENDDO
78        ENDDO        ENDDO
79    
80    #ifdef ALLOW_NONHYDROSTATIC
81          IF (nonHydrostatic) THEN
82           DO j=1,sNy
83            DO i=1,sNx
84             cg3d_b(i,j,k,bi,bj) = pf(i+1,j)-pf(i,j)
85            ENDDO
86           ENDDO
87          ENDIF
88    #endif
89    
90        DO j=jMin,jMax        DO j=jMin,jMax
91         DO i=iMin,iMax         DO i=iMin,iMax
92          pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom          pf(i,j) = yA(i,j)*gVNm1(i,j,k,bi,bj) / deltatmom
# Line 82  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM Line 95  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM
95    
96        DO j=1,sNy        DO j=1,sNy
97         DO i=1,sNx         DO i=1,sNx
98          cg2d_b(i,j,bi,bj)   =  cg2d_b(i,j,bi,bj) +          cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) +
99       &   pf(i,j+1)-pf(i,j)       &   pf(i,j+1)-pf(i,j)
100         ENDDO         ENDDO
101        ENDDO        ENDDO
102    
103    #ifdef ALLOW_NONHYDROSTATIC
104          IF (nonHydrostatic) THEN
105           DO j=1,sNy
106            DO i=1,sNx
107             cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +
108         &    pf(i,j+1)-pf(i,j)
109            ENDDO
110           ENDDO
111          ENDIF
112    #endif
113    
114        RETURN        RETURN
115        END        END
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22