FLUKA: problem in running pemf

[Fan Lei <flei@scs.dera.gov.uk>]

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.