217 |
ENDIF |
ENDIF |
218 |
#endif |
#endif |
219 |
|
|
|
DO bj=myByLo(myThid),myByHi(myThid) |
|
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
|
|
CALL THSICE_AVE( |
|
|
I bi,bj, myTime, myIter, myThid ) |
|
|
ENDDO |
|
|
ENDDO |
|
|
|
|
220 |
C add a small piece of code to check AddFluid implementation: |
C add a small piece of code to check AddFluid implementation: |
221 |
c#include "thsice_test_addfluid.h" |
c#include "thsice_test_addfluid.h" |
222 |
|
|
248 |
|
|
249 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
250 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
251 |
|
C-- Cumulate time-averaged fields and also fill-up flux diagnostics |
252 |
|
C (if not done in THSICE_DO_ADVECT call) |
253 |
|
#ifdef OLD_THSICE_CALL_SEQUENCE |
254 |
|
IF ( .TRUE. ) THEN |
255 |
|
#else /* OLD_THSICE_CALL_SEQUENCE */ |
256 |
|
IF ( thSIceAdvScheme.LE.0 ) THEN |
257 |
|
#endif /* OLD_THSICE_CALL_SEQUENCE */ |
258 |
|
CALL THSICE_AVE( |
259 |
|
I bi,bj, myTime, myIter, myThid ) |
260 |
|
ENDIF |
261 |
C-- note: If useSEAICE=.true., the stress is computed in seaice_model, |
C-- note: If useSEAICE=.true., the stress is computed in seaice_model, |
262 |
C-- and stressReduction is always set to zero |
C-- and stressReduction is always set to zero |
263 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |