RE: lattice and magnetic field

From: Stefan Roesler <sroesler_at_mail.cern.ch>
Date: Wed, 4 Aug 2010 14:37:28 +0200

Hi Makis

GEON2L has four arguments, not three. Blindly adding a fourth integer
variable removes the crash. I append below the explanation of the
parameters of GEON2L and GEON2R.

Cheers
Stefan

       ENTRY GEON2L ( LATNAM, NLATT, IRTLAT, IERR )
*
*----------------------------------------------------------------------*
* *
* GEOmetry: Name 2 Lattice *
* *
* Version november 89 by Alfredo Ferrari *
* INFN - Milan *
* *
* Last change on 05-mar-06 by Alfredo Ferrari *
* *
* *
* This is a subroutine of Fluka200x for combinatorial geometries *
* *
* *
* Input variable: *
* Latnam = lattice name *
* *
* Output variables: *
* Nlatt = lattice number *
* Irtlat = # of the (possibile) rotation associated *
* Ierr = error code (0 on success, 1 on failure) *
* *

       ENTRY GEON2R ( REGNAM, NREG, IERR )
*
*----------------------------------------------------------------------*
* *
* GEOmetry: Name 2 Region *
* *
* Version november 89 by Alfredo Ferrari *
* INFN - Milan *
* *
* Last change on 12-mar-04 by Alfredo Ferrari *
* *
* *
* This is a subroutine of Fluka200x for combinatorial geometries *
* *
* *
* Input variable: *
* Regnam = region name *
* *
* Output variables: *
* Nreg = region number *
* Ierr = error code (0 on success, 1 on failure) *
* *

On Wed, 4 Aug 2010, Chrysostomos Valderanis wrote:

> They were in the previous mail. I attach them again.
>
> Thank you,
> Makis
>
>
> ________________________________
>
> From: Stefan Roesler
> Sent: Tue 03/08/2010 22:35
> To: Chrysostomos Valderanis
> Cc: fluka-discuss_at_fluka.org
> Subject: RE: lattice and magnetic field
>
>
>
> Hi Makis
>
> Please post input and user routines so that we can check.
>
> Cheers
> Stefan
>
>
> On Tue, 3 Aug 2010, Chrysostomos Valderanis wrote:
>
> > So. just to report my investigation on the problem. It seems that calls
> > to GEON2L do not work as I expected. If I remove all the GEON2L calls
> > and hardcode their outcome into the magfld.f file it works, at least
> > partially. If I include even one call to GEON2L the program crashes.
> > I really would like to avoid hardcoding all the numbers in the routine.
> > I some on could explain to me why GEON2L doesn't work the way I expect
> > it to work I would be grateful.
> >
> > Thank you,
> > Makis
> >
> > ________________________________
> >
> > From: owner-fluka-discuss_at_mi.infn.it on behalf of Chrysostomos
> > Valderanis
> > Sent: Sat 24/07/2010 16:28
> > To: msantana_at_slac.stanford.edu; fluka-discuss_at_fluka.org
> > Subject: RE: lattice and magnetic field
> >
> >
> > 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
> >
> >
> >
> >
> > I am trying to build a geometry with lattice and magnetic fields and I
> > am f=
> > acing some problems.
> >
> > First question. My reading of the online courses says that the magnetic
> > field should only be set for the prototype region. Since I have more
> > than one magnetic field, I am using a routine (magfld.f). But what
> > should be written in this file?
> > -Should I find the prototype regions with GEON2R and set the magnetic
> > fields for them?
> > -Should I find the lattice regions with GEON2R and set the magnetic
> > field for the prototype regions they correspond to?
> > -Should I find the lattice regions and set the magnetic field as it
> > should be in the replica area?
> > For the moment I am using option2 but it doesn't seem to work. Anyway I
> > wou, ld prefer to calculate
> > the magnetic field just for the prototypes as it is a much simpler
> > expression.
> >
> > Second question. I really would like to be able to set two replicas of
> > the same prototype to different magnetic fields. Is this possible? How?
> >
> > Thank you in advance,
> > Makis
> >
> >
>
> --
> ___________________________________
>
> Stefan Roesler
> CERN, DGS-RP
> CH-1211 Geneva 23
> Switzerland
>
> Phone: +41-22-7679891
> Fax: +41-22-7669639
> E-mail: Stefan.Roesler_at_cern.ch
>
>
>

-- 
___________________________________
Stefan Roesler
CERN, DGS-RP
CH-1211 Geneva 23
Switzerland
Phone:  +41-22-7679891
Fax:    +41-22-7669639
E-mail: Stefan.Roesler_at_cern.ch
Received on Thu Aug 05 2010 - 21:04:44 CEST

This archive was generated by hypermail 2.2.0 : Thu Aug 05 2010 - 21:04:44 CEST