diff options
Diffstat (limited to 'src/macro/DXDB_guts.h')
| -rw-r--r-- | src/macro/DXDB_guts.h | 16 |
1 files changed, 11 insertions, 5 deletions
diff --git a/src/macro/DXDB_guts.h b/src/macro/DXDB_guts.h index c47e81d..45a773c 100644 --- a/src/macro/DXDB_guts.h +++ b/src/macro/DXDB_guts.h @@ -17,11 +17,17 @@ #ifdef FCODE - DXDB_OO2DX = 1.0D0/(2.0D0*DXDB_DX) - - DXDB_DXDBX = DXDB_OO2DX*(DXDB_BX_IP - DXDB_BX_IM) - DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM) - DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DXDB_DXDBX = ADM_DX_2(betax,i,j,k) + DXDB_DXDBY = ADM_DX_2(betay,i,j,k) + DXDB_DXDBZ = ADM_DX_2(betaz,i,j,k) + else + DXDB_DXDBX = ADM_DX_4(betax,i,j,k) + DXDB_DXDBY = ADM_DX_4(betay,i,j,k) + DXDB_DXDBZ = ADM_DX_4(betaz,i,j,k) + end if #endif |