7 REAL(KIND=8),
PARAMETER::
wjump_t26 = 0.032d0*(1.0d0)
15 REAL(KIND=8) :: x,x0,x1
17 REAL(KIND=8) :: a0,a1,a2,a3
23 a0 = x1**2*(3*x0-x1)/(x0-x1)**3;
24 a1 = -6.0*x0*x1/(x0-x1)**3;
25 a2 = (3.0*(x0+x1))/(x0-x1)**3;
28 vv = a0+a1*x+a2*x*x + a3*x*x*x
34 REAL(KIND=8) :: x,x0,x1
36 REAL(KIND=8) :: a0,a1,a2,a3
42 a0 = x1**2*(3*x0-x1)/(x0-x1)**3;
43 a1 = -6.0*x0*x1/(x0-x1)**3;
44 a2 = (3.0*(x0+x1))/(x0-x1)**3;
47 vv = a1+2.d0*a2*x + 3.d0*a3*x*x
52 REAL(KIND=8) :: x,x0,x1
64 REAL(KIND=8) :: x,x0,x1
85 ELSE IF ( z.LE.
z1_t26 )
THEN 100 REAL(KIND=8),
DIMENSION(2) :: vv
108 ELSE IF ( z.LE.
z1_t26 )
THEN 128 REAL(KIND=8),
DIMENSION(ne-nb+1) :: vv
130 REAL(KIND=8),
DIMENSION(2,ne-nb+1),
OPTIONAL :: pts
131 INTEGER,
DIMENSION(ne-nb+1),
OPTIONAL :: pts_ids
132 REAL(KIND=8),
DIMENSION(ne-nb+1) :: r,z
135 IF(
PRESENT(pts) .AND.
PRESENT(pts_ids) )
THEN 143 DO n = 1, ne - nb + 1
154 REAL(KIND=8),
DIMENSION(2) :: vv
155 REAL(KIND=8),
DIMENSION(2) :: pt
156 INTEGER,
DIMENSION(1) :: pt_id
158 REAL(KIND=8),
DIMENSION(2) :: tmp
real(kind=8) function, dimension(ne-nb+1) mu_bar_in_fourier_space_anal_t26(H_mesh, nb, ne, pts, pts_ids)
real(kind=8), parameter zm_t26
real(kind=8) function dsmooth_jump_up_t26(x, x0, x1)
real(kind=8) function smooth_jump_down_t26(x, x0, x1)
integer, parameter test_mode_t26
real(kind=8), parameter wjump_t26
real(kind=8) function smooth_jump_up_t26(x, x0, x1)
real(kind=8), parameter mu_disk_t26
real(kind=8), parameter alpha_t26
real(kind=8), parameter z1_t26
real(kind=8) function mu_func_t26(r, z)
real(kind=8) function, dimension(2) grad_mu_bar_in_fourier_space_anal_t26(pt, pt_id)
real(kind=8) function, dimension(2) dmu_func_t26(r, z)
real(kind=8) function dsmooth_jump_down_t26(x, x0, x1)