diff options
Diffstat (limited to 'src/macro/DZDK_guts.h')
-rw-r--r-- | src/macro/DZDK_guts.h | 26 |
1 files changed, 17 insertions, 9 deletions
diff --git a/src/macro/DZDK_guts.h b/src/macro/DZDK_guts.h index 1e205ab..cc8e5e6 100644 --- a/src/macro/DZDK_guts.h +++ b/src/macro/DZDK_guts.h @@ -13,15 +13,23 @@ #ifdef FCODE - DZDK_OO2DZ = 1D0/(2D0*DZDK_DZ) - - DZDK_DZDKXX = DZDK_OO2DZ*(DZDK_KXX_KP - DZDK_KXX_KM) - DZDK_DZDKXY = DZDK_OO2DZ*(DZDK_KXY_KP - DZDK_KXY_KM) - DZDK_DZDKXZ = DZDK_OO2DZ*(DZDK_KXZ_KP - DZDK_KXZ_KM) - DZDK_DZDKYY = DZDK_OO2DZ*(DZDK_KYY_KP - DZDK_KYY_KM) - DZDK_DZDKYZ = DZDK_OO2DZ*(DZDK_KYZ_KP - DZDK_KYZ_KM) - DZDK_DZDKZZ = DZDK_OO2DZ*(DZDK_KZZ_KP - DZDK_KZZ_KM) - +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DZDK_DZDKXX = ADM_DZ_2(kxx,i,j,k) + DZDK_DZDKXY = ADM_DZ_2(kxy,i,j,k) + DZDK_DZDKXZ = ADM_DZ_2(kxz,i,j,k) + DZDK_DZDKYY = ADM_DZ_2(kyy,i,j,k) + DZDK_DZDKYZ = ADM_DZ_2(kyz,i,j,k) + DZDK_DZDKZZ = ADM_DZ_2(kzz,i,j,k) + else + DZDK_DZDKXX = ADM_DZ_4(kxx,i,j,k) + DZDK_DZDKXY = ADM_DZ_4(kxy,i,j,k) + DZDK_DZDKXZ = ADM_DZ_4(kxz,i,j,k) + DZDK_DZDKYY = ADM_DZ_4(kyy,i,j,k) + DZDK_DZDKYZ = ADM_DZ_4(kyz,i,j,k) + DZDK_DZDKZZ = ADM_DZ_4(kzz,i,j,k) + end if #endif #ifdef CCODE |