Commit caf340ad authored by frey_m's avatar frey_m

add AMR test 9

parent 207a0941
This source diff could not be displayed because it is too large. You can view the blob instead.
Option, ECHO=FALSE;
Option, PSDUMPFREQ=4500;
Option, PSDUMPFRAME=BUNCH_MEAN;
Option, SPTDUMPFREQ = 1000000;
Option, ENABLEHDF5=TRUE;
Option, PSDUMPEACHTURN=TRUE;
Option, REPARTFREQ=5;
Option, ECHO=FALSE;
Option, STATDUMPFREQ=1;
Option, CZERO=FALSE;
Option, VERSION=20000;
Option, SCsolveFreq=2;
Option, AMR=TRUE;
Option, AMR_REGRID_FREQ=10;
Option, AMR_YT_DUMP_FREQ=1000000;
Option, MEMORYDUMP=TRUE;
Option, TELL=TRUE;
Title,string="OPAL-cycl: the first turn 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 =4;
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=32, AMR_MAXGRIDY=32,
AMR_MAXGRIDZ=32, AMR_BFX=16, AMR_BFY=16, AMR_BFZ=16,
AMR_REFX=2, AMR_REFY=2, AMR_REFZ=2, AMR_DOMAIN_RATIO={1.0, 0.75, 0.75},
AMR_TAGGING=CHARGE_DENSITY, AMR_DENSITY=1.0e-9;
beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=32*32*32, BCURRENT=2.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=2, MB_BINNING=GAMMA, MB_ETA=0.25;
endtrack;
Stop;
This source diff could not be displayed because it is too large. You can view the blob instead.
#!/bin/bash
cp AmrRingCyclotron.h5.bak AmrRingCyclotron.h5
cp AmrRingCyclotron.lbal.bak AmrRingCyclotron.lbal
cp AmrRingCyclotron.grid.bak AmrRingCyclotron.grid
cp AmrRingCyclotron.mem.bak AmrRingCyclotron.mem
cp AmrRingCyclotron-bunch-0000.smb.bak AmrRingCyclotron-bunch-0000.smb
cp AmrRingCyclotron-bunch-0001.smb.bak AmrRingCyclotron-bunch-0000.smb
mpirun -np 4 $OPAL_EXE_PATH/opal AmrRingCyclotron.in --restart -1 --info 2 --warn 0 "$@" 2>&1
This source diff could not be displayed because it is too large. You can view the blob instead.
"Neighboring bunch AMR restart test where we track two bunches by one turn. It uses the ML Poisson solver of AMReX and two AMR levels."
stat "rms_x" last 1E-7 AmrRingCyclotron-bunch-0000.smb # first bunch
stat "rms_y" last 1E-7 AmrRingCyclotron-bunch-0000.smb
stat "rms_z" last 1E-7 AmrRingCyclotron-bunch-0000.smb
stat "emit_x" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "emit_y" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "emit_z" last 1E-8 AmrRingCyclotron-bunch-0000.smb
stat "halo_x" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "halo_y" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "halo_z" last 1E-5 AmrRingCyclotron-bunch-0000.smb
stat "rms_x" last 1E-7 AmrRingCyclotron-bunch-0001.smb # second bunch
stat "rms_y" last 1E-7 AmrRingCyclotron-bunch-0001.smb
stat "rms_z" last 1E-7 AmrRingCyclotron-bunch-0001.smb
stat "emit_x" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "emit_y" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "emit_z" last 1E-8 AmrRingCyclotron-bunch-0001.smb
stat "halo_x" last 1E-5 AmrRingCyclotron-bunch-0001.smb
stat "halo_y" last 1E-5 AmrRingCyclotron-bunch-0001.smb
stat "halo_z" last 1E-5 AmrRingCyclotron-bunch-0001.smb
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
7b04ed94f4e45a299460d8eddf54fc45 AmrRingCyclotron-bunch-0000.smb
81b69f94ff2e6fde8f2b95ba5a010184 AmrRingCyclotron-bunch-0001.smb
5f9d4fa8a4e2c56f4807d4400feaf1ae AmrRingCyclotron.out
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment