353 |
C-- Calculate gradients of potential density for isoneutral |
C-- Calculate gradients of potential density for isoneutral |
354 |
C slope terms (e.g. GM/Redi tensor or IVDC diffusivity) |
C slope terms (e.g. GM/Redi tensor or IVDC diffusivity) |
355 |
IF ( useGMRedi .OR. (k.GT.1 .AND. ivdc_kappa.NE.0.) |
IF ( useGMRedi .OR. (k.GT.1 .AND. ivdc_kappa.NE.0.) |
356 |
& .OR. usePLUME .OR. doDiagsRho.GE.1 ) THEN |
& .OR. useSaltPlume .OR. doDiagsRho.GE.1 ) THEN |
357 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
358 |
IF ( debugLevel .GE. debLevB ) |
IF ( debugLevel .GE. debLevB ) |
359 |
& CALL DEBUG_CALL('FIND_RHO',myThid) |
& CALL DEBUG_CALL('FIND_RHO',myThid) |
432 |
ENDIF |
ENDIF |
433 |
|
|
434 |
#ifdef ALLOW_SALT_PLUME |
#ifdef ALLOW_SALT_PLUME |
435 |
CALL CALC_SALT_PLUME_DEPTH( rhoK, sigmaR, |
IF ( useSaltPlume ) THEN |
436 |
|
CALL SALT_PLUME_CALC_DEPTH( rhoK, sigmaR, |
437 |
& bi, bj, myTime, myIter, myThid ) |
& bi, bj, myTime, myIter, myThid ) |
438 |
|
ENDIF |
439 |
#endif |
#endif |
440 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
441 |
IF ( doDiagsRho.GE.1 ) THEN |
IF ( doDiagsRho.GE.1 ) THEN |