Yes, you are right. I thought that I had said it sometimes in one of the ma=
ny emails this thread has developed in.
If I haven't, I apologize
Alberto
________________________________________
From: Nicholas Bolibruch [nboli_at_acanac.net]
Sent: Monday, May 31, 2010 1:58 PM
To: Fasso, Alberto
Subject: Re: Synchrotron radiation
Dear Alberto,
In regards to the Source.f file you posted...
I don't think you mentioned that the SOURCE card needs WHAT(1) set to 1
for the OPEN statement to work for io stream 81. I spent enough time
today trying to debug this problem, was even running it in GDB trying to
figure out where I was getting a Floating Exception from ...
-- Nick Bolibruch Controls & Instrumentation Development Canadian Light Source Inc. Alberto Fasso' wrote: > Dear Emre, > > you have not specified what "some errors" are. But your BEAM card and you= r > source file are assuming that your maximum photon energy is 2 MeV (0.002 = GeV). > The spectrum in the sr.spectrum file is in MeV, which are translated to G= eV > inside the source routine: > DO 2 I =3D 1, 630 > READ(81,*,END=3D999) EN(I), SPEC(I) > EN(I) =3D EN(I)*1.D-3 <---- > 2 CONTINUE > On the other hand, the spectrum in the sr.spectrum file has a maximum of = 1995 > MeV (1.995 GeV), which is very unlikely. And the sampling is done between > 1.D-6 (1 keV, corresponding to the first line of the spectrum) and 1.9952= 6D-4 > (about 200 keV): > ENERGY =3D 1.D-6 + ( 1.99526D-4 - 1.D-6) * RAN01 > > To conclude: > 1) check your spectrum file: I guess it should go from 1 keV to 200 keV, = and > not from 1 keV to 2 GeV. > 2) Your BEAM card is not wrong, but if the maximum is really 200 keV, set= ting > it to 2 MeV would make your cross section tabulations less precise > without need > > Kind regards, > > Alberto > > On Mon, 31 May 2010, Emre EROGLU wrote: > > >> Dear Alberto, >> >> Many thanks for your mail. But still i have some errors. I didn't find w= here is it. I'm >> sending my .inp file, sr_source.f (modified from bl1source.f) and sr.spe= ctrum (generated >> from srplot.f). And also are there any mistakes in my BEAM card? >> Thanks. >> >> Best Regards, >> Emre >> >> >> >> On Thu, 27 May 2010 06:35:42 -0700 (PDT), Alberto Fasso' wrote >> >>> I apologize. On April 28, answering a question by Emre Eroglu, I sent >>> an example of source routine sampling from a synchrotron radiation spec= trum >>> provided as an external file. See: >>> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/2879.html >>> >>> I did not realize at that time that it was an old routine, incompatible= with >>> the latest FLUKA versions. I send here in attachment a corrected versio= n. >>> >>> Sorry! >>> >>> Alberto >>> >>> On Wed, 26 May 2010, Emre EROGLU wrote: >>> >>> >>>> Dear Alberto, >>>> >>>> I try to compile bl1source.f but i have some errors. what can i do for= that? >>>> >>>> Regards, >>>> Emre >>>> >>>> >>>> -e g77 -o bl1source.o -c -O3 -g -mcpu=3Dpentiumpro -fexpensive-optimiz= ations >>>> -funroll-loops -fstrength-reduce -Wall -fno-automatic -fno-silent -m32 >>>> -ffortran-bounds-check -I/home/eeroglu/fluka/flukapro bl1source.f >>>> `-mcpu=3D' is deprecated. Use `-mtune=3D' or '-march=3D' instead. >>>> source: >>>> bl1source.f: In subroutine `source': >>>> bl1source.f:11: >>>> INCLUDE '(AACOLL)' >>>> ^ >>>> Unable to open INCLUDE file `(AACOLL)' at (^) >>>> bl1source.f:12: >>>> INCLUDE '(BEAM)' >>>> ^ >>>> Unable to open INCLUDE file `(BEAM)' at (^) >>>> bl1source.f:13: >>>> INCLUDE '(EPISOR)' >>>> ^ >>>> Unable to open INCLUDE file `(EPISOR)' at (^) >>>> bl1source.f:17: >>>> INCLUDE '(STACK)' >>>> ^ >>>> Unable to open INCLUDE file `(STACK)' at (^) >>>> bl1source.f:18: >>>> INCLUDE '(STARS)' >>>> ^ >>>> Unable to open INCLUDE file `(STARS)' at (^) >>>> bl1source.f:39: >>>> LUSSRC =3D .TRUE. >>>> 1 2 >>>> Type disagreement between expressions at (1) and (2) >>>> bl1source.f:78: >>>> WT (LSTACK) =3D WEIGHT >>>> ^ >>>> Invalid declaration of or reference to symbol `wt' at (^) [initially s= een at (^)] >>>> bl1source.f:90: >>>> ILO (LSTACK) =3D IJHION >>>> ^ >>>> Invalid declaration of or reference to symbol `ilo' at (^) [initially = seen at (^)] >>>> bl1source.f:102: >>>> LO (LSTACK) =3D 1 >>>> ^ >>>> Invalid declaration of or reference to symbol `lo' at (^) [initially s= een at (^)] >>>> bl1source.f:104: >>>> LOUSE (LSTACK) =3D 0 >>>> ^ >>>> Invalid declaration of or reference to symbol `louse' at (^) [initiall= y seen at >>>> (^)] >>>> bl1source.f:107: >>>> SPAREK (ISPR,LSTACK) =3D ZERZER >>>> ^ >>>> Invalid declaration of or reference to symbol `sparek' at (^) [initial= ly seen at (^)] >>>> bl1source.f:111: >>>> ISPARK (ISPR,LSTACK) =3D 0 >>>> ^ >>>> Invalid declaration of or reference to symbol `ispark' at (^) [initial= ly seen at (^)] >>>> bl1source.f:116: >>>> NUMPAR (LSTACK) =3D NPARMA >>>> ^ >>>> Invalid declaration of or reference to symbol `numpar' at (^) [initial= ly seen at (^)] >>>> bl1source.f:117: >>>> NEVENT (LSTACK) =3D 0 >>>> ^ >>>> Invalid declaration of or reference to symbol `nevent' at (^) [initial= ly seen at (^)] >>>> bl1source.f:118: >>>> DFNEAR (LSTACK) =3D +ZERZER >>>> ^ >>>> Invalid declaration of or reference to symbol `dfnear' at (^) [initial= ly seen at (^)] >>>> bl1source.f:121: >>>> AGESTK (LSTACK) =3D +ZERZER >>>> ^ >>>> Invalid declaration of or reference to symbol `agestk' at (^) [initial= ly seen at (^)] >>>> bl1source.f:122: >>>> AKNSHR (LSTACK) =3D -TWOTWO >>>> ^ >>>> Invalid declaration of or reference to symbol `aknshr' at (^) [initial= ly seen at (^)] >>>> bl1source.f:124: >>>> IGROUP (LSTACK) =3D 0 >>>> ^ >>>> Invalid declaration of or reference to symbol `igroup' at (^) [initial= ly seen at (^)] >>>> bl1source.f:126: >>>> TKE(LSTACK) =3D ENERGY >>>> ^ >>>> Invalid declaration of or reference to symbol `tke' at (^) [initially = seen at (^)] >>>> bl1source.f:128: >>>> PMOM (LSTACK) =3D ENERGY >>>> ^ >>>> Invalid declaration of or reference to symbol `pmom' at (^) [initially= seen at (^)] >>>> bl1source.f:132: >>>> TX (LSTACK) =3D TINX >>>> ^ >>>> Invalid declaration of or reference to symbol `tx' at (^) [initially s= een at (^) >>>> ] >>>> bl1source.f:133: >>>> TY (LSTACK) =3D TINY >>>> ^ >>>> Invalid declaration of or reference to symbol `ty' at (^) [initially s= een at (^)] >>>> bl1source.f:134: >>>> TZ (LSTACK) =3D TINZ >>>> ^ >>>> Invalid declaration of or reference to symbol `tz' at (^) [initially s= een at (^)] >>>> bl1source.f:137: >>>> TXPOL (LSTACK) =3D -TWOTWO >>>> ^ >>>> Invalid declaration of or reference to symbol `txpol' at (^) [initiall= y seen at (^)] >>>> bl1source.f:138: >>>> TYPOL (LSTACK) =3D +ZERZER >>>> ^ >>>> Invalid declaration of or reference to symbol `typol' at (^) [initiall= y seen at (^)] >>>> bl1source.f:139: >>>> TZPOL (LSTACK) =3D +ZERZER >>>> ^ >>>> Invalid declaration of or reference to symbol `tzpol' at (^) [initiall= y seen at (^)] >>>> bl1source.f:141: >>>> XA (LSTACK) =3D XINA >>>> ^ >>>> Invalid declaration of or reference to symbol `xa' at (^) [initially s= een at (^)] >>>> bl1source.f:142: >>>> YA (LSTACK) =3D YINA >>>> ^ >>>> Invalid declaration of or reference to symbol `ya' at (^) [initially s= een at (^)] >>>> bl1source.f:143: >>>> ZA (LSTACK) =3D ZINA >>>> ^ >>>> Invalid declaration of or reference to symbol `za' at (^) [initially s= een at (^)] >>>> bl1source.f:163: >>>> NLATTC (LSTACK) =3D MLATTC >>>> ^ >>>> Invalid declaration of or reference to symbol `nlattc' at (^) [initial= ly seen at (^)] >>>> In file included from bl1source.f:11: >>>> In file included from bl1source.f:17: >>>> bl1source.f:160: warning: unused variable 'idisc' >>>> bl1source.f:143: warning: unused variable 'zina' >>>> bl1source.f:142: warning: unused variable 'yina' >>>> bl1source.f:141: warning: unused variable 'xina' >>>> bl1source.f:134: warning: unused variable 'tinz' >>>> bl1source.f:133: warning: unused variable 'tiny' >>>> bl1source.f:132: warning: unused variable 'tinx' >>>> bl1source.f:79: warning: unused variable 'weipri' >>>> bl1source.f:39: warning: unused variable 'lussrc' >>>> bl1source.f:21: warning: unused variable 'line' >>>> >>>> >>>> >>>> >>>> >>>> On Wed, 28 Apr 2010 01:42:18 -0700 (PDT), Alberto Fasso' wrote >>>> >>>>> Dear Emre, >>>>> >>>>> as you know, FLUKA for the time being cannot simulate the production = of >>>>> synchrotron radiation. A source routine, therefore, must sample from = a >>>>> spectrum produced by a separate program (see >>>>> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1183.html= ) >>>>> I attach here an example, where the spectrum is provided as an >>>>> external file >>>>> >>>>> Alberto >>>>> >>>>> On Tue, 27 Apr 2010, emreeroglu_at_uludag.edu.tr wrote: >>>>> >>>>> >>>>>> Hi all, >>>>>> >>>>>> Does anyone have source.f which is modified for synchrotron radiatio= n? >>>>>> >>>>>> Thanks in advance. >>>>>> >>>>>> Regards, >>>>>> Emre EROGLU >>>>>> >>>> Emre EROGLU >>>> Uludag University >>>> High Energy Physcics >>>> >>>> >>>> >>>> >>> -- >>> Alberto Fasso` >>> SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025 >>> Phone: (1 650) 926 4762 Fax: (1 650) 926 3569 >>> fasso_at_slac.stanford.edu >>> >> Emre EROGLU >> Uludag University >> High Energy Physcics >> >> >> > >=Received on Tue Jun 01 2010 - 10:59:09 CEST
This archive was generated by hypermail 2.2.0 : Tue Jun 01 2010 - 10:59:10 CEST