104 |
CALL PACKAGES_INIT_VARIABLES( myThid ) |
CALL PACKAGES_INIT_VARIABLES( myThid ) |
105 |
|
|
106 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
107 |
C-- Compute the surface level thickness, function of eta_n-1 |
C-- Compute the surface level thickness <-- function of etaH(n) |
108 |
C and modify hFac(C,W,S) accordingly : |
C and modify hFac(C,W,S) accordingly : |
109 |
IF (nonlinFreeSurf.GT.0) THEN |
IF (nonlinFreeSurf.GT.0) THEN |
110 |
CALL CALC_SURF_DR(etaNm1, startTime, nIter0, myThid ) |
CALL CALC_SURF_DR(etaH, startTime, nIter0, myThid ) |
111 |
CALL UPDATE_SURF_DR( startTime, nIter0, myThid ) |
CALL UPDATE_SURF_DR( startTime, nIter0, myThid ) |
112 |
ENDIF |
ENDIF |
113 |
C- update also CG2D matrix (and preconditioner) |
C- update also CG2D matrix (and preconditioner) |
130 |
ENDIF |
ENDIF |
131 |
#endif /* ALLOW_TIMEAVE */ |
#endif /* ALLOW_TIMEAVE */ |
132 |
|
|
|
#ifdef EXACT_CONSERV |
|
|
IF (exactConserv) THEN |
|
|
C-- Compute the initial Barotropic Flow Divergence : |
|
|
DO bj=myByLo(myThid),myByHi(myThid) |
|
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
|
|
CALL CALC_EXACT_ETA( bi,bj, uVel,vVel, |
|
|
I startTime, nIter0, myThid ) |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
|
#endif /* EXACT_CONSERV */ |
|
|
|
|
133 |
RETURN |
RETURN |
134 |
END |
END |