Dear Riya Dey,
in all the events reported with a warning in the portion of the the .err
file I can see, you can observe that the normalization of the cosines tx,
ty, tz, is not exactly 1 (not with the required precision: for example, for
the first warning event you have TXYZ: 1.00000004550694).
I gave a look at your source file and I saw that, to avoid a previous
error addressed by Paola, you read all three cosines from an external file.
Paola suggested you to add in this case a normalization "a posteriori"...
exactly to avoid the problems you are experiencing. This step is still
missing in your source:
* Cosines (tx,ty,tz)
TXFLK (NPFLKA) = UU(LINE)
TYFLK (NPFLKA) = VV(LINE)
TZFLK (NPFLKA) = WW(LINE)
Therefore: add the normalization by dividing by (UU**2 + VV**2 +ZZ**2)
I don't know if this is the only problem, but hope it helps
Regards,
Anna
Am Fri, 4 Sep 2020 16:25:21 +0530 schrieb "Riya Dey" :
>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:03To: riyadey_at_barc.gov.in; paola.sala_at_mi.infn.it; fluka-discuss_at_fluka.orgSubject: 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:
>>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
>>>>Date: Mon, August 17, 2020 2:51 pm
>>>>To: fluka-discuss_at_fluka.org
>>>>Cc: 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
>>>>>>Date: Fri, August 14, 2020 4:28 pm
>>>>>>To: 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
------------------------------------------------------------------------------------
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
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
- application/x-pkcs7-signature attachment: smime.p7s
Received on Fri Sep 04 2020 - 17:53:58 CEST