[Nek5000-users] turbChannel: history points and instantaneous velocity field
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Feb 20 12:13:27 CST 2018
Make sure you file is called <casename>.his
On 20 Feb 2018, at 16:55, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:
Hello,
Thank you for the reply. The istep loop shall be fine. I tried some other lines in the “if..endif” loop here and they all worked.
I also commented the if loop and I had the same error information. I really do not know why.
Regards,
Emily
_______________________________________________________________________
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 <mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180220/ddb8b695/attachment.html>
More information about the Nek5000-users
mailing list