diff options
Diffstat (limited to 'src/macro/DYDCG_guts.h')
-rw-r--r-- | src/macro/DYDCG_guts.h | 26 |
1 files changed, 17 insertions, 9 deletions
diff --git a/src/macro/DYDCG_guts.h b/src/macro/DYDCG_guts.h index 01bced2..3c328d7 100644 --- a/src/macro/DYDCG_guts.h +++ b/src/macro/DYDCG_guts.h @@ -21,15 +21,23 @@ #ifdef FCODE - DYDCG_OO2DY = 1D0/(2D0*DYDCG_DY) - - DYDCG_DYDCGXX = DYDCG_OO2DY*(DYDCG_GXX_JP - DYDCG_GXX_JM) - DYDCG_DYDCGXY = DYDCG_OO2DY*(DYDCG_GXY_JP - DYDCG_GXY_JM) - DYDCG_DYDCGXZ = DYDCG_OO2DY*(DYDCG_GXZ_JP - DYDCG_GXZ_JM) - DYDCG_DYDCGYY = DYDCG_OO2DY*(DYDCG_GYY_JP - DYDCG_GYY_JM) - DYDCG_DYDCGYZ = DYDCG_OO2DY*(DYDCG_GYZ_JP - DYDCG_GYZ_JM) - DYDCG_DYDCGZZ = DYDCG_OO2DY*(DYDCG_GZZ_JP - DYDCG_GZZ_JM) - +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DYDCG_DYDCGXX = ADM_DY_2(gxx,i,j,k) + DYDCG_DYDCGXY = ADM_DY_2(gxy,i,j,k) + DYDCG_DYDCGXZ = ADM_DY_2(gxz,i,j,k) + DYDCG_DYDCGYY = ADM_DY_2(gyy,i,j,k) + DYDCG_DYDCGYZ = ADM_DY_2(gyz,i,j,k) + DYDCG_DYDCGZZ = ADM_DY_2(gzz,i,j,k) + else + DYDCG_DYDCGXX = ADM_DY_4(gxx,i,j,k) + DYDCG_DYDCGXY = ADM_DY_4(gxy,i,j,k) + DYDCG_DYDCGXZ = ADM_DY_4(gxz,i,j,k) + DYDCG_DYDCGYY = ADM_DY_4(gyy,i,j,k) + DYDCG_DYDCGYZ = ADM_DY_4(gyz,i,j,k) + DYDCG_DYDCGZZ = ADM_DY_4(gzz,i,j,k) + end if #endif |