223 |
CALL GAD_C2_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
CALL GAD_C2_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
224 |
I dTtracerLev, rTrans, |
I dTtracerLev, rTrans, |
225 |
U b5d, c5d, d5d, |
U b5d, c5d, d5d, |
226 |
I myThid) |
I myThid ) |
227 |
ELSEIF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN |
ELSEIF ( advectionScheme.EQ.ENUM_UPWIND_1RST |
228 |
|
& .OR. advectionScheme.EQ.ENUM_DST2 ) THEN |
229 |
|
diagonalNumber = 3 |
230 |
|
CALL GAD_DST2U1_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
231 |
|
I advectionScheme, dTtracerLev, rTrans, |
232 |
|
U b5d, c5d, d5d, |
233 |
|
I myThid ) |
234 |
|
ELSEIF ( advectionScheme.EQ.ENUM_FLUX_LIMIT ) THEN |
235 |
diagonalNumber = 3 |
diagonalNumber = 3 |
236 |
CALL GAD_FLUXLIMIT_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
CALL GAD_FLUXLIMIT_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
237 |
I dTtracerLev, rTrans, localTijk, |
I dTtracerLev, rTrans, localTijk, |
238 |
U b5d, c5d, d5d, |
U b5d, c5d, d5d, |
239 |
I myThid) |
I myThid ) |
240 |
ELSEIF (advectionScheme.EQ.ENUM_UPWIND_3RD .OR. |
ELSEIF ( advectionScheme.EQ.ENUM_UPWIND_3RD |
241 |
& advectionScheme.EQ.ENUM_CENTERED_4TH) THEN |
& .OR. advectionScheme.EQ.ENUM_CENTERED_4TH |
242 |
|
& .OR. advectionScheme.EQ.ENUM_DST3 ) THEN |
243 |
diagonalNumber = 5 |
diagonalNumber = 5 |
244 |
CALL GAD_U3C4_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
CALL GAD_U3C4_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
245 |
I advectionScheme, dTtracerLev, rTrans, |
I advectionScheme, dTtracerLev, rTrans, |
246 |
U a5d, b5d, c5d, d5d, e5d, |
U a5d, b5d, c5d, d5d, e5d, |
247 |
I myThid) |
I myThid ) |
248 |
|
ELSEIF ( advectionScheme.EQ.ENUM_DST3_FLUX_LIMIT ) THEN |
249 |
|
diagonalNumber = 5 |
250 |
|
CALL GAD_DST3FL_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax, |
251 |
|
I dTtracerLev, rTrans, localTijk, |
252 |
|
U a5d, b5d, c5d, d5d, e5d, |
253 |
|
I myThid ) |
254 |
ELSE |
ELSE |
255 |
STOP 'GAD_IMPLICIT_R: Adv.Scheme in Impl form not yet coded' |
STOP 'GAD_IMPLICIT_R: Adv.Scheme in Impl form not yet coded' |
256 |
ENDIF |
ENDIF |