OPTION, PSDUMPFREQ = 50;   // 6d data written every 10 time steps (h5).
OPTION, STATDUMPFREQ = 100;  // Beam Stats written every 10 time steps (stat).
OPTION, BOUNDPDESTROYFQ = 10; // Delete lost particles, if any       
OPTION, AUTOPHASE = 0;
OPTION, VERSION = 10900;

//----------------------------------------------------------------------------
//Global Parameters
//----------------------------------------------------------------------------

REAL n_particles=1E4;      // Number of particles in simulation. (Overridden if distibution imported)

// Initial Momentum (Set very small when starting simulation from cathode emission)
REAL P0=0.01;    //inital z momentum



//----------------------------------------------------------------------------
// Injector Accelerating Cavities
//----------------------------------------------------------------------------

REAL hwl = 0.052464; //half-wavelength of cavity



L1:     TRAVELINGWAVE, L = 3.01, VOLT = 10.0, ELEMEDGE = hwl + 0.5,
        NUMCELLS = 84, MODE = 1/3,
        FAST = False, 
        FMAPFN = "L1_Cavity1.T7", 
        FREQ = 2856.0, 
        LAG = PI/4.; // LAG (phase) is converted to radians\



fit_point_beam: MONITOR, OUTFN = "fit_point_bunch.h5", 
    ELEMEDGE = L1->ELEMEDGE + L1->L + 0.5;



PCG2BCE: Line = (L1,fit_point_beam); // Ends shortly after BC1


// Cold Gaussian distribution for testing purposes
gaussian: DISTRIBUTION, TYPE = GAUSS,
        SIGMAX = 0.00076696498249,
        SIGMAPX = 0.0013038404918,
        SIGMAY = 0.00076696498249,
        SIGMAPY = 0.0013038404918,
        SIGMAZ = 0.000050,
        SIGMAPZ = 0.0,   
        EMITTED = False,
        INPUTMOUNITS = NONE,      
        WRITETOFILE = True; 
 
//----------------------------------------------------------------------------
// Define Field solvers
//----------------------------------------------------------------------------

FS_SC: Fieldsolver, FSTYPE = None, // None or FFT
            MX = 32, MY = 32, MT = 32,
            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 = 2856.0, BCURRENT = 0.0, CHARGE = -1;

//----------------------------------------------------------------------------
// Initialize Simulation
//----------------------------------------------------------------------------


TRACK, LINE = PCG2BCE, BEAM = BEAM1, MAXSTEPS = 2500000, 
    DT = {10.0e-13}, ZSTOP={5.0};

RUN, METHOD = "PARALLEL-T", BEAM = BEAM1, 
    FIELDSOLVER = FS_SC, DISTRIBUTION = gaussian;
ENDTRACK;

Stop;
Quit;