4 REAL (KIND=8),
PARAMETER,
PUBLIC :: ratio_mu_T22 = 50.d0
5 REAL (KIND=8),
PUBLIC :: b_factor_T22 = (2**6) * (ratio_mu_T22-1.d0)/(ratio_mu_T22+1.d0)
6 INTEGER,
PUBLIC :: mode_mu_T22 = 4
12 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: r,
z
13 REAL(KIND=8),
DIMENSION(SIZE(r)) :: vv
14 vv = b_factor_t22*(r*(1-r)*(
z**2-1))**3
20 REAL(KIND=8),
INTENT(IN):: r,
z
22 vv = 3 * b_factor_t22 * (
z**2-1)**3 * (r*(1-r))**2 * (1-2*r)
28 REAL(KIND=8),
INTENT(IN):: r,
z
30 vv = 3*b_factor_t22*(r*(1-r))**3*(
z**2-1)**2*(2*
z)
41 REAL(KIND=8),
DIMENSION(ne-nb+1) :: vv
43 REAL(KIND=8),
DIMENSION(2,ne-nb+1),
OPTIONAL :: pts
44 INTEGER,
DIMENSION(ne-nb+1),
OPTIONAL :: pts_ids
45 REAL(KIND=8),
DIMENSION(ne-nb+1) :: r,
z
47 IF( present(pts) .AND. present(pts_ids) )
THEN
64 REAL(KIND=8),
DIMENSION(2) :: pt,vv
65 INTEGER,
DIMENSION(1) :: pt_id
66 REAL(KIND=8),
DIMENSION(1) :: tmp,r,
z
73 IF (tmp(1) .GE. 0.d0 )
THEN
92 REAL(KIND=8),
DIMENSION(:) :: angles
95 REAL(KIND=8),
DIMENSION(nb_angles,ne-nb+1) :: vv
99 vv(ang,:) = 1/(1+
f_test_t22(h_mesh%rr(1,nb:ne),h_mesh%rr(2,nb:ne))*cos(mode_mu_t22*angles(ang)))
real(kind=8) function, dimension(ne-nb+1) mu_bar_in_fourier_space_anal_t22(H_mesh, nb, ne, pts, pts_ids)
real(kind=8) function dfdr_test_t22(r, z)
real(kind=8) function, dimension(size(r)) f_test_t22(r, z)
real(kind=8) function, dimension(2) grad_mu_bar_in_fourier_space_anal_t22(pt, pt_id)
real(kind=8) function dfdz_test_t22(r, z)
real(kind=8) function, dimension(nb_angles, ne-nb+1) mu_in_real_space_anal_t22(H_mesh, angles, nb_angles, nb, ne)
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