[Nek5000-users] Problem with multiple file inut
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jun 3 07:17:17 CDT 2013
For the default/.f format you do not need to set p66 & 67 yourself but I would start with a printout of your p66 and p67 together with the filename you want to read.
Aleks
----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Monday, June 3, 2013 7:11:25 AM
Subject: Re: [Nek5000-users] Problem with multiple file inut
Aleks, what does this mean in my particular case.
Which name should the files have and which value should I prescribe for p66 and p67.
Thanks a lot, Joerg.
Am 03.06.2013 um 14:05 schrieb nek5000-users at lists.mcs.anl.gov:
> If I am not mistaken p66 and p67 are reset to be equal to 6 later in the code for .f/default format.
>
> The printed filename though has the empty base and extension .fld as if p67 is 4...
>
> Aleks
>
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Monday, June 3, 2013 7:00:15 AM
> Subject: Re: [Nek5000-users] Problem with multiple file inut
>
> Hi Aleks,
>
> p66 and p67 is set to 0.000E+00
>
> Thanks
>
> Am 03.06.2013 um 13:52 schrieb nek5000-users at lists.mcs.anl.gov:
>
>> Hi Joerg,
>>
>> Most likely -- for some reason, the reading routines get a filename ".fld":
>>
>> 0 1 OPEN: .fld
>> .fld
>>
>> so you might also want to check a value of parameter(67)
>>
>> Best.
>> Aleks
>>
>> ----- Original Message -----
>> From: nek5000-users at lists.mcs.anl.gov
>> To: nek5000-users at lists.mcs.anl.gov
>> Sent: Monday, June 3, 2013 5:46:23 AM
>> Subject: [Nek5000-users] Problem with multiple file inut
>>
>> Hi Neks,
>>
>> we want to run a post processing job analyzing further things of our full production job.
>> For output in single file things work well with the following patch of code in the usr file of
>> the analysis program:
>>
>> IF (istep.eq.1) THEN
>>
>> do k=1,xx
>> WRITE(fldfle,'("G1_640.f",I5.5)')k
>> initc(1) = trim(fldfle)
>> call restart(1)
>>
>> …… do your analysis……..
>>
>> enddo
>>
>> ENDIF
>>
>> It does the analysis for xx snapshots G1_640.f00001 to G1_640.f000xx that we generated
>> before.
>>
>> Bigger jobs requires multiple file outputs, e.g. here each snapshot goes into 4 separate files
>> G1_640.f0000x to G1_643.f0000x
>>
>> I thought that the corresponding extension should do it:
>>
>>
>> IF (istep.eq.1) THEN
>>
>> do k=1,xx
>> k1=4+k
>> WRITE(fldfle,'("G1_64?.f",I5.5)')k1
>> initc(1) = trim(fldfle)
>> call restart(4)
>>
>> …… do your analysis……..
>>
>> enddo
>>
>> ENDIF
>>
>>
>> the patch of log file (see below) indicates that the program reads G1_640.f00005 to G1_643.f00005 correctly and
>> gets stuck afterwards when trying to read the four files G1_64*.f00006:
>>
>>
>>
>> New CG1-tolerance (RINIT*epsm) = 0.874462285180112059E-17 0.607013350432080170E-25
>> 1 1 Helmholtz VELZ F: 8.7446E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 2 Helmholtz VELZ F: 4.5200E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 3 Helmholtz VELZ F: 3.1000E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 4 Helmholtz VELZ F: 2.4871E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 5 Helmholtz VELZ F: 2.0362E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 6 Helmholtz VELZ F: 1.7312E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 7 Helmholtz VELZ F: 1.4629E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 8 Helmholtz VELZ F: 1.2836E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 9 Helmholtz VELZ F: 1.1274E-05 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 10 Helmholtz VELZ F: 9.7394E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 11 Helmholtz VELZ F: 8.3244E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 12 Helmholtz VELZ F: 6.8927E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 13 Helmholtz VELZ F: 5.5410E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 14 Helmholtz VELZ F: 4.4505E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 15 Helmholtz VELZ F: 3.5729E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 16 Helmholtz VELZ F: 2.8812E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 17 Helmholtz VELZ F: 2.3406E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 18 Helmholtz VELZ F: 1.9136E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 19 Helmholtz VELZ F: 1.5656E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 20 Helmholtz VELZ F: 1.2877E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 21 Helmholtz VELZ F: 1.0656E-06 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 22 Helmholtz VELZ F: 8.8191E-07 1.0000E-06 8.3666E-05 2.0000E+02
>> 1 Hmholtz VELZ: 21 8.8191E-07 8.7446E-05 1.0000E-06
>> L1/L2 DIV(V) : 2.0122E-21 5.8478E-05
>> L1/L2 QTL : 0.0000E+00 0.0000E+00
>> L1/L2 DIV(V)-QTL: 2.0122E-21 5.8478E-05
>> 1 5.0000E-03 2.8596E+00 Fluid done
>> filt amp .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0000 .0500
>> filt trn 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 .9500
>> 1 5.0000E-03 -1.3839E-88 1.0000E+00 -3.9134E-06 3.9717E-06 uz_t_mx
>> Reading checkpoint data
>> 0 1 OPEN: G1_640.f00005
>> byte swap: F 6.543210030 -0.2931277218E+36
>> 12288 1 OPEN: G1_643.f00005
>> 4096 1 OPEN: G1_641.f00005
>> 8192 1 OPEN: G1_642.f00005
>>
>> 1 2.2370E+02 done :: Read checkpoint data
>> avg data-throughput = -80.6MBps
>> io-nodes = 4
>>
>> 0 1 OPEN: .fld
>> .fld
>> byte_read() :: fopen failure2!
>>
>>
>> Parameter 65 in the rea file is set to 4. Could it be that the initc is not correctly set?
>>
>> Thanks a lot, Joerg.
>>
>> _______________________________________________
>> 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
>
> _______________________________________________
> 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
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
More information about the Nek5000-users
mailing list