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, |
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 |
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 |