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

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

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180220/e0b6653a/attachment-0001.html>