Commit d7eed4b6 authored by Pedro Calvo Portela's avatar Pedro Calvo Portela

Merge branch 'GasStripping-regression-test-into-master' into 'master'

Gas stripping regression test into master

See merge request !8
parents 5130d1f2 654f61b8
......@@ -62,6 +62,8 @@ RegressionTests/Elliptic-1/Elliptic-1.h5
RegressionTests/Elliptic-1/MELCO.h5
RegressionTests/Envelope-Tracker-Phase1/Envelope-Tracker-Phase1.h5
RegressionTests/ExternalFieldTest/ExternalFieldTest.h5
RegressionTests/GasStripping/GasStripping-trackOrbit.dat
RegressionTests/GasStripping/GasStripping.h5
RegressionTests/HKick-Test-1/HKick-Test-1.h5
RegressionTests/HKick-Test-2/HKick-Test-2.h5
RegressionTests/Kick-Test-1/Kick-Test-1.h5
......
Title,string="OPAL-cycl: Gas stripping in a drift";
OPTION, VERSION=20000;
OPTION, TRACE=FALSE;
OPTION, PSDUMPFREQ = 10000;
OPTION, STATDUMPFREQ = 1;
OPTION, SPTDUMPFREQ = 10000;
OPTION, REPARTFREQ = 1;
OPTION, REBINFREQ = 1;
OPTION, BOUNDPDESTROYFQ =1;
OPTION, ENABLEHDF5=FALSE;
OPTION, ASCIIDUMP=FALSE;
REAL q_e = -1.60217733e-19;
REAL nstep = 1000;
REAL turns = 0.723;
REAL npart= 100000;
REAL Edes = 1E-5;
REAL frequency = 1.0;
REAL phase = 0.0*DEGRAD;
REAL gamma = (Edes+HMMASS)/HMMASS;
REAL beta = sqrt(1-(1/gamma^2));
REAL P0 = gamma*beta*HMMASS;
CyclDrift: Cyclotron, TYPE="BANDRF", CYHARMON=1, RINIT=0.0, PHIINIT=0.0, ZINIT=0.0, PRINIT=0.0,
MAXZ=+50, MINZ=-50, MINR=0, MAXR=1100, SYMMETRY=1,
FMAPFN="./BField_0T_const.dat";
bstp_phys:PARTICLEMATTERINTERACTION, TYPE="BEAMSTRIPPING";
bstp: BEAMSTRIPPING, PRESSURE=1.0e-6, TEMPERATURE=300, GAS="AIR", stop=TRUE, PARTICLEMATTERINTERACTION=bstp_phys;
C1: LINE = (CyclDrift, bstp);
DIST1:DISTRIBUTION, TYPE=GAUSS,
SIGMAX= 1.0e-03, SIGMAPX=0.0, CORRX=0.0,
SIGMAY= 1.0e-03, SIGMAPY=0.0, CORRY=0.0,
SIGMAT= 1.0e-03,
SIGMAPT=0.0, CORRT=0.0;
Fs1:FIELDSOLVER, FSTYPE=NONE, MX=8, MY=8, MT=8,
PARFFTX=TRUE, PARFFTY=TRUE, PARFFTT=TRUE,
ITSOLVER=BiCGSTAB, INTERPL=CONSTANT, TOL=1.0e-6, MAXITERS=1000, PRECMODE=STD,
BCFFTX=open, BCFFTY=open, BCFFTT=open,
BBOXINCR=1.0, GREENSF=INTEGRATED;
BEAM1: BEAM, PARTICLE=HMINUS, MASS=HMMASS, CHARGE=-1, BFREQ=frequency, PC=P0, NPART=npart, BCURRENT=npart*q_e*frequency*1E6;
SELECT, LINE=C1;
TRACK, LINE=C1, BEAM=BEAM1, MAXSTEPS=nstep*turns, STEPSPERTURN=nstep, TIMEINTEGRATOR="RK-4";
RUN, METHOD="CYCLOTRON-T", BEAM=BEAM1, FIELDSOLVER=Fs1, DISTRIBUTION=Dist1, TURNS=1;
ENDTRACK;
Stop;
#!/bin/bash
mpirun -np 4 $OPAL_EXE_PATH/opal GasStripping.in --info 2 --warn 0 2>&1
"Drift tracking with beam stripping"
stat "numParticles" last 80
This diff is collapsed.
762e650bb90d5aecb9d240623750430c GasStripping.lbal
This diff is collapsed.
5bc065c717f4afb5aaf3b671b0552d66 GasStripping.out
This diff is collapsed.
f81fd39b462395248fe5f194fd2eb6f2 GasStripping.stat
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