93 |
# endif |
# endif |
94 |
# ifdef ALLOW_PTRACERS |
# ifdef ALLOW_PTRACERS |
95 |
# include "PTRACERS_SIZE.h" |
# include "PTRACERS_SIZE.h" |
96 |
# include "PTRACERS.h" |
# include "PTRACERS_FIELDS.h" |
97 |
# endif |
# endif |
98 |
# ifdef ALLOW_OBCS |
# ifdef ALLOW_OBCS |
99 |
# include "OBCS.h" |
# include "OBCS.h" |
317 |
guDissip(i,j) = 0. _d 0 |
guDissip(i,j) = 0. _d 0 |
318 |
gvDissip(i,j) = 0. _d 0 |
gvDissip(i,j) = 0. _d 0 |
319 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
|
cph( |
|
|
c-- need some re-initialisation here to break dependencies |
|
|
cph) |
|
320 |
# ifdef NONLIN_FRSURF |
# ifdef NONLIN_FRSURF |
321 |
# ifndef DISABLE_RSTAR_CODE |
# ifndef DISABLE_RSTAR_CODE |
322 |
dWtransC(i,j,bi,bj) = 0. _d 0 |
dWtransC(i,j,bi,bj) = 0. _d 0 |
323 |
dWtransU(i,j,bi,bj) = 0. _d 0 |
dWtransU(i,j,bi,bj) = 0. _d 0 |
324 |
dWtransV(i,j,bi,bj) = 0. _d 0 |
dWtransV(i,j,bi,bj) = 0. _d 0 |
325 |
# endif |
# endif |
326 |
# endif /* NONLIN_FRSURF */ |
# endif |
327 |
#endif /* ALLOW_AUTODIFF_TAMC */ |
#endif |
328 |
ENDDO |
ENDDO |
329 |
ENDDO |
ENDDO |
330 |
|
|
422 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
423 |
CADJ STORE gvnm1(:,:,k,bi,bj) |
CADJ STORE gvnm1(:,:,k,bi,bj) |
424 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
|
# ifndef DISABLE_RSTAR_CODE |
|
|
CADJ STORE dwtransc(:,:,bi,bj) |
|
|
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
|
|
CADJ STORE dwtransu(:,:,bi,bj) |
|
|
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
|
|
CADJ STORE dwtransv(:,:,bi,bj) |
|
|
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
|
|
# endif |
|
425 |
# ifdef ALLOW_CD_CODE |
# ifdef ALLOW_CD_CODE |
426 |
CADJ STORE unm1(:,:,k,bi,bj) |
CADJ STORE unm1(:,:,k,bi,bj) |
427 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
462 |
C-- Calculate accelerations in the momentum equations (gU, gV, ...) |
C-- Calculate accelerations in the momentum equations (gU, gV, ...) |
463 |
C and step forward storing the result in gU, gV, etc... |
C and step forward storing the result in gU, gV, etc... |
464 |
IF ( momStepping ) THEN |
IF ( momStepping ) THEN |
465 |
IF (.NOT. vectorInvariantMomentum) THEN |
#ifdef ALLOW_AUTODIFF_TAMC |
466 |
#ifdef ALLOW_MOM_FLUXFORM |
# ifdef NONLIN_FRSURF |
467 |
C |
# ifndef DISABLE_RSTAR_CODE |
468 |
# ifdef ALLOW_AUTODIFF_TAMC |
CADJ STORE dWtransC(:,:,bi,bj) |
|
# ifdef NONLIN_FRSURF |
|
|
# ifndef DISABLE_RSTAR_CODE |
|
|
CADJ STORE dwtransc(:,:,bi,bj) |
|
469 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
470 |
CADJ STORE dwtransu(:,:,bi,bj) |
CADJ STORE dWtransU(:,:,bi,bj) |
471 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
472 |
CADJ STORE dwtransv(:,:,bi,bj) |
CADJ STORE dWtransV(:,:,bi,bj) |
473 |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
CADJ & = comlev1_bibj_k, key=kkey, byte=isbyte |
|
# endif |
|
474 |
# endif |
# endif |
475 |
# endif /* ALLOW_AUTODIFF_TAMC */ |
# endif |
476 |
|
#endif |
477 |
|
IF (.NOT. vectorInvariantMomentum) THEN |
478 |
|
#ifdef ALLOW_MOM_FLUXFORM |
479 |
C |
C |
480 |
CALL MOM_FLUXFORM( |
CALL MOM_FLUXFORM( |
481 |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |