Hi again,
I implemented your suggestion. To be fair, I implemented the way I
understood it. But now fluka crashes. It crashes somewhere in the
initialization phase. I tried my input file without magnetic field and
it works. I also tried to comment out different sections of the magfld
routine but it makes no difference. Fluka crashes to the same point. I
scan the out file for signs of problems but I couldn't find any. It
crashes to the point where it tries to give a summary of the lattice
regions.
Any hint abount what is going on is very welcomed.
Thank you,
Makis
________________________________
From: owner-fluka-discuss_at_mi.infn.it on behalf of Santana, Mario
Sent: Fri 23/07/2010 19:20
To: Chrysostomos Valderanis; fluka-discuss_at_fluka.org
Subject: RE: lattice and magnetic field
Hi Makis,
What you are asking to do is possible and it is actually quite common.
You need the following instruction in your magfld.f routine
CALL GEOL2N(MLATTC,LATNAM,IRTLAT,Ierr)
And then a code that looks like in the example below:
[...]
IF (LATNAM.eq."QU01L ".or.LATNAM.eq."QU03L ".or.
$ LATNAM.eq."QU05L ".or.LATNAM.eq."QU07L ".or.
$ LATNAM.eq."QU09L ".or.LATNAM.eq."QU11L ".or.
$ LATNAM.eq."QU13L ".or.LATNAM.eq."QU15L ".or.
$ LATNAM.eq."QU17L ".or.LATNAM.eq."QU19L ".or.
$ LATNAM.eq."QU21L ".or.LATNAM.eq."QU23L ".or.
$ LATNAM.eq."QU25L ".or.LATNAM.eq."QU27L ".or.
$ LATNAM.eq."QU29L ".or.LATNAM.eq."QU31L ".or.
$ LATNAM.eq."QU33L ") THEN
ck =3D 3.85D-01 * 7.8/7.32 ! effective length=3D7.8 cm
ELSE IF (LATNAM.eq."QU02L ".or.LATNAM.eq."QU04L ".or.
$ LATNAM.eq."QU06L ".or.LATNAM.eq."QU08L ".or.
$ LATNAM.eq."QU10L ".or.LATNAM.eq."QU12L ".or.
$ LATNAM.eq."QU14L ".or.LATNAM.eq."QU16L ".or.
$ LATNAM.eq."QU18L ".or.LATNAM.eq."QU20L ".or.
$ LATNAM.eq."QU22L ".or.LATNAM.eq."QU24L ".or.
$ LATNAM.eq."QU26L ".or.LATNAM.eq."QU28L ".or.
$ LATNAM.eq."QU30L ".or.LATNAM.eq."QU32L ") THEN
ck =3D -3.85D-01 * 7.8/7.32 ! effective length=3D7.8 cm
END IF
[...]
Hope this helps,
Mario
This archive was generated by hypermail 2.2.0 : Sun Jul 25 2010 - 18:07:20 CEST