534 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
535 |
|
|
536 |
C-- Integrate hydrostatic balance for phiHyd with BC of phiHyd(z=0)=0 |
C-- Integrate hydrostatic balance for phiHyd with BC of phiHyd(z=0)=0 |
537 |
IF ( implicitIntGravWave ) THEN |
CALL CALC_PHI_HYD( |
|
CALL CALC_PHI_HYD( |
|
|
I bi,bj,iMin,iMax,jMin,jMax,k, |
|
|
I gT, gS, |
|
|
U phiHydF, |
|
|
O phiHydC, dPhiHydX, dPhiHydY, |
|
|
I myTime, myIter, myThid ) |
|
|
ELSE |
|
|
CALL CALC_PHI_HYD( |
|
538 |
I bi,bj,iMin,iMax,jMin,jMax,k, |
I bi,bj,iMin,iMax,jMin,jMax,k, |
539 |
I theta, salt, |
I theta, salt, |
540 |
U phiHydF, |
U phiHydF, |
541 |
O phiHydC, dPhiHydX, dPhiHydY, |
O phiHydC, dPhiHydX, dPhiHydY, |
542 |
I myTime, myIter, myThid ) |
I myTime, myIter, myThid ) |
|
ENDIF |
|
543 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
544 |
IF ( dPhiHydDiagIsOn ) THEN |
IF ( dPhiHydDiagIsOn ) THEN |
545 |
tmpFac = -1. _d 0 |
tmpFac = -1. _d 0 |
649 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |
650 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
651 |
I -1, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj), |
I -1, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj), |
652 |
U gU, |
U gU(1-OLx,1-OLy,1,bi,bj), |
653 |
I myThid ) |
I myThid ) |
654 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
655 |
CADJ STORE gV(:,:,:,bi,bj) = comlev1_bibj , key=idynkey, byte=isbyte |
CADJ STORE gV(:,:,:,bi,bj) = comlev1_bibj , key=idynkey, byte=isbyte |
657 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |
658 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
659 |
I -2, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj), |
I -2, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj), |
660 |
U gV, |
U gV(1-OLx,1-OLy,1,bi,bj), |
661 |
I myThid ) |
I myThid ) |
662 |
ENDIF |
ENDIF |
663 |
|
|
681 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |
682 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
683 |
I 0, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj), |
I 0, KappaRU, recip_hFacW(1-OLx,1-OLy,1,bi,bj), |
684 |
U vVelD, |
U vVelD(1-OLx,1-OLy,1,bi,bj), |
685 |
I myThid ) |
I myThid ) |
686 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
687 |
CADJ STORE uVelD(:,:,:,bi,bj) = comlev1_bibj , key=idynkey, byte=isbyte |
CADJ STORE uVelD(:,:,:,bi,bj) = comlev1_bibj , key=idynkey, byte=isbyte |
689 |
CALL IMPLDIFF( |
CALL IMPLDIFF( |
690 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |
691 |
I 0, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj), |
I 0, KappaRV, recip_hFacS(1-OLx,1-OLy,1,bi,bj), |
692 |
U uVelD, |
U uVelD(1-OLx,1-OLy,1,bi,bj), |
693 |
I myThid ) |
I myThid ) |
694 |
ENDIF |
ENDIF |
695 |
#endif /* ALLOW_CD_CODE */ |
#endif /* ALLOW_CD_CODE */ |