479 |
C recomputation. It *is* differentiable, if you need it. |
C recomputation. It *is* differentiable, if you need it. |
480 |
C Edit GAD_OPTIONS.h and #define DISABLE_MULTIDIM_ADVECTION to |
C Edit GAD_OPTIONS.h and #define DISABLE_MULTIDIM_ADVECTION to |
481 |
C disable this section of code. |
C disable this section of code. |
482 |
IF (tempAdvection .AND. tempMultiDimAdvec) THEN |
IF (tempMultiDimAdvec) THEN |
483 |
CALL GAD_ADVECTION(bi,bj,tempAdvScheme,GAD_TEMPERATURE, |
CALL GAD_ADVECTION(bi,bj,tempAdvScheme,GAD_TEMPERATURE, |
484 |
U theta,gT, |
U theta,gT, |
485 |
I myTime,myIter,myThid) |
I myTime,myIter,myThid) |
486 |
ENDIF |
ENDIF |
487 |
IF (saltAdvection .AND. saltMultiDimAdvec) THEN |
IF (saltMultiDimAdvec) THEN |
488 |
CALL GAD_ADVECTION(bi,bj,saltAdvScheme,GAD_SALINITY, |
CALL GAD_ADVECTION(bi,bj,saltAdvScheme,GAD_SALINITY, |
489 |
U salt,gS, |
U salt,gS, |
490 |
I myTime,myIter,myThid) |
I myTime,myIter,myThid) |
705 |
IF ( useAIM ) THEN |
IF ( useAIM ) THEN |
706 |
CALL AIM_AIM2DYN_EXCHANGES( myTime, myThid ) |
CALL AIM_AIM2DYN_EXCHANGES( myTime, myThid ) |
707 |
ENDIF |
ENDIF |
|
_EXCH_XYZ_R8(gT,myThid) |
|
|
_EXCH_XYZ_R8(gS,myThid) |
|
|
#else |
|
|
IF (staggerTimeStep.AND.useCubedSphereExchange) THEN |
|
|
_EXCH_XYZ_R8(gT,myThid) |
|
|
_EXCH_XYZ_R8(gS,myThid) |
|
|
ENDIF |
|
708 |
#endif /* ALLOW_AIM */ |
#endif /* ALLOW_AIM */ |
709 |
|
IF ( staggerTimeStep ) THEN |
710 |
|
IF ( useAIM .OR. useCubedSphereExchange ) THEN |
711 |
|
IF (tempStepping) _EXCH_XYZ_R8(gT,myThid) |
712 |
|
IF (saltStepping) _EXCH_XYZ_R8(gS,myThid) |
713 |
|
ELSEIF ( useGMRedi .AND. Oly.LT.4 ) THEN |
714 |
|
c .AND. GM_AdvForm .AND. .NOT.GM_AdvSeparate ) THEN |
715 |
|
IF (tempMultiDimAdvec) _EXCH_XYZ_R8(gT,myThid) |
716 |
|
IF (saltMultiDimAdvec) _EXCH_XYZ_R8(gS,myThid) |
717 |
|
ENDIF |
718 |
|
ENDIF |
719 |
|
|
720 |
#ifndef DISABLE_DEBUGMODE |
#ifndef DISABLE_DEBUGMODE |
721 |
If (debugMode) THEN |
If (debugMode) THEN |