15 |
I myIter,myTime,myThid ) |
I myIter,myTime,myThid ) |
16 |
|
|
17 |
C !DESCRIPTION: |
C !DESCRIPTION: |
18 |
C Calculates tendancy for passive tracers and integrates forward |
C Calculates tendency for passive tracers and integrates forward |
19 |
C in time. |
C in time. |
20 |
|
|
21 |
C !USES: =============================================================== |
C !USES: =============================================================== |
63 |
#ifdef ALLOW_PTRACERS |
#ifdef ALLOW_PTRACERS |
64 |
|
|
65 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
66 |
C i,j,k,bi,bj,iTracer :: loop indices |
C iTracer :: tracer index |
67 |
C iMin,iMax,jMin,jMax :: loop ranges |
C iMin,iMax,jMin,jMax :: loop ranges |
68 |
C kUp,kDown :: toggle indices for even/odd level fluxes |
C kUp,kDown :: toggle indices for even/odd level fluxes |
69 |
C km1 :: =min(1,k-1) |
C km1 :: =min(1,k-1) |
70 |
C GAD_TR :: passive tracer id (GAD_TR1+iTracer-1) |
C GAD_TR :: passive tracer id (GAD_TR1+iTracer-1) |
71 |
INTEGER i,j,iTracer |
INTEGER iTracer |
72 |
INTEGER iMin,iMax,jMin,jMax |
INTEGER iMin,iMax,jMin,jMax |
73 |
INTEGER kUp,kDown,km1 |
INTEGER kUp,kDown,km1 |
74 |
INTEGER GAD_TR |
INTEGER GAD_TR |
103 |
I PTRACERS_diffKh(iTracer), |
I PTRACERS_diffKh(iTracer), |
104 |
I PTRACERS_diffK4(iTracer), |
I PTRACERS_diffK4(iTracer), |
105 |
I KappaRtr(1-Olx,1-Oly,iTracer), |
I KappaRtr(1-Olx,1-Oly,iTracer), |
106 |
|
I gPtrNM1(1-Olx,1-Oly,1,1,1,iTracer), |
107 |
I pTracer(1-Olx,1-Oly,1,1,1,iTracer), |
I pTracer(1-Olx,1-Oly,1,1,1,iTracer), |
108 |
I GAD_TR, |
I GAD_TR, |
109 |
I PTRACERS_advScheme(iTracer), |
I PTRACERS_advScheme(iTracer), |
121 |
I surfaceForcingPtr(1-Olx,1-Oly,1,1,iTracer), |
I surfaceForcingPtr(1-Olx,1-Oly,1,1,iTracer), |
122 |
I myIter,myTime,myThid) |
I myIter,myTime,myThid) |
123 |
|
|
124 |
C If using Adams-Bashforth II, then extrapolate tendancies |
C If using Adams-Bashforth II, then extrapolate tendencies |
125 |
C gPtr is now the tracer tendency for explicit advection/diffusion |
C gPtr is now the tracer tendency for explicit advection/diffusion |
126 |
IF ( PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND |
IF ( PTRACERS_advScheme(iTracer).EQ.ENUM_CENTERED_2ND |
127 |
& .OR.PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD |
& .OR.PTRACERS_advScheme(iTracer).EQ.ENUM_UPWIND_3RD |