[Nek5000-users] SFD error message from p66

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

Hi Adam,

thank you very much.

I have solved the problem that I posted here about the SFD code and now 
it seems to work. However, I'm still doing some tests to validate the 
code. It will take me some more time. If you agree, I think it would be 
nice to discuss in private about these two versions of the SFD code. I 
will let you know how this code works and we could compare this version 
with your code.

Regarding the error that I posted, the problem was not related with 
param 66, but with the dimension of the variables.

In the correct version of the code:

ntot=nx1*ny1*nz1*nelv
call rzero(vx,ntot)

However, in the wrong version I was defining ntot as:

ntot=nx1*ny1*nz1*nz1*nelv

The wrong definition of the dimensions of the velocity vector was 
leading to the error that I posted.

Regards
SL

El 30-06-2014 09:51, nek5000-users at lists.mcs.anl.gov escribió:
> Hi SL
> 
> We have implemented SFD some time ago and we can share the code if you
> wish. It is probably good to have single tested version of the code
> than develop number of variations. We've already shared the code with
> Elia.
> Regards
> Adam
> 
> 
> On 27/06/14 17:28, nek5000-users at lists.mcs.anl.gov wrote:
>> Dear Nek's,
>> 
>> I'm trying to program a Selective Frequency Damping following the 
>> article of Akervik et. al (Physics of Fluids, 2006) and some Nek5000 
>> posts.
>> 
>> I'm able to compile the code, however, when I try to run my 
>> simulation, I find the following error message:
>> 
>>         0  0.0000E+00 Write checkpoint:
>> 
>>         0  0.0000E+00 OPEN:
>> At line 335 of file User/nek5_svn/trunk/nek/prepost.f (unit = 24, file 
>> = '')
>> Fortran runtime error: File '' does not exist
>> 
>> It seems to be something related to param 66. I've tried with p66=6, 
>> 0, 1, -1, but I always get the same error.
>> 
>> Do you know what happens?
>> 
>> Thanks
>> SL
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users