alps_var Module

Root of input file name.

Directory of input file name.

Number of dispersion solutions under consideration.

Number of dispersion solutions found in frequency map scan.

Choice of: (T) searching for roots over a map in complex frequency space, via map_read; (F) input (nroots) guesses for solutions, via solution_read

Write or suppress output to screen.

Output unit for error file.

Total number of processors.

Number of local processor.

T if iproc=0.

Integer error flag.

Perpendicular wavenumber, normalized by inverse reference inertial length, $k_\perp d_p$.

Parallel wavenumber, normalized by inverse reference inertial length, $k_\parallel d_p$.

Alfven Velocity, normalized to speed of light, $v_{Ap}/c$.

Calculate Bessel functions until the maximum is less than this value.

Number of plasma components.

Number of relativistic plasma components.

Dispersion solutions, ranging from 1 to numroots.

Maximum number of solutions.

Smallest $\omega_{\textrm{r}}/\Omega_p$ value for complex map search.

Largest $\omega_{\textrm{r}}/\Omega_p$ value for complex map search.

Smallest $\gamma/\Omega_p$ value for complex map search.

Largest $\gamma/\Omega_p$ value for complex map search.

Linear (F) or Log (T) spacing for $\omega_{\textrm{r}}/\Omega_p$ map search.

Linear (F) or Log (T) spacing for $\gamma/\Omega_p$ map search.

Number of $\omega_{\textrm{r}}/\Omega_p$ points in frequency grid.

Number of $\gamma/\Omega_p$ points in frequency grid.

Number of perpendicular momentum space grid points, $N_\perp$.

Number of parallel momentum space grid points, $N_\parallel$.

Number of grid points in relativitic $\Gamma=\sqrt{1+\frac{p_\perp^2+p_\parallel^2}{m_j^2c^2}}$, $N_\Gamma$ (Eqn. 3.14).

Number of grid points in dimensionless paralell momentum $\bar{p}_\parallel = p_\parallel/m_j c$, $N_{\bar{p}_\parallel}$.

Number of parallel momentum steps distant from the resonant momentum included in the numerical calculation of Eqn 3.5, $M_I$.

Threshold for analytical principal-value integration for evaluating Eqn 3.6 and 3.7, $t_{\textrm{lim}}$.

The species number on which this process is working.

Selection of root finding method. 0: secant 1: rtsec 2: secant_osc: an improvement to reduce oscillatory loops.

Maximum number of iterations in secant method.

Minimum threshold for secant method.

Size of bounding region for secant method.

Size of allowable difference between roots.

The ratio of a circle's circumference to its diameter.

Name of input files for distributions.

Background distribution function array $f_{0j}$; (1:nspec,0:nperp,0:npar).

Relativistic background distribution function array $f_{0j}$; (1:nspec,0:ngamma,0:npparbar).

Perpendicular and parallel derivatives of $f_{0j}$; (1:nspec,0:nperp,0:npar,1:2), with $\partial_{p_\perp} f_{0j}$ in index 1, and $\partial_{p_\parallel} f_{0j}$ in index 2.

Derivatives of $f_{0j}$; (1:nspec,0:nperp,0:npar,1:2), with $\partial_{\Gamma} f_{0j}$ in index 1 and $\partial_{\bar{p}_\parallel} f_{0j}$ in index 2.

Momentum Space Array for $f_{0j}$; (1:nspec,0:nperp,0:npar,1:2) with $p_\perp/m_p v_A$ in index 1 and $p_\parallel/m_p v_A$ in index 2.

Current density, (0:nspec) Zeroth index is sum over all species

density, (0:nspec) Zeroth index is sum over all species

Relativistic momentum space array of $\Gamma$; (1:nspec,0:ngamma,0:npparbar).

Relativistic momentum space array of $\bar{p}_\parallel$; (1:nspec,0:ngamma,0:npparbar).

number of n values to sum over, (1:nspec).

Lower and Upper limits for n values for iproc to sum over.

Ratio of species density to reference $n_j/n_p$, (1:nspec).

Ratio of species charge to reference $q_j/q_p$, (1:nspec).

Ratio of species mass to reference $m_j/m_p$, (1:nspec).

Use relativistic treatment; (1:nspec).

Wave Equation Tensor (1:3,1:3).

Susceptibility Tensor (1:nspec,1:3,1:3).

Susceptibility Tensor (1:nspec,1:3,1:3,0:1).

Array of Bessel functions; (nlim(1):nlim(2)+1,0:nperp).

Selection of method for Hybrid-Analytical Continutation, (1:nspec) 0) Use the analytic function. 1) Use the 'n_fits' functions described with 'fit_type' 2) Use a polynomial basis representation.

Polynomials of order (0,poly_order) evaluatated at (0,npar) abscissa points, (1:nspec,0:npar,0:poly_order) poly_order is taken to be the maximum order across all species.

Fit Coefficients for Polynomials from Generarl Linear Least Squares Method (1:nspec,0:nperp,0:poly_order) poly_order is taken to be the maximum order across all species. The fit is taken in one dimension at each perpendicular momentum value.

Selection of Orthogonal Basis Function, (1:nspec) 1) Chebyshev Polynomials Other Polynomials can be added in future releases.

Limit on the Evaluation of the Polynomial Representation, (1:nspec)

Selection of Maximum Order of Orthogonal Basis Function, (1:nspec)

Number of fitted functions, (1:nspec)

Type of analytic function to be fit, (1:nspec,1:maxval(nfits)); 1) Maxwellian, 2) Kappa, 3) Juettner with $p_\perp,p_\parallel$, 4) Juettner with $\Gamma,\bar{p}_\parallel$, constant $\bar{p}_\parallel$, 5) Juettner with $p_\perp,p_\parallel$; variable $\bar{p}_\parallel$, 6) Bi-Moyal distribution.

Maximum number of fitting iterations.

Inital Levenberg-Marquardt damping parameter.

Adjustment factor for Levenberg-Marquardt damping parameter.

Convergence for Levenberg-Marquardt fit.

Fit parameters, (1:nspec,0:nperp,4,maxval(n_fits)).

This parameter, $y$ in Eqn. B1, compensates for the strong $p_\perp$ dependence of $u_1$, making the fit more reliable.

If true, output fitted functions to ASCII file for each species.

If true, after map search, determine minima and refine solutions.

How many steps should be used to integrate around the resonance, $M_P$, used for integrating near poles (see section 3.1).

Select case for scans; 1) consecutive scans along input paths in wavevector space, 2) double scans of two selected parameters.

Number of wavevector scans. Must be set to 2 for scan_option=2; Must be 1 or larger for scan_option=1. 0 turns off wavevector scans.

Use logarithmic fitting, (1:nspec).

Use bi-Maxwellian/cold calculation from NHDS, (1:nspec).

Maximum number of n for NHDS bi-Maxwellian calculation, (1:nspec).

Bessel-zero for NHDS bi-Maxwellian calculation (1:nspec).

Species beta $\beta_{\parallel,j}$ for NHDS bi-Maxwellian calculation (1:nspec). If bMbetas=0.d0, then this species is treated with the cold-plasma susceptibility.

Species temperature anisotropy $T_{\perp,j}/T_{\parallel,j}$ for NHDS bi-Maxwellian calculation.

Species drift momentum for NHDS bi-Maxwellian/cold calculation, in units of $m_p v_A$ (1:nspec).

Scan parameters for each wavevector scan. Read in from scan_read.

Previous value of $k_\perp$.

Previous value of $k_\parallel$.

Current value of $k_\perp$.

Current value of $k_\parallel$.