[Nek5000-users] Hpts subroutine

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

Hello,

This email is in regards to the hpts subroutine that was added to the repo
513 and subsequent versions.
A part of the code snippet from postpro.f ( subroutine hpts ) is pasted
below :

<code>

      call copy(wrk(1,1,1,1,2),vy,ntot)
      if(if3d) call copy(wrk(1,1,1,1,2),vz,ntot)

</code>

could the second line be " if(if3d) call copy(wrk(1,1,1,1,3),vz,ntot) " else
vz would then be overwritten over vy ?

Along the same lines, I would like to know what one would have to change to
output the pressure values as well.

Thanks
Shriram
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100518/acdfd9f5/attachment.html>