RE: lattice and magnetic field

From: Chrysostomos Valderanis <Chrysostomos.Valderanis_at_cern.ch>
Date: Sat, 24 Jul 2010 16:28:13 +0200

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

Received on Sun Jul 25 2010 - 18:07:20 CEST

This archive was generated by hypermail 2.2.0 : Sun Jul 25 2010 - 18:07:20 CEST