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 |
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 |