155 |
C-- Compute the surface level thickness <-- function of etaH(n) |
C-- Compute the surface level thickness <-- function of etaH(n) |
156 |
C and modify hFac(C,W,S) accordingly : |
C and modify hFac(C,W,S) accordingly : |
157 |
IF (nonlinFreeSurf.GT.0) THEN |
IF (nonlinFreeSurf.GT.0) THEN |
158 |
CALL CALC_SURF_DR(etaH, startTime, nIter0, myThid ) |
CALL CALC_SURF_DR(etaH, startTime, -1 , myThid ) |
159 |
CALL UPDATE_SURF_DR( startTime, nIter0, myThid ) |
CALL UPDATE_SURF_DR( startTime, nIter0, myThid ) |
160 |
ENDIF |
ENDIF |
161 |
C- update also CG2D matrix (and preconditioner) |
C- update also CG2D matrix (and preconditioner) |