Menu
▾
▴
#16 bug in TOSCA maps reading (V5.1)
It seems to me that the reading of TOSCA files has a problem : while reading the header (written out in zgoubi.res) the reading of the file pass over these first lines, that are not read again in the "true" map reading. Finally that cause a read past the end of the file.
I think I solved this problem by adding a REWIND(LUN) after the header reading, to go back to the first line of the field map file.
I attached here the modified fmapw.f
Thanks for the report and the fix. Just for reference could you attach the TOSCA file and your zgoubi.dat that give the problem (assuming it will let you, otherwise post a link)? Thanks
As an illustratiob, I put you here the zgoubi.dat file I used, as well as the beginning of the TOSCA MAP.
Test Qpole
'OBJET'
563.
6
1.0 10. 1.0 10. 0.1 0.0
0. 0. 0. 0. 0. 1.
'TOSCA'
2 2
5. 1. 1. 1.
QP584
321 61 1 1
QP584_01.data
0 0. 0. 0.
2
1.
1 0. 0. 0.
'MATRIX'
1 0
'FIN'
With a part of the QP584_01.data file :
-0.1500E+02 -0.1500E+02 0.0000E+00 0.4424E-02 0.4403E-02 -0.2492E-01
-0.1500E+02 -0.1500E+02 0.5000E+00 0.4389E-02 0.4367E-02 -0.2573E-01
-0.1500E+02 -0.1500E+02 0.1000E+01 0.4349E-02 0.4328E-02 -0.2653E-01
-0.1500E+02 -0.1500E+02 0.1500E+01 0.4308E-02 0.4286E-02 -0.2733E-01
-0.1500E+02 -0.1500E+02 0.2000E+01 0.4262E-02 0.4240E-02 -0.2812E-01
-0.1500E+02 -0.1500E+02 0.2500E+01 0.4213E-02 0.4191E-02 -0.2891E-01
-0.1500E+02 -0.1500E+02 0.3000E+01 0.4159E-02 0.4138E-02 -0.2969E-01
-0.1500E+02 -0.1500E+02 0.3500E+01 0.4102E-02 0.4081E-02 -0.3046E-01
-0.1500E+02 -0.1500E+02 0.4000E+01 0.4040E-02 0.4019E-02 -0.3122E-01
-0.1500E+02 -0.1500E+02 0.4500E+01 0.3973E-02 0.3952E-02 -0.3198E-01
-0.1500E+02 -0.1500E+02 0.5000E+01 0.3901E-02 0.3880E-02 -0.3272E-01
-0.1500E+02 -0.1500E+02 0.5500E+01 0.3822E-02 0.3802E-02 -0.3346E-01
-0.1500E+02 -0.1500E+02 0.6000E+01 0.3738E-02 0.3718E-02 -0.3418E-01
(...)
(with 361*61 lines like this)
Cheers,
Antoine
Full tosca file
Zgoubi.dat input (example with only 1 Tosca file)