SetPoints_Points_t Subroutine

  subroutine SetPoints_Points_t(this,xIn)
    !! Copy user-supplied physical coordinates into the cloud.
    !! xIn must be shape (nPoints, nDim) matching the init dimensions.
    implicit none
    class(Points_t),intent(inout) :: this
    real(prec),intent(in) :: xIn(:,:)
    ! Local
    integer :: p,d

    if(size(xIn,1) /= this%nPoints .or. size(xIn,2) /= this%nDim) then
      print*,"SELF_Points_t::SetPoints: shape mismatch"
      stop 1
    endif

    do p = 1,this%nPoints
      do d = 1,this%nDim
        this%x(p,d) = xIn(p,d)
      enddo
    enddo
    this%elements = 0
    this%nCached = 0

  endsubroutine SetPoints_Points_t