172 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
173 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
174 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
175 |
IF ( nonHydrostatic .AND. zeroPsNH ) THEN |
IF ( use3Dsolver .AND. zeroPsNH ) THEN |
176 |
DO j=1,sNy |
DO j=1,sNy |
177 |
DO i=1,sNx |
DO i=1,sNx |
178 |
ks = ksurfC(i,j,bi,bj) |
ks = ksurfC(i,j,bi,bj) |
186 |
ENDIF |
ENDIF |
187 |
ENDDO |
ENDDO |
188 |
ENDDO |
ENDDO |
189 |
ELSEIF ( nonHydrostatic ) THEN |
ELSEIF ( use3Dsolver ) THEN |
190 |
DO j=1,sNy |
DO j=1,sNy |
191 |
DO i=1,sNx |
DO i=1,sNx |
192 |
ks = ksurfC(i,j,bi,bj) |
ks = ksurfC(i,j,bi,bj) |
313 |
ENDDO |
ENDDO |
314 |
|
|
315 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
316 |
IF ( nonHydrostatic ) THEN |
IF ( use3Dsolver ) THEN |
317 |
|
|
318 |
C-- Solve for a three-dimensional pressure term (NH or IGW or both ). |
C-- Solve for a three-dimensional pressure term (NH or IGW or both ). |
319 |
C see CG3D.h for the interface to this routine. |
C see CG3D.h for the interface to this routine. |