Processing math: 100%

system


Variables

Name Type Default Description
kperp double precision None

Perpendicular wavenumber, normalized by inverse reference inertial length, kdp.

kpar double precision None

Parallel wavenumber, normalized by inverse reference inertial length, kdp.

nspec integer None

Number of plasma components.

nroots integer None

Number of dispersion solutions under consideration.

use_map logical None

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

writeOut logical .true.

Write or suppress output to screen.

nperp integer None

Number of perpendicular momentum space grid points, N.

npar integer None

Number of parallel momentum space grid points, N.

ngamma integer 100

Number of grid points in relativitic Γ=1+p2+p2m2jc2, NΓ (Eqn. 3.14).

npparbar integer 200

Number of grid points in dimensionless paralell momentum ˉp=p/mjc, Nˉp.

vA double precision None

Alfven Velocity, normalized to speed of light, vAp/c.

arrayName character(len=75) None

Name of input files for distributions.

Bessel_zero double precision 1.d-45

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

secant_method integer 2

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

numiter integer 50

Maximum number of iterations in secant method.

D_threshold double precision 1.d-5

Minimum threshold for secant method.

D_prec double precision 1.d-5

Size of bounding region for secant method.

D_gap double precision 1.d-5

Size of allowable difference between roots.

positions_principal integer 5

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

Tlim double precision 0.01d0

Threshold for analytical principal-value integration for evaluating Eqn 3.6 and 3.7, tlim.

maxsteps_fit integer 500

Maximum number of fitting iterations.

lambda_initial_fit double precision 1.d0

Inital Levenberg-Marquardt damping parameter.

lambdafac_fit double precision 1.d1

Adjustment factor for Levenberg-Marquardt damping parameter.

epsilon_fit double precision 1.d-8

Convergence for Levenberg-Marquardt fit.

fit_check logical .true.

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

determine_minima logical .true.

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

n_resonance_interval integer 100

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

scan_option integer 1

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

n_scan integer 0

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.