/[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.8 by adcroft, Mon Mar 22 15:54:03 1999 UTC revision 1.9 by adcroft, Tue May 18 17:42:01 1999 UTC
# Line 5  C $Header$ Line 5  C $Header$
5  C     /==========================================================\  C     /==========================================================\
6  C     | S/R CALC_DIV_GHAT                                        |  C     | S/R CALC_DIV_GHAT                                        |
7  C     | o Form the right hand-side of the surface pressure eqn.  |  C     | o Form the right hand-side of the surface pressure eqn.  |
8    C     |==========================================================|
9  C     \==========================================================/  C     \==========================================================/
10        SUBROUTINE CALC_DIV_GHAT(        SUBROUTINE CALC_DIV_GHAT(
11       I        bi,bj,iMin,iMax,jMin,jMax,       I        bi,bj,iMin,iMax,jMin,jMax,
# Line 27  C     == Global variables == Line 28  C     == Global variables ==
28  #endif  #endif
29    
30  C     == Routine arguments ==  C     == Routine arguments ==
 C     pH - Hydrostatic pressure  
31  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
32  C                                      results will be set.  C                                      results will be set.
33  C     kUp, kDown, kM1                - Index for upper and lower layers.  C     k                              - Index of layer.
34    C     xA, yA                         - Cell face areas
35  C     myThid - Instance number for this innvocation of CALC_MOM_RHS  C     myThid - Instance number for this innvocation of CALC_MOM_RHS
36        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
37        INTEGER K        INTEGER K
# Line 51  C     pulls div(U) + dh/dt back toward z Line 52  C     pulls div(U) + dh/dt back toward z
52  C      Initialise source term on first pass  C      Initialise source term on first pass
53         DO j=jMin,jMax         DO j=jMin,jMax
54          DO i=iMin,iMax          DO i=iMin,iMax
 C        cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj)  
 C    &    -freeSurfFac*_rA(i,j,bi,bj)*  
 C    &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom  
55           cg2d_b(i,j,bi,bj) =           cg2d_b(i,j,bi,bj) =
56       &     freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*(       &     freeSurfFac*_rA(i,j,bi,bj)*horiVertRatio*(
57       &      -cg2d_x(I,J,bi,bj)/deltaTMom/deltaTMom       &      -cg2d_x(I,J,bi,bj)/deltaTMom/deltaTMom
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22