input.txt

pic 模拟程序！面向对象

Equipotential--------------------------------
[see timefunc.txt, evaluator.txt and eqpot.txt for gotchas]

Conductor parameters.
*note: uses timefunction (Generic Boundary parameters) to set spatially
	uniform time-dependent voltage on surface.  Works as a grounded
	Conductor in the electromagnetic model.
*note: if you have multiple boundary segments which are all at the same potential
	it is important for performance that you specify them all
	in the same Equipotential group.

VanalyticF = "1" (scalar) String specifying the analytic function with which
             to scale all the other variables affecting the voltage
             (the time function, C, A, a1, a0, tdelay, trise, tfall, tpulse).
             The values "x1" and "x2", when used in VanalyticF, mean the
             physical x1 and x2 coordinates.
	

Foil-----------------------------------------[currently nonfunctional]

Generic Boundary parameters, multiple segments permitted.
nemit=0  (scalar) Number of secondary electrons emitted per incident
			electron.


PortGauss --------------------------------

PortGauss launches two linearly polarized electromagnetic
pulses in 2-D Cartesian geometry. The pulses are launched 
from a given boundary by controlling the temporal and spatial
dependence of the electric field at that boundary. The
transverse spatial profile is Gaussian along the boundary.
The temporal variation can be either trapezoidal, Gaussian, 
or a half-sine.  Pulse 0 is linearly polarized in the y 
direction, while pulse 1 is linearly polarized in the z 
direction.  The pulses are launched in accordance with paraxial
theory, such that the waist (focus) is a given distance from
the port.

The parameters of pulse 0 (linearly polarized in the 
y direction) are:

 - waveLeng_p0            = wavelength [meters] of the EM wave
 - amp_p0                 = peak wave amplitude of electric field [Volts/meter] 
 - spotSize_p0            = The beam half width [meters] at the waist.
 - focus_p0               = Displacement [meters] of the focus.
 - pulShp_p0              = pulse shape (0=trapezoid; 1=Gaussian; 2=half-sine;
                                         3=wide gaussian) 
 - pulLeng_p0             = pulse length [meters]. (half-width for Gaussian, 
				full for half-sine and trapezoid). 
 - tdelay_p0              = time delay [seconds] to start pulse. 

 - offset		= an offset in MKS units from the center of the
 			  simulation to the peak of the Gaussian.  Default
			  is zero.

The parameters of pulse 1 (linearly polarized in the z 
direction) are the same, except that "_p0" is replaced
by "_p1".  

To have only one pulse, set the amplitude of the other to
zero.

The code crashes if the value of focus is set to zero. If the focus needs
to be at the boundary, set the value of the focus equal to one or two
cell lenght inside the region (Dx or 2*Dx). 
 
For option 1, gaussian shape, at t = tdelay +- tpulse/2, the value is
exp(-1), and zero beyond that.  For option 3, the pulse is actually cut
off further out (meaning your pulse is 6x longer), so the minimum value
is 4.8e-6 w.r.t the maximum.




Diagnostics--------------------------------
    For more information read diagnostic.txt
j1 (int) x1 index for first diagnostic endpoint.
k1 (int) x2 index for first diagnostic endpoint.
j2=j1 (int) x1 index for second diagnostic endpoint.
k2=k1 (int) x2 index for second diagnostic endpoint.

Alternately, you may specify boundary locations in MKS units.
However, XOOPIC will put the diagnostic on the
nearest grid point.

A1 (scalar) x1 location for first boundary endpoint. 
A2 (scalar) x2 location for first boundary endpoint.
B1 (scalar) x1 location for second boundary endpoint.
B2 (scalar) x2 location for second boundary endpoint.

VarName (string) name of variable to be plotted.  Currently the follow variables can be plotted.
For time histories or for spatial regions:
    E1 -> Ez (RZ) or Ex (XY)
    E2 -> Er (RZ) or Ey (XY)
    E3 -> Ephi (RZ) or Ez (XY)
    B1 -> Bz (RZ) or Bx (XY)
    B2 -> Br (RZ) or By (XY)
    B3 -> Bphi (RZ) or Bz (XY)
    I1 -> Iz (RZ) or Ix (XY) (only with EM field solve)
    I2 -> Ir (RZ) or Iy (XY) (only with EM field solve)
    I3 -> Iphi (RZ) or Iz (XY) (only with EM field solve)
    intEdl1 -> Ez (RZ) or Ex (XY)
    intEdl2 -> Er (RZ) or Ey (XY)
    intEdl3 -> Ephi (RZ) or Ez (XY)
	 poynting1 -> Poynting Vector in x1 (only with EM field solve)
	 poynting2 -> Poynting Vector in x2 (only with EM field solve)
	 poynting3 -> Poynting Vector in x3 (only with EM field solve)
    rho -> charge density
	 speciesName -> rho of a given species
    phi -> potential (only with electrostatic field solve)
	 Q -> surface charge on dielectrics
    LaserSpotSize  ->  Integral of y*y* Ey*Ey / Integral of Ey * Ey over 
		     the line k1 to k2, assuming the laser spot is centered
		     at (k2-k1)/2
		     Displayed is the average over HISTMAX timesteps of this
		     measure.
    WaveDirDiagnostic -> Computes over the mesh (Ey - c*Bz)/2 and (Ey + c*Bz)/2,
                         which distinguishes left and right moving waves in
                         the system (polarized in y).

                         "polarizationEB = EzBy" in the "Diagnostic" group
                         specifies the calculation of (Ez-c*By)/2 and
                         (Ez + c*By)/2 instead. If "polarizationEB" is not
                         given or "polarizationEB = EyBz", the computed 
                         diagnostics is for (Ey - c*Bz)/2 and (Ey + c*Bz)/2.
                         
                         "psd1dFlag = 1" calculates the 1d power spectral 
                         densities for the two linear combinations of E and B
                         selected via the "polarizationEB". The 1d PSD are
                         calculated along the x axis for each value of the
                         y index. By default, the psd1dFlag is turned off, 
                         i.e. its default value is zero.

                         "windowName = Hann" specifies windowing of the data
                         before the FFT is done. If "windowName" is not 
                         given in the "Diagnostic" group, the "Blackman" window 
                         is applied to the data. The following windows are 
                         implemented: "Blackman", "Bartlett", "Hamming", 
                         "Hann", and "Welch". 

                         Here is an example of this diagnostics structure:

                         Diagnostic
                         {
                         	j1 = 0
                                j2 = Nx                 
                                k1 = 0                  
                                k2 = Ny
                                VarName = WaveDirDiagnostic
                                polarizationEB = EzBy
                                psd1dFlag = 1 // calculate the 1d power spectral density
                                windowName = Hann
                                title = Wave Energy
                                x1_Label = x
                                x2_Label = y
                                x3_Label = Wave Energy
                         }                         

                         An example in a working input file is: "inp/TI_H_WDD.inp".

   PSDFieldDiag1d -> This diagnostic computes the power spectral density 
                     of em field components. 
                     The diagnostic calculates 1d FFT transforms along
                     the first spatial coordinate of the computational grid and
                     for each fixed value of the second spatial coordinate of the
                     grid. This choice is based on the assumption that a 
                     laser pulse always propagates along the first 
                     spatial direction. 

                     The resulting diagnostic displays a surface in 3D. The first
                     coordinate of the diagnostic is the second coordinate of the 
                     computational grid. The second coordinate of the diagnostic is
                     the wave number of the Fourier transform. 

                     Example of how to enable this diagnostic in the input file is
                     the following structure:

                     Diagnostic
                     {
                     	j1 = 0
                        j2 = Nx
                        k1 = 0
                        k2 = Ny
                        VarName = PSDFieldDiag1d
                        windowName = Blackman
                        title = 1D Power Spectral Density
                        x1_Label = y
                        x2_Label = kx
                        x3_Label = 1d PSD
                        fieldName = E
                        fieldComponentLabel = 3
                     }   

                     Users can again specify a window function 
                     for the data to be Fourier transformed. 
                     The example above uses the Blackman window.
                     The labels of the diagnostic plot axes as 
                     well as the the plot title
                     can also be specified via the input file. 

                     In the example above, the power spectral
                     density of the third component of the  
                     electric field is calculated, "Ez".

                     The calculation of this diagnostic is done in 
                     the PSDFieldDiag1d class. This class also 
                     contains one more diagnostic that is not
                     currently available via the input file but
                     can be enabled by setting the flag
                     "flagChirp" to true in the "psdFieldDiag1d.cpp"
                     file and recompiling with "HAVE_FFT" defined.

                     This diagnostic calculates the frequency 
                     spectrogram of field component through the
                     middle of the pulse, defined by index of the
                     second coordinate of the computational box
                     set to its half value. The diagnostic
                     selects 64 points at a time along the first
                     coordinate of the computational box and 
                     calculates the PSD of the specified field.
                     For the value in the middle of this 64 
                     points interval, the PSD is plotted along
                     the frequencies of the FFT. A translation of
                     16 points long the pulse is made and the
                     procedure is repeated. It gives the frequency
                     content of the pulse along its lengths.

                     This is an experimental diagnostic and it is
                     disabled by default. Once it is tested 
                     extensively, it will be available through
                     the input file. The intervals for the FFT
                     and the skip will be parameters in the input file
                     and not restricted to 64 and 16 since they
                     depend on the size of the box. 

                     This diagnostic is particularly good for detecting
                     chirped signals. 
                     
                     An example in a working input file is: "inp/TI_H_WDD.inp".

   PSDFieldDiag2d -> This diagnostic computes the power spectral density 
                     of em field components. 

                     The diagnostic calculates 2d FFT transforms along
                     both spatial coordinate of the computational grid. 
                     
                     Here is an example of what is needed in an input file to
                     enable it, assuming the code is compiled with the FFTW
                     library enabled and linked.

                     Diagnostic
                     {
                     	j1 = 0
                        j2 = Nx
                        k1 = 0
                        k2 = Ny
                        VarName = PSDFieldDiag2d
                        windowName = None
                        title = 2D Power Spectral Density
                        x1_Label = kx
                        x2_Label = ky
                        x3_Label = 2d PSD
                        fieldName = E
                        fieldComponentLabel = 3
                     }
                     
                     An example in a working input file is: "inp/TI_H_WDD.inp".
                         
For spatial regions
    JdotE -> JdotE 
x1_Label=x1 (string) x1 Label of the XGrafix plot
x2_Label=x2 (string) x2 Label of the XGrafix plot
x3_Label=x3 (string) x3 Label of the XGrafix plot
title=not_named (string) Title of XGrafix window. 
**************NOTE**************
    DON'T HAVE TWO TITLES THE SAME.  XGrafix DOESN'T LIKE IT.
nfft (int) number of data points for fft, must be power of 2[currently nonfunctional]
HistMax=64 (int) maximum length of history array
save=0 (int)  Flag: 1 saves the diagnostic data in the dumpfile
				Flag: 0 restarts when restarted from a dumpfile
Comb=0 (int) Every Comb'th value is left when HistMax is reached
    **Note if Comb=0 the history is a local history.**
Ave=0 (int) Averaged over Ave data points when adding to history array
integral=NULL (string) one of: line (variable dotted into dl), flux (variable dotted 
              into dS), sum (simple summation).