[Nek5000-users] hpts and pressure

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

Hello Neks,

I am trying to output pressure data at various points using the hpts routine.

The values that are output for the first time step look ok, but all
subsequent time steps seem to give garbage (i.e. the pressures given
in hpts.out are way different that what is given by visualizing
pressure in VisIt).  I know previously that this was an issue for
someone and it was thought that passive scalars may be messing the
hpts output up, but in my case I don't output temperature or pressure.

I looked at the hpts routine and noticed that it uses pm1 when it is
stacking the field variables into the working array.  Since I am using
a Pn-Pn-2 formulation for my case, do I need to use the mappr
subroutine in my usrchk to map the pressure from mesh 2 to mesh 1
before calling hpts()?

I'm assuming the issue arises because I am calling hpts more
frequently than I am outputing field files, and thus perhaps pm1 isn't
properly being updated (or perhaps it may have already been
destroyed).

Thanks,

-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do
nothing" -- Edmund Burke