VERB_code_2.3
PSD Class Reference

Phase Space Density class. More...

#include <PSD.h>

Public Member Functions

void Create_Initial_PSD (Parameters_structure::PSD parameters, Grid &grid, BoundaryCondition L_UpperBoundaryCondition)
 Create inital PSD (Steady State) More...
 
void Interpolate (PSD &otherPSD, Parameters_structure::Interpolation interpolationParameters, Grid &oldGrid, Grid &newGrid, Matrix2D< double > newGrid_pc_lowerBoundaryCondition, Matrix2D< double > newGrid_pc_upperBoundaryCondition)
 Interpolation function. More...
 
void Load_initial_f (GridElement &L, GridElement &pc, GridElement &alpha, double tau, double Kp, double min_f=1.e-99, string J_L7_function="J_L7", double fb_out=1, double fb_in=0, bool using_ip_90=true)
 
void Load_initial_f_two_zone (GridElement &L, GridElement &pc, GridElement &alpha, double tau, double tauLpp, double Kp, double min_psd=1.e-99, string J_L7_function="J_L7", double fb_out=1, double fb_in=0)
 two zone initial profile (kckim) More...
 
void Load_initial_f_from_outer_L (GridElement &L, GridElement &pc, GridElement &alpha, double tau, double Kp, Matrix2D< double > L_UpperBoundaryCondition, double min_psd=1.e-99, double fb_out=1, double fb_in=0)
 Calculate initial PSD from steady state using boundary conditions. Load initial PSD from outer L. More...
 
void Load_initial_f_2d (GridElement &L, GridElement &pc, GridElement &alpha, const char *filename)
 
void Load_initial_f_file (GridElement &L, GridElement &pc, GridElement &alpha, const char *filename, bool with_grid)
 
void Load_initial_f_maxwell (GridElement &L, GridElement &pc, GridElement &alpha)
 
void DiffusionMixTermExplicit (double dt, double Lpp, DiffusionCoefficient &Dpca, DiffusionCoefficient &DpcaLpp, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > pc_lowerBoundaryCondition, Matrix2D< double > pc_upperBoundaryCondition, Matrix2D< double > alpha_lowerBoundaryCondition, Matrix2D< double > alpha_upperBoundaryCondition, string pc_lowerBoundaryCondition_calculationType, string pc_upperBoundaryCondition_calculationType, string alpha_lowerBoundaryCondition_calculationType, string alpha_upperBoundaryCondition_calculationType)
 
void Diffusion_alpha (double dt, double Lpp, DiffusionCoefficient &Daa, DiffusionCoefficient &DaaLpp, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > alpha_lowerBoundaryCondition, Matrix2D< double > alpha_upperBoundaryCondition, string alpha_lowerBoundaryCondition_calculationType, string alpha_upperBoundaryCondition_calculationType)
 
void Diffusion_pc (double dt, double Lpp, DiffusionCoefficient &Dpcpc, DiffusionCoefficient &DpcpcLpp, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > pc_lowerBoundaryCondition, Matrix2D< double > pc_upperBoundaryCondition, string pc_lowerBoundaryCondition_calculationType, string pc_upperBoundaryCondition_calculationType)
 
void Diffusion_L (double dt, double Lpp, DiffusionCoefficient &DLL, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > lowerBoundaryCondition, Matrix2D< double > upperBoundaryCondition, string lowerBoundaryCondition_calculationType, string upperBoundaryCondition_calculationType, double tau, double tauLpp)
 
void Diffusion_pc_alpha (double dt, double Lpp, DiffusionCoefficient &Dpcpc, DiffusionCoefficient &DpcpcLpp, DiffusionCoefficient &Daa, DiffusionCoefficient &DaaLpp, DiffusionCoefficient &Dpca, DiffusionCoefficient &DpcaLpp, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > pc_lowerBoundaryCondition, Matrix2D< double > pc_upperBoundaryCondition, Matrix2D< double > alpha_lowerBoundaryCondition, Matrix2D< double > alpha_upperBoundaryCondition, string pc_lowerBoundaryCondition_calculationType, string pc_upperBoundaryCondition_calculationType, string alpha_lowerBoundaryCondition_calculationType, string alpha_upperBoundaryCondition_calculationType)
 
void Diffusion_pc_alpha_KC (double dt, double Lpp, DiffusionCoefficient &Dpcpc, DiffusionCoefficient &DpcpcLpp, DiffusionCoefficient &Daa, DiffusionCoefficient &DaaLpp, DiffusionCoefficient &Dpca, DiffusionCoefficient &DpcaLpp, GridElement &L, GridElement &pc, GridElement &alpha, Matrix3D< double > Jacobian, Matrix2D< double > pc_lowerBoundaryCondition, Matrix2D< double > pc_upperBoundaryCondition, Matrix2D< double > alpha_lowerBoundaryCondition, Matrix2D< double > alpha_upperBoundaryCondition, string pc_lowerBoundaryCondition_calculationType, string pc_upperBoundaryCondition_calculationType, string alpha_lowerBoundaryCondition_calculationType, string alpha_upperBoundaryCondition_calculationType)
 
void SourcesAndLosses (Parameters_structure parameters, GridElement &L, GridElement &pc, GridElement &alpha, AdditionalSourcesAndLosses &SL, double dt, double Lpp, double tau, double tauLpp, double Kp)
 
void Output_without_grid (double time)
 

Public Attributes

Matrix3D< double > arr
 array of PSD values
 
Parameters_structure::PSD PSD_parameters
 Parameters from ini file.
 
ofstream output_without_grid_file
 output stream for PSD
 
CalculationMatrix matr_A
 calculation matrices A, B, and C
 
CalculationMatrix matr_B
 calculation matrices A, B, and C
 
CalculationMatrix matr_C
 calculation matrices A, B, and C
 

Detailed Description

Phase Space Density class.

Do diffusions and store the result in Matrix3D. Do all possible operations with PSD (like output, loading, etc).

Definition at line 48 of file PSD.h.

Member Function Documentation

void PSD::Create_Initial_PSD ( Parameters_structure::PSD  parameters,
Grid grid,
BoundaryCondition  L_UpperBoundaryCondition 
)

Create inital PSD (Steady State)

Parameters
parameters- interpolation parameters structure
&grid- grid
L_UpperBoundaryCondition- upper L boundary values, in case we need to create initial PSD as a steady state radial diffusion solution

Initializing: storing parameters, loading initial values, making boundary conditions, initializing output parameters Simply, it is a creation of the object.

< il - L; im - pc; ia - alpha

moved to parameters.cpp if (parameters.initial_PSD_tauSteadyState <= 1.e-99) parameters.initial_PSD_tauSteadyState = 4.0/parameters.initial_PSD_Kp0;

moved to parameters.cpp if (parameters.initial_PSD_tauSteadyState <= 1.e-99) parameters.initial_PSD_tauSteadyState = 4.0/parameters.initial_PSD_Kp0;

moved to parameters.cpp if (parameters.initial_PSD_tauSteadyState <= 1.e-99) parameters.initial_PSD_tauSteadyState = 4.0/parameters.initial_PSD_Kp0;

Definition at line 47 of file PSD.cpp.

void PSD::Diffusion_alpha ( double  dt,
double  Lpp,
DiffusionCoefficient Daa,
DiffusionCoefficient DaaLpp,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  alpha_lowerBoundaryCondition,
Matrix2D< double >  alpha_upperBoundaryCondition,
string  alpha_lowerBoundaryCondition_calculationType,
string  alpha_upperBoundaryCondition_calculationType 
)
Todo:
lot's of repeated code within diffusion functions. maybe find a way to condense them

Pitch angle diffusion calculation function. Takes diffusion coefficients, boundary conditions and grid elements as input. It calculates model matrix by MakeModelMatrix_3D for 3D problem: A_3D * PSD_3D(t+1) = B_3D * PSD_3D(t) + C_3D then splits it for 1D problems: A * PSD_1D(t+1) = B * PSD_1D(t) + C and solve by tridiag method.

Parameters
dt- time step
Lpp- plasma pause location
&Daa- L diffusion coefficient
&DaaLpp- L diffusion coefficient
&L- grid element L
&pc- grid elelent pc
&alpha- grid element alpha
&Jacobian- jacobian
&alpha_lowerBoundaryCondition- lower boundary condition
&alpha_upperBoundaryCondition- upper boundary condition
&alpha_lowerBoundaryCondition_calculationType- lower boundary condition type (on value/on derivative)
&alpha_upperBoundaryCondition_calculationType- upper boundary condition type (on value/on derivative)

Definition at line 599 of file PSD.cpp.

void PSD::Diffusion_L ( double  dt,
double  Lpp,
DiffusionCoefficient DLL,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  L_lowerBoundaryCondition,
Matrix2D< double >  L_upperBoundaryCondition,
string  L_lowerBoundaryCondition_calculationType,
string  L_upperBoundaryCondition_calculationType,
double  tau,
double  tauLpp 
)

Radial diffusion calculation function. Takes diffusion coefficients, boundary conditions and grid elements as parameters. It caltulates model matrix by makeMatrix and then solve it by SolveMatrix method.

Parameters
dt- time step
Lpp- plasma pause location
&DLL- L diffusion coefficient
&L- grid element L
&pc- grid elelent pc
&alpha- grid element alpha
&Jacobian- jacobian
&L_lowerBoundaryCondition- lower boundary condition on L
&L_upperBoundaryCondition- upper boundary condition on L
&L_lowerBoundaryCondition_calculationType- lower boundary condition on L calculation type
&L_upperBoundaryCondition_calculationType- upper boundary condition on L calculation type
tau- life time upper location of the plasma pause
tauLpp- life time lower location of the plasma pause

Definition at line 895 of file PSD.cpp.

void PSD::Diffusion_pc ( double  dt,
double  Lpp,
DiffusionCoefficient Dpcpc,
DiffusionCoefficient DpcpcLpp,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  pc_lowerBoundaryCondition,
Matrix2D< double >  pc_upperBoundaryCondition,
string  pc_lowerBoundaryCondition_calculationType,
string  pc_upperBoundaryCondition_calculationType 
)

Energy diffusion calculation function. Takes diffusion coefficients, boundary conditions and grid elements as input. It calculates model matrix by MakeModelMatrix_3D for 3D problem: A_3D * PSD_3D(t+1) = B_3D * PSD_3D(t) + C_3D then splits it for 1D problems: A * PSD_1D(t+1) = B * PSD_1D(t) + C and solve by tridiag method.

Parameters
dt- time step
Lpp- plasma pause location
&Dpcpc- L diffusion coefficient
&DpcpcLpp- L diffusion coefficient with plasma pause location
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha
&Jacobian- jacobian
&pc_lowerBoundaryCondition- lower boundary condition
&pc_upperBoundaryCondition- upper boundary condition
&pc_lowerBoundaryCondition_calculationType- lower boundary condition type (on value/on derivative)
&pc_upperBoundaryCondition_calculationType- upper boundary condition type (on value/on derivative)

Definition at line 747 of file PSD.cpp.

void PSD::Diffusion_pc_alpha ( double  dt,
double  Lpp,
DiffusionCoefficient Dpcpc,
DiffusionCoefficient DpcpcLpp,
DiffusionCoefficient Daa,
DiffusionCoefficient DaaLpp,
DiffusionCoefficient Dpca,
DiffusionCoefficient DpcaLpp,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  pc_lowerBoundaryCondition,
Matrix2D< double >  pc_upperBoundaryCondition,
Matrix2D< double >  alpha_lowerBoundaryCondition,
Matrix2D< double >  alpha_upperBoundaryCondition,
string  pc_lowerBoundaryCondition_calculationType,
string  pc_upperBoundaryCondition_calculationType,
string  alpha_lowerBoundaryCondition_calculationType,
string  alpha_upperBoundaryCondition_calculationType 
)

Radial diffusion calculation function. Takes diffusion coefficients, boundary conditions and grid elements as parameters. It caltulates model matrix by makeMatrix and then solve it by SolveMatrix method.

Parameters
dt- time step
Lpp- plasma pause location
&Dpcpc- diffusion coefficient
&DpcpcLpp- diffusion coefficient
&Daa- diffusion coefficient
&DaaLpp- diffusion coefficient
&Dpca- diffusion coefficient
&DpcaLpp- diffusion coefficient
&L- grid element L
&pc- grid elelent pc
&alpha- grid element alpha
&Jacobian- jacobian
&pc_lowerBoundaryCondition- lower boundary condition
&pc_upperBoundaryCondition- upper boundary condition
&alpha_lowerBoundaryCondition- lower boundary condition
&alpha_upperBoundaryCondition- upper boundary condition
&pc_lowerBoundaryCondition_calculationType- lower boundary condition calculation type
&pc_upperBoundaryCondition_calculationType- upper boundary condition calculation type
&alpha_lowerBoundaryCondition_calculationType- lower boundary condition calculation type
&alpha_upperBoundaryCondition_calculationType- upper boundary condition calculation type

Definition at line 1053 of file PSD.cpp.

void PSD::Diffusion_pc_alpha_KC ( double  dt,
double  Lpp,
DiffusionCoefficient Dpcpc,
DiffusionCoefficient DpcpcLpp,
DiffusionCoefficient Daa,
DiffusionCoefficient DaaLpp,
DiffusionCoefficient Dpca,
DiffusionCoefficient DpcaLpp,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  pc_lowerBoundaryCondition,
Matrix2D< double >  pc_upperBoundaryCondition,
Matrix2D< double >  alpha_lowerBoundaryCondition,
Matrix2D< double >  alpha_upperBoundaryCondition,
string  pc_lowerBoundaryCondition_calculationType,
string  pc_upperBoundaryCondition_calculationType,
string  alpha_lowerBoundaryCondition_calculationType,
string  alpha_upperBoundaryCondition_calculationType 
)

Radial diffusion calculation function with Kyungchan's method. Takes diffusion coefficients, boundary conditions and grid elements as parameters. It caltulates model matrix by makeMatrix and then solve it by SolveMatrix method.

Todo:
Literally the only difference between KC version and original version is calling one custom method MakeModelMatrix_3D_KC. Find a way to compress two functions together.
Parameters
dt- time step
Lpp- plasma pause location
&Dpcpc- diffusion coefficient
&DpcpcLpp- diffusion coefficient
&Daa- diffusion coefficient
&DaaLpp- diffusion coefficient
&Dpca- diffusion coefficient
&DpcaLpp- diffusion coefficient
&L- grid element L
&pc- grid elelent pc
&alpha- grid element alpha
&Jacobian- jacobian
&pc_lowerBoundaryCondition- lower boundary condition
&pc_upperBoundaryCondition- upper boundary condition
&alpha_lowerBoundaryCondition- lower boundary condition
&alpha_upperBoundaryCondition- upper boundary condition
&pc_lowerBoundaryCondition_calculationType- lower boundary condition calculation type
&pc_upperBoundaryCondition_calculationType- upper boundary condition calculation type
&alpha_lowerBoundaryCondition_calculationType- lower boundary condition calculation type
&alpha_upperBoundaryCondition_calculationType- upper boundary condition calculation type

Definition at line 1269 of file PSD.cpp.

void PSD::DiffusionMixTermExplicit ( double  dt,
double  Lpp,
DiffusionCoefficient Dpca,
DiffusionCoefficient DpcaLpp,
GridElement L,
GridElement pc,
GridElement alpha,
Matrix3D< double >  Jacobian,
Matrix2D< double >  pc_lowerBoundaryCondition,
Matrix2D< double >  pc_upperBoundaryCondition,
Matrix2D< double >  alpha_lowerBoundaryCondition,
Matrix2D< double >  alpha_upperBoundaryCondition,
string  pc_lowerBoundaryCondition_calculationType,
string  pc_upperBoundaryCondition_calculationType,
string  alpha_lowerBoundaryCondition_calculationType,
string  alpha_upperBoundaryCondition_calculationType 
)

Mixed terms calculation of energy and pitch angle by explicit method. Uses split approximation for calculation.

Parameters
dt- time step
Lpp- plasma pause location
&Dpca- pc-alpha diffusion coefficient
&DpcaLpp- pc-alpha diffusion coefficient under plasma pause location
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha
&Jacobian- jacobian
&pc_lowerBoundaryCondition- lower boundary condition on pc
&pc_upperBoundaryCondition- upper boundary condition on pc
&alpha_lowerBoundaryCondition- lower boundary condition on alpha
&alpha_upperBoundaryCondition- upper boundary condition on alpha
pc_lowerBoundaryCondition_calculationType- pc lower boundary condition type,
pc_upperBoundaryCondition_calculationType- pc upper boundary condition type,
alpha_lowerBoundaryCondition_calculationType- alpha lower boundary condition type,
alpha_upperBoundaryCondition_calculationType- alpha upper boundary condition type.

AM_Split_LR

\( LHS = G(GridElement) * Dpca * \frac {PSD[x]-PSD[x-1]}{GridElement[x]-GridElement[x-1]} \)

\( RHS = G(GridElement) * Dpca * \frac {PSD[x+1]-PSD[x]}{GridElement[x+1]-GridElement[x]} \) for VF::G() and x is iteration index and GridElement is pc or alpha

\( Lmar = \frac{1.0}{{pc}^2} * \frac{RHS[x] - RHS[x-1]}{pc[x] - pc[x-1]} \) \( Lmal = \frac{1.0}{{pc}^2} * \frac{LHS[x+1] - LHS[x]}{pc[x+1] - pc[x]} \)

\( Lamr = \frac{1.0}{G(\alpha)} * \frac{RHS[x]-RHS[x-1]}{\alpha[x] - \alpha[x-1]} \) \( Laml = \frac{1.0}{G(\alpha)} * \frac{LHS[x+1]-LHS[x]}{\alpha[x+1] - \alpha[x]} \)

\( PSD += \frac{dt}{2} * (Lmar + Lamr + Lmal + Laml) \)

AM_Split_C

\( PSD\_m\_r = Dpca[m+1] * Jacobian[m+1] * \frac{PSD[m+1][a+1] - PSD[m+1][a-1]}{\alpha[m+1][a+1] - \alpha[m+1][a-1]} \) \( PSD\_m\_l = Dpca[m-1] * Jacobian[m-1] * \frac{PSD[m-1][a+1] - PSD[m-1][a-1]}{\alpha[m-1][a+1] - \alpha[m-1][a-1]} \) \( Lma = \frac{1.0}{Jacobian} * \frac{PSD\_m\_r - PSD\_m\_l}{pc[m+1] - pc[m-1]} \)

\( PSD\_a\_r = Dpca[a+1] * Jacobian[a+1] * \frac{PSD[m+1][a+1] - PSD[m-1][a+1]}{pc[m+1][a+1] - pc[m-1][a+1]} \) \( PSD\_a\_l = Dpca[a-1] * Jacobian[a-1] * \frac{PSD[m+1][a-1] - PSD[m-1][a-1]}{pc[m+1][a-1] - pc[m-1][a-1]} \) \( Lam = \frac{1.0}{Jacobian} * \frac{PSD\_a\_r - PSD\_a\_l}{\alpha[a+1] - \alpha[a-1]} \)

\( PSD += dt * (Lma + Lam) \)

Definition at line 414 of file PSD.cpp.

void PSD::Interpolate ( PSD otherPSD,
Parameters_structure::Interpolation  interpolationParameters,
Grid oldGrid,
Grid newGrid,
Matrix2D< double >  newGrid_pc_lowerBoundaryCondition,
Matrix2D< double >  newGrid_pc_upperBoundaryCondition 
)

Interpolation function.

It does interpolation... for log(function) or for just function, depends on parameters.

Parameters
&otherPSD- PSD to interpolate to the values from
interpolationParameters- parameters for interpolation
&oldGrid- old grid
&newGrid- new grid
newGrid_pc_lowerBoundaryCondition- lower energy boundary VALUE
newGrid_pc_upperBoundaryCondition- upper energy boundary VALUE
Todo:
If grid types are the same, it's not necessary that grids are the same

Definition at line 1557 of file PSD.cpp.

void PSD::Load_initial_f ( GridElement L,
GridElement pc,
GridElement alpha,
double  tau,
double  Kp,
double  min_f = 1.e-99,
string  J_L7_function = "J_L7",
double  fb_out = 1,
double  fb_in = 0,
bool  using_ia_90 = true 
)

Calculate initial PSD from steady state.

Parameters
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha
tau- life time
Kp- Kp value
min_f- minimum of function (should be some positive value)
J_L7_function- parameter, J_L7 function. Can be J_L7 or J_L7_corrected
fb_out- value at the outer L boundary
fb_in- value at the inner L boundary
using_ia_90- flag indicating if using ia_90 value for calculating L values

Definition at line 1790 of file PSD.cpp.

void PSD::Load_initial_f_2d ( GridElement L,
GridElement pc,
GridElement alpha,
const char *  filename 
)

Load initial PSD from 2d-file (for 2d calculations).

Parameters
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha
*filename- file name
Todo:
not exactly sure what this does

Definition at line 2227 of file PSD.cpp.

void PSD::Load_initial_f_file ( GridElement L,
GridElement epc,
GridElement alpha,
const char *  filename,
bool  withGrid 
)

Load initial PSD from file

Parameters
&L- grid element L
&epc- grid element pc
&alpha- grid element alpha
filename- filename
withGrid- indicate if the file contains grid values

Definition at line 1742 of file PSD.cpp.

void PSD::Load_initial_f_from_outer_L ( GridElement L,
GridElement pc,
GridElement alpha,
double  tau,
double  Kp,
Matrix2D< double >  L_UpperBoundaryCondition,
double  min_f = 1.e-99,
double  fb_out = 1,
double  fb_in = 0 
)

Calculate initial PSD from steady state using boundary conditions. Load initial PSD from outer L.

Parameters
&L- radial distance
&pc- MeV
&alpha- pitch angle
tau- life time
Kp- Kp index value
L_UpperBoundaryCondition- L's upper boundary condition
min_f- minimum value of vector
fb_out- value of outer boundary
fb_in- value of inner boundary

Definition at line 2183 of file PSD.cpp.

void PSD::Load_initial_f_maxwell ( GridElement L,
GridElement pc,
GridElement alpha 
)

Load initial PSD from maxwell (for 2d calculations).

Parameters
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha

Definition at line 2335 of file PSD.cpp.

void PSD::Load_initial_f_two_zone ( GridElement L,
GridElement pc,
GridElement alpha,
double  tau,
double  tauLpp,
double  Kp,
double  min_f = 1.e-99,
string  J_L7_function = "J_L7",
double  fb_out = 1,
double  fb_in = 0 
)

two zone initial profile (kckim)

Calculate initial PSD from steady state (kckim). Load initial PSD with two zone

Parameters
&L- grid element L
&pc- grid element pc
&alpha- grid element alpha
tau- lifetime
tauLpp- lifetime inside plasmasphere
Kp- scale diffusion
min_f- minimum value possible
J_L7_function- name of boundary flux function at L=7
fb_out- value of outer boundary
fb_in- value of inner boundary

Definition at line 1926 of file PSD.cpp.

void PSD::Output_without_grid ( double  time)

PSD output.

Parameters
time- time

Definition at line 1720 of file PSD.cpp.

void PSD::SourcesAndLosses ( Parameters_structure  parameters,
GridElement L,
GridElement pc,
GridElement alpha,
AdditionalSourcesAndLosses SL,
double  dt,
double  Lpp,
double  tau,
double  tauLpp,
double  Kp 
)

Scaling initial PSD if changed upper L-boundary (not 7Re) Sources and losses term from the FP eq.

Parameters
parameters- parameters to be used
&L- L-grid array
&pc- pc-grid array
&alpha- alpha-grid array
&SL- sources/losses 3D matrix
dt- time step
Lpp- plasmapause location
tau- lifetime outside of the plasmapause
tauLpp- lifetime inside of the plasmapause
Kp- Kp index

Definition at line 319 of file PSD.cpp.


The documentation for this class was generated from the following files: