Gas stripping regression test into master

.gitignore

@@ -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

RegressionTests/GasStripping/GasStripping.in

+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;

RegressionTests/GasStripping/GasStripping.local

+#!/bin/bash
+
+mpirun -np 4 $OPAL_EXE_PATH/opal GasStripping.in --info 2 --warn 0 2>&1

RegressionTests/GasStripping/GasStripping.rt

+"Drift tracking with beam stripping"
+stat "numParticles" last 80

RegressionTests/GasStripping/reference/GasStripping.lbal.md5

+762e650bb90d5aecb9d240623750430c  GasStripping.lbal

RegressionTests/GasStripping/reference/GasStripping.out.md5

+5bc065c717f4afb5aaf3b671b0552d66  GasStripping.out

RegressionTests/GasStripping/reference/GasStripping.stat.md5

+f81fd39b462395248fe5f194fd2eb6f2  GasStripping.stat