/[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.6 by cnh, Fri Nov 6 22:44:44 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  C     /==========================================================\  C     /==========================================================\
6  C     | S/R CALC_DIV_GHAT                                        |  C     | S/R CALC_DIV_GHAT                                        |
# 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.6
  ViewVC Help
Powered by ViewVC 1.1.22