Re: Re: [fluka-discuss]: problems of two-steps method
Dear Andrea,
Thanks for your reply. I ran it as your suggestions, unfortunately, it has some errors " STOP TOO MANY ERRORS IN GEOMETRY: STOP", but I see no geometry errors in FLAIR, so I am very puzzled with the error, the geometry is not complex, so where are errors from?
Thanks again.
Best regards
Yang
> -----原始邮件-----
> 发件人: "Andrea Fontana" <andrea.fontana_at_pv.infn.it>
> 发送时间: 2018-11-02 19:16:20 (星期五)
> 收件人: "YANG Tao" <yangt_at_ihep.ac.cn>, fluka-discuss_at_fluka.org
> 抄送:
> 主题: Re: [fluka-discuss]: problems of two-steps method
>
> Dear Yang,
> I ran your example and I have found a couple of FORTRAN errors
> in the file I/O instructions:
>
> - in file 1.f in the OPEN instruction, you have put 'NEW' as status,
> but this will create a new file for each WRITE statement (the
> function is called at *each* boundary crossing): simply write
> 'UNKNOWN' or do not specify the status.
> With this correction, at step 1 you write the Data1.dat file
> (named 1001_Data1.dat with Fluka naming conventions).
>
> - in file source.f, the problem is due to the fact that UNIT=88
> should be opened in the 2.inp file with the Fluka card OPEN
> and not directly in the FORTRAN code.
> Moreover, you read 20000 lines, but your file has only 10000:
> therefore you hit the end of file and this generates the EOF error.
> To be on the safe side, I have put a protection in the READ
> statement as:
>
> READ (88,102,END=101)
> ......
> 101 CONTINUE
>
> so that the code knows how to handle this situation.
> At this point, step 2 also works, but the I have seen tracking
> errors, seemingly due to geometry problems...
>
> I leave it to you now, but I am available for further interactions
> if needed.
>
> In attachment, a revised version of your files.
>
> Kind regards,
> Andrea
>
> Il 02/11/2018 08:16, YANG Tao ha scritto:
> >
> > Hi, everyone!
> >
> >
> > Recently, I simulate using so called two-steps method to simulate the
> > activation problem. File 1.input and 1.f use FLUSCW.F +USRBDX to
> > calculate the generated particles, and 2.f use the 1st step results as
> > the source. But when I run the 2nd step, I get an error（seen in
> > 2.log）, "At line 73 of file source.f (unit = 88, file = '../try.dat')
> > Fortran runtime error: End of file", I try many ways but cannot solve
> > it. could anyone give a suggestion, any suggestion is appreciate!
> >
> > Best regards!
> >
> > Yang
> >
>
> --
> ========================================================================
> Dr. Andrea Fontana tel: +39 0382 987991
> Istituto Nazionale fax: +39 0382 423241
> di Fisica Nucleare
> Sezione di Pavia e-mail: andrea.fontana_at_pv.infn.it
> Via Bassi 6 web : www.pv.infn.it/~fontana
> 27100 PAVIA, Italy
> ========================================================================
>
(Ƨځ맲r*'~&x%Zm 0nߖi'iw
Received on Sat Nov 03 2018 - 05:51:15 CET
This archive was generated by hypermail 2.3.0
: Sat Nov 03 2018 - 05:51:23 CET