diff options
Diffstat (limited to 'src/macro/LIEK_guts.h')
| -rw-r--r-- | src/macro/LIEK_guts.h | 290 |
1 files changed, 74 insertions, 216 deletions
diff --git a/src/macro/LIEK_guts.h b/src/macro/LIEK_guts.h index b18b3f9..bdb6141 100644 --- a/src/macro/LIEK_guts.h +++ b/src/macro/LIEK_guts.h @@ -22,6 +22,8 @@ #ifdef FCODE +#include "ADM_Derivative.h" + /* Advection. */ LIEK_LKXX = 0.0D0 @@ -34,96 +36,47 @@ if (admmacros_advectionx .eq. 0) then /* CENTER X */ LIEK_LKXX = LIEK_LKXX + DXDK_DXDKXX*LIEK_BX - LIEK_LKYY = LIEK_LKYY + DXDK_DXDKYY*LIEK_BX - LIEK_LKZZ = LIEK_LKZZ + DXDK_DXDKZZ*LIEK_BX - LIEK_LKXY = LIEK_LKXY + DXDK_DXDKXY*LIEK_BX - LIEK_LKXZ = LIEK_LKXZ + DXDK_DXDKXZ*LIEK_BX - LIEK_LKYZ = LIEK_LKYZ + DXDK_DXDKYZ*LIEK_BX else if (admmacros_advectionx .eq. 1) then /* UPWIND1 X */ - LIEK_LKXX = LIEK_LKXX + LIEK_BX \ - *(kxx(i+1,j,k) - kxx(i,j,k))/dx - - LIEK_LKYY = LIEK_LKYY + LIEK_BX \ - *(kyy(i+1,j,k) - kyy(i,j,k))/dx - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BX \ - *(kzz(i+1,j,k) - kzz(i,j,k))/dx - - LIEK_LKXY = LIEK_LKXY + LIEK_BX \ - *(kxy(i+1,j,k) - kxy(i,j,k))/dx - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BX \ - *(kxz(i+1,j,k) - kxz(i,j,k))/dx - - LIEK_LKYZ = LIEK_LKYZ + LIEK_BX \ - *(kyz(i+1,j,k) - kyz(i,j,k))/dx + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DX_P1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DX_P1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DX_P1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DX_P1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DX_P1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DX_P1(kyz,i,j,k) else if (admmacros_advectionx .eq. -1) then /* UPWIND1 X */ - LIEK_LKXX = LIEK_LKXX + LIEK_BX \ - *(kxx(i,j,k) - kxx(i-1,j,k))/dx - - LIEK_LKYY = LIEK_LKYY + LIEK_BX \ - *(kyy(i,j,k) - kyy(i-1,j,k))/dx - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BX \ - *(kzz(i,j,k) - kzz(i-1,j,k))/dx - - LIEK_LKXY = LIEK_LKXY + LIEK_BX \ - *(kxy(i,j,k) - kxy(i-1,j,k))/dx - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BX \ - *(kxz(i,j,k) - kxz(i-1,j,k))/dx - - LIEK_LKYZ = LIEK_LKYZ + LIEK_BX \ - *(kyz(i,j,k) - kyz(i-1,j,k))/dx + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DX_M1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DX_M1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DX_M1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DX_M1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DX_M1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DX_M1(kyz,i,j,k) else if (admmacros_advectionx .eq. 2) then /* UPWIND2 X */ - LIEK_LKXX = LIEK_LKXX - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i+1,j,k) + kxx(i+2,j,k)) - - LIEK_LKYY = LIEK_LKYY - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i+1,j,k) + kyy(i+2,j,k)) - - LIEK_LKZZ = LIEK_LKZZ - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i+1,j,k) + kzz(i+2,j,k)) - - LIEK_LKXY = LIEK_LKXY - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i+1,j,k) + kxy(i+2,j,k)) - - LIEK_LKXZ = LIEK_LKXZ - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i+1,j,k) + kxz(i+2,j,k)) - - LIEK_LKYZ = LIEK_LKYZ - 0.5D0*LIEK_BX/dx \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i+1,j,k) + kyz(i+2,j,k)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DX_P2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DX_P2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DX_P2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DX_P2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DX_P2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DX_P2(kyz,i,j,k) else if (admmacros_advectionx .eq. -2) then /* UPWIND2 X */ - LIEK_LKXX = LIEK_LKXX + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i-1,j,k) + kxx(i-2,j,k)) - - LIEK_LKYY = LIEK_LKYY + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i-1,j,k) + kyy(i-2,j,k)) - - LIEK_LKZZ = LIEK_LKZZ + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i-1,j,k) + kzz(i-2,j,k)) - - LIEK_LKXY = LIEK_LKXY + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i-1,j,k) + kxy(i-2,j,k)) - - LIEK_LKXZ = LIEK_LKXZ + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i-1,j,k) + kxz(i-2,j,k)) - - LIEK_LKYZ = LIEK_LKYZ + 0.5D0*LIEK_BX/dx \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i-1,j,k) + kyz(i-2,j,k)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DX_M2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DX_M2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DX_M2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DX_M2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DX_M2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DX_M2(kyz,i,j,k) end if @@ -131,96 +84,47 @@ if (admmacros_advectiony .eq. 0) then /* CENTER Y */ LIEK_LKXX = LIEK_LKXX + DYDK_DYDKXX*LIEK_BY - LIEK_LKYY = LIEK_LKYY + DYDK_DYDKYY*LIEK_BY - LIEK_LKZZ = LIEK_LKZZ + DYDK_DYDKZZ*LIEK_BY - LIEK_LKXY = LIEK_LKXY + DYDK_DYDKXY*LIEK_BY - LIEK_LKXZ = LIEK_LKXZ + DYDK_DYDKXZ*LIEK_BY - LIEK_LKYZ = LIEK_LKYZ + DYDK_DYDKYZ*LIEK_BY else if (admmacros_advectiony .eq. 1) then /* UPWIND1 Y */ - LIEK_LKXX = LIEK_LKXX + LIEK_BY \ - *(kxx(i,j+1,k) - kxx(i,j,k))/dy - - LIEK_LKYY = LIEK_LKYY + LIEK_BY \ - *(kyy(i,j+1,k) - kyy(i,j,k))/dy - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BY \ - *(kzz(i,j+1,k) - kzz(i,j,k))/dy - - LIEK_LKXY = LIEK_LKXY + LIEK_BY \ - *(kxy(i,j+1,k) - kxy(i,j,k))/dy - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BY \ - *(kxz(i,j+1,k) - kxz(i,j,k))/dy - - LIEK_LKYZ = LIEK_LKYZ + LIEK_BY \ - *(kyz(i,j+1,k) - kyz(i,j,k))/dy + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_P1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_P1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_P1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_P1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_P1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_P1(kyz,i,j,k) else if (admmacros_advectiony .eq. -1) then /* UPWIND1 Y */ - LIEK_LKXX = LIEK_LKXX + LIEK_BY \ - *(kxx(i,j,k) - kxx(i,j-1,k))/dy - - LIEK_LKYY = LIEK_LKYY + LIEK_BY \ - *(kyy(i,j,k) - kyy(i,j-1,k))/dy - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BY \ - *(kzz(i,j,k) - kzz(i,j-1,k))/dy - - LIEK_LKXY = LIEK_LKXY + LIEK_BY \ - *(kxy(i,j,k) - kxy(i,j-1,k))/dy - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BY \ - *(kxz(i,j,k) - kxz(i,j-1,k))/dy - - LIEK_LKYZ = LIEK_LKYZ + LIEK_BY \ - *(kyz(i,j,k) - kyz(i,j-1,k))/dy + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_M1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_M1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_M1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_M1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_M1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_M1(kyz,i,j,k) else if (admmacros_advectiony .eq. 2) then /* UPWIND2 Y */ - LIEK_LKXX = LIEK_LKXX - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i,j+1,k) + kxx(i,j+2,k)) - - LIEK_LKYY = LIEK_LKYY - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i,j+1,k) + kyy(i,j+2,k)) - - LIEK_LKZZ = LIEK_LKZZ - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i,j+1,k) + kzz(i,j+2,k)) - - LIEK_LKXY = LIEK_LKXY - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i,j+1,k) + kxy(i,j+2,k)) - - LIEK_LKXZ = LIEK_LKXZ - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i,j+1,k) + kxz(i,j+2,k)) - - LIEK_LKYZ = LIEK_LKYZ - 0.5D0*LIEK_BY/dy \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i,j+1,k) + kyz(i,j+2,k)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_P2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_P2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_P2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_P2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_P2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_P2(kyz,i,j,k) else if (admmacros_advectiony .eq. -2) then /* UPWIND2 Y */ - LIEK_LKXX = LIEK_LKXX + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i,j-1,k) + kxx(i,j-2,k)) - - LIEK_LKYY = LIEK_LKYY + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i,j-1,k) + kyy(i,j-2,k)) - - LIEK_LKZZ = LIEK_LKZZ + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i,j-1,k) + kzz(i,j-2,k)) - - LIEK_LKXY = LIEK_LKXY + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i,j-1,k) + kxy(i,j-2,k)) - - LIEK_LKXZ = LIEK_LKXZ + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i,j-1,k) + kxz(i,j-2,k)) - - LIEK_LKYZ = LIEK_LKYZ + 0.5D0*LIEK_BY/dy \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i,j-1,k) + kyz(i,j-2,k)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_M2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_M2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_M2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_M2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_M2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_M2(kyz,i,j,k) end if @@ -228,96 +132,50 @@ if (admmacros_advectionz .eq. 0) then /* CENTER Z */ LIEK_LKXX = LIEK_LKXX + DZDK_DZDKXX*LIEK_BZ - LIEK_LKYY = LIEK_LKYY + DZDK_DZDKYY*LIEK_BZ - LIEK_LKZZ = LIEK_LKZZ + DZDK_DZDKZZ*LIEK_BZ - LIEK_LKXY = LIEK_LKXY + DZDK_DZDKXY*LIEK_BZ - LIEK_LKXZ = LIEK_LKXZ + DZDK_DZDKXZ*LIEK_BZ - LIEK_LKYZ = LIEK_LKYZ + DZDK_DZDKYZ*LIEK_BZ else if (admmacros_advectionz .eq. 1) then /* UPWIND1 Z */ - LIEK_LKXX = LIEK_LKXX + LIEK_BZ \ - *(kxx(i,j,k+1) - kxx(i,j,k))/dz - - LIEK_LKYY = LIEK_LKYY + LIEK_BZ \ - *(kyy(i,j,k+1) - kyy(i,j,k))/dz - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BZ \ - *(kzz(i,j,k+1) - kzz(i,j,k))/dz - - LIEK_LKXY = LIEK_LKXY + LIEK_BZ \ - *(kxy(i,j,k+1) - kxy(i,j,k))/dz - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BZ \ - *(kxz(i,j,k+1) - kxz(i,j,k))/dz - - LIEK_LKYZ = LIEK_LKYZ + LIEK_BZ \ - *(kyz(i,j,k+1) - kyz(i,j,k))/dz + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_P1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_P1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_P1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_P1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_P1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_P1(kyz,i,j,k) else if (admmacros_advectionz .eq. -1) then /* UPWIND1 Z */ - LIEK_LKXX = LIEK_LKXX + LIEK_BZ \ - *(kxx(i,j,k) - kxx(i,j,k-1))/dz - - LIEK_LKYY = LIEK_LKYY + LIEK_BZ \ - *(kyy(i,j,k) - kyy(i,j,k-1))/dz - - LIEK_LKZZ = LIEK_LKZZ + LIEK_BZ \ - *(kzz(i,j,k) - kzz(i,j,k-1))/dz - - LIEK_LKXY = LIEK_LKXY + LIEK_BZ \ - *(kxy(i,j,k) - kxy(i,j,k-1))/dz - - LIEK_LKXZ = LIEK_LKXZ + LIEK_BZ \ - *(kxz(i,j,k) - kxz(i,j,k-1))/dz + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_M1(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_M1(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_M1(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_M1(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_M1(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_M1(kyz,i,j,k) else if (admmacros_advectionz .eq. 2) then /* UPWIND2 Z */ - LIEK_LKXX = LIEK_LKXX - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i,j,k+1) + kxx(i,j,k+2)) - - LIEK_LKYY = LIEK_LKYY - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i,j,k+1) + kyy(i,j,k+2)) - - LIEK_LKZZ = LIEK_LKZZ - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i,j,k+1) + kzz(i,j,k+2)) - - LIEK_LKXY = LIEK_LKXY - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i,j,k+1) + kxy(i,j,k+2)) - - LIEK_LKXZ = LIEK_LKXZ - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i,j,k+1) + kxz(i,j,k+2)) - - LIEK_LKYZ = LIEK_LKYZ - 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i,j,k+1) + kyz(i,j,k+2)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_P2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_P2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_P2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_P2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_P2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_P2(kyz,i,j,k) else if (admmacros_advectionz .eq. -2) then /* UPWIND2 Z */ - LIEK_LKXX = LIEK_LKXX + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxx(i,j,k) - 4.0D0*kxx(i,j,k-1) + kxx(i,j,k-2)) - - LIEK_LKYY = LIEK_LKYY + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kyy(i,j,k) - 4.0D0*kyy(i,j,k-1) + kyy(i,j,k-2)) + LIEK_LKXX = LIEK_LKXX + ADM_ADV_DY_M2(kxx,i,j,k) + LIEK_LKYY = LIEK_LKYY + ADM_ADV_DY_M2(kyy,i,j,k) + LIEK_LKZZ = LIEK_LKZZ + ADM_ADV_DY_M2(kzz,i,j,k) + LIEK_LKXY = LIEK_LKXY + ADM_ADV_DY_M2(kxy,i,j,k) + LIEK_LKXZ = LIEK_LKXZ + ADM_ADV_DY_M2(kxz,i,j,k) + LIEK_LKYZ = LIEK_LKYZ + ADM_ADV_DY_M2(kyz,i,j,k) - LIEK_LKZZ = LIEK_LKZZ + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kzz(i,j,k) - 4.0D0*kzz(i,j,k-1) + kzz(i,j,k-2)) - - LIEK_LKXY = LIEK_LKXY + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxy(i,j,k) - 4.0D0*kxy(i,j,k-1) + kxy(i,j,k-2)) - - LIEK_LKXZ = LIEK_LKXZ + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kxz(i,j,k) - 4.0D0*kxz(i,j,k-1) + kxz(i,j,k-2)) - - LIEK_LKYZ = LIEK_LKYZ + 0.5D0*LIEK_BZ/dz \ - *(3.0D0*kyz(i,j,k) - 4.0D0*kyz(i,j,k-1) + kyz(i,j,k-2)) end if - /* Extra terms in the Lie derivative. */ LIEK_LKXX = LIEK_LKXX + 2.0D0*(DXDB_DXDBX*LIEK_KXX \ |