The coefficient correspond to the maximum of the imaginary part of the
modified wave number to be 2e-6 and the integrated error to be 1.54e-11.
The subroutine was renamed to to OCMP9P, as pentadiagonal.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The coefficient correspond to the maximum of the imaginary part of the
modified wave number to be 2e-6 and the integrated error to be 1.38e-11.
The subroutine was renamed to to OCMP7P, as pentadiagonal.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The coefficient correspond to the maximum of the imaginary part of the
modified wave number to be 2e-6 and the integrated error to be 6.4e-12.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
There two variants selected using parameter 'ocmp_scheme_mode':
1) 'stable' - which guarantees the scheme numerical instability by
enforcing the stable dissipation errors
2) 'minimum' - which enforced the dissipation errors to be smaller than
2e-5
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
Restore the original Suresh & Huynh limiter and add modification by
Ahn & Lee for compact schemes. Additionally, pass the variable
positivity flag and for positive variables modifie the lower limit
using the minimum value obtained from the cubic spline interpolation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This is an improvement in the case a job restarted with a higher
number of processes. The reading of data blocks is spread among new
processes, so it better balances memory.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The thikness of the current sheet is set to the square root of the
resistivity, if it is not set explictely in the parameter file.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This is the Alfvén speed along the X direction defining the velocity
scale, time scale, and the Lundquist number.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
Parameter 'zeta' controls if the thermal pressure (zeta=0.0) or
guide field (zeta=1.0) should contribute to make the total
pressure uniform.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
By default this parameter is set to zero and the behavior is the same as
before. However, if this parameter is set to a value larger than zero,
it determines the driving time scale for the Ornstein–Uhlenbeck method,
and the injection power for the Alvelius method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
