Re: FLUKA: problem in running pemf
You are right, you can only handle one material at a time.
>From the FLUKA manual, chapter 13:
"PEMF can only process one material at the time, thus producing only
single-material data sets. However, several such data sets can be appended
into a single file in order to be accessed by FLUKA"
In practice:
prepare two separate input files al.peg and u.peg, create two pemf files,
and concatenate them to make a single file:
cat u.pemf al.pemf > ual.pemf
Alberto
On Mon, 26 Mar 2001, Fan Lei wrote:
> Dear fluka team,
>
> I have enountered some problem in running the pemf on Linux machine. I
> created a simple
> "ual.peg" file in fluka/pemf
> -------------
> ELEM
> &INP
> IRAYL=1
> RHO=2.6989
> &END
> ALUMINUM AL
> AL
> ENER
> $INP AE=0.516,UE=100000001.,AP=0.001,UP=100000000. $END
> PWLF
> $INP NALE=300,NALG=400,NALR=100 $END
> DECK
> $INP $END
> ELEM
> &INP
> IRAYL=1
> RHO=18.95
> &END
> URANIUM U
> U
> ENER
> $INP AE=0.516,UE=100000001.,AP=0.001,UP=100000000. $END
> PWLF
> $INP NALE=300,NALG=400,NALR=100 $END
> DECK
> $INP $END
> ----------------------------------
>
> and tried
>
> % ./rpemf ual
>
> The running of the script seems to have completed ok and no ual.pemf file
> created. a inspection of the ual.lis file showes the fiest ELEM is treated
> ok but
> it ran into problem with the second element:
>
> ---------------------------
> [flei@holly:pemf]# tail -50 ual.lis
> 1
>
> ************************************************************
> * *
> * OPT = DECK *
> * *
> ************************************************************
>
>
> 1
>
> ************************************************************
> * *
> * OPT = ELEM *
> * *
> ************************************************************
>
>
> MEDIUM=URANIUM U ,STERNHEIMER ID=URANIUM
> U
> ATOMIC SYMBOLS ARE: U
>
>
>
> IN SUBROUTINE MIX:
>
>
> NUMBER OF ELEMENTS = 1, DENSITY= 18.9500 (GM/CC) (ACTUAL
> DENSITY=RHOx 1.00000 )
> 0 I Z(I) WA(I) PZ(I) RHOZ(I)
> INDEX PERIODIC ATOMIC PROPORTION PROPORTION
> NUMBER WEIGHT BY NUMBER BY WEIGHT
>
> 1 92.0000 238.030 1.00000 238.030
> Z VARIABLES--WM,ZC,ZT,ZA,ZB,ZAB,ZF,ZG,ZP
> 2.380300E+02 9.200000E+01 8.577018E+03 4.473559E+04 -1.292782E+04
> 3.180777E+04 3.387412E+03 -1.507263E+00
> ZV,ZU,ZS,ZE,ZX,RLC
> -1.902203E+00 -3.647037E-01 8.556000E+03 -2.579228E+04 7.858256E+03
> 3.166309E-01
> (I,XSI,ZZX,FZC,FCOUL,ZZ,I=1,NE)
> 1 1.228453E+00 8.577018E+03 4.507195E-01 3.949405E-01 8.556000E+03
> ZTILDE, ACCRIT, ZBREMS, FCOULZ, GMOLIE
> 9.100000E+01 1.526144E-04 9.261219E+01 3.949405E-01 3.739131E-02
> ZSQLPM, ZLNLPM, ES1LPM, E0LPM
> 8.577018E+03 3.708488E+00 6.009640E-04 1.946491E+07
> EDEN= 4.4107857E+24
>
> BLCCRA = 1.50540
>
> BLCC,XCC,TEFF0,XR0= 2.83115E+04 5.81745E+00 7.06561E-04 3.28030E-01
>
>
> STOPPED IN SPINIT: INCORRECT USER-OVERRIDE OF SSB-DATA
> [flei@holly:pemf]#
>
> ------------------------------------------------------
>
> I have switched the order of U and AL in the peg file and the result was the
> same.
>
> It seems "rpemf" cannot handle more than one element in the .peg file. Is
> this right? Or I
> must have done something wrong in the "ual.peg" file!
>
>
> Best Regards.
>
> Fan
>
> -------------------------------------------------------------------
> Dr. Fan Lei
> Space Department, DERA
> G081, A8 Building
> Farnborough, GU14 0LX
> U.K.
>
> Tel . 44-1252-393904
> Fax. 44-1252-396330
> email: flei@scs.dera.gov.uk
>
>
>
>
> The information contained in this e-mail and any subsequent correspondence
>is private and is intended solely for the intended recipient(s). For those
>other than the intended recipient any disclosure, copying, distribution, or
>any action taken or omitted to be taken in reliance on such information is
>prohibited and may be unlawful.
>
--
Alberto Fasso'
CERN-EP, CH-1211 Geneve 23 (Switzerland)
Fax: (41 22) 767 9480
Alberto.Fasso@cern.ch
Partial thread listing: