sample_points Subroutine

public subroutine sample_points(vals, points, a, b, c, n_points, n_f, n_k, d, grid)

vals(1:3, p) = this rank's part of the trilinear interpolation of the DIR_X VERT fields a, b and c at points(:, p), one thread per point; see field_sample_points. grid holds the local VERT dims, global VERT dims, rank offsets and periodic flags (1/0), three entries each, then the padded y-block count of the DIR_X layout.

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(out) :: vals(3,n_points)
real(kind=dp), intent(in) :: points(3,n_points)
real(kind=dp), intent(in) :: a(SZ,n_f,n_k)
real(kind=dp), intent(in) :: b(SZ,n_f,n_k)
real(kind=dp), intent(in) :: c(SZ,n_f,n_k)
integer, intent(in) :: n_points
integer, intent(in) :: n_f
integer, intent(in) :: n_k
real(kind=dp), intent(in) :: d(3)
integer, intent(in) :: grid(13)