69 delta, c_dyn, test_filter, mij, lij, num, den)
70 logical,
intent(in) :: if_ext
71 real(kind=
rp),
intent(in) :: t
72 integer,
intent(in) :: tstep
73 type(
coef_t),
intent(in) :: coef
74 type(
field_t),
intent(inout) :: nut
75 type(
field_t),
intent(in) :: delta
76 type(
field_t),
intent(inout) :: c_dyn
78 type(
field_t),
intent(inout) :: mij(6), lij(6)
79 type(
field_t),
intent(inout) :: num, den
81 type(
field_t),
pointer :: u, v, w
84 type(
field_t),
pointer :: s11, s22, s33, s12, s13, s23, s_abs
85 real(kind=
rp) :: alpha
86 integer :: temp_indices(7)
89 if (tstep .eq. 1)
then
95 if (if_ext .eqv. .true.)
then
114 call strain_rate(s11%x, s22%x, s33%x, s12%x, s13%x, s23%x, u, v, w, coef)
116 call coef%gs_h%op(s11%x, s11%dof%size(),
gs_op_add)
117 call coef%gs_h%op(s22%x, s11%dof%size(),
gs_op_add)
118 call coef%gs_h%op(s33%x, s11%dof%size(),
gs_op_add)
119 call coef%gs_h%op(s12%x, s11%dof%size(),
gs_op_add)
120 call coef%gs_h%op(s13%x, s11%dof%size(),
gs_op_add)
121 call coef%gs_h%op(s23%x, s11%dof%size(),
gs_op_add)
124 s12%x_d, s13%x_d, s23%x_d, &
129 s11, s22, s33, s12, s13, s23, &
130 s_abs, test_filter, delta, alpha, &
134 call device_col2(nut%x_d, coef%mult_d, nut%dof%size())
148 type(
field_t),
intent(inout) :: lij(6)
149 type(
field_t),
pointer,
intent(in) :: u, v, w
151 integer,
intent(in) :: n
154 type(
field_t),
pointer :: fu, fv, fw, fuu, fvv, fww, fuv, fuw, fvw
155 integer :: temp_indices(9)
168 call test_filter%apply(fu, u)
169 call test_filter%apply(fv, v)
170 call test_filter%apply(fw, w)
175 lij(4)%x_d, lij(5)%x_d, lij(6)%x_d, &
176 u%x_d, v%x_d, w%x_d, &
177 fu%x_d, fv%x_d, fw%x_d, &
178 fuu%x_d, fvv%x_d, fww%x_d, &
179 fuv%x_d, fuw%x_d, fvw%x_d, n)
182 call test_filter%apply(fuu, fuu)
183 call test_filter%apply(fvv, fvv)
184 call test_filter%apply(fww, fww)
185 call test_filter%apply(fuv, fuv)
186 call test_filter%apply(fuw, fuw)
187 call test_filter%apply(fvw, fvw)
193 lij(4)%x_d, lij(5)%x_d, lij(6)%x_d, &
194 fuu%x_d, fvv%x_d, fww%x_d, &
195 fuv%x_d, fuw%x_d, fvw%x_d, n)
220 s11, s22, s33, s12, s13, s23, &
221 s_abs, test_filter, delta, alpha, &
223 type(
field_t),
intent(inout) :: nut, c_dyn
224 type(
field_t),
intent(inout) :: num, den
225 type(
field_t),
intent(in) :: lij(6)
226 type(
field_t),
intent(inout) :: mij(6)
227 type(
field_t),
intent(inout) :: s11, s22, s33, s12, s13, s23, s_abs
229 type(
field_t),
intent(in) :: delta
230 real(kind=
rp),
intent(in) :: alpha
231 type(
coef_t),
intent(in) :: coef
232 integer,
intent(in) :: n
234 real(kind=
rp) :: delta_ratio2
235 integer :: temp_indices(13)
236 type(
field_t),
pointer :: fs11, fs22, fs33, fs12, fs13, fs23, fs_abs, &
237 fsabss11, fsabss22, fsabss33, &
238 fsabss12, fsabss13, fsabss23
240 delta_ratio2 = ((test_filter%nx-1.0_rp)/(test_filter%nt-1.0_rp))**2
259 call test_filter%apply(fs_abs, s_abs)
260 call test_filter%apply(fs11, s11)
261 call test_filter%apply(fs22, s22)
262 call test_filter%apply(fs33, s33)
263 call test_filter%apply(fs12, s12)
264 call test_filter%apply(fs13, s13)
265 call test_filter%apply(fs23, s23)
270 mij(4)%x_d, mij(5)%x_d, mij(6)%x_d, &
271 s_abs%x_d, s11%x_d, s22%x_d, s33%x_d, &
272 s12%x_d, s13%x_d, s23%x_d, &
273 fs_abs%x_d, fs11%x_d, fs22%x_d, fs33%x_d, &
274 fs12%x_d, fs13%x_d, fs23%x_d, &
275 fsabss11%x_d, fsabss22%x_d, fsabss33%x_d, &
276 fsabss12%x_d, fsabss13%x_d, fsabss23%x_d, &
280 call test_filter%apply(fsabss11, fsabss11)
281 call test_filter%apply(fsabss22, fsabss22)
282 call test_filter%apply(fsabss33, fsabss33)
283 call test_filter%apply(fsabss12, fsabss12)
284 call test_filter%apply(fsabss13, fsabss13)
285 call test_filter%apply(fsabss23, fsabss23)
289 mij(4)%x_d, mij(5)%x_d, mij(6)%x_d, &
290 lij(1)%x_d, lij(2)%x_d, lij(3)%x_d, &
291 lij(4)%x_d, lij(5)%x_d, lij(6)%x_d, &
292 fsabss11%x_d, fsabss22%x_d, fsabss33%x_d, &
293 fsabss12%x_d, fsabss13%x_d, fsabss23%x_d, &
294 num%x_d, den%x_d, c_dyn%x_d, delta%x_d, &
295 s_abs%x_d, nut%x_d, alpha, coef%mult_d, n)
subroutine, public device_mij_compute_part1(m11_d, m22_d, m33_d, m12_d, m13_d, m23_d, s_abs_d, s11_d, s22_d, s33_d, s12_d, s13_d, s23_d, fs_abs_d, fs11_d, fs22_d, fs33_d, fs12_d, fs13_d, fs23_d, fsabss11_d, fsabss22_d, fsabss33_d, fsabss12_d, fsabss13_d, fsabss23_d, delta_ratio2, n)
part 1 of the computing of the mij field
subroutine, public device_mij_nut_compute_part2(m11_d, m22_d, m33_d, m12_d, m13_d, m23_d, l11_d, l22_d, l33_d, l12_d, l13_d, l23_d, fsabss11_d, fsabss22_d, fsabss33_d, fsabss12_d, fsabss13_d, fsabss23_d, num_d, den_d, c_dyn_d, delta_d, s_abs_d, nut_d, alpha, mult_d, n)
part 1 of the computing of the mij field