254 |
ENDIF |
ENDIF |
255 |
#endif /* ALLOW_SHELFICE */ |
#endif /* ALLOW_SHELFICE */ |
256 |
|
|
257 |
|
#ifdef ALLOW_ICEFRONT |
258 |
|
IF ( useICEFRONT .AND. fluidIsWater ) THEN |
259 |
|
#ifdef ALLOW_DEBUG |
260 |
|
IF ( debugLevel .GE. debLevB ) |
261 |
|
& CALL DEBUG_CALL('ICEFRONT_THERMODYNAMICS',myThid) |
262 |
|
#endif |
263 |
|
C compute temperature and (virtual) salt flux at the |
264 |
|
C ice-front ocean interface |
265 |
|
CALL TIMER_START('ICEFRONT_THERMODYNAMICS [DO_OCEANIC_PHYS]', |
266 |
|
& myThid) |
267 |
|
CALL ICEFRONT_THERMODYNAMICS( myTime, myIter, myThid ) |
268 |
|
CALL TIMER_STOP( 'ICEFRONT_THERMODYNAMICS [DO_OCEANIC_PHYS]', |
269 |
|
& myThid) |
270 |
|
ENDIF |
271 |
|
#endif /* ALLOW_ICEFRONT */ |
272 |
|
|
273 |
C-- Freeze water at the surface |
C-- Freeze water at the surface |
274 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
275 |
CADJ STORE theta = comlev1, key = ikey_dynamics, |
CADJ STORE theta = comlev1, key = ikey_dynamics, |