563 |
C-- Calculate gradients of potential density for isoneutral |
C-- Calculate gradients of potential density for isoneutral |
564 |
C slope terms (e.g. GM/Redi tensor or IVDC diffusivity) |
C slope terms (e.g. GM/Redi tensor or IVDC diffusivity) |
565 |
IF ( useGMRedi .OR. (k.GT.1 .AND. ivdc_kappa.NE.0.) |
IF ( useGMRedi .OR. (k.GT.1 .AND. ivdc_kappa.NE.0.) |
566 |
|
& .OR. usePP81 .OR. useMY82 .OR. useGGL90 |
567 |
& .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN |
& .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN |
568 |
IF (k.GT.1) THEN |
IF (k.GT.1) THEN |
569 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
743 |
IF (debugMode) CALL DEBUG_CALL('PP81_CALC',myThid) |
IF (debugMode) CALL DEBUG_CALL('PP81_CALC',myThid) |
744 |
#endif |
#endif |
745 |
CALL PP81_CALC( |
CALL PP81_CALC( |
746 |
I bi, bj, myTime, myThid ) |
I bi, bj, sigmaR, myTime, myIter, myThid ) |
747 |
ENDIF |
ENDIF |
748 |
#endif /* ALLOW_PP81 */ |
#endif /* ALLOW_PP81 */ |
749 |
|
|
754 |
IF (debugMode) CALL DEBUG_CALL('MY82_CALC',myThid) |
IF (debugMode) CALL DEBUG_CALL('MY82_CALC',myThid) |
755 |
#endif |
#endif |
756 |
CALL MY82_CALC( |
CALL MY82_CALC( |
757 |
I bi, bj, myTime, myThid ) |
I bi, bj, sigmaR, myTime, myIter, myThid ) |
758 |
ENDIF |
ENDIF |
759 |
#endif /* ALLOW_MY82 */ |
#endif /* ALLOW_MY82 */ |
760 |
|
|
770 |
#endif |
#endif |
771 |
CALL TIMER_START('GGL90_CALC [DO_OCEANIC_PHYS]', myThid) |
CALL TIMER_START('GGL90_CALC [DO_OCEANIC_PHYS]', myThid) |
772 |
CALL GGL90_CALC( |
CALL GGL90_CALC( |
773 |
I bi, bj, myTime, myIter, myThid ) |
I bi, bj, sigmaR, myTime, myIter, myThid ) |
774 |
CALL TIMER_STOP ('GGL90_CALC [DO_OCEANIC_PHYS]', myThid) |
CALL TIMER_STOP ('GGL90_CALC [DO_OCEANIC_PHYS]', myThid) |
775 |
ENDIF |
ENDIF |
776 |
#endif /* ALLOW_GGL90 */ |
#endif /* ALLOW_GGL90 */ |