[Nek5000-users] Problem with multiple file inut

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Jun 3 10:13:51 CDT 2013


Dear Aleks, thanks I will try this out. Joerg.

Am 03.06.2013 um 17:07 schrieb nek5000-users at lists.mcs.anl.gov:

> Sorry, Joerg,
> 
> I missed the part of your email where you said
> 
> "
> Bigger jobs requires multiple file outputs, e.g. here each snapshot goes into 4 separate files G1_640.f0000x to G1_643.f0000x
> "
> 
> so as I understand now, you case works fine while restarting from a single file, i.e. G1_640.f00005, and your question is how to restart from a field written to 4 files G1_64[0-3].f00005
> 
> If my understand is (finely :) correct then you need to
> 
>      character*13 fnames
> c
> c               1234567890123
> c
>      fnames = 'G1_64?.f00005'
> 
>      call chcopy (initc(1),fnames,13)
>      call restart(1)
> 
> No need to change p67
> 
> Aleks
> 
> 
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Monday, June 3, 2013 9:11:49 AM
> Subject: Re: [Nek5000-users] Problem with multiple file inut
> 
> 
> Hi Joerg,
> 
> Sorry - just saw this email chain - which makes my previous
> response senseless.
> 
> What you're doing looks correct to me.
> 
> Note that p65 has no role in input, only output (why?).
> 
> We can start digging into this.
> 
> Paul
> 
> 
> 
> 
> 
> On Mon, 3 Jun 2013, nek5000-users at lists.mcs.anl.gov wrote:
> 
>> 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



More information about the Nek5000-users mailing list