diff options
Diffstat (limited to 'src/macro/DZDB_guts.h')
-rw-r--r-- | src/macro/DZDB_guts.h | 76 |
1 files changed, 4 insertions, 72 deletions
diff --git a/src/macro/DZDB_guts.h b/src/macro/DZDB_guts.h index 75b708e..5797db7 100644 --- a/src/macro/DZDB_guts.h +++ b/src/macro/DZDB_guts.h @@ -17,79 +17,11 @@ #ifdef FCODE - DZDB_OO2DZ = 1D0/(2D0*DZDB_DZ) - - if (admmacros_advectionz .eq. 0) then /* CENTER */ - - DZDB_DZDBX = DZDB_OO2DZ*(DZDB_BX_KP - DZDB_BX_KM) - - else if (admmacros_advectionz .eq. 1) then /* UPWIND1 */ - - DZDB_DZDBX = (DZDB_BX_KP - DZDB_BX)/DZDB_DZ - - else if (admmacros_advectionz .eq. -1) then /* UPWIND1 */ - - DZDB_DZDBX = (DZDB_BX - DZDB_BX_KM)/DZDB_DZ - - else if (admmacros_advectionz .eq. 2) then /* UPWIND2 */ - - DZDB_DZDBX = - 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BX - 4.0D0*DZDB_BX_KP + DZDB_BX_KPP) - - else if (admmacros_advectionz .eq. -2) then /* UPWIND2 */ - - DZDB_DZDBX = 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BX - 4.0D0*DZDB_BX_KM + DZDB_BX_KMM) - - end if - - if (admmacros_advectionz .eq. 0) then /* CENTER */ - - DZDB_DZDBY = DZDB_OO2DZ*(DZDB_BY_KP - DZDB_BY_KM) - - else if (admmacros_advectionz .eq. 1) then /* UPWIND1 */ - - DZDB_DZDBY = (DZDB_BY_KP - DZDB_BY)/DZDB_DZ - - else if (admmacros_advectionz .eq. -1) then /* UPWIND1 */ - - DZDB_DZDBY = (DZDB_BY - DZDB_BY_KM)/DZDB_DZ - - else if (admmacros_advectionz .eq. 2) then /* UPWIND2 */ - - DZDB_DZDBY = - 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BY - 4.0D0*DZDB_BY_KP + DZDB_BY_KPP) - - else if (admmacros_advectionz .eq. -2) then /* UPWIND2 */ - - DZDB_DZDBY = 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BY - 4.0D0*DZDB_BY_KM + DZDB_BY_KMM) - - end if - - if (admmacros_advectionz .eq. 0) then /* CENTER */ - - DZDB_DZDBZ = DZDB_OO2DZ*(DZDB_BZ_KP - DZDB_BZ_KM) - - else if (admmacros_advectionz .eq. 1) then /* UPWIND1 */ - - DZDB_DZDBZ = (DZDB_BZ_KP - DZDB_BZ)/DZDB_DZ - - else if (admmacros_advectionz .eq. -1) then /* UPWIND1 */ - - DZDB_DZDBZ = (DZDB_BZ - DZDB_BZ_KM)/DZDB_DZ - - else if (admmacros_advectionz .eq. 2) then /* UPWIND2 */ - - DZDB_DZDBZ = - 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BZ - 4.0D0*DZDB_BZ_KP + DZDB_BZ_KPP) - - else if (admmacros_advectionz .eq. -2) then /* UPWIND2 */ - - DZDB_DZDBZ = 0.5D0/DZDB_DZ \ - *(3.0D0*DZDB_BZ - 4.0D0*DZDB_BZ_KM + DZDB_BZ_KMM) + DZDB_OO2DZ = 1.0D0/(2.0D0*DZDB_DZ) - end if + 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) #endif |