152 |
C Starting values of U, V, W, temp., salt. and tendency terms |
C Starting values of U, V, W, temp., salt. and tendency terms |
153 |
C are set here. Fields are either set to default or read from |
C are set here. Fields are either set to default or read from |
154 |
C stored files. |
C stored files. |
|
#ifndef ALLOW_OFFLINE |
|
155 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
156 |
IF (debugMode) CALL DEBUG_CALL('INI_FIELDS',myThid) |
IF (debugMode) CALL DEBUG_CALL('INI_FIELDS',myThid) |
157 |
#endif |
#endif |
158 |
CALL INI_FIELDS( myThid ) |
CALL INI_FIELDS( myThid ) |
|
#endif |
|
159 |
_BARRIER |
_BARRIER |
160 |
|
|
161 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
184 |
#endif |
#endif |
185 |
|
|
186 |
C-- Initial conditions are convectively adjusted (for historical reasons) |
C-- Initial conditions are convectively adjusted (for historical reasons) |
|
#ifndef ALLOW_OFFLINE |
|
187 |
IF ( startTime .EQ. baseTime .AND. cAdjFreq .NE. 0. ) THEN |
IF ( startTime .EQ. baseTime .AND. cAdjFreq .NE. 0. ) THEN |
188 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
189 |
IF (debugMode) CALL DEBUG_CALL('CONVECTIVE_ADJUSTMENT_INI',myThid) |
IF (debugMode) CALL DEBUG_CALL('CONVECTIVE_ADJUSTMENT_INI',myThid) |
202 |
ENDDO |
ENDDO |
203 |
ENDDO |
ENDDO |
204 |
_BARRIER |
_BARRIER |
205 |
END IF |
ENDIF |
|
#endif |
|
206 |
|
|
207 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
208 |
C-- Compute the surface level thickness <-- function of etaH(n) |
C-- Compute the surface level thickness <-- function of etaH(n) |
229 |
ENDIF |
ENDIF |
230 |
#endif /* NONLIN_FRSURF */ |
#endif /* NONLIN_FRSURF */ |
231 |
|
|
|
#ifndef ALLOW_OFFLINE |
|
232 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
233 |
IF (debugMode) CALL DEBUG_CALL('INTEGR_CONTINUITY',myThid) |
IF (debugMode) CALL DEBUG_CALL('INTEGR_CONTINUITY',myThid) |
234 |
#endif |
#endif |
238 |
C-- Integrate continuity vertically for vertical velocity |
C-- Integrate continuity vertically for vertical velocity |
239 |
CALL INTEGR_CONTINUITY( bi, bj, uVel, vVel, |
CALL INTEGR_CONTINUITY( bi, bj, uVel, vVel, |
240 |
I startTime, nIter0, myThid ) |
I startTime, nIter0, myThid ) |
241 |
|
|
242 |
ENDDO |
ENDDO |
243 |
ENDDO |
ENDDO |
|
#endif |
|
244 |
|
|
245 |
#ifdef EXACT_CONSERV |
#ifdef EXACT_CONSERV |
246 |
IF ( exactConserv ) THEN |
IF ( exactConserv ) THEN |
290 |
IF (taveFreq.GT.0.) THEN |
IF (taveFreq.GT.0.) THEN |
291 |
CALL TIMEAVE_STATVARS(startTime, nIter0, bi, bj, myThid) |
CALL TIMEAVE_STATVARS(startTime, nIter0, bi, bj, myThid) |
292 |
ENDIF |
ENDIF |
293 |
cswdptr -- add --- |
|
294 |
#ifdef ALLOW_PTRACERS |
#ifdef ALLOW_PTRACERS |
295 |
che3 needs an IF ( usePTRACERS ) THEN |
IF ( usePTRACERS ) THEN |
296 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
297 |
IF (debugMode) CALL DEBUG_CALL('PTRACERS_STATVARS',myThid) |
IF (debugMode) CALL DEBUG_CALL('PTRACERS_STATVARS',myThid) |
298 |
#endif |
#endif |
299 |
CALL PTRACERS_STATVARS(startTime, nIter0, bi, bj, myThid) |
CALL PTRACERS_STATVARS(startTime, nIter0, bi, bj, myThid) |
300 |
#endif |
ENDIF |
301 |
cswdptr -- end add --- |
#endif /* ALLOW_PTRACERS */ |
302 |
|
|
303 |
C-- end bi,bj loop. |
C-- end bi,bj loop. |
304 |
ENDDO |
ENDDO |
305 |
ENDDO |
ENDDO |
306 |
#endif /* ALLOW_TIMEAVE */ |
#endif /* ALLOW_TIMEAVE */ |
307 |
|
|
|
C AMM |
|
308 |
#ifdef ALLOW_GRIDALT |
#ifdef ALLOW_GRIDALT |
309 |
if (useGRIDALT) then |
if (useGRIDALT) then |
310 |
CALL TIMER_START('GRIDALT_UPDATE [INITIALISE_VARIA]',mythid) |
CALL TIMER_START('GRIDALT_UPDATE [INITIALISE_VARIA]',mythid) |
312 |
CALL TIMER_STOP ('GRIDALT_UPDATE [INITIALISE_VARIA]',mythid) |
CALL TIMER_STOP ('GRIDALT_UPDATE [INITIALISE_VARIA]',mythid) |
313 |
endif |
endif |
314 |
#endif |
#endif |
|
C AMM |
|
315 |
|
|
316 |
C-- Fill in overlap regions for wVel : |
C-- Fill in overlap regions for wVel : |
317 |
#ifndef ALLOW_OFFLINE |
_EXCH_XYZ_R8(wVel,myThid) |
|
_EXCH_XYZ_R8(wVel,myThid) |
|
|
#endif |
|
318 |
|
|
319 |
C-- Finally summarise the model state |
C-- Finally summarise the model state |
320 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |