Commit eefb63c3 authored by frey_m's avatar frey_m

add AMR test 10

parent caf340ad
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=AMR_MG, 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,
AMR_MG_REBALANCE=TRUE, AMR_MG_NSWEEPS=12, ITSOLVER=SA, AMR_MG_NORM=LINF;
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 AMR_MG Poisson solver based on Trilinos packages 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
d9ed3f591d5eb0f611bd12c1991f40bc AmrRingCyclotron-bunch-0000.smb
This source diff could not be displayed because it is too large. You can view the blob instead.
8d3416a2f9454786da21e434753ac815 AmrRingCyclotron-bunch-0001.smb
76478ffdc24c524dae6331565fc579a3 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