531 |
#else |
#else |
532 |
Kgm_tmp = GM_isopycK |
Kgm_tmp = GM_isopycK |
533 |
#endif |
#endif |
534 |
#ifdef ALLOW_KAPGM_CONTROL |
#if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL)) |
535 |
& + GM_skewflx*kapgm(i,j,k,bi,bj) |
& + GM_skewflx*kapgm(i,j,k,bi,bj) |
536 |
#else |
#else |
537 |
& + GM_skewflx*GM_background_K |
& + GM_skewflx*GM_background_K |
715 |
#else |
#else |
716 |
& ( GM_isopycK |
& ( GM_isopycK |
717 |
#endif |
#endif |
718 |
#ifdef ALLOW_KAPGM_CONTROL |
#if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL)) |
719 |
& - GM_skewflx*kapgm(i,j,k,bi,bj) |
& - GM_skewflx*kapgm(i,j,k,bi,bj) |
720 |
#else |
#else |
721 |
& - GM_skewflx*GM_background_K |
& - GM_skewflx*GM_background_K |
776 |
|
|
777 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
778 |
C-- 3rd k loop : compute Tensor Coeff. at V point |
C-- 3rd k loop : compute Tensor Coeff. at V point |
779 |
|
|
780 |
#ifdef ALLOW_KPP |
#ifdef ALLOW_KPP |
781 |
IF ( useKPP ) THEN |
IF ( useKPP ) THEN |
782 |
DO j=2-Oly,sNy+Oly |
DO j=2-Oly,sNy+Oly |
928 |
#else |
#else |
929 |
& ( GM_isopycK |
& ( GM_isopycK |
930 |
#endif |
#endif |
931 |
#ifdef ALLOW_KAPGM_CONTROL |
#if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL)) |
932 |
& - GM_skewflx*kapgm(i,j,k,bi,bj) |
& - GM_skewflx*kapgm(i,j,k,bi,bj) |
933 |
#else |
#else |
934 |
& - GM_skewflx*GM_background_K |
& - GM_skewflx*GM_background_K |