82 |
cg2dNorm = myNorm |
cg2dNorm = myNorm |
83 |
_BEGIN_MASTER( myThid ) |
_BEGIN_MASTER( myThid ) |
84 |
CcnhDebugStarts |
CcnhDebugStarts |
85 |
WRITE(msgBuf,'(A,F)') '// CG2D normalisation factor = ', cg2dNorm |
WRITE(msgBuf,'(A,E40.25)') '// CG2D normalisation factor = ', cg2dNorm |
86 |
CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
87 |
WRITE(msgBuf,*) ' ' |
WRITE(msgBuf,*) ' ' |
88 |
CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
CALL PRINT_MESSAGE( msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
133 |
aC = -( |
aC = -( |
134 |
& aW2d(I,J,bi,bj) + aW2d(I+1,J ,bi,bj) |
& aW2d(I,J,bi,bj) + aW2d(I+1,J ,bi,bj) |
135 |
& +aS2d(I,J,bi,bj) + aS2D(I ,J+1,bi,bj) |
& +aS2d(I,J,bi,bj) + aS2D(I ,J+1,bi,bj) |
136 |
& +freeSurfFac*myNorm* Gravity*rhoConst* |
& +freeSurfFac*myNorm* horiVertRatio* |
137 |
& rA(I,J,bi,bj)/deltaTMom/deltaTMom |
& rA(I,J,bi,bj)/deltaTMom/deltaTMom |
138 |
& ) |
& ) |
139 |
aCs = -( |
aCs = -( |
140 |
& aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj) |
& aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj) |
141 |
& +aS2d(I,J-1,bi,bj) + aS2d(I ,J ,bi,bj) |
& +aS2d(I,J-1,bi,bj) + aS2d(I ,J ,bi,bj) |
142 |
& +freeSurfFac*myNorm* Gravity*rhoConst* |
& +freeSurfFac*myNorm* horiVertRatio* |
143 |
& rA(I,J-1,bi,bj)/deltaTMom/deltaTMom |
& rA(I,J-1,bi,bj)/deltaTMom/deltaTMom |
144 |
& ) |
& ) |
145 |
aCw = -( |
aCw = -( |
146 |
& aW2d(I-1,J,bi,bj) + aW2d(I ,J ,bi,bj) |
& aW2d(I-1,J,bi,bj) + aW2d(I ,J ,bi,bj) |
147 |
& +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj) |
& +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj) |
148 |
& +freeSurfFac*myNorm* Gravity*rhoConst* |
& +freeSurfFac*myNorm* horiVertRatio* |
149 |
& rA(I-1,J,bi,bj)/deltaTMom/deltaTMom |
& rA(I-1,J,bi,bj)/deltaTMom/deltaTMom |
150 |
& ) |
& ) |
151 |
IF ( aC .EQ. 0. ) THEN |
IF ( aC .EQ. 0. ) THEN |