diff options
Diffstat (limited to 'src/macro/DYDB_guts.h')
-rw-r--r-- | src/macro/DYDB_guts.h | 76 |
1 files changed, 4 insertions, 72 deletions
diff --git a/src/macro/DYDB_guts.h b/src/macro/DYDB_guts.h index 43f5b91..54cba6c 100644 --- a/src/macro/DYDB_guts.h +++ b/src/macro/DYDB_guts.h @@ -17,79 +17,11 @@ #ifdef FCODE - DYDB_OO2DY = 1D0/(2D0*DYDB_DY) - - if (admmacros_advectiony .eq. 0) then /* CENTER */ - - DYDB_DYDBX = DYDB_OO2DY*(DYDB_BX_JP - DYDB_BX_JM) - - else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */ - - DYDB_DYDBX = (DYDB_BX_JP - DYDB_BX)/DYDB_DY - - else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */ - - DYDB_DYDBX = (DYDB_BX - DYDB_BX_JM)/DYDB_DY - - else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */ - - DYDB_DYDBX = - 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BX - 4.0D0*DYDB_BX_JP + DYDB_BX_JPP) - - else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */ - - DYDB_DYDBX = 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BX - 4.0D0*DYDB_BX_JM + DYDB_BX_JMM) - - end if - - if (admmacros_advectiony .eq. 0) then /* CENTER */ - - DYDB_DYDBY = DYDB_OO2DY*(DYDB_BY_JP - DYDB_BY_JM) - - else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */ - - DYDB_DYDBY = (DYDB_BY_JP - DYDB_BY)/DYDB_DY - - else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */ - - DYDB_DYDBY = (DYDB_BY - DYDB_BY_JM)/DYDB_DY - - else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */ - - DYDB_DYDBY = - 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BY - 4.0D0*DYDB_BY_JP + DYDB_BY_JPP) - - else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */ - - DYDB_DYDBY = 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BY - 4.0D0*DYDB_BY_JM + DYDB_BY_JMM) - - end if - - if (admmacros_advectiony .eq. 0) then /* CENTER */ - - DYDB_DYDBZ = DYDB_OO2DY*(DYDB_BZ_JP - DYDB_BZ_JM) - - else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */ - - DYDB_DYDBZ = (DYDB_BZ_JP - DYDB_BZ)/DYDB_DY - - else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */ - - DYDB_DYDBZ = (DYDB_BZ - DYDB_BZ_JM)/DYDB_DY - - else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */ - - DYDB_DYDBZ = - 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BZ - 4.0D0*DYDB_BZ_JP + DYDB_BZ_JPP) - - else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */ - - DYDB_DYDBZ = 0.5D0/DYDB_DY \ - *(3.0D0*DYDB_BZ - 4.0D0*DYDB_BZ_JM + DYDB_BZ_JMM) + DYDB_OO2DY = 1.0D0/(2.0D0*DYDB_DY) - end if + 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) #endif |