...
 
Commits (314)
......@@ -149,19 +149,6 @@ RegressionTests/EGunCTF3-2/reference/EGunCTF3-2.out -text
RegressionTests/EGunCTF3-2/reference/EGunCTF3-2.out.md5 -text
RegressionTests/EGunCTF3-2/reference/EGunCTF3-2.stat -text
RegressionTests/EGunCTF3-2/reference/EGunCTF3-2.stat.md5 -text
RegressionTests/Envelope-Tracker-Phase1/Envelope-Tracker-Phase1.in -text
RegressionTests/Envelope-Tracker-Phase1/Envelope-Tracker-Phase1.local -text
RegressionTests/Envelope-Tracker-Phase1/Envelope-Tracker-Phase1.rt -text
RegressionTests/Envelope-Tracker-Phase1/Envelope-Tracker-Phase1.sge -text
RegressionTests/Envelope-Tracker-Phase1/FIND1-MSOL10.T7 -text
RegressionTests/Envelope-Tracker-Phase1/FINSS-RGUN.dat -text
RegressionTests/Envelope-Tracker-Phase1/FinPhase1.phases -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.lbal -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.lbal.md5 -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.out -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.out.md5 -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.stat -text
RegressionTests/Envelope-Tracker-Phase1/reference/Envelope-Tracker-Phase1.stat.md5 -text
RegressionTests/ExternalFieldTest/CTF3_Ez_ASTRA.opal -text
RegressionTests/ExternalFieldTest/ExternalFieldTest.in -text
RegressionTests/ExternalFieldTest/ExternalFieldTest.local -text
......
# use >git clean -f to get rid of all untracked files
.keytab
*~
*.pyc
*.png
*_DesignPath.dat
*_ElementPositions.py
*_ElementPositions.sdds
*_ElementPositions.txt
*_Monitors.stat
*_3D.opal
*_DesignPath.dat
*_ElementPositions.py
*_ElementPositions.sdds
*_ElementPositions.txt
*_Monitors.stat
*_3D.opal
# For OPAL-maps
*.dispersion
*.map
RegressionTests/*/*-RT.o
RegressionTests/*/*.lbal
RegressionTests/*/*.out
RegressionTests/*/*.stat
RegressionTests/*/AMF*.h5
RegressionTests/*/KMA*.h5
RegressionTests/*/KMB*.h5
RegressionTests/*/M1*.h5
RegressionTests/*/MM*.h5
RegressionTests/*/MMF*.h5
RegressionTests/*/QMA*.h5
RegressionTests/*/QMF*.h5
RegressionTests/*/data
RegressionTests/*/errormsg.txt
RegressionTests/*/test.dat
RegressionTests/*/timing.dat
RegressionTests/results.xml
RegressionTests/results/
RegressionTests/timing.dat
RegressionTests/AWAGun-1/AWAGun-1.h5
RegressionTests/AWAGun-TrackBack-1/AWAGun-TrackBack-1.h5
RegressionTests/BC1-1_fixme/BC1-1.h5
RegressionTests/BeamLine-1/BeamLine-1.h5
RegressionTests/BeamLine-2/BeamLine-2.h5
RegressionTests/CSRBendDrift/CSRBendDrift.h5
RegressionTests/DCGun-Cornell-1/DCGun-Cornell-1.h5
RegressionTests/DCGun-Cornell-2/DCGun-Cornell-2.h5
RegressionTests/Degrader-1/Degrader-1.h5
RegressionTests/Degrader-1/M*.h5
RegressionTests/Degrader-2/Degrader-2.h5
RegressionTests/Distribution-Binomial-1/Distribution-Binomial-1.h5
RegressionTests/Distribution-Gauss-1/Distribution-Gauss-1.h5
RegressionTests/Distribution-Gauss-2/Distribution-Gauss-2.h5
RegressionTests/EGunCTF3-1/EGunCTF3-1.h5
RegressionTests/EGunCTF3-2/EGunCTF3-2.h5
RegressionTests/Elliptic-1/ELCO.h5
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/LaserEmission-1/LaserEmission-1.h5
RegressionTests/MAP-Circular/MAP-Circular.h5
RegressionTests/MAP-FODO/MAP-FODO.h5
RegressionTests/MAP-FODODipole/MAP-FODODipole.h5
RegressionTests/MAP-Gantry2/MAP-Gantry2.h5
RegressionTests/OpalRingTest/OpalRingTest-trackOrbit.dat
RegressionTests/OpalRingTest/OpalRingTest.h5
RegressionTests/OpalRingTest/PROBE1.hist
RegressionTests/OpalRingTest/PROBE1.loss
RegressionTests/OpalRingTest/PROBE1.loss
RegressionTests/OpalRingTest/PROBE1.peaks
RegressionTests/OpalRingTest/RINGDEF.loss
RegressionTests/OpalRingTest/RINGDEF.loss
RegressionTests/PROSCAN-1/PROSCAN-1.h5
RegressionTests/PROSCAN-2/PROSCAN-2.h5
RegressionTests/PROSCAN-3/PROSCAN-3.h5
RegressionTests/PROSCAN-4/PROSCAN-4.h5
RegressionTests/PSIGUN-1/PSIGUN-1.h5
RegressionTests/ParallelPlateFurman1/ParallelPlateFurman1.SurfaceLoss.h5
RegressionTests/ParallelPlateFurman1/ParallelPlateFurman1.h5
RegressionTests/ParallelPlateVaughan1/ParallelPlateVaughan1.h5
RegressionTests/PerfectDiode/PerfectDiode.h5
RegressionTests/Quad-Simple-Test-1/Quad-Simple-Test-1.h5
RegressionTests/Quad-Simple-Test-1/Quad-Simple-Test-1.lbal
RegressionTests/RBend-Simple-Test-1/RBend-Simple-Test-1.h5
RegressionTests/RBend-Simple-Test-2/RBend-Simple-Test-2.h5
RegressionTests/RBend-Simple-Test-3/RBend-Simple-Test-3.h5
RegressionTests/RBend-Simple-Test-4/RBend-Simple-Test-4.h5
RegressionTests/RBend-Simple-Test-4/RBend-Simple-Test-4.h5
RegressionTests/RestartTest-1/RestartTest-1.h5
RegressionTests/RestartTest-2/RestartTest-2-trackOrbit.dat
RegressionTests/RestartTest-2/RestartTest-2.h5
RegressionTests/RestartTest-3/RestartTest-3.h5
RegressionTests/RestartTest-4/RestartTest-4.h5
RegressionTests/RestartTest-5/RestartTest-5.h5
RegressionTests/RestartTest-6/RestartTest-6.h5
RegressionTests/RestartTest-6/RestartTest-6.h5
RegressionTests/RestartTest-6/RestartTest-6.h5.ref
RegressionTests/RestartTest-6/RestartTest-6.lbal.ref
RegressionTests/RestartTest-6/RestartTest-6.stat.md5
RegressionTests/RestartTest-6/RestartTest-6.stat.ref
RegressionTests/RestartTest-7/RestartTest-7.h5
RegressionTests/RestartTest-8/RestartTest-8.h5
RegressionTests/RingCyclotron/RingCyclotron-trackOrbit.dat
RegressionTests/RingCyclotron/RingCyclotron.h5
RegressionTests/RingCyclotronMTS/RingCyclotronMTS-trackOrbit.dat
RegressionTests/RingCyclotronMTS/RingCyclotronMTS.h5
RegressionTests/RingCyclotronMatched/RingCyclotronMatched-trackOrbit.dat
RegressionTests/RingCyclotronMatched/RingCyclotronMatched.h5
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle-Angle0.dat
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle-Angle1.dat
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle-Angle2.dat
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle-afterEachTurn.dat
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle-trackOrbit.dat
RegressionTests/RingCyclotronSingleParticle/RingCyclotronSingleParticle.h5
RegressionTests/SAAMG-Test-1/SAAMG-Test-1.h5
RegressionTests/SAAMG-Test-1/timing_MX16_MY16_MZ16_nProc1_recB0_numB0_nLHS1
RegressionTests/SBend-Simple-Test-1/SBend-Simple-Test-1.h5
RegressionTests/SBend-Simple-Test-2/SBend-Simple-Test-2.h5
RegressionTests/SBend-Simple-Test-3/SBend-Simple-Test-3.h5
RegressionTests/SBend-Simple-Test-4/SBend-Simple-Test-4.h5
RegressionTests/Scan-1/Scan-1.h5
RegressionTests/Slit-1/MSLIT1.h5
RegressionTests/Slit-1/SLIT1.h5
RegressionTests/Slit-1/Slit-1.h5
RegressionTests/Slit-2/MSLIT1.h5
RegressionTests/Slit-2/SLIT1.h5
RegressionTests/Slit-2/Slit-2.h5
RegressionTests/VKick-Test-1/VKick-Test-1.h5
RegressionTests/opalimpact-1/opalimpact-1.h5
RegressionTests/VFFA-1/RINGPROBE*
RegressionTests/VFFA-1/*.dat
RegressionTests/VFFA-1/*.h5
OPTION, PSDUMPFREQ = 300000; // 6d data written every 300000 time steps (h5).
OPTION, STATDUMPFREQ = 10; // Beam Stats written every 10 time steps (stat).
OPTION, BOUNDPDESTROYFQ=10; // Delete lost particles, if any
OPTION, AUTOPHASE=4; // Autophase is on, and phase of max energy
// gain will be found automatically for cavities
Option, VERSION=10900;
Title, string="AWA Photoinjector";
//----------------------------------------------------------------------------
//Global Parameters
REAL rf_freq = 1.3e3; //RF frequency. (Hz)
REAL n_particles = 1E4; //Number of particles in simulation.
REAL beam_bunch_charge = 1e-9; //Charge of bunch. (C)
//Initial Momentum Calculation
REAL Edes = 1.4e-9; //initial energy in GeV
REAL gamma = (Edes+EMASS)/EMASS;
REAL beta = sqrt(1-(1/gamma^2));
REAL P0 = gamma*beta*EMASS; //inital z momentum
//Printing initial energy and momentum to terminal output.
value , {Edes, P0, OPALVERSION};
//----------------------------------------------------------------------------
// RF photoinjector (Gun)
//
// Cavity/RF field.
//
// L: physical element length (real in m). Length (of field map) (m).
// VOLT: field scaling factor (real). RF field magnitude (MV/m).
// FMAPFN: file name, or full path name to file (string)
// ELEMEDGE: physical start of the element on the floor (m)
// TYPE: specifies "STANDING", "TRAVELLING" or "SINGLE GAP" structure
// FREQ: RF frequency of cavity (real in MHz). Resonance frequency.
// LAG: cavity phase (radians)
//
REAL gun_inj_phase = 0.0; // Rf injection phase (in degrees) with restpect
// to max phase. i.e. 0.0 inj phase corresponds
// to running on crest (max energy).
GUNSOURCE: SOURCE, ELEMEDGE = 0.0;
GUN: RFCavity, L = 0.2927, VOLT = 60.0, ELEMEDGE = 0.0, TYPE = "STANDING",
FMAPFN = "DriveGun.T7", FREQ = rf_freq, LAG = (gun_inj_phase*Pi)/180.0;
//LAG (phase) is converted to radians
//----------------------------------------------------------------------------
// Solenoids
//
// L: Physcial element length (m)
// ELEMEDGE: Physcial start of element (m)
// KS: Solenoid strength (Rad/m)
// FMAPFM: Field file (string)
// Note: OPAL scales the field file based on the max magnetic
// field value in the file, not Bz on axis. The max field
// value is normalized to 1 [T], and scaled with KS.
// i.e. The max value in the BF_559 file = 0.162544398 [T].
// Therefore, setting KS = 0.162544398 runs the magnet at max current.
REAL KSBF = 0.162544398;
if (OPALVERSION>10500)
KSBF = KSBF/1.3528;
// Buck and Focusing Solenoid
BF: Solenoid, L = 0.5, ELEMEDGE=0.0, KS = KSBF,
FMAPFN = "BF_550.T7";
// Matching Solenoid
// Max value in the M_440 file = 1.973966
// Conversion from Current (I) to Scaling factor:
// Scaling Factor (KS) = (I / Imax) * 1.973966
// Where Imax for this Matching solenoid = 440 [A]
REAL I = 273; //[A]
REAL SF = (I/440)*1.973966;
if (OPALVERSION>10500)
SF = SF/3.2306;
M: Solenoid, L = 0.5, ELEMEDGE=0.0, KS = SF,
FMAPFN = "M_440.T7";
value,{KSBF,SF};
//----------------------------------------------------------------------------
//Drift after gun.
DR1: DRIFT, L = 10.0, ELEMEDGE = 0.5;
//----------------------------------------------------------------------------
// DEFINE BEAM LINE
GS: Line = (GUNSOURCE, GUN, BF, M);
// Complete accelerator
DRIVE: Line = (GS, DR1);
//----------------------------------------------------------------------------
// INITIAL DISTRIBUTION
//
// Flattop distribution.
// SIGMAX/Y: RMS radius of transverse beam size (Laser radius in m).
// TRISE/FALL: Rise time/fall time in longitudinal direction (s).
// TPULSEFWHM: FWHM in longitudinal direction (s).
// CUTOFFLONG: Longitudinal cuttoff in units of sigma.
// NBIN: Number of energy bins to use during emission.
// EMISSIONSTEPS: Number of steps during emssion.
// Emission time step is adjusted to fit this number.
// EKIN: Kinetic energy of electrons at emission (eV).
// ELASER: Energy of laser (eV).
// W: Photocathode work functioin (eV).
// FE: Fermi energy of photocathode (eV).
// CATHTEMP: Operating temperature of photocathode (K).
// Note, ELASER, W, FE, and CATHTEMP are used for the NONEQUIL emission model.
// These values are not necessary when using other models.
// Note, If you want a Gaussian, in the longitudinal direction:
// TRISE/FALL = 1.6869*simgar
// Sigmar = FWHM / 2.35 (this equation is only valid for 1.3 GHz)
Dist: DISTRIBUTION, TYPE = FLATTOP,
SIGMAX = 0.00075,
SIGMAY = 0.00075,
TRISE = 6.0e-12,
TFALL = 6.0e-12,
TPULSEFWHM = 20.0e-12,
CUTOFFLONG = 4.0,
NBIN = 9,
EMISSIONSTEPS = 100,
EMISSIONMODEL = ASTRA,
EKIN = 0.55,
EMITTED = True,
WRITETOFILE = True; //Saves the distribution to a text file
// Note on emission time step: FWHM pulse width divided by emission
// steps gives the time step for the emissions process.
// i.e 20.0e-12 / 100 gives a time step of 2e-13 (s) during emission.
// This is the not the same as the time step used in rest of the file.
//----------------------------------------------------------------------------
// Define Field solvers
// The mesh sizes should be a factor of 2
// for most efficient space charge calculation.
FS_SC: Fieldsolver, FSTYPE = FFT,
MX = 16, MY = 16, MT = 16, // SC grid size is 32^3
PARFFTX = false,
PARFFTY = false,
PARFFTT = true, // parallel in the z direction only
BCFFTX = open,
BCFFTY = open,
BCFFTT = open,
BBOXINCR = 1,
GREENSF = INTEGRATED;
//----------------------------------------------------------------------------
// Electron Beam Definition
BEAM1: BEAM, PARTICLE = ELECTRON, pc = P0, NPART = n_particles,
BFREQ = rf_freq, BCURRENT = beam_bunch_charge * rf_freq * 1e6, CHARGE = -1;
//----------------------------------------------------------------------------
// Simulate the beamline using TRACK and RUN.
// Note, different time steps are set based on the z location in the beam line.
// In the case below, 1.0e-13 is used for 0.0 to 0.4 m,
// and 3.0e-12 is used from 0.4 to 5 m.
TRACK, LINE = DRIVE, BEAM = BEAM1, MAXSTEPS = 19000,
DT = {1.0e-13, 3.0e-12}, ZSTOP={0.4, 2.0};
RUN, METHOD = "PARALLEL-T", BEAM = BEAM1,
FIELDSOLVER = FS_SC, DISTRIBUTION = Dist;
ENDTRACK;
Stop;
Quit;
\ No newline at end of file
#!/bin/bash
mpirun -np 4 "$OPAL_EXE_PATH/opal" AWAGun-1.in --info 2 "$@" 2>&1
"AWAGun-1"
stat "rms_x" last 1E-5 #this is a comment
stat "rms_y" last 1E-5 #this is a comment
stat "rms_s" last 1E-5 #this is a comment
stat "emit_x" last 1E-5 #this is a comment
stat "emit_y" last 1E-5 #this is a comment
stat "emit_s" last 1E-5 #this is a comment
stat "energy" last 1E-5 #this is a comment
#!/bin/bash
#$ -cwd
#$ -j y
#$ -pe orte 1
#$ -N BeamLine-1-RT
#$ -v LD_LIBRARY_PATH,OPAL_EXE_PATH,OPENMPI,REG_TEST_DIR
MACHINE_FILE=$TMPDIR/machinefile
awk '/^merlin/ {print $1" slots="$2}' $PE_HOSTFILE > $MACHINE_FILE
cp $MACHINE_FILE machinefile.last
cd $REG_TEST_DIR
OPAL="$OPAL_EXE_PATH/opal AWAGun-1.in --commlib mpi --info 3 --warn 0 2>&1"
CMD="$OPENMPI/bin/mpirun -x LD_LIBRARY_PATH -machinefile $MACHINE_FILE -np $NSLOTS $OPAL "
$CMD
\ No newline at end of file
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
57cdf40cbfaff54774103f708f3c6820 AWAGun-1.lbal
This source diff could not be displayed because it is too large. You can view the blob instead.
6e870b6499af82692a857a6b81d62566 AWAGun-1.out
This diff is collapsed.
1e9d085cc86605a51cdd6203cddd17e2 AWAGun-1.stat
OPTION, PSDUMPFREQ = 300000; // 6d data written every 300000 time steps (h5).
OPTION, STATDUMPFREQ = 10; // Beam Stats written every 10 time steps (stat).
OPTION, BOUNDPDESTROYFQ=10; // Delete lost particles, if any
OPTION, AUTOPHASE=4; // Autophase is on, and phase of max energy
// gain will be found automatically for cavities
Option, VERSION=10900;
Title, string="AWA Photoinjector";
//----------------------------------------------------------------------------
//Global Parameters
REAL rf_freq = 1.3e3; //RF frequency. (Hz)
REAL n_particles = 1E4; //Number of particles in simulation.
REAL beam_bunch_charge = 1e-9; //Charge of bunch. (C)
//Initial Momentum Calculation
REAL Edes = 1.4e-9; //initial energy in GeV
REAL gamma = (Edes+EMASS)/EMASS;
REAL beta = sqrt(1-(1/gamma^2));
REAL P0 = gamma*beta*EMASS; //inital z momentum
//Printing initial energy and momentum to terminal output.
value , {Edes, P0, OPALVERSION};
//----------------------------------------------------------------------------
// RF photoinjector (Gun)
//
// Cavity/RF field.
//
// L: physical element length (real in m). Length (of field map) (m).
// VOLT: field scaling factor (real). RF field magnitude (MV/m).
// FMAPFN: file name, or full path name to file (string)
// ELEMEDGE: physical start of the element on the floor (m)
// TYPE: specifies "STANDING", "TRAVELLING" or "SINGLE GAP" structure
// FREQ: RF frequency of cavity (real in MHz). Resonance frequency.
// LAG: cavity phase (radians)
//
REAL gun_inj_phase = 0.0; // Rf injection phase (in degrees) with restpect
// to max phase. i.e. 0.0 inj phase corresponds
// to running on crest (max energy).
GUNSOURCE: SOURCE, ELEMEDGE = 0.0;
GUN: RFCavity, L = 0.2927, VOLT = 60.0, ELEMEDGE = 0.0, TYPE = "STANDING",
FMAPFN = "DriveGun.T7", FREQ = rf_freq, LAG = (gun_inj_phase*Pi)/180.0;
//LAG (phase) is converted to radians
//----------------------------------------------------------------------------
// Solenoids
//
// L: Physcial element length (m)
// ELEMEDGE: Physcial start of element (m)
// KS: Solenoid strength (Rad/m)
// FMAPFM: Field file (string)
// Note: OPAL scales the field file based on the max magnetic
// field value in the file, not Bz on axis. The max field
// value is normalized to 1 [T], and scaled with KS.
// i.e. The max value in the BF_559 file = 0.162544398 [T].
// Therefore, setting KS = 0.162544398 runs the magnet at max current.
REAL KSBF = 0.162544398;
if (OPALVERSION>10500)
KSBF = KSBF/1.3528;
// Buck and Focusing Solenoid
BF: Solenoid, L = 0.5, ELEMEDGE=0.0, KS = KSBF,
FMAPFN = "BF_550.T7";
// Matching Solenoid
// Max value in the M_440 file = 1.973966
// Conversion from Current (I) to Scaling factor:
// Scaling Factor (KS) = (I / Imax) * 1.973966
// Where Imax for this Matching solenoid = 440 [A]
REAL I = 273; //[A]
REAL SF = (I/440)*1.973966;
if (OPALVERSION>10500)
SF = SF/3.2306;
M: Solenoid, L = 0.5, ELEMEDGE=0.0, KS = SF,
FMAPFN = "M_440.T7";
value,{KSBF,SF};
//----------------------------------------------------------------------------
//Drift after gun.
DR1: DRIFT, L = 10.0, ELEMEDGE = 0.5;
//----------------------------------------------------------------------------
// DEFINE BEAM LINE
GS: Line = (GUNSOURCE, GUN, BF, M);
// Complete accelerator
DRIVE: Line = (GS, DR1);
//----------------------------------------------------------------------------
// INITIAL DISTRIBUTION
//
// Flattop distribution.
// SIGMAX/Y: RMS radius of transverse beam size (Laser radius in m).
// TRISE/FALL: Rise time/fall time in longitudinal direction (s).
// TPULSEFWHM: FWHM in longitudinal direction (s).
// CUTOFFLONG: Longitudinal cuttoff in units of sigma.
// NBIN: Number of energy bins to use during emission.
// EMISSIONSTEPS: Number of steps during emssion.
// Emission time step is adjusted to fit this number.
// EKIN: Kinetic energy of electrons at emission (eV).
// ELASER: Energy of laser (eV).
// W: Photocathode work functioin (eV).
// FE: Fermi energy of photocathode (eV).
// CATHTEMP: Operating temperature of photocathode (K).
// Note, ELASER, W, FE, and CATHTEMP are used for the NONEQUIL emission model.
// These values are not necessary when using other models.
// Note, If you want a Gaussian, in the longitudinal direction:
// TRISE/FALL = 1.6869*simgar
// Sigmar = FWHM / 2.35 (this equation is only valid for 1.3 GHz)
Dist: DISTRIBUTION, TYPE = FROMFILE,
FNAME = "DistAt2m.dat";
// Note on emission time step: FWHM pulse width divided by emission
// steps gives the time step for the emissions process.
// i.e 20.0e-12 / 100 gives a time step of 2e-13 (s) during emission.
// This is the not the same as the time step used in rest of the file.
//----------------------------------------------------------------------------
// Define Field solvers
// The mesh sizes should be a factor of 2
// for most efficient space charge calculation.
FS_SC: Fieldsolver, FSTYPE = FFT,
MX = 16, MY = 16, MT = 16, // SC grid size is 32^3
PARFFTX = false,
PARFFTY = false,
PARFFTT = true, // parallel in the z direction only
BCFFTX = open,
BCFFTY = open,
BCFFTT = open,
BBOXINCR = 1,
GREENSF = INTEGRATED;
//----------------------------------------------------------------------------
// Electron Beam Definition
BEAM1: BEAM, PARTICLE = ELECTRON, pc = P0, NPART = n_particles,
BFREQ = rf_freq, BCURRENT = beam_bunch_charge * rf_freq * 1e6, CHARGE = -1;
//----------------------------------------------------------------------------
// Simulate the beamline using TRACK and RUN.
// Note, different time steps are set based on the z location in the beam line.
// In the case below, 1.0e-13 is used for 0.0 to 0.4 m,
// and 3.0e-12 is used from 0.4 to 5 m.
TRACK, LINE = DRIVE, BEAM = BEAM1, MAXSTEPS = {17200, 1800},
DT = {1.0e-13, 3.0e-12}, ZSTART = 2.0, ZSTOP = {0.0, 0.4};
RUN, METHOD = "PARALLEL-T", BEAM = BEAM1,
FIELDSOLVER = FS_SC, DISTRIBUTION = Dist, TRACKBACK = TRUE;
ENDTRACK;
Stop;
Quit;
\ No newline at end of file
#!/bin/bash
mpirun -np 4 $OPAL_EXE_PATH/opal AWAGun-TrackBack-1.in --info 2 "$@" 2>&1
"Track AWAGun-1 backward"
stat "rms_x" avg 1E-5 #this is a comment
stat "rms_y" avg 1E-5 #this is a comment
stat "rms_s" avg 1E-5 #this is a comment
stat "emit_x" avg 1E-5 #this is a comment
stat "emit_y" avg 1E-5 #this is a comment
stat "emit_s" avg 1E-5 #this is a comment
stat "energy" avg 1E-5 #this is a comment
......@@ -2,14 +2,14 @@
#$ -cwd
#$ -j y
#$ -pe orte 1
#$ -N Envelope-Tracker-Phase1-RT
#$ -v LD_LIBRARY_PATH,OPENMPI,OPAL_EXE_PATH,REG_TEST_DIR
#$ -N BeamLine-1-RT
#$ -v LD_LIBRARY_PATH,OPAL_EXE_PATH,OPENMPI,REG_TEST_DIR
MACHINE_FILE=$TMPDIR/machinefile
awk '/^merlin/ {print $1" slots="$2}' $PE_HOSTFILE > $MACHINE_FILE
awk '/^merlin/ {print $1" slots="$2}' $PE_HOSTFILE > $MACHINE_FILE
cp $MACHINE_FILE machinefile.last
cd $REG_TEST_DIR
OPAL="$OPAL_EXE_PATH/opal Envelope-Tracker-Phase1.in --commlib mpi --info 3 --warn 0 "
OPAL="$OPAL_EXE_PATH/opal AWAGun-TrackBack-1.in --commlib mpi --info 3 --warn 0 2>&1"
CMD="$OPENMPI/bin/mpirun -x LD_LIBRARY_PATH -machinefile $MACHINE_FILE -np $NSLOTS $OPAL "
$CMD
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
328817d5c7f70644001d89b2716e1aa7 AWAGun-TrackBack-1.lbal
d5e1ffaa311a109ef9feff9db8b8ffc0 AWAGun-TrackBack-1.out
7e759719016b954d999c7c108ec96a32 AWAGun-TrackBack-1.stat