362 |
& .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN |
& .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN |
363 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
364 |
IF ( debugLevel .GE. debLevB ) |
IF ( debugLevel .GE. debLevB ) |
365 |
& CALL DEBUG_CALL('FIND_RHO',myThid) |
& CALL DEBUG_CALL('FIND_RHO_2D',myThid) |
366 |
#endif |
#endif |
367 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
368 |
CADJ STORE theta(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE theta(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
369 |
CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
370 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
371 |
CALL FIND_RHO( |
CALL FIND_RHO_2D( |
372 |
I bi, bj, iMin, iMax, jMin, jMax, k, k, |
I iMin, iMax, jMin, jMax, k, |
373 |
I theta, salt, |
I theta(1-OLx,1-OLy,k,bi,bj), |
374 |
O rhoK, |
I salt (1-OLx,1-OLy,k,bi,bj), |
375 |
I myThid ) |
O rhoK, |
376 |
|
I k, bi, bj, myThid ) |
377 |
|
|
378 |
IF (k.GT.1) THEN |
IF (k.GT.1) THEN |
379 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
380 |
CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
381 |
CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte |
382 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
383 |
CALL FIND_RHO( |
CALL FIND_RHO_2D( |
384 |
I bi, bj, iMin, iMax, jMin, jMax, k-1, k, |
I iMin, iMax, jMin, jMax, k, |
385 |
I theta, salt, |
I theta(1-OLx,1-OLy,k-1,bi,bj), |
386 |
O rhoKm1, |
I salt (1-OLx,1-OLy,k-1,bi,bj), |
387 |
I myThid ) |
O rhoKm1, |
388 |
|
I k-1, bi, bj, myThid ) |
389 |
ENDIF |
ENDIF |
390 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
391 |
IF ( debugLevel .GE. debLevB ) |
IF ( debugLevel .GE. debLevB ) |