02.05.2014 Views

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

PENELOPE 2003 - OECD Nuclear Energy Agency

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

198 Chapter 6. Structure and operation of the code system<br />

NBTL ... number of bins for track-length distributions of primary particles. Limits<br />

of the interval where these distributions are tallied.<br />

-- Defaults: NBTL=100, TLMIN=0, TLMAX=2.5*RANGE(EPMAX,KPARP,1)<br />

ABSEN ... indicates a material M for which we require the code to tally the distribution<br />

of absorbed energy (up to three different materials can be selected, a separate line<br />

for each).<br />

-- Default: off<br />

DOSE2D ... the program will tally 2D, depth-radius, dose and deposited charge distributions<br />

in the body KL, KC (i.e. the cylinder KC of layer KL). The numbers NZ<br />

and NR of z- and r-bins have to be specified by the user, they must be between 1<br />

and 100. Up to three different bodies can be selected, a DOSE2D line for each body.<br />

-- Default: off<br />

Note: The labels KL, KC that identify a given cylinder are defined by the ordering<br />

in the input geometry list. These labels are written on the output geometry<br />

report.<br />

IFORCE ... activates forcing of interactions of type ICOL of particles KPAR in cylinder<br />

KC of layer KL. FORCE is the forcing factor and WLOW, WHIG are the limits of the<br />

weight window where interaction forcing is active.<br />

-- Default: no interaction forcing<br />

RESUME ... the program will read the dump file named filename1.ext (18 characters)<br />

and resume the simulation from the point where it was left. Use this option<br />

very, very carefully. Make sure that the input data file is fully consistent with the<br />

one used to generate the dump file.<br />

-- Default: off<br />

DUMPTO ... generate a dump file named filename2.ext (18 characters) after completing<br />

the simulation run. This allows resuming the simulation to improve statistics.<br />

-- Default: off<br />

NSIMSH ... desired number of simulated showers. Notice that NTOT is an INTEGER*4<br />

value and, hence, it cannot exceed 2 31 − 1.<br />

-- Default: NTOT=2147 million<br />

RSEED ... seeds of the random number generator.<br />

-- Default: ISEED1=12345; ISEED2=54321<br />

TIME ... allotted simulation time, in sec.<br />

-- Default: ITIME=100<br />

The program PENCYL (as well as the other example MAIN programs) is aborted when<br />

an incorrect input datum is found. The conflicting quantity usually appears in the last<br />

line of the output file. If the trouble is with arrays having dimensions smaller than<br />

required, the program indicates how the problem can be solved (this usually requires<br />

editing the source file, be careful).<br />

The example of input file given below belongs to the PENCYL file set included in the<br />

distribution package. It corresponds to the simulation of a narrow photon beam with

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!