diff options
Diffstat (limited to 'src/macro/DZDB_guts.h')
-rw-r--r-- | src/macro/DZDB_guts.h | 16 |
1 files changed, 10 insertions, 6 deletions
diff --git a/src/macro/DZDB_guts.h b/src/macro/DZDB_guts.h index 5797db7..4b01986 100644 --- a/src/macro/DZDB_guts.h +++ b/src/macro/DZDB_guts.h @@ -17,12 +17,16 @@ #ifdef FCODE - DZDB_OO2DZ = 1.0D0/(2.0D0*DZDB_DZ) - - DZDB_DZDBX = DZDB_OO2DZ*(DZDB_BX_KP - DZDB_BX_KM) - DZDB_DZDBY = DZDB_OO2DZ*(DZDB_BY_KP - DZDB_BY_KM) - DZDB_DZDBZ = DZDB_OO2DZ*(DZDB_BZ_KP - DZDB_BZ_KM) - +#include "ADM_Derivative.h" + if (spatial_order.eq.2) then + DZDB_DZDBX = ADM_DZ_2(betax,i,j,k) + DZDB_DZDBY = ADM_DZ_2(betay,i,j,k) + DZDB_DZDBZ = ADM_DZ_2(betaz,i,j,k) + else + DZDB_DZDBX = ADM_DZ_4(betax,i,j,k) + DZDB_DZDBY = ADM_DZ_4(betay,i,j,k) + DZDB_DZDBZ = ADM_DZ_4(betaz,i,j,k) + end if #endif #ifdef CCODE |