215 |
& *implicitNHPress*implicDiv2DFlow |
& *implicitNHPress*implicDiv2DFlow |
216 |
tmpSurf = tmpSurf / (1. _d 0 + tmpSurf ) |
tmpSurf = tmpSurf / (1. _d 0 + tmpSurf ) |
217 |
ENDIF |
ENDIF |
218 |
ks = ksurfC(i,j,bi,bj) |
ks = kSurfC(i,j,bi,bj) |
219 |
IF ( ks.LE.Nr ) THEN |
IF ( ks.LE.Nr ) THEN |
220 |
aC3d(i,j,ks,bi,bj) = aC3d(i,j,ks,bi,bj) |
aC3d(i,j,ks,bi,bj) = aC3d(i,j,ks,bi,bj) |
221 |
& - freeSurfFac*recip_Bo(i,j,bi,bj) |
& - freeSurfFac*recip_Bo(i,j,bi,bj) |
319 |
_EXCH_XYZ_RS(zML, myThid) |
_EXCH_XYZ_RS(zML, myThid) |
320 |
_EXCH_XYZ_RS(zMU, myThid) |
_EXCH_XYZ_RS(zMU, myThid) |
321 |
|
|
322 |
IF ( debugLevel .GE. debLevB ) THEN |
IF ( debugLevel .GE. debLevC ) THEN |
323 |
CALL WRITE_FLD_XYZ_RS( 'zMC',' ',zMC, 0, myThid ) |
CALL WRITE_FLD_XYZ_RS( 'zMC',' ',zMC, 0, myThid ) |
324 |
CALL WRITE_FLD_XYZ_RS( 'zML',' ',zML, 0, myThid ) |
CALL WRITE_FLD_XYZ_RS( 'zML',' ',zML, 0, myThid ) |
325 |
CALL WRITE_FLD_XYZ_RS( 'zMU',' ',zMU, 0, myThid ) |
CALL WRITE_FLD_XYZ_RS( 'zMU',' ',zMU, 0, myThid ) |