Option, TFS=FALSE;
Option, ECHO=FALSE;
Option, PSDUMPFREQ=1;

Title,string="ISIS Test";

Edes=0.000296;
gamma=(Edes+HMMASS)/HMMASS;
beta=sqrt(1-(1/gamma^2));
gambet=gamma*beta;
P0 = gamma*beta*HMMASS;
brho = (HMMASS*1.0e9*gambet) / CLIGHT;


Q1: EQUADRUPOLE, L=0.095504, V=4000.0, ELEMEDGE=0.850087, APERTURE={0.0254};
Q2: EQUADRUPOLE, L=0.120900, V=-4000.0, ELEMEDGE=1.013307, APERTURE={0.0254};
Q3: EQUADRUPOLE, L=0.120900, V=0.0, ELEMEDGE=1.380993, APERTURE={0.0254};

l1: Line = (Q1,Q2,Q3);

Dist1:DISTRIBUTION, DISTRIBUTION=fromfile, FNAME="init_coords_fixed.dat";

Fs1:FIELDSOLVER, FSTYPE=fft, MX=32, MY=32, MT=2048, 
                 PARFFTX=true, PARFFTY=true, PARFFTT=false,
                 BCFFTX=open, BCFFTY=open, BCFFTT=open, 
                 BBOXINCR=1, GREENSF=INTEGRATED;

beam1: BEAM, PARTICLE=HMINUS, pc=P0, CHARGE=-1, SPACECHARGE=false, NPART=50001, BCURRENT=583e-6, BFREQ=23.06e6;

Select, Line=l1;

track,line=l1, beam=beam1, MAXSTEPS=100000, DT=1.0e-10, ZSTOP=1.53;
 run, method = "PARALLEL-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1;
endtrack;
Stop;
Quit;