aboutsummaryrefslogtreecommitdiff
path: root/src/macro/DZDB_guts.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/macro/DZDB_guts.h')
-rw-r--r--src/macro/DZDB_guts.h76
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