3 |
|
|
4 |
#include "PACKAGES_CONFIG.h" |
#include "PACKAGES_CONFIG.h" |
5 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
6 |
|
#ifdef ALLOW_AUTODIFF |
7 |
|
# include "AUTODIFF_OPTIONS.h" |
8 |
|
#endif |
9 |
#ifdef ALLOW_MOM_COMMON |
#ifdef ALLOW_MOM_COMMON |
10 |
# include "MOM_COMMON_OPTIONS.h" |
# include "MOM_COMMON_OPTIONS.h" |
11 |
#endif |
#endif |
92 |
#ifdef ALLOW_CD_CODE |
#ifdef ALLOW_CD_CODE |
93 |
# include "CD_CODE_VARS.h" |
# include "CD_CODE_VARS.h" |
94 |
#endif |
#endif |
95 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF |
96 |
# include "tamc.h" |
# include "tamc.h" |
97 |
# include "tamc_keys.h" |
# include "tamc_keys.h" |
98 |
# include "FFIELDS.h" |
# include "FFIELDS.h" |
114 |
# ifdef ALLOW_MOM_FLUXFORM |
# ifdef ALLOW_MOM_FLUXFORM |
115 |
# include "MOM_FLUXFORM.h" |
# include "MOM_FLUXFORM.h" |
116 |
# endif |
# endif |
117 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF */ |
118 |
|
|
119 |
C !CALLING SEQUENCE: |
C !CALLING SEQUENCE: |
120 |
C DYNAMICS() |
C DYNAMICS() |
486 |
CADJ STORE totPhiHyd (:,:,k,bi,bj) |
CADJ STORE totPhiHyd (:,:,k,bi,bj) |
487 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
488 |
CADJ STORE phiHydLow (:,:,bi,bj) |
CADJ STORE phiHydLow (:,:,bi,bj) |
489 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
490 |
CADJ STORE theta (:,:,k,bi,bj) |
CADJ STORE theta (:,:,k,bi,bj) |
491 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
492 |
CADJ STORE salt (:,:,k,bi,bj) |
CADJ STORE salt (:,:,k,bi,bj) |
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 */ |