[Nek5000-users] Problem with multiple file inut

Mon Jun 3 10:07:41 CDT 2013

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