diff options
Diffstat (limited to 'src/macro/DXDB_guts.h')
-rw-r--r-- | src/macro/DXDB_guts.h | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/src/macro/DXDB_guts.h b/src/macro/DXDB_guts.h index 3b1531b..295da7b 100644 --- a/src/macro/DXDB_guts.h +++ b/src/macro/DXDB_guts.h @@ -43,50 +43,50 @@ end if - if (admmacros_advectiony .eq. 0) then /* CENTER */ + if (admmacros_advectionx .eq. 0) then /* CENTER */ DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM) - else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */ + else if (admmacros_advectionx .eq. 1) then /* UPWIND1 */ DXDB_DXDBY = (DXDB_BY_IP - DXDB_BY)/DXDB_DX - else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */ + else if (admmacros_advectionx .eq. -1) then /* UPWIND1 */ DXDB_DXDBY = (DXDB_BY - DXDB_BY_IM)/DXDB_DX - else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */ + else if (admmacros_advectionx .eq. 2) then /* UPWIND2 */ DXDB_DXDBY = - 0.5D0/DXDB_DX \ *(3.0D0*DXDB_BY - 4.0D0*DXDB_BY_IP + DXDB_BY_IPP) - else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */ + else if (admmacros_advectionx .eq. -2) then /* UPWIND2 */ DXDB_DXDBY = 0.5D0/DXDB_DX \ *(3.0D0*DXDB_BY - 4.0D0*DXDB_BY_IM + DXDB_BY_IMM) end if - if (admmacros_advectionz .eq. 0) then /* CENTER */ + if (admmacros_advectionx .eq. 0) then /* CENTER */ DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM) - else if (admmacros_advectionz .eq. 1) then /* UPWIND1 */ + else if (admmacros_advectionx .eq. 1) then /* UPWIND1 */ DXDB_DXDBZ = (DXDB_BZ_IP - DXDB_BZ)/DXDB_DX - else if (admmacros_advectionz .eq. -1) then /* UPWIND1 */ + else if (admmacros_advectionx .eq. -1) then /* UPWIND1 */ DXDB_DXDBZ = (DXDB_BZ - DXDB_BZ_IM)/DXDB_DX - else if (admmacros_advectionz .eq. 2) then /* UPWIND2 */ + else if (admmacros_advectionx .eq. 2) then /* UPWIND2 */ DXDB_DXDBZ = - 0.5D0/DXDB_DX \ *(3.0D0*DXDB_BZ - 4.0D0*DXDB_BZ_IP + DXDB_BZ_IPP) - else if (admmacros_advectionz .eq. -2) then /* UPWIND2 */ + else if (admmacros_advectionx .eq. -2) then /* UPWIND2 */ - DXDB_DXDBY = 0.5D0/DXDB_DX \ + DXDB_DXDBZ = 0.5D0/DXDB_DX \ *(3.0D0*DXDB_BZ - 4.0D0*DXDB_BZ_IM + DXDB_BZ_IMM) end if |