CE    PROGRAM CASCADE							 CA   4
CE    PROGRAM CASCADE !(INPUT,OUTPUT,TAPE10,TAPE11,TAPE12,TAPE13,TAPE14) CA   5
CE									 CA   6
CE    EVAPORATION CODE BASED ON THE STATISTICAL THEORY			 CA   7
CE    FOR COMPOUND-NUCLEUS REACTIONS (HAUSER-FESHBACH THEORY)		 CA   8
CE    A DESCRIPTION OF THE FORMALISM USED IS GIVEN IN			 CA   9
CE    NUCL. PHYS. A280 (1977) 267					 CA  10
CE    F. PUEHLHOFER,  GSI DARMSTADT					 CA  11
CE									 CA  12
CE    THIS IS THE VERSION CASCADE I   (APRIL 1980)			 CA  13
CE									 CA  14
C EFG	Modified as described below by Elspeth Garman ,N.P.L. OXFORD.May 1982.
C SS    Modified as described below by Stephen Steadman, N.B.I. June 1983.
C JJG   Modified as described below by Jens Joergen Gaardhoeje, NBI
CE									 CA  15
CE    FILE 10: MASS TABLE, GENERATED BY EBTABLE SEPARATELY		 CA  16
CE                          AS A PERMANENT FILE,READ IN MAIN		 CA  17
CE    FILE 11: TABLE OF TRANSMISSION COEFFICIENTS, GENERATED BY TLCALC	 CA  18
CE                          SEPARATELY AS A PERMANENT FILE,		 CA  19
C SS                        READ IN TLREAD AS A BLOCK FOR USE BY TSUMME  CA  20
C SS            TSUMME accesses table in COMMON/C8
CE    FILE 12: TEMPORARY STORAGE OF LEVEL DENSITY MATRICES		 CA  21
CE                          (RTABLE, GETRHO)				 CA  22
CE                          SPACE: IPSMAX * 16 KBYTE MAXIMUM		 CA  23
C EFG		Direct Access file now - see below.
CE    FILE 13: TEMPORARY STORAGE OF POPULATION MATRICES WT(64,64)	 CA  24
CE                          (GETW, STOREW)				 CA  25
CE                          SPACE: ITMAX * 16 KBYTE MAXIMUM		 CA  26
C EFG		Direct Access file now - see below.
CE    FILE 14: SCRATCH FILE FOR SUBROUTINE CLEAN, SAME SIZE AS 13	 CA  27
C EFG		NOT NEEDED IN THE PRESENT VERSION
CE									 CA  28
C EFG	**********************************************************************
C EFG	Changes to programming in this version of CASCADE
C EFG
C EFG	1/  Notes about FILE 12 & 13 (direct access). All input/output to/from
C EFG   these files conforms to the DEC VAX 11/780 standards. Files are
C EFG   defined in the DEFFIL subroutine and new versions of GETW, STOREW,
C EFG   CLEAN, GETRHO and RTABLE are provided.
C EFG	See DEFFIL for adaption of file-defining statements to other compilers.
C EFG	READ(K,REC=N) means : binary read record number N from unit number K.
C EFG	2/  All PRINT statements have been changed to WRITE(K,format number)
C SS    3/  Instead of .plt files, "spectra" files for sap are created. The 
C SS    needed coding is labelled jjg.
C EFG	3/  A file called 'CASCADE.PLT' is opened on disk after line CA1659 and
C EFG	the summed evaporation spectrum of gamma rays is written to it(after
C EFG	line CA1686). If the cutoffs are set low so that the spectra for
C EFG	individual exit channels are printed out, a plot file for each channel 
C EFG	is produced with successive version numbers,X.i.e.CASCADE.PLT;X
C SS    4/  A table of needed transmission coeff. is read in TLREAD as array TL
C EFG	**********************************************************************
C EFG	**********************************************************************
C EFG 	Notes on modifications to physics in this version of CASCADE.
C EFG	
C EFG	A/ The gamma-ray spectrum emitted from the compound nucleus is now
C EFG	  included. However, if the first gamma ray comes from the initial 
C EFG	  compound nucleus, the rest of this decay chain is not calculated.
C EFG	  (N.B.The standard version DOES NOT calculate this
C EFG	  contribution to the summed evaporation spectrum at all.)
C EFG	  Changes necessary were provided by F.Puehlhofer and are marked
C EFG	  'PUEH' in the comment column. Lines changed were CA 321,359-361,445,
C EFG	  469,480. Lines were added after CA 440 and 564.
C EFG	
C EFG	B/ A Lorentzian GDR strength function can be included in the 
C EFG	  transmission coefficients. If the input E1 strength (input card 10)
C EFG	  is set negative, an extra input line is read by the program (input 
C EFG	  card 10a) to specify the GDR strength, excitation energy and FWHM.
C EFG	  The necessary changes were supplied by J.Barrette,B.N.L. and are 
C EFG	  marked 'BN X' in the comment column. (N.B. BN X numbers are not
C EFG	  sequential as BNL have a CDC computer and many of the changes to 
C EFG	  reduce the size of common blocks and overlay CASCADE were not needed
C EFG	  for a VAX.)   
C EFG	  Lines (no. of lines in parentheses) were added by BNL after:CA 97(1),
C EFG	  CA 149(3),CA 488(9),CA 984(2),CA 1010(1).Lines CA 146,CA 152,CA 489,
C EFG	  changed. Line CA 1748 changed from 'CALL STANLV' to 'CONTINUE'
C EFG	  and whole subroutine STANLV deleted.
C EFG	
C EFG	C/ The GDR can be represented by a double Lorentzian by putting input 
C EFG	  card 10a in as: strength, energy 1, FWHM 1, FRAC2, energy 2, FWHM 2,
C EFG	  where the first entry is the total strength of the double Lorentzian
C EFG	  and FRAC2 = fraction of total strength in the second Lorentzian.
C EFG	  Lines added after BN 26 and BN 27.
C EFG	
C EFG	D/ A mistake in the BNL changes regarding the energy dependent cutoff
C EFG	  has been corrected by adding 2 lines after CA 495.
C JJG  
C JJG   E/ AN INTERPOLATION TO ZERO OF THE GAMMA TRANSMISSION COEFFFICIENTS
C JJG      FROM 22 TO 32 MEV HAS BEEN REMOVED.
C JJG   F/ INCLUDE CALLS TO NBI-TAL SPECTRUM ANALYSIS SYSTEM FOR GAMMA 
C JJG      SPECTRUM OUTPUT
C JJG   G/ SIMULATE THE GIANT QUADRUPOLE RESONANCE IN THE STATISTICAL DECAY
C JJG      BY A DOUBLE LORENTZIAN FUNCTION:
C JJG      USE CARD 10B IF E2 STRENGTH IS < 0. THE INPUT PARAMETERS ARE
C JJG      SIMILAR TO THOSE NEEDED FOR THE GIANT DIPOLE RESONANCE.
C EFG	***********************************************************************
CE								!	 CA  29
