diff options
Diffstat (limited to 'src/macro/DYDB_guts.h')
| -rw-r--r-- | src/macro/DYDB_guts.h | 16 |
1 files changed, 11 insertions, 5 deletions
diff --git a/src/macro/DYDB_guts.h b/src/macro/DYDB_guts.h index 54cba6c..bd92f57 100644 --- a/src/macro/DYDB_guts.h +++ b/src/macro/DYDB_guts.h @@ -17,11 +17,17 @@ #ifdef FCODE - DYDB_OO2DY = 1.0D0/(2.0D0*DYDB_DY) - - DYDB_DYDBX = DYDB_OO2DY*(DYDB_BX_JP - DYDB_BX_JM) - DYDB_DYDBY = DYDB_OO2DY*(DYDB_BY_JP - DYDB_BY_JM) - DYDB_DYDBZ = DYDB_OO2DY*(DYDB_BZ_JP - DYDB_BZ_JM) +#include "ADM_Derivative.h" + + if (spatial_order.eq.2) then + DYDB_DYDBX = ADM_DY_2(betax,i,j,k) + DYDB_DYDBY = ADM_DY_2(betay,i,j,k) + DYDB_DYDBZ = ADM_DY_2(betaz,i,j,k) + else + DYDB_DYDBX = ADM_DY_4(betax,i,j,k) + DYDB_DYDBY = ADM_DY_4(betay,i,j,k) + DYDB_DYDBZ = ADM_DY_4(betaz,i,j,k) + end if #endif |