193 |
fVerT(i,j,1) = 0. _d 0 |
fVerT(i,j,1) = 0. _d 0 |
194 |
fVerT(i,j,2) = 0. _d 0 |
fVerT(i,j,2) = 0. _d 0 |
195 |
rTransKp1(i,j)= 0. _d 0 |
rTransKp1(i,j)= 0. _d 0 |
196 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
197 |
|
localTij(i,j) = 0. _d 0 |
198 |
|
#endif |
199 |
ENDDO |
ENDDO |
200 |
ENDDO |
ENDDO |
201 |
|
|
317 |
|
|
318 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
319 |
C-- X direction |
C-- X direction |
320 |
|
C- Advective flux in X |
321 |
|
DO j=1-Oly,sNy+Oly |
322 |
|
DO i=1-Olx,sNx+Olx |
323 |
|
af(i,j) = 0. |
324 |
|
ENDDO |
325 |
|
ENDDO |
326 |
|
C |
327 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
328 |
|
cph-test |
329 |
|
CADJ STORE localTij(:,:) = |
330 |
|
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
331 |
|
CADJ STORE af(:,:) = |
332 |
|
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
333 |
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
334 |
|
C |
335 |
IF (calc_fluxes_X) THEN |
IF (calc_fluxes_X) THEN |
336 |
|
|
337 |
C- Do not compute fluxes if |
C- Do not compute fluxes if |
351 |
ENDIF |
ENDIF |
352 |
#endif |
#endif |
353 |
|
|
|
C- Advective flux in X |
|
|
DO j=1-Oly,sNy+Oly |
|
|
DO i=1-Olx,sNx+Olx |
|
|
af(i,j) = 0. |
|
|
ENDDO |
|
|
ENDDO |
|
|
|
|
354 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
355 |
#ifndef DISABLE_MULTIDIM_ADVECTION |
# ifndef DISABLE_MULTIDIM_ADVECTION |
356 |
CADJ STORE localTij(:,:) = |
CADJ STORE localTij(:,:) = |
357 |
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
358 |
#endif |
# endif |
359 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif /* ALLOW_AUTODIFF_TAMC */ |
360 |
|
|
361 |
IF ( advectionScheme.EQ.ENUM_UPWIND_1RST |
IF ( advectionScheme.EQ.ENUM_UPWIND_1RST |
472 |
|
|
473 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
474 |
C-- Y direction |
C-- Y direction |
475 |
|
cph-test |
476 |
|
C- Advective flux in Y |
477 |
|
DO j=1-Oly,sNy+Oly |
478 |
|
DO i=1-Olx,sNx+Olx |
479 |
|
af(i,j) = 0. |
480 |
|
ENDDO |
481 |
|
ENDDO |
482 |
|
C |
483 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
484 |
|
cph-test |
485 |
|
CADJ STORE localTij(:,:) = |
486 |
|
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
487 |
|
CADJ STORE af(:,:) = |
488 |
|
CADJ & comlev1_bibj_k_gad_pass, key=passkey, byte=isbyte |
489 |
|
#endif /* ALLOW_AUTODIFF_TAMC */ |
490 |
|
C |
491 |
IF (calc_fluxes_Y) THEN |
IF (calc_fluxes_Y) THEN |
492 |
|
|
493 |
C- Do not compute fluxes if |
C- Do not compute fluxes if |