I’m new to using FLUKA and FLAIR and I’ve been trying to just model a constant magnetic field using a modified magfld.f routine to compare with the result I get when I just change the MGNFIELD card to By=1.0 T, Bx=Bz=1.0 T (with the aim to model more complicated fields later). I’ve attached my input file and modified magfld.f routine (called magfld_dipole.f), the executable builds fine. I’m scoring dose in a water phantom after going through the constant magnet, and the result is exactly the same as when there is no magnetic field (despite the fact that it works when I just use MGNFIELD). Can anyone let me know what I’ve done wrong? I’ve downloaded the latest FLAIR and FLUKA from the fluka.cern website and I’m running it on Fedora on virtual box.
Dipole.flair (2.1 KB)
magfld_dipole.f (2.3 KB)