[Nek5000-users] Problem with multiple file inut
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jun 3 06:52:16 CDT 2013
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
More information about the Nek5000-users
mailing list