Input deck summary
------------------

First field given for each card is the format, the rest are the input
variable names.

Card 1 : 80A1 ITEXT
Card 2 : * IZP IAP IZT IAT ELAB UXOFET
Card 3 : * JP JT IP12 JCN
If JCN < 0
    Card 3A: * widwid
Card 4 : * KOPTK IPS1 IPSMAX NNX NPX NAX
Card 5 : * IZE(4) IAE(4) JE(4) IZFF DAF FFB SIERK VSADOFF
Card 6 : * CL0 DIFF SIGMCN
If sigmcn<0
    Card 6A: * TLGAUSSC TLGAUSSW tlgaussoff
Card 7 : * FTHETA KOPTLD
Card 8 : * DALDM UTR ULDM UJTR UJLDM KOPTLQ KOPTEB
Card 9 : * R0LDM DEF DEFS
If R0LDM < 0 then
    R0LDM = 0
    Card 9A: * JNOD OFFEV OFFODD TRAMO
Card 10: * XYE1 XYM1 XYE2  CJG1 CJG2 XYENH GMIN
if XYE1<0 then
    XYE1 = 0
    Card 10A: * DIPOLE DIPOLW DIPOLS
if XYENH<0 then
    Must have declared R0LDM<0
    XYENH=-XYENH
    Card 10B: * E2BANDPAR
Card 11: * WGR CGR CVCBE CVCB VK CVK
Card 12: * EXR0 EXH CJC LCO
Card 13: * KOUTW KOUTL KEVAP KGAMMA i2dpop
IF KOPTK = 1 then
    Card 13A: * ((IPSZO(IZF,IPS) IZF=1,4) IPS=IH1,IH2)
Cards 14: 2I5,4F10.3 IZ IA DA DELTA FTHETA CK
repeated until IZ = 0
Cards 15:
    Cards 15A: 2I5,F10.3 IZZ IAA EX
    if IZZ != 0 then
        Card 15A: 7(F7.3,I2,A1) (EXL(I) JL(I) PL(I) I=1,7)
        repeated until PL(I) is blank or null
    repeated until IZZ = 0

CARD 1    Title
---------------

ITEXT	Any alphanumeric title of 80 or less characters.

CARD 2    Entrance channel
--------------------------

IZP     Atomic number of projectile.

IAP     Mass number of projectile.

IZT     Atomic number of target nucleus.

IAT     Mass number of target nucleus.

ELAB    Beam energy in MeV. This is interpreted as the excitation energy
	of the compound nucleus EXCN if JCN>=0 (see card 3) and if IAP = 0.

UXOFET	???

CARD 3    Spins and parity in entrance channel
----------------------------------------------

JP      Spin of projectile (twice the value for an odd nucleus).

JT      Spin of target (twice the value for an odd nucleus).

IP12    Product of intrinsic parities of target and projectile.

        = 0  For even parity in the entrance channel.

        = 1  For odd parity in the entrance channel.

        = 2  For identical 0+ nuclei.

        < 0  For unknown parity in the entrance channel.

JCN
        < 0 : Then the population in the CN taken to be a gaussian
	centred on -JCN with a width given in card 3A.

	= 0 : Distribution calculated according to card 6.

        > 0 : Then the CN population is assumed to be in one spin JCN with a
	population cross section of 1000 mb. The spin is JCN+1/2 for odd
	mass. The excitation energy EXCN = ELAB (card 2) if IAP = 0. This
	option is used for example, if one wants to study the influence of
	the angular momentum on the CN on the decay mode. The parameters on
	card 6 are then irrelevant but should be entered nevertheless.


CARD 3A   Width of JCN distribution
-----------------------------------

This card only read if JCN < 0 on card 3.

widwid	Width of gaussian distribution of JCN

CARD 4    Structure of the decay cascade
----------------------------------------

This card specifies the structure of the so-called decay cascade, which
contains all nuclei considered as possible decay products and defines their
sequence in the various decay chains. The parameters are usually computed
internally if all values are entered as zeros.

KOPTK   Number of decay steps (for any  combination  of decays).

        = 1 The cascade structure is read in from card 13A (check subroutine
	KASKAD).

        = 0 Internal computation (= EXCN/12, but >=3).

        > 1 Internal computation with number of steps given. (check cross
	sections "last step above threshold").

IPS1    Number of decaying nuclei (only for KOPTK < 0).

IPSMAX  Number of nuclei in the cascade ( <=500 ) (only for KOPTK < 0).

NNX     Maximum excess of neutron minus proton emissions.

NPX     Same for proton emission.

NAX     Maximum number of alpha emissions.

The last three parameters may be used to trim the wings of the decay cascade
in case of high excitation.

CARD 5    4th particle decay mode, fission
------------------------------------------

IZE(4)  Atomic number of the 4th particle to be taken into account in
	addition to n-, p- and alpha emission.  Make sure that the
	corresponding transmission coefficients are available.  No
	additional decay mode if entered zero.

IAE(4)  Mass number of the 4th particle.

JE(4)   Spin of the 4th particle (twice the value for odd mass).

IPE(4)  Parity of the 4th particle.
                = 0  even.
                = 1  odd.
                Else unknown.

EXC4	Excitation energy of 4th decay particle (since 4th decay particle
	may come out in an excited state). In this case, JE(4) would also be
	the spin of the excited state of the 4th decay particle.

IZFF            = 0  No fission competition.
                <>0  With fission.

DAF	af = A/DAF = level density parameter at the saddle point (default
	value: DALDM).

FFB	Fraction of the liquid-drop fission barrier (default value: 1.0).

SIERK	=0 Use standard liquid drop fission barriers. If iz <= 19 then
	these are always used even if SIERK != 0.
	
	!= 0 Use Sierk fission barriers (recommended 

VSADOFF	This parameter is ignored, it's forced to 0 in the code

CARD 6    Angular momentum distribution in the compound nucleus
---------------------------------------------------------------

CL0	Maximum angular momentum (a very approximate value is calculated
	internally if entered zero).

DIFF   	Diffuseness (default value: 2 hbar).  NOTE: The transmission
	coefficient is given by the expression
                  TL = 1/[1+Exp{(L-CL0)/DIFF}]
        and the total fusion cross section by
                  SIGMCN = Sum over L of SIGML
        where
                  SIGML = (2L+1)*TL.

SIGMCN	Instead of CL0 the total fusion cross section (in mb) may be
	inserted.

CARD 6A   Transmission coefficient spread
-----------------------------------------

This card is only read if sigmcn < 0, and is ignored in any event (the code
has been commented out).

TLGAUSSC

TLGAUSSW

tlgaussoff

CARD 7    Level density parameters at low excitation (EX<=UTR)
--------------------------------------------------------------

FTHETA	Fraction of the rigid-sphere moment of inertia used as effective
	value for calculating the yrast line at low excitation energy
	(default value: 0.85).

KOPTLD	Option for the level density parameters DA, DELTA (and CK).

        = 0 Interpolation of values of Dilg et al., Nucl.Phys.A217, for the
	mass range A < 45.

        =2 Same for A > 45.

CARD 8    Level density parameters at high excitation (liquid drop region,
	  EX >= ULDM) ----------------------------------------------------
----------------------

DALDM	Level density parameter constant for ALDM = A/DALDM (1/MeV),
	(default value: 8 MeV)

UTR, ULDM
        Interpolation range between the low-energy region (<=UTR) and the
	liquid-drop region (>=ULDM).  Valid for DA and DELTA (default
	values: 40/A**(1/3) MeV and 80/A**(1/3) MeV, respectively).

UJTR , UJLDM
        Same for the moment of inertia, (default values: UTR and ULDM,
	respectively).

KOPTLQ  Level density parameter option.

        = 0 DALDM as above, DELTALDM corresponding to the liquid-drop g.s.

        = 1 Same parameters DA and DELTA as at low excitation.

KOPTEB  Option for liquid-drop mass formula.
                = 0  Myers-Swiatecki Lysekil.
                = 1  Myers droplet (recommended).
                = 2  Same with Wigner term.
                = 3  Groote, Hilf, Takahashi.
                = 4  Same with Wigner term
                = 5  Seeger.
                = 6  Same with Wigner term.

CARD 9    Yrast line
--------------------

R0LDM  	< 0: Use parameters in card 9A

	> 0: Radius parameter (in fm) for calculating the moment of inertia
	and the spin cutoff parameter (default value: rms-radius from Myers,
	Nucl.Phys. A204(1973)465).

DEF, DEFS
        Constants used for parameterizing the yrast line of a rotating
	liquid drop:
                 EROT = I(I+1)*hbar**2/(2J),
                 J = J0(1+DEF*L**2+DEFS*L**4),
        (default values are calculated internally to fit the Cohen, Plasil,
	Swiatecki results).

CARD 9A   Yrast line, extra parameters
--------------------------------------

JNOD	???

OFFEV	???

OFFODD	???

TRAMO	???

CARD 10   Gamma-decay widths
----------------------------

XYE1    E1 gamma-width in Weisskopf units (default value: .0001).

        NOTE: if entered negative an extra card 10a is read for GDR
	parameters.

XYM1	M1 gamma-width in Weisskopf units (default value: .03).

XYE2	E2 gamma-width in Weisskopf units (default value: 5).

	NOTE: Some versions of Cascade read in GQR parameters if XYE2<0,
	this version has no support for the GQR.

CJG1	Lower edge for spin transition range for E2 enhancement.

CJG2	Upper edge for spin transition range for E2 enhancement.

XYENH	Enhanced E2 strength in W.u. above CJG2.

GMIN    Minimum gamma-decay, GMIN*(XYE1+XYM1+XYE2). This is important only
	for lower edge of the population matrix. (default value: 1.E-6).

CARD 10a  GDR parameters
------------------------

This card only read in if XYE1 < 0 on card 10

DIPOLE	???

DIPOLW	???

DIPOLS	???

Card 10b  ???
-------------

This card only read in if XYENH < 0 on card 10, and is only meaningful if
R0LDM < 0 on card 9.

E2BANDPAR
	???

CARD 11   Cut-offs
------------------

WGR     Populations below the given value in any decaying nucleus are
	ignored, (default value: 0.003 mb/MeV.hbar). Check lost cross
	sections.

CGR     Cutoff for the decay intensity of a 1MeV.1hbar population
	element in a specific channel (default value: 0.25*WGR).

        The following four parameters are used to control the upper energy
	limit of the 64MeV.64hbar population matrices.

CVCBE   Fraction of the Coulomb barrier of an emitted charged particle
	assumed as minimum kinetic energy (default value: 0.30, if EXCN < 60
	then 0.1).

CVCB    Relative strength of the quadratic term (default value: 0.04).

VK      Minimum kinetic energy above Coulomb barrier (default value:
	EXCN/100).

CVK     Quadratic term (default value 0.05).

CARD 12   Cut-offs
------------------

EXR0    EX-range of the population matrices in the first decay steps
	(default value: 32 MeV).

EXH     Above this excitation energy a step size of 2 hbar in spin is used
	in the decaying nucleus, above 2*EXH 3 hbar steps are used (default
	value: 30 MeV).

CJC     Number of spin steps st the high-spin side of a population, where 1
	hbar steps are used (default value: ACN/20, but .GE. 4). CAUTION:
	The spin region with strongly varying fission competition should be
	covered with 1 hbar steps.

LCO     Angular momentum cutoff for all level densities (default value: CL0
	+2*DIFF).

CARD 13   Output control
------------------------

KOUTW   Lower cross section limit (in mb), above which population matrices
	are printed (default value: 100mb).

KOUT    Same for decay probabilities (default value: 1000 mb).

KEVAP	Same for evaporation spectra (default value: 1000 mb).

KGAMMA
        = 0 Do not calculate gamma-decay below particle thresholds.

        <>0 Calculate gamma-decay below particle thresholds (important shape
	of gamma spectrum at low energies and for gamma multiplicities).

        = 2 Calculate only particle and gamma spectra and population
	matrices following an initial gamma decay in the compound nucleus.

        =-2 Calculate only particle and gamma spectra and population
	matrices following an initial 4th particle decay in the compound
	nucleus.

i2dpop	=1 write out population matrix into i2dpop.dat

CARDS 13A Read in cascade structure
-----------------------------------

This card only read if KOPTK = 1 on card 4.

Input: IPSZO(IZF=1,4; IPS=1,IPS1) = assignment of numbers to the daughter
nuclei produced by n-, p-, alpha- and 4th decay (IZF=1,4) of each particle
emitting nucleus IPS = 1,IPS1 in the decay cascade. Defines length of decay
chains and order of computation. Insert 0 if decay is to be neglected.

CARDS 14 Individual level density parameters for low excitation energy
----------------------------------------------------------------------

IZ      Atomic number Z.

IA      Mass number A.

DA      Constant for calculating the level density parameter a = A/DA, (DA ~
	8 MeV).

DELTA	Zero of the thermal excitation energy for the g.s., U = EX - DELTA.

FTHETA  = J/Jrigid (like on card 7).

CK	Multiplicative constant in the level density (valid for the whole
	excitation energy range).

One card for each nucleus; arbitrary number of cards; arbitrary order.
Overwrites internally calculated standard parameters. End defined by a card
on which IZ = 0. This card is always necessary.

CARDS 15  Individual spectra
----------------------------

For each nucleus, for which low-lying levels are to be read in:

	a) One card containing IZZ,IAA,EX
	b) Cards containing EXL,JL,PL; format 7(F7.2,I2,A1).

End defined by a blank card, or an EOR card.

The input values are written, if they are needed, into a level table (for a
maximum of 99 nuclei).  If two spectra are supplied for the same or the
mirror nucleus, the second overwrites the first one. A maximum of 50 levels
can be given (with arbitrary order except for yrast levels or with
interruption except for the first position on the card).

CARD a

IZZ     Atomic number Z.

IAA     Mass number A.

EX      
	Excitation energy (in MeV) below which the analytical level density
	is erased and replaced by levels and above which the levels are
	interpreted as yrast levels (level density erased above the spin of
	the level). Yrast levels must be read in with increasing or constant
	spin at increasing excitation energy.

CARDS b   Format(7(F7.2,I2,2A)

EXL     Level excitation energy in MeV.

JL      Level spin (twice for odd nuclei).

PL      Level parity.

        = +  Positive parity.

        = -  Negative parity.

        = 0 or blank start reading cards for another nucleus (see card a).

        Any other character => undefined parity.


