[Nek5000-users] Hpts subroutine

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

Yes you're right this is a bug. I just coded that up very quickly and haven't had the time to check it. I will fix it and add support for pressure.

Stefan


On May 18, 2010, at 5:41 PM, nek5000-users at lists.mcs.anl.gov wrote:

> 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
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users