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() |
485 |
kkey = (idynkey-1)*Nr + k |
kkey = (idynkey-1)*Nr + k |
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) |
489 |
|
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 |