[Nek5000-users] Problem with multiple file inut

Mon Jun 3 09:11:49 CDT 2013

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
>