[fluka-discuss]: Trouble running fusrbv.f

From: Ryan Brosch <rbrosch_at_asu.edu>
Date: Wed, 27 Dec 2017 17:02:57 +0700

Dear Fluka experts,

I would like to score energy deposited per unit volume of a proton beam in
a water tank, binned laterally by angle. For this I am trying to use the
USRBIN 18.0 setting with fusrbv.f

While I can compile and link fusrbv.f, the routine is not executed by FLUKA
as it runs with the following script.

./fff fusrbv.f
./lfluka -o myfluka -m fluka fusrbv.o
$FLUPRO/flutil/rfluka -N0 -M1 angle -e myfluka

Right now I am using the default fusrbv.f, modified only with a print
statement to let me know it is working, as specified in the manual. The
statement does not appear, which tells me the routine is not being called.
I am using 32bit g77 (version 3.4.6 20060404 (Red Hat 3.4.6-41.fc26))
in 64bit Fedora 27 with FLUKA installed from rpm.

I would appreciate suggestions about either a different way to score
lateral dose by total scattering angle using built in scoring (USRBIN does
not appear to score this way by default) or ideas why my user routine is
not running. I've read on the CERN website that there are some known issues
linking with gfortran and user routines, but I'm not sure this applies to
me since I definitely have FLUFOR=g77 set.

I've attached my input cards. Thanks for your time and happy new year.


File Edit Options Buffers Tools
Help

*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
Secondary proton energy deposition in a water phantom
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
BEAM -0.2288 PROTON
BEAMPOS 0.0 0.0 0.0
*
GEOBEGIN
COMBNAME
  0 0 Water target inside vacuum
RPP body1 -500000.0 +5000000.0 -500000.0 +500000.0 -500000.0 +500000.0
RPP body2 -100000.0 +1000000.0 -100000.0 +100000.0 -100.0 +100000.0
RPP body3 -200.0 +200.0 -200.0 +200.0 0.0 +200.0
END
* black hole
regBH1 5 +body1 -body2
* vacuum around
regVA2 5 +body2 -body3
regWA 5 +body3
END
GEOEND
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*MATERIAL
WATER
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Water target
ASSIGNMAT WATER regWA
* 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 KeV
EMFCUT -0.000010 0.000010 1.0 WATER
PROD-CUT
IONFLUCT -1.0 -1.0 WATER WATER
DISCARD 7 8
THRESHOLd 1.0 1.0
USRBIN 18.0 228.0 +65.0 10.0 10.0 50.
hadrons
USRBIN 0.0 100.0 100.0 500.0 &
STEPSIZE 0.1 0.1 regWA regWA
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
RANDOMIZE 1.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
START 10000
STOP



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

Received on Wed Dec 27 2017 - 12:49:36 CET

This archive was generated by hypermail 2.3.0 : Wed Dec 27 2017 - 12:49:55 CET