diff options
Diffstat (limited to 'src/Derivatives_2_1.F90')
| -rw-r--r-- | src/Derivatives_2_1.F90 | 36 |
1 files changed, 6 insertions, 30 deletions
diff --git a/src/Derivatives_2_1.F90 b/src/Derivatives_2_1.F90 index 788d325..4ea8593 100644 --- a/src/Derivatives_2_1.F90 +++ b/src/Derivatives_2_1.F90 @@ -4,13 +4,13 @@ subroutine deriv_gf_2_1 ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) + use All_Coeffs_mod + implicit none DECLARE_CCTK_FUNCTIONS DECLARE_CCTK_PARAMETERS - CCTK_REAL, parameter :: zero = 0.0 - integer, parameter :: wp = kind(zero) CCTK_INT, intent(IN) :: ni, nj, nk CCTK_REAL, dimension(ni,nj,nk), intent(IN) :: var CCTK_INT, intent(IN) :: dir @@ -28,10 +28,7 @@ subroutine deriv_gf_2_1 ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) logical, save :: first = .true. if ( first ) then - a(1) = 1.0_wp/2.0_wp - - q(1,1) = -1.0_wp; q(2,1) = 1.0_wp; q(3,1) = zero - q(1,2) = -1.0_wp/2.0_wp; q(2,2) = zero; q(3,2) = 1.0_wp/2.0_wp + call coeffs_1_2_1 ( a, q ) first = .false. end if @@ -112,13 +109,13 @@ end subroutine deriv_gf_2_1 subroutine up_deriv_gf_2_1 ( var, ni, nj, nk, dir, bb, gsize, delta, up, dvar ) + use All_Coeffs_mod + implicit none DECLARE_CCTK_FUNCTIONS DECLARE_CCTK_PARAMETERS - CCTK_REAL, parameter :: zero = 0.0 - integer, parameter :: wp = kind(zero) CCTK_INT, intent(IN) :: ni, nj, nk CCTK_REAL, dimension(ni,nj,nk), intent(IN) :: var CCTK_INT, intent(IN) :: dir @@ -137,28 +134,7 @@ subroutine up_deriv_gf_2_1 ( var, ni, nj, nk, dir, bb, gsize, delta, up, dvar ) logical, save :: first = .true. if ( first ) then - a1(-1) = -1.000000000000000000000000000000000000000_wp - a1(0) = 1.000000000000000000000000000000000000000_wp - a1(1) = 0.0_wp - - q1(1,1) = 0.0_wp - q1(2,1) = 0.0_wp - q1(3,1) = 0.0_wp - q1(1,2) = -1.000000000000000000000000000000000000000_wp - q1(2,2) = 1.000000000000000000000000000000000000000_wp - q1(3,2) = 0.0_wp - - a2(-1) = 0.0_wp - a2(0) = -1.000000000000000000000000000000000000000_wp - a2(1) = 1.000000000000000000000000000000000000000_wp - - q2(1,1) = -2.000000000000000000000000000000000000000_wp - q2(2,1) = 2.000000000000000000000000000000000000000_wp - q2(3,1) = 0.0_wp - q2(1,2) = 0.0_wp - q2(2,2) = -1.000000000000000000000000000000000000000_wp - q2(3,2) = 1.000000000000000000000000000000000000000_wp - + call coeffs_up_2_1 ( a1, q1, a2, q2 ) first = .false. end if |