/[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.2 by cnh, Sat Aug 22 17:51:07 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 52  C        cg2d_b(i,j,bi,bj) = cg2d_b(i,j, Line 52  C        cg2d_b(i,j,bi,bj) = cg2d_b(i,j,
52  C    &    -freeSurfFac*_rA(i,j,bi,bj)*  C    &    -freeSurfFac*_rA(i,j,bi,bj)*
53  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom  C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom
54           cg2d_b(i,j,bi,bj) =           cg2d_b(i,j,bi,bj) =
55       &    -freeSurfFac*_rA(i,j,bi,bj)*       &    -freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*
56       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom
57          ENDDO          ENDDO
58         ENDDO         ENDDO
# Line 78  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    
86    
87        RETURN        RETURN
88        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22