ICOOL   Revision History

3.05  19 Mar 07	add SOL(10) for given on-axis field;
		add HEL(5) for Fourier desciption of helical multipoles;
		add helical wiggler model WIG(5);
		modify BSOL(3) to allow 4th and 5th order;
		modify BSOL(3) to include a0 multipole;
		modify BSOL(3) to increase execution speed;
		modify space charge calculation of bunch sigmas;
		fix normalization error  in BROD;
		fix errors in WIG(4);
		fix error with z0 on beam restart;

3.01  18 Sep 06	add bent axial current rod (BROD) field model;
		add static electric field model (EFLD), J.S. Berg;
		add new phase models 5 and 6 with quadratic dependence
		 of energy on distance;
		modify ACCEL(2) to allow transverse displacement of cavity;
		modify tracking to allow particles to start with z > 0;
		modify logic to require a REFP command for beam restarts;
		modify control variable F9DP to allow 3 digit exponents;
		modify grid interpolation to give B=0 outside grid;
		update particle masses from PDG;
		update material mean excitation potentials; 
		increase maximum number of characters per line
		 for name substitution;
		assign mass(6) to Li7;
		fix bug in GO_REGION that affected polarization;
		fix bug in probability in RESTRICT_LOSS;
		fix bug in square root in MUON_DECAY;

2.96  10 Apr 06	add Fano scattering model;
		add Tollestrup scattering model;
		add simple helical wiggler model;
		add planar wiggler model;
		modify default scattering model from Bethe to Fano;
		modify SUB command to only process to first ENDS;
		modify ELMS to read ELMSCOM.TXT for paths;
		modify histogram HCPRN command to write to specified file;
		modify curvature factor logic in BSOL;
		increase number of particles in memory to 300000;
		fix error in phasemodel 4 time;
		fix error in Twiss beam correlation;

2.92  17 Nov 05	changed ACCEL model 5 to read Superfish output;
		changed HORN model 2 to allow radially stacked horns;
		fixed error in BSOL and DIPOLE for negative fields;
		fixed error in ROTATE for x and y axes;
		fixed bug in CUTV;

2.89  03 Jun 05	add SOL(9) to get on-axis field from Fourier coefficients;
		add DENP command for making variable density profile;
		add ability to specify material on exterior of a WEDGE;
		force VARSTEP=F when PHMODREF=2;
		fix bug in MAKE_BEAM that gave initial particles with z < 0;
		fix bug in writing Bz to for009;
		move FIELD call for phasemodels to get ACCEL(11) to work
   		 correctly;

2.85  01 Feb 05	add BMALT beam variable to generate particles with both signs;
		add CUTV control variable to cut on ICOOL variables;
		add dispersion correlation for initial beam;
		add new solenoid model 8 for hard-edge with adjustable end fields;
		add ELMS routines for correlated energy loss and scattering;
		modify BSOL(1) to allow independent curvature;
		modify ICOOL version number to allow character expressions;
		modify phasemodels 3 and 4 so that fixed refernce particle stepping
		 exactly matches normal particle stepping;
		modify RTUPLE so that OUTPUT information for reference and normal
		 particles are the same;
		modify calculation of energy gain for phasemodel 4;
		modify BSOL(3) to specify geometric curvature, allow 2000 points;
		modify HELIX model 3 for multifilar multipoles;
		fix error in edge field for ACCEL(13);
		fix error in random number seeds;

2.79  01 Nov 04	add polyethlyene as a material;
		add parameter 9 to ACCEL(2) for equal-longitudinal length cavity
		 in a circular geometry;
		add normal and skew quad aperture cut;
		add new phasemodel 2 to iteratively find 0-crossing time in
		 rf cavities (Scott Berg);
		add ACCEL model 13 for open-cell RF cavities
		 with edge focusing (Scott Berg);
		add parameter 7 to DIP(3) for negative bend, radial sector;
		add control variable DIAGREF;
		modify field superposition by calling CELL fields before REG fields.
		 this removes conflict between BSOL(4) and ACCEL(11);
		modify rotation for EDGE:DIP3 model;
		fix error in RF diagnostics file;
		remove extraneous random generation code from REG_SUMMARY;
		
2.77  07 Jul 04	add scaling function &SCL to name substitution;
		add new control variable SUMMARY to print out
		 region summary tables;
		add Wang pion production from proton interactions;
		modify ICOOL region logic to no longer use internal direct 
		 access file for007;
		modify name substitution logic to allow & in titles;
		modify REPEAT loops to allow REF2;
		modify BSOL logic for increased speed;
		modify INTERACTION logic to include all hadrons;
		modify NEIGHBOR logic, only allow 100 regions;
		modify RKICK logic, only allow 100 kicks;
	 	modify terminal listing to every 250 particles;
		delete phasemodel=2 logic;
		fix error using GRID with REF2;
		fix incorrect error message when just running with the production
		 region;
		fix bug in muon decay angular distribution;
		fix error in DATA statement location in BBLOCS;
		fix bug in HELIX (4) "segmentation fault";

2.72  08 Apr 04	modify HELIX models 1,2 and 4;
		modify BLOCK models to allow option of series algorithm;
		change substitution character from \ to & for UNIX compatibility;
		modified definition of internal variables 21-27;
		modify phasemodel 4 to allow x,y to change in ACCEL regions;
		change BZFLDPRD to control variable;
		change PZCORR emittance variable to IPZCOR;
		print out total good particles at end of run;
		modify SHEET(5) to allow absolute coordinates in grid;
		increase maximum length of substitution value string;
		modify stepping logic to give integral # of equal steps
		 for fixed stepping;
		change TROD command to ROD model;
		change TSOL command to SOL model;
		modify EDGE(FACE) to make horizontal focusing optional; 
		fix bug calling ZHIST for second reference particle;
		fix bug in SHEET moving grid;

2.67  01 Dec 03	add name substitutions for for001.dat;
		add random variable initialization for for001.dat;
		add interaction control variable FASTDECAY;
		add new BLOCK algorithm;
		modify RKICK command for scaling, random azimuth,
		 solenoidal kicks;
		modify DISP and TILT to allow uniform random azimuthal
		 directions;
		modify muon decay routines from 3-body phase space
		 to unpolarized V-A distributions;
		modify DIP(3) to allow user control of curvature;
		modify BLOOP and BSHEET routines;
		modify 2D emittance calculations to use canonical values;
		fix error in for003 reference particle input;
		fix error in restricted Landau distribution;
		fix error in SCATTER that gave wrong dE/dx in gas-filled
		 RF cavities;

2.62  18 Sep 03 add control variable F9DP to set # of decimal places for
                 floating point variables in for009;
		add pseudocommand RKICK to generate kicks from random magnetic
		 field errors;
		add parameter 12 to ACCEL(10) to set phase using average of 2 
		 reference particles;
		add new parameter when reading in field grids to
                 longitudinally shift values stored in memory;
		modify 2D and 3D magnetic field grids generated by ICOOL to
                 have same format as input grids;
		modify format of DVAR and ROTATE commands to specify whether
                 change applies to reference particle;
		modify x' and y' definition in DVAR;
		modify DIP(3) to allow general field index n;
		modify FORCERP to conserve total momentum;
		fix error in Rutherford scattering for diffuse medium;
		fix error in HELIX tracking parameters;
		fix error in BACKGROUND initialization (U. Bravar);
		fix error in Lz correlation for internal beam generation;
		fix error in arclength calculations;

2.57  29 Apr 03 add EDGE command for fringe field kicks from hard-edge fields;
		new subroutine HRDEND for multipole end kicks (S. Berg);
                SKICK and FACE commands replaced by EDGE;
		add DVAR command to change particle variables by fixed amount;
		DELE command replaced by DVAR;
		modify DIP model 4 to eliminate higher multipoles in field;
		modify STUS field to allow unformatted grids, spline coefficient
		 input, and user control of curvature;
		modify BSOL model 4 multipole fields;
                modify GRID to allow predefining 3D grids;
		fix problem with ACCEL (11) and STUS cell fields;
		fix problem with NaN written to for009;

2.50  02 Jan 03 add DIP model (4) with rotated pole faces;
                add g curvature factor to BSOL(4);
                add gas H and HE materials;
                add DENS pseudocommand;
                add PHANTOM control variable;
                add p key input during execution;
                add Non-Isoceles Absorber;
                modify BSOL(4) to allow 200 Fourier coefficients;
                modify logic to allow REFP in repeats and cells (Fukui);
                modify TRANSPORT to check for unphysical momentum (Berg);
                modify BSHEET to match RFOFO code;
                modify BSOL(4) to change sign of h,g for By < 0
                 and q < 0;
                modify DIP, HDIP, QUAD, SQUA & SEX to agree with MC207
                 expansion;
                modify DIP, HDIP to change sign of h,g;

2.45  29 Aug 02 rewrote BSOL model 4 for speed;
		added Bx data to BSOL model 4;
		modify FSAV logic to work with two reference particles;
		add higher order terms to DIP model 1;
		modify TRANSPORT to eliminate "length" parameter;
		modified DERIV code (Berg);
		add control variable FORCERP to force reference particle onto
                axis;
2.40  17 May 02 add RESET and DELE commands;  
		add ACCEL model 11 for straight pillbox in dipole region;
		add ACCEL model 12 for sector RF cavity with rectangular
		 cross section (Berg);
		modified Balbekov beam correlation;
		modified REFP logic to not reset transverse position or
		 momentum;
		modified ACCEL model 9 to eliminate scaling of field
		 strength with radius;
		modify random number generator for double precision;
		fix error in DIP model 4;
		fix error in HORN model 1;
		fix error in ellipical APERTURE command (Bracker);
2.33  4 Mar 02  add ability to put blank lines and comments anywhere in
		 for001.dat;
		add unique IPTYP codes for neutrinos;
		changed radiation length of liquid hydrogen to 866 cm;
		changed density of LIH to 0.78 g/cm^3;
		changed radiation length of LiH to 102 cm;
		fix initialization error in RUTH_SCAT (Fukui);
		fix error (x,y not defined) in BSOL model 3 (C. Wang);
		fix error in SKICK command for -ve charge particles;
		modified spacing in print out;
		fix logic error in PHASEMODEL (Bracker);
		fix error in ACCEL model > 5;
		fix error in HELIX model 3;
		fix initial field error in OUTPUT (Penn);
2.32  9 Jan 02  add GRID command to allow multiple, simultaneous r-z grids;
		add BEGS command to control the amount of region description
		 that is repeated using the NSECTIONS control variable;
		add liquid helium (LHE) as defined material type;
		modify definition of Balbekov correlation;
		modify r-z grid interpolation for speed (Penn);
		modify elliptic integral routines for speed (Penn);
		replaced rf phase with arclength in FOR009.DAT;
		fix test for KICKER in GO_REGION;
		fix error in DIP model 4;
		fix error in phasemodel = 2;
		fix error in uniform beam distribution (Fukui);
2.30  5 Dec 01  fix instances of single precision constant arguments
		 in calls to double precision subroutines;
2.29  3 Dec 01  modify BSOL, model 4 to include sextupole components
		 and longitudinally-varying coefficients;
		add test on Boris stepping in curved regions;
		modify contents of RF diagnostic file;
		define all real variables as double precision;
		fix error in azimuthal distribution of scattered particles
		 for steeply inclined tracks (Penn);
		fix error in SKICK command;
		fix error in NEIGHBOR command;
		fix error in uniform beam distribution (Fukui);
2.27  30 Oct 01 modify DIP, model 1 to include hard-edged higher
		 multipoles;
		add DIP, model 4 for hard-edge combined function dipole;
		add test on 3rd order contribution to BSOL, model 3;
		modify fitting in BSOL, model 3 for increased accuracy;
		new BSOL, model 4 with Fourier component input and
		 5th order fields;
		add SKICK command for fringe field at end of hard-edge 
		 solenoid;
		new beam correlation 10 for Balbekov version of amplitude-
		 momentum correlation;
		for BLOCK and COIL the current scale factor has been moved
		 to FOR001 from the external files;
		add dispersion calculations to covariance print out;
		add new variable DISCORR to enable corrections to 2D
		 emittances for dispersion;
		fix bugs in TWISS initial beam correlations;
		add tests in scattering routines to ensure algorithms
		 are in range of validity;
		replace LEN_TRIM calls with LASTNB (S. Berg);
		fix region summary print out for negative file numbers;
2.24  10 Sep 01 add HORN magnetic field model;
		add FACE psuedocommand for rectangular dipoles;
		add Boris push integrator option (Stoltz,Penn);
		add new SCALESTEP variable to modify all step
		 sizes in a problem simultaneously (G. Penn);
		enhance error message print out to log file;
		modify ICOOL to be compatible with OPTICOOL (Bracker,Kirk);
		fix errors in SOLENOID, model 6;
2.23 25 Jul 01  add beam envelope calcualtion mode of running ICOOL (G. Penn); 
		add two new kicker models;
		add ability to save neutrino production data;
		add scale factors to bent solenoid model 3;
		add energy loss model 4 with dE = const * dz; 
		add lead as material type;
		add solenoid model with user-supplied r-z map;
		modify bent solenoid fields to use C. Wang 3rd order expansions;
		modify dipole fields to use C. Wang 3rd order expansions;
		correct bug in BESSI1;
		fix data type error in GO_REGION;
2.20  5 Apr 01  add 2nd reference particle; 
		add ACCEL model 10 for adiabatic buncher;
		add space charge model 4 for single bunch in a train;
		add diagnostic print out for space charge kicks;
		add ability to insert comment lines in region defintion;
		modify argument checking for I0,I1,K0,K1 Bessel functions;
		modify stepping argument checking for transverse position,
		 momentum;
		modify DZCUT logic in SHEET (Penn);
		fix unit 7 file size error;
		fix background field errors;
2.17 23 Jan 01  add RING geometry type;
		add SHEET, model 4 for "moving" grid;
		add TILT and DISP pseudoregions;
		add DECLEV=5 to suppress muon decays and sum accumulated 
		 decay length;
		modify SHEET, ... to allow -ve file # => opposite polarity;
		increase radial grid size to 100;
		fix interpolation logic in BENTSOL;
		fix sign error in BLOCK;
2.15 20 Nov 00  add axi-symmetric radial wedge(Penn);
		add beam correlation type 9 (equal time in solenoid);
		modify LBNL emittance routine(Penn);
		modify WIGGLER routine(Penn);
		modify summary table to give sheet file numbers;
		modify ECALC9 postprocessor(Penn);
		fix real*8 declaration problem in HELICITY, BL_...;
		fix problem with HELIX end field attenuation;
		fix problem with Bessel functions for small arguments;
2.13 13 Oct 00  fix error in 2nd order constraint in BSOL, DIP, QUAD;
		fix error in reusing files in BSOL model 3;
		fix sign error in initial beam angular momentum correlation;
2.12 9 Oct 00   add HELIX field models;    
		add horizontal dipole element;
		add skew quad element;
		add BSOL model 3 for user input of fields on axis;
		add ACCEL model 9 for sector shaped rf cavity;
		add LINTERACT variable for pion inelastic nuclear 
		 interactions;
		modify tracking to take into account curvature in y-s plane;
		modify calculation of electric field in curved regions;
		modify BACKGROUND to add BENTBKG parameter;
		modify STUS input to add field magnitude scaling factor;
		modify DIP, QUAD and SEX input parameters;
2.08 28 Jul 00  add helical wiggler field (Penn);
		add thick solenoidal current block (Gallardo);
		add new DECLEV models with no mu->e decays;
		add new diagnostic file for ACCEL model 8;
		add SOLENOID model 5 using block conductor;
		modify COIL and SHEET input format to include title field
		 and current scaling parameter;
		modify TYPEing for CRAY portability;
		fix bug in ACCEL model 6,8 reset parameter;
		fix bug in depolarization calculation for elastic scattering;
2.06 11 May 00  add ACCEL model 8, induction linac model with waveform 
		 from user file;
		add add parameter (10) to COIL and SHEET to control
		 calculation of field grid;
		improve spin tracking;
		fix error in calculation of depolarization from elastic
		 scattering;
2.05 2 May 00   add DUMMY pseudocommand (for OUTPUT);
		add ACCEL model 7 for induction linac;
		add new spin flip probability for Coulomb and
		 elastic scattering (Gallardo);
		modify phasemodel 3 and 4 to be read from REFP command;
		modify CELL input commands to keep field tag and its
		 parameters together;
		modify time limit default to be very long;
		modify internal time diagnostics to use ct;
		modify ACCEL model 4 to use parameters 4-7 for phase;
		modify ACCEL model 6 so that voltage pulse is non-zero before
		 pulse "start" time;
		modify LBNL emittance calculation (Penn);
		improve accuracy in muon spin tracking (Fukui);
		fix bug that generated 2 ref par header lines in for004.dat;
2.03 2 Feb 00   modified induction linac to use microseconds and allow
		 up to 8 voltage coefficients (R. Scrivens);
		enhanced run time status message (R. Scrivens);
		add hemispherical absorber end region;
		fix bug in spin tracking (lost initial polarization
		 after 1st step;
		changed spin tracking to use reference particle coordinate
		 system and Jackson spin transformations.
		 (Y. Fukui, J. Gallardo);
		added more spin diagnotic print out;
		change logic so GOODTRACK=F resets IPFLG on beam input
		 and allows program to process failed tracks;
		fix bug in generating beam output file for failed tracks;
		fix bug for sheet routine RJ failures;
2.02  7 Dec 99  added induction linac model;
		improved accuracy of spin tracking;
		reformated region summary and other log file print out;
		added spin information to beam input/output files;
		new control variable PRNMAX;
		fix bug with # of sections > 1 giving multiple PRODs;
		fix bug with delta production in wedge absorbers;
2.01 22 Nov 99  added calculation of emittance and related quantities
		 according to the LBNL prescription;
		added 2D Mafia accelerator model;
		added mercury as allowed material;
		increased number of allowed z-histories to 20;
		increased default accuracy in stopping at end of region;
		modify for009.dat output format;
		modify routine PZ_AMP_CORR to include weights;
		fix bug preventing the use of 100 emittance and covariance
		planes;
		fix rtuplen logic to allow output at every step;
		fix error in computing energy loss with delta rays;
2.00 26 Oct 99  increase accuracy of dE/dx calculation using D.P. variables;
		use higher accuracy tracking equations in straight regions;
1.99 25 Oct 99  fix error in calculating longitudinal emittance;
		fix error in space charge model 3;
		fix logic error that caused erroneous (-36) errors;
		fix error in generating particle decays;
		fix mu polarization calculation;
		reformat for009 header;
		add weighted emittance calculation( C.Kim, Y.Fukui);
		add mean value to covariance output;
		increase allowed # emit, covar planes, radial grid points;
1.98 8 Oct 99   fix bug in OUTPUT1 logic;
		list every 25th particle on the screen;
1.97 29 Sep 99  control variable RFPHASE now allows user to input 
		rf freq and gradient from file;
		new control variable OUTPUT1 for output of production info;
		new control variable SPINTRK to control spin tracking;
		3-body phase space decays implemented;
		mu polarization in pi decay implemented;
		6-D emittance calculations combined together in new table;
		fix bug using 2nd REFP in problem
		fix bug in ZHIST output;
		fix bugs in for009.dat OUTPUT file;
		fix bug in for004.dat beam info file;
1.96 26 Aug 99  include patch for UNIX problem reading records from file 7;
		add ability to produce all decay products at 90 deg in
		parent particle rest frame;
		add polarization info to OUTPUT file;
		add TREFMEAN to rf diagnostic output;
		fix error in time distribution for ifail=-1;
		fix error in header for for004.dat for izfile=-1;
		fix error in dipole initialization;
1.95 3 Aug 99   fix error in space charge model 2; 
		add space charge model 3;
1.94 2 Jul 99   fix logic error in OUT_STATS;
		add transverse free space charge model (2)
1.93 1 Jun 99   fix output pseudocommand logic error
1.92 17 May 99  decrease MXZG to 200
1.91 14 May 99  interpolate h and h' in background fields
1.90 13 May 99  save curvature in table for background fields;
		fix error in calculating curvature for background;
1.89 7 May 99   add static field, user input for background field;
		replace BFIELD parameter ZCUTBKG with 
		ZMINBKG and ZMAXBKG (INPUT change);
		increased allowed background grid size to (30,30,400);
		fixed logic errors in setting up background field;
1.88 12 Mar 99  add R-history diagnostics (INPUT change);
		add FFCR control variable for print control;
		add dTANH model for Li rod;
		modify input beam correlations (INPUT change); 
		add initial angular momentum correlation;
		add initial RF bucket correlations;
		add initial Twiss parameter correaltions;
		modify Z-history input to use Ymin,Ymax and always
		 23 bins (INPUT change);
		stop using 2 sets of plane numbering;
		removed redundant print out of 3rd order moments;
		fix error in calculation of canonical momentum 
		 correction to emittance;
1.87 7 Feb 99   add control variable PZMINTRK to stop tracking very low
		momentum tracks;
		increased background field grid size;
		fix error in background field calculation;
		fix error in getting emittances after TRANSPORT region;
		fix calculation of 6D emittance;
		fix error in getting output from OUTPUT command;
		fix error in phasemodel=3 for reference particle from
		external beam file;
1.86 11 Dec 98  allow reference particle to be redefined;
		add variable BZFLDPRD to define Bz at production;
		add field level offset parm to solenoid, model 1;
		make equations sensitive to sign of particle charge;
		changes accelerator, model 3 to travelling wave linac;
		fix bug in covariance calculation after 1st requested plane;
		fix bug in solenoid, model 4 fields;
1.85 13 Nov 98  increase mxzgr -> 5000, mxpar -> 50000;
		fix bug not permitting W & LIH materials;
1.84 28 Oct 98  allow covariance calculation at arbitrary planes;
		add calculation of LBL emittance;
		refine default parameter values;
		fix bug in emittance calculation when sigma=0;
		fix bug in decay daughter particle distributions;
1.83 21 Sep 98  add delta ray interactions;
		add multiple scattering fudge factors;
		modify step size algorithm to include material medium;
		fix decay routine;
		fix problem with random number re-seeding;
1.82 11 May 98  add space charge model 1 (image charge in can);
		fix logic error in ACCEL, model=4;
1.81 27 Apr 98  fix logic errors in FILL_ZHIST, OUT_STATS,
		SIMULATE npart; sheet grid interp=3;
		add warning when grid point is on current sheet;
1.80 8 Apr 98   add RF diagnostic output file;
		add phasemodel=5 RF phase input file;
		modify OUTPUT to include electric field, phase;
		fix 2-body decay kinematics;
1.79 3 Apr 98   fix program restart with reference particle
		(this changes FOR003,4.DAT file formats)
1.78 23 Mar 98  add Pz - transverse amplitude correction to long emit
1.77 16 Mar 98  fix background field initialization problem
1.76 11 Mar 98  fix h=0 error in long bent solenoids;
		fix for007 errors from OUTPUT in cell + repeat loops;
		fix bsheet model 1,2 sign error;
		add Fe and Cu materials;
		add background field grid output to file;
1.75 26 Feb 98  fix error in BSHEET;
		add {coil,sheet} field grid output to file;
		fix normalization in covariance print out;
1.74 18 Feb 98  add 4 sigma tail cut for emittance;
		add canonical momentum corection in solenoid;
		add rtuple control variable for OUTPUT;
		remove OUTPUT from z-history;
		force fixed steps if multiple radial regions;
1.73 5 Feb 98   add covariance of final particle distribution;
		use common/RZGRID/ for both COIL and SHEET;
		add user TRANSPORT matrix input;
		implement background fields;
1.72 30 Jan 98  add asymmetric, polynomial shaped wedge;
		add BUNCHCUT variable;
		add calculation of 3rd and 4th moments;
		force fixed step size in wedge regions
1.71 28 Jan 98  fix longitudinal emittance normalization;
		add OUTPUT via z-history;
		add weight field in input and output files;
		change OUTPUT format;
		increase COIL, SHEET radial grid points;
		allow use of phase0 in phasemodel=3
1.70 26 Jan 98  fix angle effect on step size in material
1.69 15 Jan 98  add quad,cubic interpolation to COIL,SHEET grids;
		add constant offset parameter for solenoid;
		fix problem with energy loss when scattering on;
		fix sextupole constraint in dipole field;
		fix problem with ROTATE
1.68 22 Dec 97  add constant offset parameter for bent solenoid;
		increase # of allowed sheets
1.67 5 Dec 97   added 2nd order bent solenoid; stopped tracking
		phasemodel=3 reference particle
1.66 28 Nov 97  added current sheet fields
1.65 23 Nov 97  removed limit on NPART; revised phasemodel=3 logic
1.64 21 Nov 97  azimuthally adjustable wedge; add electric field
		terms to accelerator equations of motion
1.63 14 Nov 97  changed region-particle loop order; background field
1.62 11 Nov 97  change to spatial steps in input;
		use accelerator coordinates exclusively
1.61 30 Oct 97  increase COIL z grid to 1000, # coils to 1000
1.60 30 Oct 97  add APERTURE, OUTPUT commands, GRADREF variable;
		made code F77-compatible; corrected SCATTER
1.59 27 Oct 97  add phasemodel=4
1.58 23 Oct 97  add liquid hydrogen and other materials
1.57 8 Oct 97   increased array sizes for grid in COIL
1.56 7 Oct 97   ROTATE command; title field
1.55 6 Oct 97   improved portability; ntuple option; neighbor region field
1.54 24 Sep 97  added polynomial shaped absorber wedge
1.53 24 Sep 97  added PHASEMODEL including ability to set cavity phase using
		reference particle
1.52 23 Sep 97  added bent solenoid; modified spatial stepping in DERIV
		to use external field; increased # of field parameters to 15
1.51 2 Sep 97   added REPEAT; fixed over stepping problem; made TSTEP a region
		parameter; no scattering or straggling for reference
		particle
1.49 26 Aug 97  add BETAPERP control; fix bug in wedge logic
1.48 13 Aug 97  generalized material input; Moliere for compounds
1.47 6 Aug 97   add ref particle logic; control on RN seed
1.46 18 Jul 97  fix bug in wedge energy loss
1.45 17 Jul 97  fix bug in TE011 field components  
1.44 30 May 97  allow control of radial grid size [COIL]
1.43 6 May 97   allow forcing phase0 to given value; Hist limit raised to 20
1.42 2 May 97   allow multiple coil grids in problem, repeated use of grid, 
		CELLFLIP variable
1.41 29 Apr 97  S.W. deflection cavity, LiH, increase coil grid size   
1.40 8 Apr 97   changed coordinate system for TM011 in ACCEL
1.39 28 Mar 97  Fixed bugs in DERIV, SETUP, OUT_ZHIST,
		    IN_{C,L,R}WEDGE, DIPOLE
1.38 17 Mar 97  added CWED; changed FOFO Bz to sine function
1.37 10 Mar 97  add TM011 model to ACCEL
1.35 27 Feb 97  change stsp,pol arrays to real*8
1.34 18 Feb 97  add new diagnostic variables
1.33 4 Feb 97   add single loop coil model
1.32 30 Jan 97  add TANH model for quad and sextupole
1.31 23 Jan 97  add constraints for dipole TANH model
1.30 21 Jan 97  add wedge absorber 
1.29 16 Jan 97  use accelerator coordinates for dipole tracking
1.28 26 Dec 96  separate lambda and ELEN parameters for TANH solenoid
1.27 23 Dec 96  add constant Ez model to ACCEL
1.26 20 Dec 96  add order variable to TANH solenoid model
1.25 8 Dec 96   add jcel, kcel accelerator phase logic
1.24 4 Dec 96   add DPHASE0 variable to ACCEL
1.23 7 Nov 96   include spin tracking
1.22 28 Oct 96  generalize beam correlations; include Palmer amp
1.21 8 Oct  96  add arclength variable
1.20 1 Oct  96  add z-history print out
1.19 17 Sep 96  fix long emit; add traveling wave accelerator
1.18 16 Sep 96  set B=0 beyond grid in COIL,SOLENOID
1.17 4 Sep  96  add z history; Lz variable
1.16 21 Aug 96  fix coil grid over cell
1.15 15 Aug 96  add grid model to COIL, SOLENOID
1.14 8 Aug  96  fix multi-beam error; phi,pphi,pr plots
1.13 6 Aug  96  add sine z + Bessel r model to SOLENOID