4 REAL (KIND=8),
PARAMETER,
PRIVATE :: ratio_mu_T17 = 50.d0
5 REAL (KIND=8),
PRIVATE :: b_factor_T17_anal = (2**6) * (1.d0 - 1/ratio_mu_T17)
13 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: r,
z
14 REAL(KIND=8),
DIMENSION(SIZE(r)) :: vv
15 vv = b_factor_t17_anal*(r*(1-r)*(
z**2-1))**3
21 REAL(KIND=8),
INTENT(IN):: r,
z
23 vv = 3 * b_factor_t17_anal * (
z**2-1)**3 * (r*(1-r))**2 * (1-2*r)
29 REAL(KIND=8),
INTENT(IN):: r,
z
31 vv = 3*b_factor_t17_anal*(r*(1-r))**3*(
z**2-1)**2*(2*
z)
42 REAL(KIND=8),
DIMENSION(ne-nb+1) :: vv
44 REAL(KIND=8),
DIMENSION(2,ne-nb+1),
OPTIONAL :: pts
45 INTEGER,
DIMENSION(ne-nb+1),
OPTIONAL :: pts_ids
46 REAL(KIND=8),
DIMENSION(ne-nb+1) :: r,
z
48 IF( present(pts) .AND. present(pts_ids) )
THEN
64 REAL(KIND=8),
DIMENSION(2) :: vv
65 REAL(KIND=8),
DIMENSION(2) :: pt
66 INTEGER,
DIMENSION(1) :: pt_id
67 REAL(KIND=8),
DIMENSION(1) :: r,
z, tmp
real(kind=8) function, dimension(size(r)) f_test_t17(r, z)
real(kind=8) function, dimension(2) grad_mu_bar_in_fourier_space_anal_t17(pt, pt_id)
real(kind=8) function, dimension(ne-nb+1) mu_bar_in_fourier_space_anal_t17(H_mesh, nb, ne, pts, pts_ids)
real(kind=8) function dfdr_test_t17(r, z)
real(kind=8) function dfdz_test_t17(r, z)
section doc_intro_frame_work_num_app Numerical approximation subsection doc_intro_fram_work_num_app_Fourier_FEM Fourier Finite element representation The SFEMaNS code uses a hybrid Fourier Finite element formulation The Fourier decomposition allows to approximate the problem’s solutions for each Fourier mode modulo nonlinear terms that are made explicit The variables are then approximated on a meridian section of the domain with a finite element method The numerical approximation of a function f $f f is written in the following generic z