# RE: lattice and magnetic field

From: Chrysostomos Valderanis <Chrysostomos.Valderanis_at_cern.ch>
Date: Wed, 4 Aug 2010 15:01:35 +0200

-----_=_NextPart_001_01CB33D5.2B5DA7EF
Content-Type: text/plain; charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

Hi Stefan,

To support myself about that blatant omission I have to say that when
I google for this routine the link I ended up is the =
http://www.fluka.org/content/course/NEA/lectures/Lattice.pdf
which clearly sugests in page 10 that GEON2L gets three arguments.
The same is true for the presentation at
http://www.fluka.org/content/course/NEA/lectures/UserRoutines.pdf
where in page 14 there is a suggestion that similar routines exist for
the lattice geometry.

Thank you again,
Makis
________________________________

From: Stefan Roesler
Sent: Wed 04/08/2010 14:37
To: Chrysostomos Valderanis
Cc: fluka-discuss_at_fluka.org
Subject: RE: lattice and magnetic field

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 =3D lattice name =
*
* *
* Output variables: *
* Nlatt =3D lattice number =
*
* Irtlat =3D # of the (possibile) rotation associated =
*
* Ierr =3D 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 =3D region name =
*
* *
* Output variables: *
* Nreg =3D region number =
*
* Ierr =3D 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.
>=20
> Thank you,
> Makis
>=20
>
> ________________________________
>
> 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 =3D3D 3.85D-01 * 7.8/7.32 ! effective length=3D3D7.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 =3D3D -3.85D-01 * 7.8/7.32 ! effective length=3D3D7.8 =
cm
> > END IF
> >
> > [...]
> >
> > Hope this helps,
> >
> > Mario
> >
> >
> >
> >
> > I am trying to build a geometry with lattice and magnetic fields and =
I
> > am f=3D
> > 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
------_=_NextPart_001_01CB33D5.2B5DA7EF
Content-Type: text/html; charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<HTML dir=3Dltr><HEAD><TITLE>RE: lattice and magnetic field</TITLE>=0A=
<META content=3D"text/html; charset=3Dunicode" http-equiv=3DContent-Type>=0A=
<BODY>=0A=
<DIV dir=3Dltr><FONT color=3D#000000 size=3D2 face=3DArial>Hi =
Stefan,</FONT></DIV>=0A=
<DIV dir=3Dltr><FONT size=3D2 face=3DArial></FONT>&nbsp;</DIV>=0A=
<DIV dir=3Dltr><FONT size=3D2 face=3DArial>Thank you very much for your =
<DIV dir=3Dltr><FONT size=3D2 face=3DArial>To support myself about that =
blatant omission I have to&nbsp;say that&nbsp;&nbsp;when I google for =
this routine the link I ended up is the&nbsp;<A =
href=3D"http://www.fluka.org/content/course/NEA/lectures/Lattice.pdf">htt=
p://www.fluka.org/content/course/NEA/lectures/Lattice.pdf</A>  =
&nbsp;</FONT>which clearly sugests in page 10 that GEON2L gets three =
arguments.<BR>The same is true for the presentation at<A =
href=3D"http://www.fluka.org/content/course/NEA/lectures/UserRoutines.pdf=
">http://www.fluka.org/content/course/NEA/lectures/UserRoutines.pdf</A>&n=
bsp; where in page 14 there is a suggestion that similar routines exist =
for the lattice geometry.</DIV>=0A=
<DIV dir=3Dltr>&nbsp;</DIV>=0A=
<DIV dir=3Dltr>&nbsp;</DIV>=0A=
<DIV dir=3Dltr>Thank you again,</DIV>=0A=
<DIV dir=3Dltr>Makis</DIV>=0A=
<DIV dir=3Dltr>&nbsp;</DIV>=0A=
<DIV dir=3Dltr>=0A=
<HR tabIndex=3D-1>=0A=
</DIV>=0A=
<DIV dir=3Dltr><FONT size=3D2 face=3DTahoma><B>From:</B>  Stefan =
Roesler<BR><B>Sent:</B>  Wed 04/08/2010 14:37<BR><B>To:</B>  Chrysostomos =
Valderanis<BR><B>Cc:</B>  fluka-discuss_at_fluka.org<BR><B>Subject:</B>  RE: =
lattice and magnetic field<BR></FONT><BR></DIV></DIV>=0A=
<DIV>=0A=
<P><FONT size=3D2>Hi Makis<BR><BR>GEON2L has four arguments, not three. =
Blindly adding a fourth integer<BR>variable removes the crash. I append =
below the explanation of the<BR>parameters of GEON2L and =
GEON2R.<BR><BR>Cheers<BR>Stefan<BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
ENTRY GEON2L ( LATNAM, NLATT, IRTLAT, IERR =
)<BR>*<BR>*--------------------------------------------------------------=
--------*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; GEOmetry: Name 2 =
Lattice&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Version november =
89&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; by&nbsp;&nbsp; Alfredo =
Ferrari&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp; INFN - =
Milan&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Last change on 05-mar-06&nbsp; =
by&nbsp;&nbsp; Alfredo =
Ferrari&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; This is a subroutine of Fluka200x for =
combinatorial geometries&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Input =
variable:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp; *<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Latnam&nbsp;&nbsp; =3D lattice =
name&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Output =
variables:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Nlatt&nbsp;&nbsp;&nbsp; =3D lattice =
number&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp; *<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Irtlat&nbsp;&nbsp; =
=3D # of the (possibile) rotation =
associated&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Ierr&nbsp;&nbsp;&nbsp;&nbsp; =3D error code (0 on success, 1 on =
failure)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR><BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ENTRY GEON2R ( REGNAM, NREG, =
IERR =
)<BR>*<BR>*--------------------------------------------------------------=
--------*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; GEOmetry: Name 2 =
Region&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Version november =
89&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; by&nbsp;&nbsp; Alfredo =
Ferrari&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp; INFN - =
Milan&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Last change on 12-mar-04&nbsp; =
by&nbsp;&nbsp; Alfredo =
Ferrari&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; This is a subroutine of Fluka200x for =
combinatorial geometries&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Input =
variable:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp; *<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Regnam&nbsp;&nbsp; =3D region =
name&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp; Output =
variables:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Nreg&nbsp;&nbsp;&nbsp;&nbsp; =3D region =
number&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp; *<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
Ierr&nbsp;&nbsp;&nbsp;&nbsp; =3D error code (0 on success, 1 on =
failure)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
; =
*<BR>*&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
p;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
*<BR><BR><BR><BR>On Wed, 4 Aug 2010, Chrysostomos Valderanis =
wrote:<BR><BR>&gt; They were in the previous mail. I attach them =
again.<BR>&gt;&nbsp;<BR>&gt; Thank you,<BR>&gt; =
Makis<BR>&gt;&nbsp;<BR>&gt;<BR>&gt; =
________________________________<BR>&gt;<BR>&gt; From: Stefan =
Roesler<BR>&gt; Sent: Tue 03/08/2010 22:35<BR>&gt; To: Chrysostomos =
Valderanis<BR>&gt; Cc: fluka-discuss_at_fluka.org<BR>&gt; Subject: RE: =
lattice and magnetic field<BR>&gt;<BR>&gt;<BR>&gt;<BR>&gt; Hi =
Makis<BR>&gt;<BR>&gt; Please post input and user routines so that we can =
check.<BR>&gt;<BR>&gt; Cheers<BR>&gt; Stefan<BR>&gt;<BR>&gt;<BR>&gt; On =
Tue, 3 Aug 2010, Chrysostomos Valderanis wrote:<BR>&gt;<BR>&gt;&gt; So. =
just to report my investigation on the problem. It seems that =
calls<BR>&gt;&gt; to GEON2L do not work as I expected. If I remove all =
the GEON2L calls<BR>&gt;&gt; and hardcode their outcome into the =
magfld.f file it works, at least<BR>&gt;&gt; partially. If I include =
even one call to GEON2L the program crashes.<BR>&gt;&gt; I really would =
like to avoid hardcoding all the numbers in the routine.<BR>&gt;&gt; I =
some on could explain to me why GEON2L doesn't work the way I =
expect<BR>&gt;&gt; it to work I would be grateful.<BR>&gt;&gt;<BR>&gt; =
&gt; Thank you,<BR>&gt;&gt; Makis<BR>&gt;&gt;<BR>&gt;&gt; =
________________________________<BR>&gt;&gt;<BR>&gt;&gt; From: =
owner-fluka-discuss_at_mi.infn.it on behalf of Chrysostomos<BR>&gt;&gt; =
Valderanis<BR>&gt;&gt; Sent: Sat 24/07/2010 16:28<BR>&gt;&gt; To: =
msantana_at_slac.stanford.edu; fluka-discuss_at_fluka.org<BR>&gt;&gt; =
Subject: RE: lattice and magnetic field<BR>&gt;&gt;<BR>&gt; =
&gt;<BR>&gt;&gt; Hi again,<BR>&gt;&gt;<BR>&gt;&gt; I implemented your =
suggestion. To be fair, I implemented the way I<BR>&gt;&gt; understood =
it. But now fluka crashes. It crashes somewhere in the<BR>&gt;&gt; =
initialization phase. I tried my input file without magnetic field =
and<BR>&gt;&gt; it works. I also tried to comment out different =
sections of the magfld<BR>&gt;&gt; routine but it makes no difference. =
Fluka crashes to the same point. I<BR>&gt;&gt; scan the out file for =
signs of problems but I couldn't find any. It<BR>&gt;&gt; crashes to =
the point where it tries to give a summary of the lattice<BR>&gt;&gt; =
regions.<BR>&gt;&gt; Any hint abount what is going on is very =
welcomed.<BR>&gt;&gt;<BR>&gt;&gt; Thank you,<BR>&gt;&gt; =
Makis<BR>&gt;&gt;<BR>&gt;&gt;<BR>&gt;&gt; =
________________________________<BR>&gt;&gt;<BR>&gt;&gt; From: =
owner-fluka-discuss_at_mi.infn.it on behalf of Santana, Mario<BR>&gt;&gt; =
Sent: Fri 23/07/2010 19:20<BR>&gt;&gt; To: Chrysostomos Valderanis; =
fluka-discuss_at_fluka.org<BR>&gt;&gt; Subject: RE: lattice and magnetic =
field<BR>&gt;&gt;<BR>&gt;&gt;<BR>&gt;&gt;<BR>&gt;&gt; Hi =
Makis,<BR>&gt;&gt;<BR>&gt;&gt; What you are asking to do is possible =
and it is actually quite common.<BR>&gt;&gt;<BR>&gt;&gt; You need the =
following instruction in your magfld.f routine<BR>&gt;&gt;<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CALL =
GEOL2N(MLATTC,LATNAM,IRTLAT,Ierr)<BR>&gt;&gt;<BR>&gt;&gt; And then a =
code that looks like in the example below:<BR>&gt;&gt;<BR>&gt;&gt; =
[...]<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
IF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (LATNAM.eq."QU01L&nbsp; =
".or.LATNAM.eq."QU03L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU05L&nbsp; ".or.LATNAM.eq."QU07L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU09L&nbsp; ".or.LATNAM.eq."QU11L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU13L&nbsp; ".or.LATNAM.eq."QU15L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU17L&nbsp; ".or.LATNAM.eq."QU19L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU21L&nbsp; ".or.LATNAM.eq."QU23L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU25L&nbsp; ".or.LATNAM.eq."QU27L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU29L&nbsp; ".or.LATNAM.eq."QU31L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU33L&nbsp; ") THEN<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
ck =3D3D&nbsp; 3.85D-01 * 7.8/7.32&nbsp; ! effective length=3D3D7.8 =
cm<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ELSE IF =
(LATNAM.eq."QU02L&nbsp; ".or.LATNAM.eq."QU04L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU06L&nbsp; ".or.LATNAM.eq."QU08L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU10L&nbsp; ".or.LATNAM.eq."QU12L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU14L&nbsp; ".or.LATNAM.eq."QU16L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU18L&nbsp; ".or.LATNAM.eq."QU20L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU22L&nbsp; ".or.LATNAM.eq."QU24L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU26L&nbsp; ".or.LATNAM.eq."QU28L&nbsp; ".or.<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
\$&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
LATNAM.eq."QU30L&nbsp; ".or.LATNAM.eq."QU32L&nbsp; ") THEN<BR>&gt; =
&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
ck =3D3D -3.85D-01 * 7.8/7.32&nbsp; ! effective length=3D3D7.8 =
cm<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; END =
IF<BR>&gt;&gt;<BR>&gt;&gt; [...]<BR>&gt;&gt;<BR>&gt;&gt; Hope this =
helps,<BR>&gt;&gt;<BR>&gt;&gt; Mario<BR>&gt;&gt;<BR>&gt;&gt;<BR>&gt; =
&gt;<BR>&gt;&gt;<BR>&gt;&gt; I am trying to build a geometry with =
lattice and magnetic fields and I<BR>&gt;&gt; am f=3D<BR>&gt;&gt; =
acing some problems.<BR>&gt;&gt;<BR>&gt;&gt; First question. My =
reading of the online courses says that the magnetic<BR>&gt;&gt; field =
should only be set for the prototype region. Since I have more<BR>&gt; =
&gt; than one magnetic field, I am using a routine (magfld.f). But =
what<BR>&gt;&gt; should be written in this file?<BR>&gt;&gt; -Should I =
find the prototype regions with GEON2R and set the magnetic<BR>&gt;&gt; =
fields for them?<BR>&gt;&gt; -Should I find the lattice regions with =
GEON2R&nbsp; and set the magnetic<BR>&gt;&gt; field for the prototype =
regions they correspond to?<BR>&gt;&gt; -Should I find the lattice =
regions and set the magnetic field as it<BR>&gt;&gt; should be in the =
replica area?<BR>&gt;&gt; For the moment I am using option2 but it =
doesn't seem to work. Anyway I<BR>&gt;&gt; wou, ld prefer to =
calculate<BR>&gt;&gt; the magnetic field just for the prototypes as it =
is a much simpler<BR>&gt;&gt; expression.<BR>&gt;&gt;<BR>&gt;&gt; =
Second question. I really would like to be able to set two replicas =
of<BR>&gt;&gt; the same prototype to different magnetic fields. Is this =
possible? How?<BR>&gt;&gt;<BR>&gt;&gt; Thank you in advance,<BR>&gt; =
&gt; Makis<BR>&gt;&gt;<BR>&gt;&gt;<BR>&gt;<BR>&gt; --<BR>&gt; =
___________________________________<BR>&gt;<BR>&gt; Stefan =
Roesler<BR>&gt; CERN, DGS-RP<BR>&gt; CH-1211 Geneva 23<BR>&gt; =
Switzerland<BR>&gt;<BR>&gt; Phone:&nbsp; +41-22-7679891<BR>&gt; =
Fax:&nbsp;&nbsp;&nbsp; +41-22-7669639<BR>&gt; E-mail: =
Stefan.Roesler_at_cern.ch<BR>&gt;<BR>&gt;<BR>&gt;<BR><BR>--<BR>_____________=
______________________<BR><BR>Stefan Roesler<BR>CERN, DGS-RP<BR>CH-1211 =
Geneva 23<BR>Switzerland<BR><BR>Phone:&nbsp; =
+41-22-7679891<BR>Fax:&nbsp;&nbsp;&nbsp; +41-22-7669639<BR>E-mail: =
Stefan.Roesler_at_cern.ch<BR></FONT></P></DIV></BODY></HTML>
------_=_NextPart_001_01CB33D5.2B5DA7EF--
```
Received on Thu Aug 05 2010 - 20:05:55 CEST

This archive was generated by hypermail 2.2.0 : Thu Aug 05 2010 - 20:06:25 CEST