_______________________________________________________________________ BEAM3D INPUT DATA _______________________________________________________________________ CN1 BSCOIL : Source magnetic field data filename. CN2 EEND : Extraction electrode data filename. CN3 EPULL : Puller electrode data filename. CN4 BFCOIL : Focussing solenoid data filename. CN5 9999 : IR--seed for random generation of : initial particle positions, : IR=1, manually input all positions CN6 400 : NP- number of rays, (NPmax=400) CN7 30*1,20*2,20*3,20*4,20*5,10*6,10*7,250*1,20*2 : N(i)*Q(i) (No. rays of same ion : species with charge state Q(i)) CN8 30*3,20*4,20*5,20*6,20*7,10*2,10*1,250*1,20*1 : N(i)*I(i) (N(i) No. rays, I(i) : Current (uA) carried by each ray), : Total Current = sum (N(i)*I(i)) CN9 130*14,270*4 : N(j)*A(j) (No. rays of the same ion : mass number A(j)) CN10 130*2,270*3 : N(j)*Tt(j) (No. rays with initial : thermal transverse energy Tt eV/per : charge state) CN11 130*1.5,270*2 : N(j)*Tl(j) (No. rays with initial : thermal longitudinal energy : (Tl eV/per charge state) CN12 .1,50,120 : DZ, NZ, M. DZ (mm) is the stepsize of : calculation, NZ is the number of steps : to take before a data point will be : written into the output data file : OPTF. The Z interval between writes : is then NZ*DZ. : (M+1) is the total data points written : to OPTF, and M*NZ*DZ is the total ray : tracing length. CN13 0.004,.00508,.0327,.076,.156,.3,.7 : EXR,PUR,EXSP,RPIP,EXEFL,EXBFL,FSFL (in M) : EXR - extraction aperture radius; : PUR - puller aperture radius; : EXSP - the spacing between the extraction : aperture and the puller aperture; : RPIP - the radius of the transport beam pipe; : EXEFL - Electric field data length in : the Z direction for files : EEND and EPULL; : EXBFL - Source magnetic field data : length (Z direction) for BSCOIL; : FSFL - Focussing solenoid field data : length (Z direction) for BFCOIL. CN14 0.,1,.05,2.2,.5 : FNEU1, NF, ZN1, ZN2, FNEU_NF : FNEU1 - degree of neutralization for : reducing the strength of the : space charge force over the whole : calculation. For FNEU1=1, no space : charge force will be calculated, : while FNEU = .2, 80% of the total : space charge force is used in : the ray tracing calculation. : NF = 1 indicates that in the region : [ZN1, ZN2], the level of space charge : force neutralization will be : defined by FNEU_NF. CN15 10000,0,1 : V_ext,V_pul. (Volts) : V_ext - extraction electrode voltage. : V_pul - puller electrode voltage : bpercent- a scaling factor for : the source magnetic field. This : parameter is intended to eliminate : the regeneration of the source : magnetic fields when a simple : scaling is desired. CN16 72.,.7544,.015,1.9 : SOL_CURI(A),SOL_LOC(M),OSLIT(M),OSLT_LOC(M) : SOL_CURI - the focussing solenoid current; : SOL_LOC - the location of the focussing : solenoid. It is the distance between the : extraction aperture (Z=0) and the point : where the focussing solenoid data is : taken to be effective. : OSLIT - the radius of the beam defining slit : at the object point of the analyzing dipole. : OSLT_LOC - the distance between the extraction : aperture and the object slit. ________________________________________________________________________ BEAM3D_READ INPUT DATA ________________________________________________________________________ CN1 m2a.GEO : Extraction electrodes, focussing solenoid and : object slit geometry data file. CN2 400 : NP CN3 .7544,.015 : SOl_LOC, OSLIT. CN4 .95,2,5 : FRAC, No. of ellipse plot, data point distance : interval. FRAC is the fitting factor for : the least square ellipse fit, FRACmax=1. CN5 .392,.592 : Locations (M) that ellipses will be generated and plotted. CN6 131,380 : Rays used to generate the ellipses. ________________________________________________________________________ ** Format of Source Magnetic Field Date File * (Br(r) and Bz(r) in Tesla, Z in M) Z(1),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8),Br(9),Br(10),Br(11) Z(1),Br(12),Br(13),Br(14),Br(15),Br(16),Br(17),Br(18),Br(19) Z(1),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8),Bz(9),Bz(10),Bz(11) Z(1),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16),Bz(17),Bz(18),Bz(19) ...... Z(4),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8),Br(9),Br(10),Br(11) Z(4),Br(12),Br(13),Br(14),Br(15),Br(16),Br(17),Br(18),Br(19) Z(4),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8),Bz(9),Bz(10),Bz(11) Z(4),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16),Bz(17),Bz(18),Bz(19) ...... Z(n),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8),Br(9),Br(10),Br(11) Z(n),Br(12),Br(13),Br(14),Br(15),Br(16),Br(17),Br(18),Br(19) Z(n),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8),Bz(9),Bz(10),Bz(11) Z(n),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16),Bz(17),Bz(18),Bz(19) * Note: Z(1)=-.003 M Z(4)=0.0 M Z(n)=EXBFL DZ=1 mm, DR=.1 mm for i from 1 to 11, DR=.5 mm for i from 12 to 19. *** In the POISSON calculation, the magnetic field direction should point th e first stage (opposit to the down stream beam line direction) direction, that means the Bz is positive and Br is negative. In BEAM_3D, the magnetic field is assumed point from source to the focussing solenoid, and the BEAM_3D rouine reverses the magnetic field direction accordingly. See the following sketches. ** Format of End Electrode Electric Field Data File * (Er(r) and Ez(r) in V/M, Z matches the Z in Source Magnetic * Field Data File) Z(1),Er(1),Er(2),Er(3),Er(4),Er(5),Er(6),Er(7),Er(8),Er(9) Z(1),Er(10),Er(11),Er(12),Er(13),Er(14),Er(15),Er(16),Er(17),Er(18) Z(1),Ez(1),Ez(2),Ez(3),Ez(4),Ez(5),Ez(6),Ez(7),Ez(8),Ez(9) Z(1),Ez(10),Ez(11),Ez(12),Ez(13),Ez(14),Ez(15),Ez(16),Ez(17),Ez(18) ...... Z(4),Er(1),Er(2),Er(3),Er(4),Er(5),Er(6),Er(7),Er(8),Er(9) Z(4),Er(10),Er(11),Er(12),Er(13),Er(14),Er(15),Er(16),Er(17),Er(18) Z(4),Ez(1),Ez(2),Ez(3),Ez(4),Ez(5),Ez(6),Ez(7),Ez(8),Ez(9) Z(4),Ez(10),Ez(11),Ez(12),Ez(13),Ez(14),Ez(15),Ez(16),Ez(17),Ez(18) ...... Z(m),Er(1),Er(2),Er(3),Er(4),Er(5),Er(6),Er(7),Er(8),Er(9) Z(m),Er(10),Er(11),Er(12),Er(13),Er(14),Er(15),Er(16),Er(17),Er(18) Z(m),Ez(1),Ez(2),Ez(3),Ez(4),Ez(5),Ez(6),Ez(7),Ez(8),Ez(9) Z(m),Ez(10),Ez(11),Ez(12),Ez(13),Ez(14),Ez(15),Ez(16),Ez(17),Ez(18) Note: Z(1)=-.003 M Z(4)=0.0 M Z(m)=EXEFL DZ=1 mm, DR=1.27 mm for i from 1 to 18. *** In the POISSON calculation, 10000 V voltage should be applied to the end electrode (Vpuller = 0) in order to match the scale factor inside the BEAM_3D routine. ** The Electric Field Data File for the Puller Electrode has the same Format as the End Electrode Data File, except -100 V is applied on the puller electrode (the end electrode Vend = 0) in order to match the scale in BEAM_3D. ** Format of the Focussing Solenoid Magnetic Field Date File * (Br(r) and Bz(r) in Tesla, Z in M) Z(1),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8) Z(1),Br(9),Br(10),Br(11),Br(12),Br(13),Br(14),Br(15),Br(16) Z(1),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8) Z(1),Bz(9),Bz(10),Bz(11),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16) ...... Z(i),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8) Z(i),Br(9),Br(10),Br(11),Br(12),Br(13),Br(14),Br(15),Br(16) Z(i),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8) Z(i),Bz(9),Bz(10),Bz(11),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16) ...... Z(k),Br(1),Br(2),Br(3),Br(4),Br(5),Br(6),Br(7),Br(8) Z(k),Br(9),Br(10),Br(11),Br(12),Br(13),Br(14),Br(15),Br(16) Z(k),Bz(1),Bz(2),Bz(3),Bz(4),Bz(5),Bz(6),Bz(7),Bz(8) Z(k),Bz(9),Bz(10),Bz(11),Bz(12),Bz(13),Bz(14),Bz(15),Bz(16) * Note: Here Z(n)s (in M) are in the POISSON calculation coordinates, Z(n)<0 for n0 for n>i. Z(k)-Z(i)=FSFL. BEAM_3D will converts these Z(n)s into the calculation coordinates Z when the calculation is carried out. Illustraction of the POISSON calculation coordinates and the definition of FSFL can be seen in the following sketch. ** Geometry Data File Definition and Format: * (All lengths in mm) * Format N,Z(1)....Z(N),R(1)....R(N) * N is the number of points that a part (end plate or puller or focussing solenoid or object slit) will be drawn in the output plot. Z(i)s are in the BEAM_3D calculation coordinates except for the focussing solenoid, that is Z=0 is the location of the edge of the extraction aperture in the puller side. * Note Zs: Since the focussing solenoid position can be varied in BEAM_3D calculation, special attention should be paid to this solenoid coordinates Zs. In the plot routine, the focussing solenoid will be plot at Z=SOL_LOC+Zs to match its real position. * Data Sample 6, 9.525,9.525,0.,-3.175,4.445,4.445,50.,21.59,4.,4.,21.59,50 (End plate) 9, 32.7,37.785,51.76,89.86,89.86,69.54,69.54,59.38,32.7 (Puller) 5.08,5.08,19.05,19.05,22.85,22.85,31.496,31.496,5.08 8, 101.29,133.29,155.26,155.26,75.89,75.89,101.29,101.29 (Ground Ring) 20.32,20.32,24.765,29.845,29.845,26.035,26.035,20.32 4, 0.,0.,428.88,428.88,150.,75.,75.,150. (Focussing Solenoid, note Zs) 4, 1900.,1900.,1905.,1905.,75.,0.,0.,75. (Object Slit) *** Plot Routine The Plot routine