diff options
Diffstat (limited to 'src/macro/DZDCG_guts.h')
-rw-r--r-- | src/macro/DZDCG_guts.h | 25 |
1 files changed, 16 insertions, 9 deletions
diff --git a/src/macro/DZDCG_guts.h b/src/macro/DZDCG_guts.h index 82f3a09..08664d0 100644 --- a/src/macro/DZDCG_guts.h +++ b/src/macro/DZDCG_guts.h @@ -21,15 +21,22 @@ #ifdef FCODE - DZDCG_OO2DZ = 1D0/(2D0*DZDCG_DZ) - - DZDCG_DZDCGXX = DZDCG_OO2DZ*(DZDCG_GXX_KP - DZDCG_GXX_KM) - DZDCG_DZDCGXY = DZDCG_OO2DZ*(DZDCG_GXY_KP - DZDCG_GXY_KM) - DZDCG_DZDCGXZ = DZDCG_OO2DZ*(DZDCG_GXZ_KP - DZDCG_GXZ_KM) - DZDCG_DZDCGYY = DZDCG_OO2DZ*(DZDCG_GYY_KP - DZDCG_GYY_KM) - DZDCG_DZDCGYZ = DZDCG_OO2DZ*(DZDCG_GYZ_KP - DZDCG_GYZ_KM) - DZDCG_DZDCGZZ = DZDCG_OO2DZ*(DZDCG_GZZ_KP - DZDCG_GZZ_KM) - +#include "ADM_Derivative.h" + if (spatial_order.eq.2) then + DZDCG_DZDCGXX = ADM_DZ_2(gxx,i,j,k) + DZDCG_DZDCGXY = ADM_DZ_2(gxy,i,j,k) + DZDCG_DZDCGXZ = ADM_DZ_2(gxz,i,j,k) + DZDCG_DZDCGYY = ADM_DZ_2(gyy,i,j,k) + DZDCG_DZDCGYZ = ADM_DZ_2(gyz,i,j,k) + DZDCG_DZDCGZZ = ADM_DZ_2(gzz,i,j,k) + else + DZDCG_DZDCGXX = ADM_DZ_4(gxx,i,j,k) + DZDCG_DZDCGXY = ADM_DZ_4(gxy,i,j,k) + DZDCG_DZDCGXZ = ADM_DZ_4(gxz,i,j,k) + DZDCG_DZDCGYY = ADM_DZ_4(gyy,i,j,k) + DZDCG_DZDCGYZ = ADM_DZ_4(gyz,i,j,k) + DZDCG_DZDCGZZ = ADM_DZ_4(gzz,i,j,k) + end if #endif #ifdef CCODE |