Dear Alessia,
first a general comment. You say that you are a Fluka beginner and that
you run some basic examples.
A beginner should first learn the basics, and mgdraw is not part of the
basics. It needs to be used only in complex situations that a beginner is
not aware of. Even in the Fluka courses for beginners, students are not
taught how to use mgdraw and other user routines. This stuff is covered
only in the advanced courses.
By the way, Fluka has been designed so that it can be used even without
any knowledge of programming, except in very exceptional cases. If you
insist on trying the user routines of course you can, but you are on your
own: you should not expect much help from this discussion list.
Having said that, in this particular case it is clear where the error
comes from.
I invite you to have a look at the manual. The whole description of mgdraw
can be found in 3 places: 1) the description of command USERDUMP, 2) Chap. 11,
3) Chap 13, sec. 13.2.13.
Even looking only at the description of USERDUMP, you read:
WHAT(2) : if the default version of MGDRAW is used, number of the
unformatted output unit.
Values of WHAT(2) < 21 must be avoided because of possible
conflicts with FLUKA pre-defined units.
and you have set WHAT(2) = 2, which is < 21:
USERDUMP 101.0 2.0 dump
(in fact, 2 is the unit number for the output of the random number
generator).
By the way, why did you write "dump" at the position of WHAT(5)? Not very
important, since according to the manual, WHAT(5...6) = not used. But you
must have had something in mind...
Alberto
On Wed, 13 May 2015, Saggio Alessia wrote:
> Hi fluka users,
> I'm a Fluka beginner and I'm trying to run some basic examples.
>
> I've downloaded the example.inp from fluka website, and when I run it with
>
> $FLUPRO/flutil/rfluka -N0 -M1 example.inp
>
> all works fine.
>
> Then I modified the mgdraw.f routine, adding a few lines of code (I guess I should use the card USERDUMP in my example.inp). Now I type:
>
> $FLUPRO/flutil/fff mgdraw_modified.f
> $FLUPRO/flutil/lfluka -o myfluka -m fluka mgdraw_modified.o
>
> (and up to here all is fine)
>
> $FLUPRO/flutil/rfluka -e myfluka -N0 -M1 example.inp
>
> and I get:
> ======================= Running FLUKA for cycle # 1 =======================
>
> Removing links
> Removing temporary files
> Saving output and random number seed
> No ranexample002 generated!
> Saving additional files generated
> No additional files to move
> End of FLUKA run
>
>
> In practise, I get this error only when adding the USERDUMP card in the .inp file!
> Do you have any idea?
>
> Cheers, Alessia
>
>
>
>
> Here is the example.inp (the same as the original, I've only added the line with USERDUMP):
>
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> TITLE
> Charged pion fluence inside and around a proton-irradiated Be target
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BEAM 50.E+00 PROTON
> BEAMPOS 0.0 0.0 -50.0
> *
> GEOBEGIN COMBNAME
> 0 0 A simple Be target inside vacuum
> RPP body1 -5000000.0 +5000000.0 -5000000.0 +5000000.0 -5000000.0 +5000000.0
> RPP body2 -1000000.0 +1000000.0 -1000000.0 +1000000.0 -100.0 +1000000.0
> RPP body3 -10.0 +10.0 -10.0 +10.0 0.0 +5.0
> * plane to separate the upstream and downstream part of the target
> XYP body4 2.5
> END
> * black hole
> regBH1 5 +body1 -body2
> * vacuum around
> regVA2 5 +body2 -body3
> * Be target 1st half
> regBE3 5 +body3 +body4
> * Be target 2nd half
> regBE4 5 +body3 -body4
> END
> GEOEND
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> MATERIAL 4.0 9.0122 1.848 5.0 BERYLLIU
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * Be target, 1st and 2nd half
> ASSIGNMAT BERYLLIU regBE3 regBE4
> * External Black Hole
> ASSIGNMAT BLCKHOLE regBH1
> * Vacuum
> ASSIGNMAT VACUUM regVA2
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * e+e- and gamma production threshold set at 10 MeV
> EMFCUT -0.010 0.010 1.0 BERYLLIU PROD-CUT
> * score in each region energy deposition and stars produced by primaries
> SCORE ENERGY BEAMPART
> USERDUMP 101.0 2.0 dump
> * Boundary crossing fluence in the middle of the target (log intervals, one-way)
> USRBDX 99.0 PIONS+- -47.0 regBE3 regBE4 400. piFluenUD
> USRBDX +50.0 +50.0 0.0 10.0 &
> * Boundary crossing current in the middle of the target (log intervals, one-way)
> USRBDX -1.0 PIONS+- -47.0 regBE3 regBE4 400. piCurrUD
> USRBDX +50.0 +50.0 0.0 10.0 &
> * Tracklength fluence inside the target, Upstream part and Downstream part
> * Logarithmic energy intervals
> USRTRACK -1.0 PIONS+- -48.0 regBE3 1000.0 20. piFluenU
> USRTRACK 50.0 0.001 &
> USRTRACK -1.0 PIONS+- -49.0 regBE4 1000.0 20. piFluenD
> USRTRACK 50.0 0.001 &
> * Cartesian binning of the pion fluence inside and around the target
> USRBIN 10.0 PIONS+- -50.0 50.0 50.0 50. piFluBin
> USRBIN -50.0 -50.0 -10.0 100.0 100.0 60.0 &
> * Cartesian binning of the deposited energy inside the target
> USRBIN 10.0 ENERGY -51.0 10.0 10.0 5. Edeposit
> USRBIN -10.0 -10.0 0.0 20.0 20.0 5.0 &
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> RANDOMIZE 1.0
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> START 1000.0
> STOP
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
>
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Wed May 13 2015 - 15:56:03 CEST