[Nek5000-users] hpts issue: not dumping anything

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

Hi Michael,

Somehow your buffered I/O device does not get flushed when the nek5000
process ends! Typically the OS triggers a close of all open file
handles which will flush the I/O buffer in case of a write access.

I guess there is an easy way to fix this. Just flush the considered
unit (file) explicitly every n timesteps in userchk(). Note: This is a
quick and dirty workaround.

Stefan


On 3/28/11, nek5000-users at lists.mcs.anl.gov
<nek5000-users at lists.mcs.anl.gov> wrote:
>
> Hi Michael,
>
> You would need to close it only after your
> run is done?  Otherwise, it closes and (I guess)
> reopens after each call and then you for sure
> wouldn't get anything.
>
> It's hard for me to check w/o being on precisely
> the same platform.
>
> Unfortunately I'm on travel right now - otherwise
> we could perhaps sort it out through skype.
>
> Paul
>
> On Sun, 27 Mar 2011, nek5000-users at lists.mcs.anl.gov wrote:
>
>>
>>
>> Hi Paul,
>>
>>
>>
>> Do you mean add to the hpts routine within postpro.f an explicit close( )
>> command after writing to hpts.out?  I see what you mean
>>
>> that nothing actually shows up unless the run is closed, in which case all
>> of my data now appears.
>>
>>
>>
>> - Michael
>>
>>
>> ----- Original Message -----
>> From: nek5000-users at lists.mcs.anl.gov
>> To: "Nekton User List" <nek5000-users at lists.mcs.anl.gov>
>> Sent: Tuesday, March 15, 2011 1:54:08 AM GMT -06:00 US/Canada Central
>> Subject: Re: [Nek5000-users] hpts issue: not dumping anything
>>
>> Hi Michael, I've had some experiences on IBM machines in the past where
>> the file never appeared unless it was explicitly closed. Surprisingly,
>> this was the case even if I properly exitted. I've not looked at the hpts
>> source, but this may be the issue. Can you try closing the file prior to
>> exitting? Obviously this is a pain in a situation where you want to let
>> your run keep going until terminated by the end of your allocated block,
>> but I wouldn't know anyway around it. There may be a compiler option to
>> force it to not buffer the I/O so parsimoniously. Paul On Mon, 14 Mar
>> 2011, nek5000-users at lists.mcs.anl.gov wrote: > > > Hi, > > > > I am having
>> an issue with hpts.  The issue is on a machine with the operating system
>> AIX 5.3.  The compilers > > are mpxlf for Fortran, and mpcc for  the C
>> compiler. The history point routine works on all other machines I > > have
>> used, except for this one. > > > > The simulation runs fine, saying that
>> it is dumping history points to hpts.ou
>  t in the output log.  However, nothing is ever actually written > > to
> hpts.out, not even the header. > > > > I am not sure what the issue is. I
> suspect something with the compilers or operating system as this is the only
> difference > > between this case, and all other cases.   Any thoughts? > > >
>> - Michael M.
>> _______________________________________________ Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users