54 |
_RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly) |
_RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly) |
55 |
_RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly) |
_RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly) |
56 |
_RL firstResidual,lastResidual |
_RL firstResidual,lastResidual |
57 |
_RL tmpFac, convertVol2Mass |
_RL tmpFac |
58 |
INTEGER numIters |
INTEGER numIters |
59 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
60 |
CEOP |
CEOP |
61 |
|
|
|
IF ( buoyancyRelation .eq. 'OCEANICP' ) THEN |
|
|
convertVol2Mass = gravity*rhoConstFresh |
|
|
ELSE |
|
|
convertVol2Mass = 1. _d 0 |
|
|
ENDIF |
|
|
|
|
62 |
C-- Save previous solution & Initialise Vector solution and source term : |
C-- Save previous solution & Initialise Vector solution and source term : |
63 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
64 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
72 |
ENDDO |
ENDDO |
73 |
ENDDO |
ENDDO |
74 |
IF (useRealFreshWaterFlux) THEN |
IF (useRealFreshWaterFlux) THEN |
75 |
tmpFac = freeSurfFac*convertVol2Mass |
tmpFac = freeSurfFac*convertEmP2rUnit |
76 |
IF (exactConserv) |
IF (exactConserv) |
77 |
& tmpFac = freeSurfFac*convertVol2Mass*implicDiv2DFlow |
& tmpFac = freeSurfFac*convertEmP2rUnit*implicDiv2DFlow |
78 |
DO j=1,sNy |
DO j=1,sNy |
79 |
DO i=1,sNx |
DO i=1,sNx |
80 |
cg2d_b(i,j,bi,bj) = |
cg2d_b(i,j,bi,bj) = |