101 |
C-- Initialize variable data for packages |
C-- Initialize variable data for packages |
102 |
CALL PACKAGES_INIT_VARIABLES( myThid ) |
CALL PACKAGES_INIT_VARIABLES( myThid ) |
103 |
|
|
104 |
|
#ifdef NONLIN_FRSURF |
105 |
|
C-- Compute the surface level thickness, function of eta_n-1 |
106 |
|
C and modify hFac(C,W,S) accordingly : |
107 |
|
IF (nonlinFreeSurf.GT.0) THEN |
108 |
|
CALL CALC_SURF_DR(etaNm1, startTime, nIter0, myThid ) |
109 |
|
CALL UPDATE_SURF_DR( startTime, nIter0, myThid ) |
110 |
|
ENDIF |
111 |
|
C- update also CG2D matrix (and preconditioner) |
112 |
|
IF ( nonlinFreeSurf.GT.2) THEN |
113 |
|
CALL UPDATE_CG2D( startTime, nIter0, myThid ) |
114 |
|
ENDIF |
115 |
|
#endif |
116 |
|
|
117 |
C-- Finally summarise the model state |
C-- Finally summarise the model state |
118 |
CALL STATE_SUMMARY( myThid ) |
CALL STATE_SUMMARY( myThid ) |
119 |
|
|