/[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.4 by cnh, Tue Sep 8 01:37:49 1998 UTC revision 1.5 by cnh, Wed Oct 28 03:11:36 1998 UTC
# Line 39  C     == Local variables == Line 39  C     == Local variables ==
39        INTEGER i,j        INTEGER i,j
40        _RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41    
 CcnhDEbugStarts  
       INTEGER i1, j1  
 CcnhDebugEnds  
   
42  C--   Pressure equation source term  C--   Pressure equation source term
43  C     Term is the vertical integral of the divergence of the  C     Term is the vertical integral of the divergence of the
44  C     time tendency terms along with a relaxation term that  C     time tendency terms along with a relaxation term that
# Line 82  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM Line 78  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTM
78    
79        DO j=1,sNy        DO j=1,sNy
80         DO i=1,sNx         DO i=1,sNx
81          cg2d_b(i,j,bi,bj) =  cg2d_b(i,j,bi,bj) +          cg2d_b(i,j,bi,bj)   =  cg2d_b(i,j,bi,bj) +
82       &   pf(i,j+1)-pf(i,j)       &   pf(i,j+1)-pf(i,j)
83         ENDDO         ENDDO
84        ENDDO        ENDDO
85    
 CcnhDebugSTarts  
 C     IF ( K .EQ. 10 .OR. K .EQ. 1 ) THEN  
 C      i1 = 54  
 C      j1 = 26  
 C      WRITE(0,*) ' @ I = ', i1, ' , J = ', j1, ' K = ', K  
 C      WRITE(0,*) '  cg2d_b = ', cg2d_b(i1,j1,bi,bj)  
 C     ENDIF  
 CcnhDebugEnds  
   
86    
87        RETURN        RETURN
88        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22