diff options
Diffstat (limited to 'src/macro/DXDK_guts.h')
| -rw-r--r-- | src/macro/DXDK_guts.h | 25 |
1 files changed, 17 insertions, 8 deletions
diff --git a/src/macro/DXDK_guts.h b/src/macro/DXDK_guts.h index 6a6df3f..fd8af92 100644 --- a/src/macro/DXDK_guts.h +++ b/src/macro/DXDK_guts.h @@ -13,14 +13,23 @@ #ifdef FCODE - DXDK_OO2DX = 1D0/(2D0*DXDK_DX) - - DXDK_DXDKXX = DXDK_OO2DX*(DXDK_KXX_IP - DXDK_KXX_IM) - DXDK_DXDKXY = DXDK_OO2DX*(DXDK_KXY_IP - DXDK_KXY_IM) - DXDK_DXDKXZ = DXDK_OO2DX*(DXDK_KXZ_IP - DXDK_KXZ_IM) - DXDK_DXDKYY = DXDK_OO2DX*(DXDK_KYY_IP - DXDK_KYY_IM) - DXDK_DXDKYZ = DXDK_OO2DX*(DXDK_KYZ_IP - DXDK_KYZ_IM) - DXDK_DXDKZZ = DXDK_OO2DX*(DXDK_KZZ_IP - DXDK_KZZ_IM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DXDK_DXDKXX = ADM_DX_2(kxx,i,j,k) + DXDK_DXDKXY = ADM_DX_2(kxy,i,j,k) + DXDK_DXDKXZ = ADM_DX_2(kxz,i,j,k) + DXDK_DXDKYY = ADM_DX_2(kyy,i,j,k) + DXDK_DXDKYZ = ADM_DX_2(kyz,i,j,k) + DXDK_DXDKZZ = ADM_DX_2(kzz,i,j,k) + else + DXDK_DXDKXX = ADM_DX_4(kxx,i,j,k) + DXDK_DXDKXY = ADM_DX_4(kxy,i,j,k) + DXDK_DXDKXZ = ADM_DX_4(kxz,i,j,k) + DXDK_DXDKYY = ADM_DX_4(kyy,i,j,k) + DXDK_DXDKYZ = ADM_DX_4(kyz,i,j,k) + DXDK_DXDKZZ = ADM_DX_4(kzz,i,j,k) + end if #endif |