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

Use the tensor in Stix's (10-57) (true) or Multiply by kperp^2 d_ref^2 (false)

Write or suppress output to screen.

Enable optional startup profiling through the /system/ namelist. Defaults to .false. for compatibility with existing input files. startup_tick and startup_tock measure elapsed wall time using MPI_Wtime, not CPU time. Normal ALPS startup enables profiling only on rank 0; this flag is not broadcast. startup_report prints only on rank 0, independently of writeOut.

Relativistic cache mode: off (default), auto, or refresh; used only on rank 0. Set in the optional &relativistic_params group, read only for relativistic runs.

Existing directory for native content-addressed cache entries and temporary files. Shared by relativistic species; configured only through &relativistic_params.

Accumulated local startup wall times in seconds, initially zero. Repeated startup_tock calls add to the selected phase, so the normal rank-0 report sums work across species. No MPI reduction is performed.

Phase 10 includes MPI distribution and harmonic/Bessel setup, but excludes MPI initialization and namelist input. It overlaps phases 1--9 and must not be added to them. Phase 14 overlaps phases 3--8 and 11--13; phase 10 also includes cache work. startup_report does not reset these accumulators.

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 intervals in relativistic $\Gamma=\sqrt{1+\frac{p_\perp^2+p_\parallel^2}{m_j^2c^2}}$, $N_\Gamma$ (Eqn. 3.14), giving array bounds 0:ngamma. Shared &relativistic_params setting; default 100, minimum 2 for relativistic runs.

Number of intervals in dimensionless parallel momentum $\bar{p}_\parallel = p_\parallel/m_j c$, $N_{\bar{p}_\parallel}$, giving bounds 0:npparbar. Shared &relativistic_params setting; default 200, minimum 4 for relativistic runs.

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.

Tolerance for secant method=1, alps_fn

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)

Selection of derivative type for Orthogonal Basis Function, (1:nspec) 0: Finite Difference 1: Analytic Evaluation

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$.