diff options
Diffstat (limited to 'src/macro/DA_guts.h')
-rw-r--r-- | src/macro/DA_guts.h | 18 |
1 files changed, 11 insertions, 7 deletions
diff --git a/src/macro/DA_guts.h b/src/macro/DA_guts.h index 2b3ea63..cceb5c8 100644 --- a/src/macro/DA_guts.h +++ b/src/macro/DA_guts.h @@ -12,13 +12,17 @@ #ifdef FCODE - DA_OO2DX = 1D0/(2D0*DA_DX) - DA_OO2DY = 1D0/(2D0*DA_DY) - DA_OO2DZ = 1D0/(2D0*DA_DZ) - - DA_DXDA = DA_OO2DX*(DA_A_IP - DA_A_IM) - DA_DYDA = DA_OO2DY*(DA_A_JP - DA_A_JM) - DA_DZDA = DA_OO2DZ*(DA_A_KP - DA_A_KM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DA_DXDA = ADM_DX_2(alp,i,j,k) + DA_DYDA = ADM_DY_2(alp,i,j,k) + DA_DZDA = ADM_DZ_2(alp,i,j,k) + else + DA_DXDA = ADM_DX_4(alp,i,j,k) + DA_DYDA = ADM_DY_4(alp,i,j,k) + DA_DZDA = ADM_DZ_4(alp,i,j,k) + end if #endif |