diff options
Diffstat (limited to 'src/macro/DDA_guts.h')
-rw-r--r-- | src/macro/DDA_guts.h | 31 |
1 files changed, 17 insertions, 14 deletions
diff --git a/src/macro/DDA_guts.h b/src/macro/DDA_guts.h index de4c849..7f115f0 100644 --- a/src/macro/DDA_guts.h +++ b/src/macro/DDA_guts.h @@ -12,20 +12,23 @@ #ifdef FCODE - DDA_OODX2 = 1D0/(DDA_DX*DDA_DX) - DDA_OODY2 = 1D0/(DDA_DY*DDA_DY) - DDA_OODZ2 = 1D0/(DDA_DZ*DDA_DZ) - DDA_OO4DXDY = 1D0/(4D0*DDA_DX*DDA_DY) - DDA_OO4DXDZ = 1D0/(4D0*DDA_DX*DDA_DZ) - DDA_OO4DYDZ = 1D0/(4D0*DDA_DY*DDA_DZ) - - DDA_DXXDA = DDA_OODX2*(DDA_A_IP - 2D0*DDA_A + DDA_A_IM) - DDA_DYYDA = DDA_OODY2*(DDA_A_JP - 2D0*DDA_A + DDA_A_JM) - DDA_DZZDA = DDA_OODZ2*(DDA_A_KP - 2D0*DDA_A + DDA_A_KM) - - DDA_DXYDA = DDA_OO4DXDY*(DDA_A_IPJP-DDA_A_IPJM-DDA_A_IMJP+DDA_A_IMJM) - DDA_DXZDA = DDA_OO4DXDZ*(DDA_A_IPKP-DDA_A_IPKM-DDA_A_IMKP+DDA_A_IMKM) - DDA_DYZDA = DDA_OO4DYDZ*(DDA_A_JPKP-DDA_A_JPKM-DDA_A_JMKP+DDA_A_JMKM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DDA_DXXDA = ADM_DXX_2(alp,i,j,k) + DDA_DXYDA = ADM_DXY_2(alp,i,j,k) + DDA_DXZDA = ADM_DXZ_2(alp,i,j,k) + DDA_DYYDA = ADM_DYY_2(alp,i,j,k) + DDA_DYZDA = ADM_DYZ_2(alp,i,j,k) + DDA_DZZDA = ADM_DZZ_2(alp,i,j,k) + else + DDA_DXXDA = ADM_DXX_4(alp,i,j,k) + DDA_DXYDA = ADM_DXY_4(alp,i,j,k) + DDA_DXZDA = ADM_DXZ_4(alp,i,j,k) + DDA_DYYDA = ADM_DYY_4(alp,i,j,k) + DDA_DYZDA = ADM_DYZ_4(alp,i,j,k) + DDA_DZZDA = ADM_DZZ_4(alp,i,j,k) + end if #endif |