Dear Ryan,
the fusrbv routine is not called because your USRBIN is wrong, in
particular the value 100 you input as step [NOT number of bins] both for
the region and lattice sets makes no sense, since this way you selected
only the region #0 (that does not exist), in the lattice #0 (i.e. outside
any possible replica, which is indeed your case since you have no
replica). The only element where you can expect a not null dose is region
#3 for lattice #0, so your scoring intervals have to be defined
accordingly.
In order to get a lateral dose profile, I'd rather adopt the built-in
cylindrical USRBIN around the beam axis, where each bin is characterized
by a (r,z) coordinate pair that can be then associated to a polar angle.
Best wishes
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. +41 22 7678962
fax +41 22 7668854
On Wed, 27 Dec 2017, Ryan Brosch wrote:
> 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 ToolsHelp
>
> *...+....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 Fri Jan 05 2018 - 19:03:59 CET