ICOOL Revision History 3.20 22 Jul 09 add KICK(3); add NPSKIP control variable; add d, He3, Li7 as defined particles; add stainless steel as defined material; modify for009 header for compatibility with HistoRoot; modify material interaction routines to take particle charge into account; modify dE/dx for compounds; modify max number of SuperFish axial grid points; modify DENS command to refer to reference (not current) density value; move LBNL emittance calculation to EMIT command; delete old HELIX(5); fix error in KICK for t=0; fix error in setting up SuperFish grid; fix error in CUTV command; fix error in DENS logic; fix error in DENP logic; fix error using WEDGE with two materials; fix error in WEDGE for steps that partially enter wedge; fix error in old HELIX(6), rename as new HELIX(5); 3.17 20 Jan 09 modify ACCEL(5) and (11) to allow SuperFish cavities in dipole regions; modify APERATURE to allow off-center rectangular cuts; modify multipole expansions for orders 1-3 and constant h; modify DIP, HDIP, QUAD, SQUA and SEx commands so model=1 uses 2nd order; modify some IFLD and IINT routines to use double precision constants (S. Berg); fix errors in HELIX(6); fix initialization error in DECAY; 3.13 15 Jul 08 add HELIX model 6 including torsion; removed &RAN command; modified BSOL(3) logic; modified straggling default to LDRAY=true; fix errors in HELIX; fix errors in SAMCS implementation; 3.10 26 Feb 08 add TAPER command; add BLOCK(4) model for tapering; add BSOL(5) model for tapering; add SAMCS model for scattering and straggling; modify ipnum incrementing with decays; modify for007 print out to 5 decimal places; modify name substitution to allow & ... / namelist format; fix problem of * welcome print out with termout; fix problem in SOL(1); fix Lz problem in LOADVAR; 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