[Nek5000-users] Interpolating data onto an evenly spaced grid

[nek5000-users at lists.mcs.anl.gov]

I want to  interpolate my nek5000 data fields (*.fxxxxx) onto an evenly
spaced grid (and globally evenly spaced so g2gi does not help here since my
elements are not all of the same size). It seems like hpts is set up for
this, where I read in a data file, call hpts (the data in hpts.in is my
evenly spaced grid), and output the data to hpts.out. However, I have some
pretty big grids, so I need to use the parallelization feature (i.e., set
lhis to npts/np, here npts = total number of interpolation points and np =
number of processors). However, then the data output to hpts.out is jumbled
up. If I instead set lhis to the total number of points, the output data is
fine. I suspect there is a bug somewhere in  hpts such that it mixes up the
data.

Any suggestions?

Cheers,

Janet

-- 
Janet Scheel
Assistant Professor of Physics
Occidental College
1600 Campus Road, M21
Los Angeles, CA 90041
(323) 259-2777
jscheel at oxy.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131113/c859980d/attachment.html>