46 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
47 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
48 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
49 |
|
#ifdef INCLUDE_CD_CODE |
50 |
etaNm1(i,j,bi,bj) = etaN(i,j,bi,bj) |
etaNm1(i,j,bi,bj) = etaN(i,j,bi,bj) |
51 |
|
#endif |
52 |
cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj) |
cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj) |
53 |
cg2d_b(i,j,bi,bj) = 0. |
cg2d_b(i,j,bi,bj) = 0. |
54 |
#ifdef USE_NATURAL_BCS |
#ifdef USE_NATURAL_BCS |
94 |
& -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom |
& -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom |
95 |
& *( etaN(i,j,bi,bj) |
& *( etaN(i,j,bi,bj) |
96 |
& +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity ) |
& +phi_nh(i,j,1,bi,bj)*horiVertRatio/gravity ) |
|
C-jmc |
|
|
c & -freeSurfFac*_rA(i,j,bi,bj)*recip_Bo(i,j,bi,bj) |
|
|
c & *( cg2d_x(i,j,bi,bj) + phi_nh(i,j,1,bi,bj) ) |
|
|
c & /deltaTMom/deltaTMom |
|
|
C-jmc |
|
97 |
ENDDO |
ENDDO |
98 |
ENDDO |
ENDDO |
99 |
#else |
#else |
100 |
DO j=1,sNy |
IF ( exactConserv ) THEN |
101 |
DO i=1,sNx |
c IF (nonlinFreeSurf.GT.0) THEN |
102 |
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) |
DO j=1,sNy |
103 |
& -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom |
DO i=1,sNx |
104 |
& * etaN(i,j,bi,bj) |
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) |
105 |
|
& -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom |
106 |
|
& * etaH(i,j,bi,bj) |
107 |
|
ENDDO |
108 |
ENDDO |
ENDDO |
109 |
ENDDO |
ELSE |
110 |
|
DO j=1,sNy |
111 |
|
DO i=1,sNx |
112 |
|
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) |
113 |
|
& -freeSurfFac*_rA(i,j,bi,bj)/deltaTMom/deltaTMom |
114 |
|
& * etaN(i,j,bi,bj) |
115 |
|
ENDDO |
116 |
|
ENDDO |
117 |
|
ENDIF |
118 |
#endif |
#endif |
119 |
|
|
120 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |