...
 
Commits (7)
21
0.0e+00 0.63013 1.08
5.0e-02 0.68626 1.08
1.0e-01 0.73857 1.00
1.5e-01 0.78677 9.20e-1
2.0e-01 0.83060 8.30e-1
2.5e-01 0.86979 7.35e-1
3.0e-01 0.90416 6.36e-1
3.5e-01 0.93349 5.34e-1
4.0e-01 0.95762 4.29e-1
4.5e-01 0.97643 3.21e-1
5.0e-01 0.98981 2.12e-1
5.5e-01 0.99766 1.01e-1
6.0e-01 0.99998 -9.31e-3
6.5e-01 0.99673 -1.20e-1
7.0e-01 0.98793 -2.30e-1
7.5e-01 0.97363 -3.40e-1
8.0e-01 0.95393 -4.47e-1
8.5e-01 0.92891 -5.52e-1
9.0e-01 0.89871 -6.53e-1
9.5e-01 0.86352 -7.51e-1
1.0e+00 0.82354 -7.51e-1
\ No newline at end of file
21
0.0e+00 0.109737 2.91
5.0e-02 0.258532 2.91
1.0e-01 0.401421 2.76
1.5e-01 0.535168 2.55
2.0e-01 0.656550 2.27
2.5e-01 0.762840 1.95
3.0e-01 0.851785 1.58
3.5e-01 0.921153 1.17
4.0e-01 0.969231 7.40e-1
4.5e-01 0.995244 2.92e-1
5.0e-01 0.998514 -1.65e-1
5.5e-01 0.978652 -6.20e-1
6.0e-01 0.936425 -1.05
6.5e-01 0.872880 -1.47
7.0e-01 0.789161 -1.85
7.5e-01 0.687305 -2.19
8.0e-01 0.569821 -2.48
8.5e-01 0.439238 -2.71
9.0e-01 0.298508 -2.88
9.5e-01 0.150978 -2.98
1.0e+00 0.222050 -2.98
\ No newline at end of file
Option, ECHO=FALSE;
Option, PSDUMPFREQ=100000;
Option, SPTDUMPFREQ = 10;
Option, PSDUMPEACHTURN=true;
Option, REPARTFREQ=20;
Option, ECHO=FALSE;
Option, STATDUMPFREQ=1;
Option, CZERO=FALSE;
Option, VERSION=10900;
Option, AMR=TRUE;
Option, AMR_YT_DUMP_FREQ=100000;
Option, AMR_REGRID_FREQ=5;
Option, TELL=TRUE;
Option, MEMORYDUMP=FALSE;
Title,string="OPAL-cycl: the first 11 turns acceleration in PSI 590MeV Ring";
REAL Edes=.072;
REAL gamma=(Edes+PMASS)/PMASS;
REAL beta=sqrt(1-(1/gamma^2));
REAL gambet=gamma*beta;
REAL P0 = gamma*beta*PMASS;
REAL brho = (PMASS*1.0e9*gambet) / CLIGHT;
//value,{gamma,brho,Edes,beta,gambet};
REAL phi01=139.4281;
REAL phi02=phi01+180.0;
REAL phi04=phi01;
REAL phi05=phi01+180.0;
REAL phi03=phi01+10.0;
REAL volt1st=0.9;
REAL volt3rd=0.9*4.0*0.112;
REAL turns = 11;
REAL nstep=360;
REAL frequency=50.650;
REAL frequency3=3.0*frequency;
ring: Cyclotron, TYPE="RING", CYHARMON=6, PHIINIT=0.0,
PRINIT=-0.000174, RINIT=2130.0, SYMMETRY=8.0, RFFREQ=frequency,
FMAPFN="s03av.nar";
rf0: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=35.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi01;
rf1: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=125.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi02;
rf2: RFCavity, VOLT=volt3rd, FMAPFN="Cav3.dat", TYPE="SINGLEGAP",
FREQ=frequency3,RMIN = 1900.0, RMAX = 4500.0, ANGLE=170.0, PDIS = 452.0,
GAPWIDTH = 250.0, PHI0=phi03;
rf3: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=215.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi04;
rf4: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=305.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi05;
l1: Line = (ring,rf0,rf1,rf2,rf3,rf4);
Dist1:DISTRIBUTION, TYPE=gauss,
sigmax = 2.0e-03,
sigmapx = 1.0e-7,
corrx = 0.0,
sigmay = 2.0e-03,
sigmapy = 1.0e-7,
corry = 0.0,
sigmat = 2.0e-03,
sigmapt = 3.394e-4,
corrt=0.0;
Fs1:FIELDSOLVER, FSTYPE=ML, MX=32, MY=32, MT=32,
PARFFTX=true, PARFFTY=true, PARFFTT=true,
BCFFTX=dirichlet, BCFFTY=dirichlet, BCFFTT=dirichlet,
BBOXINCR=20, AMR_MAXLEVEL=2, AMR_MAXGRIDX=16, AMR_MAXGRIDY=16, AMR_MAXGRIDY=16,
AMR_TAGGING=CHARGE_DENSITY, AMR_SCALING=0.5, AMR_DENSITY=1.0e-10;
beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=32768, BCURRENT=1.0E-3, CHARGE=1.0, BFREQ= frequency;
Select, Line=l1;
TRACK,LINE=l1, BEAM=beam1, MAXSTEPS=nstep*turns, STEPSPERTURN=nstep,TIMEINTEGRATOR="RK-4";
run, method = "CYCLOTRON-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1, MBMODE=FORCE, TURNS=3, MB_BINNING=GAMMA;
endtrack;
Stop;
This diff is collapsed.