VERB_code_2.3
|
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 | |
Phase Space Density class.
Do diffusions and store the result in Matrix3D. Do all possible operations with PSD (like output, loading, etc).
void PSD::Create_Initial_PSD | ( | Parameters_structure::PSD | parameters, |
Grid & | grid, | ||
BoundaryCondition | L_UpperBoundaryCondition | ||
) |
Create inital PSD (Steady State)
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;
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 | ||
) |
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.
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) |
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.
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 |
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.
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) |
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.
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 |
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.
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 |
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.
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) \)
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.
&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 |
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.
&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 |
void PSD::Load_initial_f_2d | ( | GridElement & | L, |
GridElement & | pc, | ||
GridElement & | alpha, | ||
const char * | filename | ||
) |
void PSD::Load_initial_f_file | ( | GridElement & | L, |
GridElement & | epc, | ||
GridElement & | alpha, | ||
const char * | filename, | ||
bool | withGrid | ||
) |
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.
&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 |
void PSD::Load_initial_f_maxwell | ( | GridElement & | L, |
GridElement & | pc, | ||
GridElement & | alpha | ||
) |
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
&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 |
void PSD::Output_without_grid | ( | double | time | ) |
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 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 |