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