diff options
Diffstat (limited to 'src/macro/DYDK_guts.h')
| -rw-r--r-- | src/macro/DYDK_guts.h | 26 |
1 files changed, 17 insertions, 9 deletions
diff --git a/src/macro/DYDK_guts.h b/src/macro/DYDK_guts.h index 26364b3..7dbfc0c 100644 --- a/src/macro/DYDK_guts.h +++ b/src/macro/DYDK_guts.h @@ -13,15 +13,23 @@ #ifdef FCODE - DYDK_OO2DY = 1D0/(2D0*DYDK_DY) - - DYDK_DYDKXX = DYDK_OO2DY*(DYDK_KXX_JP - DYDK_KXX_JM) - DYDK_DYDKXY = DYDK_OO2DY*(DYDK_KXY_JP - DYDK_KXY_JM) - DYDK_DYDKXZ = DYDK_OO2DY*(DYDK_KXZ_JP - DYDK_KXZ_JM) - DYDK_DYDKYY = DYDK_OO2DY*(DYDK_KYY_JP - DYDK_KYY_JM) - DYDK_DYDKYZ = DYDK_OO2DY*(DYDK_KYZ_JP - DYDK_KYZ_JM) - DYDK_DYDKZZ = DYDK_OO2DY*(DYDK_KZZ_JP - DYDK_KZZ_JM) - +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DYDK_DYDKXX = ADM_DY_2(kxx,i,j,k) + DYDK_DYDKXY = ADM_DY_2(kxy,i,j,k) + DYDK_DYDKXZ = ADM_DY_2(kxz,i,j,k) + DYDK_DYDKYY = ADM_DY_2(kyy,i,j,k) + DYDK_DYDKYZ = ADM_DY_2(kyz,i,j,k) + DYDK_DYDKZZ = ADM_DY_2(kzz,i,j,k) + else + DYDK_DYDKXX = ADM_DY_4(kxx,i,j,k) + DYDK_DYDKXY = ADM_DY_4(kxy,i,j,k) + DYDK_DYDKXZ = ADM_DY_4(kxz,i,j,k) + DYDK_DYDKYY = ADM_DY_4(kyy,i,j,k) + DYDK_DYDKYZ = ADM_DY_4(kyz,i,j,k) + DYDK_DYDKZZ = ADM_DY_4(kzz,i,j,k) + end if #endif |