47     integer num_xyz(3)  local
70     num_xyz(d) = ISHFT(1, sd)
82     integer num_xyz(3)  local
111       num_elem = num_elem * num_xyz(i)
112       nd(i) = num_xyz(i) * (p - 1) + 1
129         e_xyz(j) = mod(re, num_xyz(j))
130         re = re/num_xyz(j)
198     integer num_xyz(3)  local
218       nd(i) = num_xyz(i) * (p - 1) + 1
239 …      coords(scalar_size * (j - 1) + i) = ((d_1d/(p - 1)) + nodes(mod(d_1d, p - 1) + 1))/num_xyz(j)
259     integer num_elem, num_xyz(3), elem_qpts  variable
384     call getcartesianmeshsize(fe_dim, sol_degree, prob_size, num_xyz)
387     write (*, '(A16, I8)', advance='no') 'Mesh size: nx = ', num_xyz(1)
389       write (*, '(A7, I8)', advance='no') ',  ny = ', num_xyz(2)
392       write (*, '(A7, I8)', advance='no') ',  nz = ', num_xyz(3)
399 …call buildcartesianrestriction(ceed, fe_dim, num_xyz, mesh_degree, num_comp_x, mesh_size, num_qpts…
401 …call buildcartesianrestriction(ceed, fe_dim, num_xyz, sol_degree, 1, sol_size, num_qpts, sol_restr…
414     call setcartesianmeshcoords(fe_dim, num_xyz, mesh_degree, mesh_coords, exact_volume, err)
457       num_elem = num_elem * num_xyz(i)