diff options
Diffstat (limited to 'src/macro/DXDCG_guts.h')
| -rw-r--r-- | src/macro/DXDCG_guts.h | 26 |
1 files changed, 17 insertions, 9 deletions
diff --git a/src/macro/DXDCG_guts.h b/src/macro/DXDCG_guts.h index 464fdb7..b922068 100644 --- a/src/macro/DXDCG_guts.h +++ b/src/macro/DXDCG_guts.h @@ -21,15 +21,23 @@ #ifdef FCODE - DXDCG_OO2DX = 1D0/(2D0*DXDCG_DX) - - DXDCG_DXDCGXX = DXDCG_OO2DX*(DXDCG_GXX_IP - DXDCG_GXX_IM) - DXDCG_DXDCGXY = DXDCG_OO2DX*(DXDCG_GXY_IP - DXDCG_GXY_IM) - DXDCG_DXDCGXZ = DXDCG_OO2DX*(DXDCG_GXZ_IP - DXDCG_GXZ_IM) - DXDCG_DXDCGYY = DXDCG_OO2DX*(DXDCG_GYY_IP - DXDCG_GYY_IM) - DXDCG_DXDCGYZ = DXDCG_OO2DX*(DXDCG_GYZ_IP - DXDCG_GYZ_IM) - DXDCG_DXDCGZZ = DXDCG_OO2DX*(DXDCG_GZZ_IP - DXDCG_GZZ_IM) - +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DXDCG_DXDCGXX = ADM_DX_2(gxx,i,j,k) + DXDCG_DXDCGXY = ADM_DX_2(gxy,i,j,k) + DXDCG_DXDCGXZ = ADM_DX_2(gxz,i,j,k) + DXDCG_DXDCGYY = ADM_DX_2(gyy,i,j,k) + DXDCG_DXDCGYZ = ADM_DX_2(gyz,i,j,k) + DXDCG_DXDCGZZ = ADM_DX_2(gzz,i,j,k) + else + DXDCG_DXDCGXX = ADM_DX_4(gxx,i,j,k) + DXDCG_DXDCGXY = ADM_DX_4(gxy,i,j,k) + DXDCG_DXDCGXZ = ADM_DX_4(gxz,i,j,k) + DXDCG_DXDCGYY = ADM_DX_4(gyy,i,j,k) + DXDCG_DXDCGYZ = ADM_DX_4(gyz,i,j,k) + DXDCG_DXDCGZZ = ADM_DX_4(gzz,i,j,k) + end if #endif #ifdef CCODE |