From: h jalali (h_bjalali@yahoo.com)
Date: Thu Jan 04 2007 - 09:13:21 CET
Dear FLUKA users
I want to calculate energy deposition due to a carbon source with distribution of energy in the Al shell.Runnig my input was crashed by this error
================================ Running FLUKA for cycle # 1 =============================================
/home/hb/fluka/flutil/rfluka: line 311: 3770 Aborted (core dumped) ${EXE} <$INPF 2>$LOGF >$LOGF
I do not know what the problem is.Please guide me.
My source.f and input are attached.
Best Regards
Hamide
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
GLOBAL
TITLE
calculation of energy deposition in Al sell
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
BEAM -100.E+00 HEAVYION
SOURCE
*
GEOBEGIN COMBINAT
A cylindrical AL target inside vacuum
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
SPH 1 0.0 0.0 0.0 1.E06
SPH 2 0.0 0.0 0.0 8.E05
RCC 3 0.0 0.0 0.0 0.0 0.0 1000.0
400.0
RCC 4 0.0 0.0 4.0 0.0 0.0 992.0
398.0
END
* black hole
BH1 5 +1 -2
* vacuum around
VA2 5 +2 -3
* Al target(outer layer)
AL3 5 +3 -4
* Al target(inner layer)
AL4 5 +4
END
GEOEND
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 13.0 26.981 2.699 10.0 ALUMINUM
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Al target, Shell
ASSIGNMAT 10.0 3.0
* Inside of shell
ASSIGNMAT 2.0 4.0
* External Black Hole
ASSIGNMAT 1.0 1.0
* Vacuum
ASSIGNMAT 2.0 2.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* score in each region energy deposition by primaries
SCORE 208.0
*
USRBDX 99.0 208.0 -47.0 3 4 5.0E8 FluenUD
USRBDX 100.0 +50.0 10.0 10.0 &
* R-PHI-Z binning of the fluence inside and around the target
USRBIN 17.0 208.0 -48.0 400.0 1000.0 Edeposit
USRBIN 398.0 992.0 30.0 10.0 30.0 &
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
RANDOMIZE 1.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
START 100000.0
STOP
*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
SUBROUTINE SOURCE ( NOMORE )
INCLUDE '(DBLPRC)'
INCLUDE '(DIMPAR)'
INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2006 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* New source for FLUKA9x-FLUKA200x: *
* *
* Created on 07 january 1990 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 03-mar-06 by Alfredo Ferrari *
* *
* This is just an example of a possible user written source routine. *
* note that the beam card still has some meaning - in the scoring the *
* maximum momentum used in deciding the binning is taken from the *
* beam momentum. Other beam card parameters are obsolete. *
* *
*----------------------------------------------------------------------*
* Version for isotropic Heavy ion source
*
DIMENSION NORENE(0:13), ENELOW(13), ENEUP(13)
*----------------------------------------------------------------------*
INCLUDE '(BEAMCM)'
INCLUDE '(FHEAVY)'
INCLUDE '(FLKSTK)'
INCLUDE '(IOIOCM)'
INCLUDE '(LTCLCM)'
INCLUDE '(PAPROP)'
INCLUDE '(SOURCM)'
INCLUDE '(SUMCOU)'
*
LOGICAL LFIRST
*
SAVE LFIRST
DATA LFIRST / .TRUE. /
*======================================================================*
* *
* BASIC VERSION *
* *
*======================================================================*
* Heavy ion energy group
DATA ENELOW /
& 0.10162E-02, 0.11349E-02, 0.44499E-02, 0.96386E-02, 0.31153E-01,
& 0.10069E+00, 0.30375E+00, 0.98177E+00, 0.31732E+01, 0.95724E+01,
& 0.30939E+02, 0.70819E+02, 0.98629E+03/
DATA ENEUP /
& 0.10304E-02, 0.11506E-02, 0.45118E-02, 0.97726E-02, 0.31586E-01,
& 0.10209E+00, 0.30797E+00, 0.99541E+00, 0.32137E+01, 0.97054E+01,
& 0.33610E+02, 0.71803E+02, 0.10000E+03/
* Normalized spectrum
DATA NORENE / 0.D0,
& 0.14766E-02, 0.35102E-02, 0.23169E-01, 0.26918E-01, 0.35186E-01,
& 0.78632E-01, 0.25713E+00, 0.57673E+00, 0.86827E+00, 0.97527E+00,
& 0.99680E+00, 0.99952E+00, 0.10000E+01/
* +-------------------------------------------------------------------*
NOMORE = 0
* +-------------------------------------------------------------------*
* | First call initializations:
IF ( LFIRST ) THEN
* | *** The following 3 cards are mandatory ***
TKESUM = ZERZER
LFIRST = .FALSE.
LUSSRC = .TRUE.
* | *** User initialization ***
write(lunout,*)
write(lunout,*) '--->modified source.f for Carbon heavy ion'
END IF
* |
* +-------------------------------------------------------------------*
* Push one source particle to the stack. Note that you could as well
* push many but this way we reserve a maximum amount of space in the
* stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
* NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
* WTFLK (NPFLKA) = ONEONE
* WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
* +-------------------------------------------------------------------*
* Sample the energy group
XI = FLRNDM(DUMMY)
DO 500 K = 1 ,13
IF (XI .LE. NORENE(K)) THEN
ENERGY = (ENEUP(K) - ENELOW(K))/2
GO TO 600
END IF
500 CONTINUE
STOP ' Failed to sample the energy group '
600 CONTINUE
count=count+1
NPFLKA = NPFLKA + 1
WTFLK (NPFLKA) = ONEONE
WEIPRI = WEIPRI + WTFLK (NPFLKA)
* +--------------------------------------------------------------------*
* Decode Heavy ion particle ID for 100*A+Z FLUKA number
IDpart = -1206
IF (IDpart .LT. -402) THEN
IJBEAM = -2
IDpart = abs (IDpart)
IPROA = IDpart / 100
IPROZ = IDpart - 100 * IPROA
END IF
* +--------------------------------------------------------------------*
* | (Radioactive) isotope:
IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
IARES = IPROA
IZRES = IPROZ
IISRES = IPROM
CALL STISBM ( IARES, IZRES, IISRES )
IJHION = IPROZ * 1000 + IPROA
IJHION = IJHION * 100 + KXHEAV
IONID = IJHION
CALL DCDION ( IONID )
CALL SETION ( IONID )
* |
* +-------------------------------------------------------------------*
* | Heavy ion:
ELSE IF ( IJBEAM .EQ. -2 ) THEN
IJHION = IPROZ * 1000 + IPROA
IJHION = IJHION * 100 + KXHEAV
IONID = IJHION
CALL DCDION ( IONID )
CALL SETION ( IONID )
ILOFLK (NPFLKA) = IJHION
* | Flag this is prompt radiation
LRADDC (NPFLKA) = .FALSE.
* |
* +-------------------------------------------------------------------*
* | Normal hadron:
ELSE
IONID = IJBEAM
ILOFLK (NPFLKA) = IJBEAM
* | Flag this is prompt radiation
* LRADDC (NPFLKA) = .FALSE.
END IF
* |
* +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
LOFLK (NPFLKA) = 1
* User dependent flag:
LOUSE (NPFLKA) = 0
* User dependent spare variables:
DO 100 ISPR = 1, MKBMX1
SPAREK (ISPR,NPFLKA) = ZERZER
100 CONTINUE
* User dependent spare flags:
DO 200 ISPR = 1, MKBMX2
ISPARK (ISPR,NPFLKA) = 0
200 CONTINUE
* Save the track number of the stack particle:
ISPARK (MKBMX2,NPFLKA) = NPFLKA
NPARMA = NPARMA + 1
NUMPAR (NPFLKA) = NPARMA
NEVENT (NPFLKA) = 0
DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
AGESTK (NPFLKA) = +ZERZER
AKNSHR (NPFLKA) = -TWOTWO
* Group number for "low" energy neutrons, set to 0 anyway
IGROUP (NPFLKA) = 0
* Kinetic energy of the particle (GeV)
TKEFLK (NPFLKA) = ENERGY
* TKEFLK (NPFLKA) = SQRT ( PBEAM**2 + AM (IONID)**2 ) - AM (IONID)
* Particle momentum
* PMOFLK (NPFLKA) = PBEAM
PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
* + TWOTWO * AM (ILOFLK(NPFLKA)) ) )
* +------------------------------------------------------------------------*
* Isotropic angular distribution
YI = FLRNDM(DUMMY)
COSTH = ONEONE - 2.D0 * YI
IF (COSTH .GT. ONEONE) THEN
COSTH = ONEONE
ELSE IF (COSTH .LT. -ONEONE) THEN
COSTH = -ONEONE
END IF
ZI = FLRNDM(DUMMY)
ANG = TWOPIP * ZI
UBEAM = SQRT (ONEONE - COSTH * COSTH) * COS(ANG)
VBEAM = SQRT (ONEONE - COSTH * COSTH) * SIN(ANG)
WBEAM = COSTH
* +-----------------------------------------------------------------------*
* Cosines (tx,ty,tz)
TXFLK (NPFLKA) = UBEAM
TYFLK (NPFLKA) = VBEAM
TZFLK (NPFLKA) = WBEAM
* TZFLK (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
* & - TYFLK (NPFLKA)**2 )
* Polarization cosines:
TXPOL (NPFLKA) = -TWOTWO
TYPOL (NPFLKA) = +ZERZER
TZPOL (NPFLKA) = +ZERZER
* Particle coordinates
XFLK (NPFLKA) = XBEAM
YFLK (NPFLKA) = YBEAM
ZFLK (NPFLKA) = ZBEAM
XBEAM = 4E05
YBEAM = 4E05
ZBEAM = 4E05
* Calculate the total kinetic energy of the primaries: don't change
IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 )
& THEN
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
& * WTFLK (NPFLKA)
ELSE
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
END IF
RADDLY (NPFLKA) = ZERZER
* Here we ask for the region number of the hitting point.
* NREG (NPFLKA) = ...
* The following line makes the starting region search much more
* robust if particles are starting very close to a boundary:
CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
& NRGFLK(NPFLKA), IDISC )
* Do not change these cards:
CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
NLATTC (NPFLKA) = MLATTC
CMPATH (NPFLKA) = ZERZER
CALL SOEVSV
RETURN
*=== End of subroutine Source =========================================*
END
This archive was generated by hypermail 2.1.6 : Sun Jan 07 2007 - 11:24:17 CET