RE: [fluka-discuss]: [Fwd: Run time error if cycle number is greater than 2]

From: Riya Dey <riyadey_at_barc.gov.in>
Date: Fri, 4 Sep 2020 16:25:21 +0530

Dear Anna Ferrari,

 

Thank you for your quick response. Yes, I have noticed this large size of error file also. But I am unable to understand the issue. I cannot attach the whole error file since it is having a large size. I am copying first few lines from that file. Please let me know if these lines represent any error.

 

 

 

1NUMBER OF BEAM NUMBER OF BEAM APPROXIMATE NUMBER AVERAGE TIME USED TIME LEFT (RESERVED NUMBER OF STARS

PARTICLES HANDLED PARTICLES LEFT OF BEAM PARTICLES BY A BEAM PARTICLE 10000.0 SECONDS CREATED

                                            THAT CAN STILL BE FOR PRINTOUT)

                                            HANDLED

 

NEXT SEEDS: 0 0 0 0 0 0 2 3039 0 0

EMFGEO, NCOUN: 1 TXYZ: 1.00000004550694 ETEMF, ICHEMF: 1.0000000E-04 0

X,Y,Z -97.4989700000000 165.896400000000 254.519100000000

 U,V,W 0.996221500000000 6.038126000000000E-02 -6.242493000000000E-02

USTEP,NREG 23616.1146498612 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

          1 999999 999999 1.2203903E+00 1.0000000E+30 0

 NEXT SEEDS: 2 0 0 0 0 0 2 3039 0 0

EMFGEO, NCOUN: 2 TXYZ: 0.999999958941740 ETEMF, ICHEMF: 6.4911190E-05 0

X,Y,Z 19.8812300000000 188.014200000000 265.233400000000

 U,V,W 0.634992300000000 -0.719271400000000 0.281839300000000

 USTEP,NREG 272.279550938379 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 14 TXYZ: 0.999999936862250 ETEMF, ICHEMF: 9.4525240E-05 0

X,Y,Z 51.6857900000000 176.905400000000 122.327200000000

 U,V,W -0.888056900000000 -0.445434000000000 -0.113769200000000

 USTEP,NREG 5960.73348525489 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 22 TXYZ: 1.00000001753137 ETEMF, ICHEMF: 8.1174060E-05 0

X,Y,Z -34.8704300000000 128.937400000000 348.859600000000

 U,V,W 0.333192900000000 -0.547496000000000 -0.767613600000000

 USTEP,NREG 8970.92524982406 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 46 TXYZ: 1.00000005098995 ETEMF, ICHEMF: 1.0000000E-04 0

X,Y,Z 11.7637000000000 155.453800000000 74.7174300000000

 U,V,W 0.442499400000000 -0.896752400000000 5.428175000000000E-03

USTEP,NREG 6191.67121009877 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 47 TXYZ: 1.00000002428595 ETEMF, ICHEMF: 8.1562130E-05 0

X,Y,Z 100.947300000000 -117.605600000000 326.406300000000

 U,V,W -0.773137300000000 0.633731500000000 2.535992000000000E-02

USTEP,NREG 10263.0788519149 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 48 TXYZ: 1.00000004210062 ETEMF, ICHEMF: 1.0000000E-04 0

X,Y,Z 11.0643000000000 139.018700000000 343.357500000000

 U,V,W 0.821131500000000 -0.387564600000000 -0.418971100000000

 USTEP,NREG 11109.9389887756 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 60 TXYZ: 1.00000004068161 ETEMF, ICHEMF: 4.0579370E-05 0

X,Y,Z -143.052100000000 -43.6912800000000 332.767400000000

 U,V,W -9.821579000000000E-02 0.589275500000000 -0.801940200000000

 USTEP,NREG 794.070453573721 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 73 TXYZ: 0.999999988519330 ETEMF, ICHEMF: 4.6905810E-05 0

X,Y,Z -195.567300000000 -13.8361200000000 239.522100000000

 U,V,W 0.281709300000000 0.892844000000000 -0.351382200000000

 USTEP,NREG 1030.04510610702 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 85 TXYZ: 0.999999964424330 ETEMF, ICHEMF: 8.1063840E-05 0

X,Y,Z 95.0013400000000 -38.6331500000000 371.703900000000

 U,V,W -0.263048100000000 -0.742140400000000 -0.616468400000000

 USTEP,NREG 760.476787027911 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 91 TXYZ: 1.00000001459732 ETEMF, ICHEMF: 7.3936660E-05 0

X,Y,Z 134.363700000000 147.832900000000 190.417500000000

 U,V,W -0.815988400000000 0.107829000000000 -0.567922400000000

 USTEP,NREG 3198.18152153661 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F T F F F

EMFGEO, NCOUN: 106 TXYZ: 0.999999987293090 ETEMF, ICHEMF: 5.6331240E-05 0

X,Y,Z -46.0188600000000 -114.599300000000 357.319000000000

 U,V,W 0.170981200000000 0.667360200000000 -0.724841900000000

 USTEP,NREG 3144.14818144011 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 121 TXYZ: 0.999999954274850 ETEMF, ICHEMF: 5.4100980E-05 0

X,Y,Z -110.219400000000 -140.546500000000 110.009100000000

 U,V,W 0.514078000000000 -0.372841900000000 0.772471800000000

 USTEP,NREG 746.812500117736 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 123 TXYZ: 0.999999974027802 ETEMF, ICHEMF: 8.3119700E-05 0

X,Y,Z 143.073900000000 -139.446200000000 190.802100000000

 U,V,W -0.576457500000000 0.811857500000000 -9.265055000000000E-02

USTEP,NREG 4258.90171438895 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F T F F F

EMFGEO, NCOUN: 155 TXYZ: 1.00000003084133 ETEMF, ICHEMF: 9.1891330E-05 0

X,Y,Z -138.555100000000 87.5020500000000 314.655500000000

 U,V,W 0.627243800000000 0.439949200000000 -0.642658500000000

 USTEP,NREG 622.548832032971 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 168 TXYZ: 0.999999875193840 ETEMF, ICHEMF: 3.3612900E-05 0

X,Y,Z -155.805400000000 -69.9509000000000 304.074700000000

 U,V,W 0.800899800000000 0.255211800000000 -0.541688400000000

 USTEP,NREG 5309.68317891873 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F T F F F

EMFGEO, NCOUN: 181 TXYZ: 1.00000004517609 ETEMF, ICHEMF: 9.8289070E-05 0

X,Y,Z 13.4972000000000 -83.6529000000000 18.8370400000000

 U,V,W -9.944928000000000E-02 0.991597700000000 8.272779000000000E-02

USTEP,NREG 3476.39032798103 2

LAGAIN,LSNSCT,LLDA,LSTNEW,LSENSE F F F F F

EMFGEO, NCOUN: 184 TXYZ: 1.00000004818643 ETEMF, ICHEMF: 4.1525260E-05 0

X,Y,Z 187.948600000000 -56.5727000000000 238.403700000000

 U,V,W -0.589292700000000 0.796493500000000 -0.135396700000000

 

 

 

Thanks and regards,

 

Riya Dey

From: Anna Ferrari [mailto:a.ferrari_at_hzdr.de]
Sent: 04 September 2020 16:03
To: riyadey_at_barc.gov.in; paola.sala_at_mi.infn.it; fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: [Fwd: Run time error if cycle number is greater than 2]

 

Dear Riya Dey,

 

  I see in your screenshot that you are generating an enormous .err file (about 2.8 GB)... with parallel 10 subdirectories you fill your disk space in few minutes (as I see). To understand your problem please look at this .err file.

 

 Kind regards,

  Anna

 

Am Fri, 04 Sep 2020 14:17:50 +0530 schrieb riyadey_at_barc.gov.in <mailto:riyadey_at_barc.gov.in> :

Dear Paola,

 

I modified the source.f where z cosines are read from WWLINE.

 

For 1000 histories, it is running smoothly with 10 spawns, but for more than 1E+7, a run time error is occurring and a message is appearing that there is no disk space available in root. I am using a system with 40 GB disk space. I am attaching one screenshot of the size of the output files. Is this problem occurring due to the source file ? Or is there any unnecessary big file generating during the simulation?

 

 

It will be very helpful if you can guide something about this issue.

 

Thanks and regards,

 

Riya Dey

 

 

 

 

On 2020-08-21 14:51, Paola Sala wrote:

 

Hello it is very difficult to debug user routines, especially when not all the elements are available (dada file). From the log file, you can see that the problem is indeed in your source.f routine. I suspect that something goes wrong in the calculation of the z cosine, resulting in a squareroot of negative value. This can come from rounding of the other two cosines, maybe written to disk with too few significant digits. Maybe it would be better to keep all the three of them as read from the data file, and apply a normalisation a posteriori on all three. But, of course, I am not sure that this is te reason. You can again print values (only for the cycle that crashes) and/or inspect the core file with the gdb linux utility.

 

Regards Paola

Dear FLUKA experts,

 

I am re-sending this mail if anyone can help me to solve the issue.

 

 

---------------------------- Original Message ----------------------------

Subject: Run time error if cycle number is greater than 2

From: riyadey_at_barc.gov.in <mailto:riyadey_at_barc.gov.in>

Date: Mon, August 17, 2020 2:51 pm

To: fluka-discuss_at_fluka.org <mailto:fluka-discuss_at_fluka.org>

Cc: paola.sala_at_mi.infn.it <mailto:paola.sala_at_mi.infn.it>

--------------------------------------------------------------------------

 

Dear Paola,

 

Thank you for the explanation. I have modified the source file and took

print out of 50 lines for checking and it is working properly. I have used

this source file with a fluka input containing voxel information. The

source file is compiled properly. But i am getting run time error for a

history of 1E+8.

 

I tested it for 1000 histories also.For 2 cycles, there is no error. If

the cycle is greater than 2, it is showing finished with errors.

 

I have added the modified source file also and two pngs where the list of

output files for 1st cycle and 3rd cycle are shown. Also the error files,

log files for 3rd cycle are attached.

 

Can you please guide me how to solve this ??

 

 

 

 

 

Hello First of all, It would be better to place the reading of your file inside the IF-block starting with "IF (LFIRST)" , so that it is read only once. As it is now, it is re-read at each primary.

 

Then, I do not know how the txt file is formatted, but I suppose that you have one line for each (x,y,z,E,wt,u,v,w) set? If yess, the read statement is not correct: at each READ, one record will be load , I mean one line, and only one quantity will be read in from this record, the first one, and assigned to A(i). Thus, the code will try to read 1946728 LINES, not 1946728 values. Again, if my assumption is correct, you can do two things: A loop on 243341 lines, reading 8 separate vectors (which will make things much easier) : read (30,*) x(i),y(i),z(i),E(i),wt(i),u(i),v(i),w(i) A single read, without loop : read(30,*)A I suggest in both cases to add prinout lines to check that everything is working correctly. regards Paola

Dear experts,

 

I am resending this mail since in the previous mail, there was one

attachment (fluence.txt) which has large mb size and thus I received

the

following instruction

Sorry, but your message was blocked by the system with the message:

"Message too long (>10000000 chars)"

 

I hope, with out this file, still my issue can be addressed.

 

 

 

 

 

---------------------------- Original Message

----------------------------

Subject: Error generated while reading x,y,z,energy,u,v,w from a file

using source.f

From: riyadey_at_barc.gov.in <mailto:riyadey_at_barc.gov.in>

Date: Fri, August 14, 2020 4:28 pm

To: fluka-discuss_at_fluka.org <mailto:fluka-discuss_at_fluka.org>

--------------------------------------------------------------------------

 

Dear experts,

 

I am trying to read position coordinates, energy, weight and

directional

cosines from a file (see the attached fluence.txt file). Total number

of

data are 1946728 and if we consider (x,y,z,E,wt,u,v,w) as a set, then

this

number will be 243341 (=1946728/8). I want to store these data in an

array

A and subsequently I have used this A to assign XFLK (NPFLKA) and other

variables.

 

Although there is no error in the compilation, the timed out error is

appearing and in the log file the following errors are generated.

 

At line 61 of file source1.f (unit = 30, file = '../fluence.txt')

Fortran runtime error: End of file

 

Error termination. Backtrace:

#0 0x7f1b456e6aed in ???

#1 0x7f1b456e7605 in ???

#2 0x7f1b456e82ab in ???

#3 0x7f1b459280a3 in ???

#4 0x7f1b45921316 in ???

#5 0x7f1b459222b1 in ???

#6 0x55f34f7cc910 in source_

 

Can you please guide me to solve this issue ?

 

 

Paola Sala INFN Milano tel. Milano +39-0250317374 tel. CERN +41-227679148

 

 

 

Paola Sala INFN Milano tel. Milano +39-0250317374 tel. CERN +41-227679148

__________________________________________________________________________ You can manage unsubscription from this mailing list athttps://www.fluka.org/fluka.php?id=acc_info

------------------------------------------------------------------------------------

Dr. Anna Ferrari

Institute of Radiation Physics, High Energy Density Division

Helmholtz-Zentrum Dresden-Rossendorf e.V.

Bautzner Landstraße 400, D-01328 Dresden

Tel. +49 351 260 2872

http://www.hzdr.de


 

Vorstand: Prof. Dr. Sebastian M. Schmidt, Dr. Ulrich Breuer

Vereinsregister: VR 1693 beim Amtsgericht Dresden



__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Sep 04 2020 - 14:20:32 CEST

This archive was generated by hypermail 2.3.0 : Fri Sep 04 2020 - 14:21:04 CEST