4 dt, list_mode, un_m1, un, pn_m1, pn, phin_m1, phin)
7 REAL(KIND=8),
INTENT(OUT):: time
8 REAL(KIND=8),
INTENT(IN) :: dt
9 INTEGER,
DIMENSION(:),
INTENT(IN) :: list_mode
10 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: un_m1, un
11 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: pn_m1, pn, phin_m1, phin
17 REAL(KIND=8),
INTENT(OUT):: time
18 REAL(KIND=8),
INTENT(IN) :: dt
19 INTEGER,
DIMENSION(:),
INTENT(IN) :: list_mode
20 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: tempn_m1, tempn
25 dt, list_mode, level_set_m1, level_set)
28 REAL(KIND=8),
INTENT(OUT):: time
29 REAL(KIND=8),
INTENT(IN) :: dt
30 INTEGER,
DIMENSION(:),
INTENT(IN) :: list_mode
31 REAL(KIND=8),
DIMENSION(:,:,:,:),
INTENT(OUT):: level_set, level_set_m1
35 opt_density, opt_tempn) result(vv)
36 INTEGER ,
INTENT(IN) :: type
37 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
38 INTEGER ,
INTENT(IN) :: mode, i
39 REAL(KIND=8),
INTENT(IN) :: time
40 REAL(KIND=8),
INTENT(IN) :: re
41 CHARACTER(LEN=2),
INTENT(IN) :: ty
42 REAL(KIND=8),
DIMENSION(:,:,:),
OPTIONAL,
INTENT(IN) :: opt_density
43 REAL(KIND=8),
DIMENSION(:,:,:),
OPTIONAL,
INTENT(IN) :: opt_tempn
44 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
48 INTEGER ,
INTENT(IN) :: type
49 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
50 INTEGER ,
INTENT(IN) :: m
51 REAL(KIND=8),
INTENT(IN) ::
t
52 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
56 INTEGER ,
INTENT(IN) :: type
57 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
58 INTEGER ,
INTENT(IN) :: m, interface_nb
59 REAL(KIND=8),
INTENT(IN) ::
t
60 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
64 INTEGER ,
INTENT(IN) :: type
65 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
66 INTEGER,
INTENT(IN) :: m
67 REAL(KIND=8),
INTENT(IN) ::
t
68 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
72 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
73 REAL(KIND=8),
INTENT(IN) ::
t
74 REAL(KIND=8),
DIMENSION(SIZE(rr,2),6) :: vv
78 INTEGER ,
INTENT(IN) :: type
79 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
80 INTEGER ,
INTENT(IN) :: m
81 REAL(KIND=8),
INTENT(IN) ::
t
82 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
86 INTEGER,
INTENT(IN) :: type
87 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
88 INTEGER ,
INTENT(IN) :: m
89 REAL(KIND=8),
INTENT(IN) ::
t
90 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
94 INTEGER ,
INTENT(IN) :: type
95 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
96 INTEGER ,
INTENT(IN) :: m, interface_nb
97 REAL(KIND=8),
INTENT(IN) ::
t
98 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
104 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr_gauss
105 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: angles
106 INTEGER,
INTENT(IN) :: nb_angles
107 INTEGER,
INTENT(IN) :: nb, ne
108 REAL(KIND=8),
INTENT(IN) :: time
109 REAL(KIND=8),
DIMENSION(nb_angles,ne-nb+1) :: vv
115 INTEGER ,
INTENT(IN) :: type, n_start
116 INTEGER,
INTENT(IN) :: mode
117 REAL(KIND=8),
INTENT(IN) ::
t
118 REAL(KIND=8),
DIMENSION(H_Mesh%np) :: vv
121 FUNCTION sub_vexact(m, H_mesh) RESULT(vv) !Set uniquement a l'induction
124 INTEGER,
INTENT(IN) :: m
125 REAL(KIND=8),
DIMENSION(H_mesh%np,6) :: vv
129 CHARACTER(LEN=1),
INTENT(IN) :: char_h_b
130 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
131 INTEGER,
INTENT(IN) :: m
132 REAL(KIND=8),
DIMENSION(SIZE(rr,2),6) :: vv
135 FUNCTION sub_hexact(H_mesh, TYPE, rr, m, mu_H_field, t) RESULT(vv)
138 INTEGER ,
INTENT(IN) :: type
139 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
140 INTEGER ,
INTENT(IN) :: m
141 REAL(KIND=8),
INTENT(IN) ::
t
142 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: mu_h_field
143 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
147 INTEGER ,
INTENT(IN) :: type
148 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
149 INTEGER ,
INTENT(IN) :: m
150 REAL(KIND=8),
INTENT(IN) :: mu_phi,
t
151 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
155 mesh_id, opt_b_ext) result(vv)
156 INTEGER ,
INTENT(IN) :: type
157 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: rr
158 INTEGER ,
INTENT(IN) :: m
159 REAL(KIND=8),
INTENT(IN) :: mu_phi, sigma, mu_h,
t
160 INTEGER ,
INTENT(IN) :: mesh_id
161 REAL(KIND=8),
DIMENSION(6),
OPTIONAL,
INTENT(IN) :: opt_b_ext
166 INTEGER,
INTENT(IN) :: type
167 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: rr
168 INTEGER,
INTENT(IN) :: m
169 REAL(KIND=8),
INTENT(IN) :: mu_phi, sigma, mu_h,
t
174 list_mode, hn1, hn, phin1, phin)
177 REAL(KIND=8),
INTENT(OUT):: time
178 REAL(KIND=8),
INTENT(IN) :: dt
179 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: mu_h_field
180 REAL(KIND=8),
INTENT(IN) :: mu_phi
181 INTEGER,
DIMENSION(:),
INTENT(IN) :: list_mode
182 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: hn, hn1
183 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: phin, phin1
189 REAL(KIND=8),
DIMENSION(ne-nb+1) :: vv
190 INTEGER,
INTENT(IN) :: nb, ne
191 REAL(KIND=8),
DIMENSION(2,ne-nb+1),
OPTIONAL :: pts
192 INTEGER,
DIMENSION(ne-nb+1),
OPTIONAL :: pts_ids
196 REAL(KIND=8),
DIMENSION(2),
INTENT(in):: pt
197 INTEGER,
DIMENSION(1),
INTENT(in) :: pt_id
198 REAL(KIND=8),
DIMENSION(2) :: vv
204 REAL(KIND=8),
DIMENSION(:),
INTENT(IN) :: angles
205 INTEGER,
INTENT(IN) :: nb_angles
206 INTEGER,
INTENT(IN) :: nb, ne
207 REAL(KIND=8),
INTENT(IN) :: time
208 REAL(KIND=8),
DIMENSION(nb_angles,ne-nb+1) :: vv
214 REAL(KIND=8),
DIMENSION(SIZE(H_mesh%rr,2)) :: vv
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 t
real(kind=8) function sub_jexact_gauss(TYPE, rr, m, mu_phi, sigma, mu_H, t, mesh_id, opt_B_ext)
real(kind=8) function, dimension(size(rr, 2), 6) sub_h_b_quasi_static(char_h_b, rr, m)
real(kind=8) function, dimension(size(rr, 2)) sub_phiexact(TYPE, rr, m, mu_phi, t)
real(kind=8) function, dimension(size(rr, 2)) sub_source_in_ns_momentum(TYPE, rr, mode, i, time, Re, ty, opt_density, opt_tempn)
real(kind=8) function sub_eexact_gauss(TYPE, rr, m, mu_phi, sigma, mu_H, t)
real(kind=8) function, dimension(h_mesh%np, 6) sub_vexact(m, H_mesh)
real(kind=8) function, dimension(size(rr, 2)) sub_hexact(H_mesh, TYPE, rr, m, mu_H_field, t)