plasma_osc { Plasma oscillation Electrons are given a kick in the beginning by sepcifying a spatially varying initial electric field. The field solve is electromagnetic. If the electrons are non-relativistic (as is the case here), then it is faster to use the electrostatic approximation (ElectrostaticFlag = 4; see description of the input parameters in the manual file) } Variables { // General numerical parameters PI = 3.14159 // ********************************************************************** // General physical parameters // ********************************************************************** electronMassMKS = 9.1094e-31 electronCharge = -1.6022e-19 speedLight = 2.9979e8 epsilon0MKS = 8.854e-12 electronDensityMKS = 5e23 plasmaFrequency = sqrt(electronCharge*electronCharge*electronDensityMKS/(electronMassMKS*epsilon0MKS)) ds = 0.05 * speedLight / plasmaFrequency // ********************************************************************** // Grid parameters // ********************************************************************** Nx = 256 Ny = 32 dx = ds dy = dx Lx = Nx * dx Ly = Ny * dy simulationVolume = Lx * Ly numCells = Nx * Ny // // parameters for electrons // totalNumElectrons = electronDensityMKS * simulationVolume numElectronsPerCell = 100 numElectronPtcls = numElectronsPerCell * numCells electronNumRatio = totalNumElectrons / numElectronPtcls //the following will be useful for the two-stream instability problem rmsElectronSpeedMKS = 3.e+04 electronVelocityMKS = 1.e+08 peakCurrentElectrons = electronCharge * electronDensityMKS * Ly * electronVelocityMKS d = 1. / sqrt( 1./(dx*dx) + 1./(dy*dy) ) timeStep = 0.9 * d / speedLight wl=1 } Region { Grid { J = Nx x1s = 0.0 x1f = Lx n1 = 1.0 K = Ny x2s = 0.0 x2f = Ly n2 = 1.0 Geometry = 1 // 2D (x-y) slab geometry PeriodicFlagX1 = 1 // periodic in x1 (x) PeriodicFlagX2 = 1 // periodic in x2 (y) } Control { dt = timeStep // Turn off the initial Poisson solve initPoissonSolve = 0 // E1init=1e5*sin(2 * 3.14 * x1 / (.5*Lx) ) //this is the initial kick electric field, which is not used in the calculation beyond the first step. } Species { name = electrons m = electronMassMKS q = electronCharge } Species { name = ions m = electronMassMKS*1800 q = -electronCharge } // Load the right-going plasma electrons over the entire simulation region //Load //{ // speciesName = electrons // density = electronDensityMKS //if you want to play with the nonuniform density of electrons, can turn this on // analyticF = (1 + 0.1*sin(2 * 3.14 * x1 / (.5*Lx) )) * electronDensityMKS //however,then you need to turn on initPoissonSolve in the Control block // x1MinMKS = 0.0 // x1MaxMKS = Lx // x2MinMKS = 0.0 // x2MaxMKS = Ly // np2c = electronNumRatio // specify MKS units for all velocities // units = MKS // v1drift = 0*electronVelocityMKS // v1thermal = 0*rmsElectronSpeedMKS // v2thermal = 0. // v3thermal = 0. //} Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 0 x1MaxMKS = Lx/ 100* 1 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 0.500000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 1 x1MaxMKS = Lx/ 100* 2 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 1.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 2 x1MaxMKS = Lx/ 100* 3 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 2.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 3 x1MaxMKS = Lx/ 100* 4 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 3.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 4 x1MaxMKS = Lx/ 100* 5 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 4.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 5 x1MaxMKS = Lx/ 100* 6 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 5.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 6 x1MaxMKS = Lx/ 100* 7 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 6.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 7 x1MaxMKS = Lx/ 100* 8 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 7.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 8 x1MaxMKS = Lx/ 100* 9 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 8.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 9 x1MaxMKS = Lx/ 100* 10 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 9.50000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 10 x1MaxMKS = Lx/ 100* 11 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 10.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 11 x1MaxMKS = Lx/ 100* 12 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 11.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 12 x1MaxMKS = Lx/ 100* 13 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 12.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 13 x1MaxMKS = Lx/ 100* 14 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 13.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 14 x1MaxMKS = Lx/ 100* 15 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 14.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 15 x1MaxMKS = Lx/ 100* 16 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 15.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 16 x1MaxMKS = Lx/ 100* 17 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 16.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 17 x1MaxMKS = Lx/ 100* 18 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 17.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 18 x1MaxMKS = Lx/ 100* 19 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 18.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 19 x1MaxMKS = Lx/ 100* 20 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 19.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 20 x1MaxMKS = Lx/ 100* 21 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 20.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 21 x1MaxMKS = Lx/ 100* 22 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 21.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 22 x1MaxMKS = Lx/ 100* 23 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 22.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 23 x1MaxMKS = Lx/ 100* 24 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 23.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 24 x1MaxMKS = Lx/ 100* 25 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 24.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 25 x1MaxMKS = Lx/ 100* 26 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 25.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 26 x1MaxMKS = Lx/ 100* 27 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 26.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 27 x1MaxMKS = Lx/ 100* 28 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 27.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 28 x1MaxMKS = Lx/ 100* 29 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 28.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 29 x1MaxMKS = Lx/ 100* 30 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 29.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 30 x1MaxMKS = Lx/ 100* 31 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 30.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 31 x1MaxMKS = Lx/ 100* 32 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 31.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 32 x1MaxMKS = Lx/ 100* 33 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 32.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 33 x1MaxMKS = Lx/ 100* 34 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 33.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 34 x1MaxMKS = Lx/ 100* 35 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 34.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 35 x1MaxMKS = Lx/ 100* 36 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 35.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 36 x1MaxMKS = Lx/ 100* 37 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 36.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 37 x1MaxMKS = Lx/ 100* 38 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 37.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 38 x1MaxMKS = Lx/ 100* 39 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 38.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 39 x1MaxMKS = Lx/ 100* 40 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 39.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 40 x1MaxMKS = Lx/ 100* 41 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 40.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 41 x1MaxMKS = Lx/ 100* 42 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 41.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 42 x1MaxMKS = Lx/ 100* 43 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 42.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 43 x1MaxMKS = Lx/ 100* 44 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 43.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 44 x1MaxMKS = Lx/ 100* 45 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 44.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 45 x1MaxMKS = Lx/ 100* 46 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 45.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 46 x1MaxMKS = Lx/ 100* 47 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 46.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 47 x1MaxMKS = Lx/ 100* 48 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 47.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 48 x1MaxMKS = Lx/ 100* 49 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 48.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 49 x1MaxMKS = Lx/ 100* 50 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 49.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 50 x1MaxMKS = Lx/ 100* 51 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 50.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 51 x1MaxMKS = Lx/ 100* 52 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 51.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 52 x1MaxMKS = Lx/ 100* 53 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 52.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 53 x1MaxMKS = Lx/ 100* 54 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 53.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 54 x1MaxMKS = Lx/ 100* 55 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 54.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 55 x1MaxMKS = Lx/ 100* 56 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 55.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 56 x1MaxMKS = Lx/ 100* 57 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 56.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 57 x1MaxMKS = Lx/ 100* 58 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 57.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 58 x1MaxMKS = Lx/ 100* 59 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 58.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 59 x1MaxMKS = Lx/ 100* 60 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 59.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 60 x1MaxMKS = Lx/ 100* 61 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 60.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 61 x1MaxMKS = Lx/ 100* 62 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 61.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 62 x1MaxMKS = Lx/ 100* 63 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 62.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 63 x1MaxMKS = Lx/ 100* 64 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 63.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 64 x1MaxMKS = Lx/ 100* 65 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 64.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 65 x1MaxMKS = Lx/ 100* 66 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 65.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 66 x1MaxMKS = Lx/ 100* 67 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 66.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 67 x1MaxMKS = Lx/ 100* 68 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 67.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 68 x1MaxMKS = Lx/ 100* 69 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 68.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 69 x1MaxMKS = Lx/ 100* 70 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 69.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 70 x1MaxMKS = Lx/ 100* 71 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 70.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 71 x1MaxMKS = Lx/ 100* 72 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 71.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 72 x1MaxMKS = Lx/ 100* 73 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 72.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 73 x1MaxMKS = Lx/ 100* 74 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 73.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 74 x1MaxMKS = Lx/ 100* 75 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 74.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 75 x1MaxMKS = Lx/ 100* 76 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 75.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 76 x1MaxMKS = Lx/ 100* 77 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 76.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 77 x1MaxMKS = Lx/ 100* 78 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 77.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 78 x1MaxMKS = Lx/ 100* 79 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 78.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 79 x1MaxMKS = Lx/ 100* 80 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 79.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 80 x1MaxMKS = Lx/ 100* 81 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 80.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 81 x1MaxMKS = Lx/ 100* 82 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 81.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 82 x1MaxMKS = Lx/ 100* 83 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 82.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 83 x1MaxMKS = Lx/ 100* 84 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 83.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 84 x1MaxMKS = Lx/ 100* 85 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 84.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 85 x1MaxMKS = Lx/ 100* 86 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 85.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 86 x1MaxMKS = Lx/ 100* 87 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 86.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 87 x1MaxMKS = Lx/ 100* 88 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 87.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 88 x1MaxMKS = Lx/ 100* 89 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 88.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 89 x1MaxMKS = Lx/ 100* 90 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 89.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 90 x1MaxMKS = Lx/ 100* 91 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 90.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 91 x1MaxMKS = Lx/ 100* 92 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 91.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 92 x1MaxMKS = Lx/ 100* 93 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 92.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 93 x1MaxMKS = Lx/ 100* 94 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 93.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 94 x1MaxMKS = Lx/ 100* 95 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 94.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 95 x1MaxMKS = Lx/ 100* 96 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 95.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 96 x1MaxMKS = Lx/ 100* 97 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 96.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 97 x1MaxMKS = Lx/ 100* 98 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 97.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 98 x1MaxMKS = Lx/ 100* 99 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 98.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } Load { speciesName = electrons density = electronDensityMKS x1MinMKS = Lx/ 100* 99 x1MaxMKS = Lx/ 100* 100 x2MinMKS = 0.0 x2MaxMKS = Ly np2c = electronNumRatio units = MKS v1drift = electronVelocityMKS*sin(2*PI*wl/ 100* 99.5000) v1thermal = rmsElectronSpeedMKS v2thermal = 0. v3thermal = 0. } }