1 |
|
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
5 |
|
|
97 |
C | |
C | |
98 |
C |-- CALC_PHI_HYD |
C |-- CALC_PHI_HYD |
99 |
C | |
C | |
100 |
|
C |-- STORE_PRESSURE |
101 |
|
C | |
102 |
C |-- MOM_FLUXFORM |
C |-- MOM_FLUXFORM |
103 |
C | |
C | |
104 |
C |-- MOM_VECINV |
C |-- MOM_VECINV |
152 |
_RL phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL phiSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
153 |
_RL KappaRU (1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr) |
_RL KappaRU (1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr) |
154 |
_RL KappaRV (1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr) |
_RL KappaRV (1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr) |
|
_RL sigmaX (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
|
|
_RL sigmaY (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
|
|
_RL sigmaR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
|
155 |
|
|
156 |
INTEGER iMin, iMax |
INTEGER iMin, iMax |
157 |
INTEGER jMin, jMax |
INTEGER jMin, jMax |
224 |
ENDDO |
ENDDO |
225 |
ENDDO |
ENDDO |
226 |
|
|
227 |
|
C-- Call to routine for calculation of |
228 |
|
C Eliassen-Palm-flux-forced U-tendency, |
229 |
|
C if desired: |
230 |
|
#ifdef INCLUDE_EP_FORCING_CODE |
231 |
|
CALL CALC_EP_FORCING(myThid) |
232 |
|
#endif |
233 |
|
|
234 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
235 |
C-- HPF directive to help TAMC |
C-- HPF directive to help TAMC |
236 |
CHPF$ INDEPENDENT |
CHPF$ INDEPENDENT |
348 |
I myThid ) |
I myThid ) |
349 |
ENDIF |
ENDIF |
350 |
|
|
351 |
|
C calculate pressure from phiHyd and store it on common block |
352 |
|
C variable pressure |
353 |
|
CALL STORE_PRESSURE( bi, bj, k, phiHyd, myThid ) |
354 |
|
|
355 |
|
|
356 |
C-- Calculate accelerations in the momentum equations (gU, gV, ...) |
C-- Calculate accelerations in the momentum equations (gU, gV, ...) |
357 |
C and step forward storing the result in gUnm1, gVnm1, etc... |
C and step forward storing the result in gUnm1, gVnm1, etc... |
358 |
IF ( momStepping ) THEN |
IF ( momStepping ) THEN |
399 |
C-- end of dynamics k loop (1:Nr) |
C-- end of dynamics k loop (1:Nr) |
400 |
ENDDO |
ENDDO |
401 |
|
|
|
|
|
|
|
|
402 |
C-- Implicit viscosity |
C-- Implicit viscosity |
403 |
IF (implicitViscosity.AND.momStepping) THEN |
IF (implicitViscosity.AND.momStepping) THEN |
404 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |