[Nek5000-users] turbChannel: history points and instantaneous velocity field

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

Hi Emily,

Did you call the hpts() subroutine without calling the istep loop?

Sincerely,

On Mon, Feb 19, 2018 at 3:41 PM <nek5000-users at lists.mcs.anl.gov> wrote:

> Dear Nek users,
>
>
>
> I am attempting to use History Points features in turbChannel example and
> unfortunately I run into an issue. In the subroutine userchck, I added:
>
>
>
> if(istep. eq. 1) then
>
>               call hpts()                       !checkpoints
>
> endif
>
>
>
> where the file is name “check.his” , the format of which follows the user
> guide:
>
> 2             !number of checkpoints
>
> 1. -0.98 1.   !x,y,z locations
>
> 2. -0.67 4.
>
>
>
> However, I received an error saying: Cannot open history file in
> subroutine hpts(). Has anyone run into this error before?
>
>
>
> Another question really bothering me is that I could not figure out how to
> write out the instantaneous velocity field (u(x,y,z),v(x,y,z),w(x,y,z)) at
> a specific time step (istep). Could anyone please help me on this as well?
>
>
>
> Regards,
>
>
>
> Emily
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-- 

*Amitvikram Dutta*

Graduate Research Assistant

Fluid Mechanics Research Lab

Multi-Physics Interaction Lab

University of Waterloo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180219/5cfcc243/attachment.html>