$ set verify$    !
$    !    This command file runs the program PERC.FOR.
$    !    It was used for: W. Bauer, Phys. Rev. C 38, 1297 (1988).
$    !
$ filenm = "test"
$ define for015 'filenm'.mom   ! = output file for moments event-by-event
$ define for016 'filenm'.his   ! = final mass and multiplicity distributions
$ define for017 'filenm'.sav   ! = intermediate distr. (can be deleted in end)
$ run perc_public
197,				TARGET MASS
1,				PROJECTILE MASS
0.82,				PBREAK
38383,				IRAND0
2,				IMPPAR (0->CONST, 1->WoodsSaxon, 2->Glauber)
1.0,				STOSS (IMP. PAR., IF IMPPAR = 0)
0.,7.,				B_UNTEN, B_OBEN (FM)
38,				ALPHA (IN DEGREE, IF IMPPAR = 0)  ***  ignore ***
5000,				NORUNS (# OF SIMULATED EVENTS)
30,				MCUT (CUTOFF MASS FOR MOMENTS)
1,197				MLOW, MHIGH (CUTOFFS MASSES FOR TOTAL MULT.)
1,				ifix (0 => fixed #of bonds, else fluctuating)
0,				io_mom (1 => write out moments, else not)
0,				io_sav (1 => intemeditae output, else not)
$ set noverify
$ !
$ !  Comments:
$ !  ---------
$ !
$ !  PBREAK     =   breaking probability for impact parameter 0.
$ !  IRAND0     =   seed for random number generator (RAN is used on a VAX)
$ !  IMPPAR     =   test flag for different dependencies of the breaking
$ !                 probability on impact parameter.  Usually the value
$ !                 2 is used for proton induced reactions.
$ !  STOSS      =   impact parameter (fm).  This is only used, if calculations
$ !                 for only one impact parameter are desired.  If IMPPAR = 1
$ !                 or 2, the impact parameter for every run is randomly 
$ !                 determined, and STOSS is ignored.
$ !  B_UNTEN,B_OBEN = window of impact parameters (fm) which are used for the
$ !                 output of distributions.  For inclusive reactions, choose:
$ !                 B_UNTEN = 0. and B_OBEN = Radius_targ + Radius_proj (fm).
$ !  ALPHA      =   specifies the reaction plane angle, and is only used
$ !                 with IMPPAR = 1
$ !  NORUNS     =   number of fragmentation events calculated.
$ !  MCUT       =   upper cutoff mass for the calculation of the moments
$ !                 of the mass distribution (should be < A_targ/3)
$ !  MLOW,MHIGH =   lower and upper mass bound for particles to be counted
$ !                 in the total multiplicity.  These two numbers can be
$ !                 chosen to compute the multiplicity of fragments in a given
$ !                 mass bin.  To compute the total multiplicity, set:
$ !                 MLOW = 1, MHIGH = A_targ
$ !  ifix       =   choice to use fixed number of broken bonds (fraction =
$ !                 breaking probability) (=> 0) or fluctuating
$ !  io_mom     =   if 1, then the logs of the moments of mass distribution
$ !                 are written out, else not
$ !  io_sav     =   if 1, then the mass distribution is saved after every 10%
$ !                 of the events have been completed, else not
$ !
$ ! Comment for heavy ion reactions:
$ !
$ ! To simulate thermalized heavy ion systems at a certain excitation energy,
$ ! choose as the target mass in the calculation the combined mass of target
$ ! and projectile in the heavy ion reaction.  The projectile mass should then
$ ! be set to 1, STOSS should be a large number (example: 100), and IMPPAR
$ ! should be set to 0.  BOBEN must in this case be larger than STOSS!