[Nek5000-users] Problem with multiple file inut

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

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.