2213 Commits

Author SHA1 Message Date
3d16d5c13d Update copyrights.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-02-04 17:35:04 -03:00
8c0e6dea29 PYTHON: Add argument 'maxlev' to function dataset().
This argument is an alternative way to reduce the resolution of the
data. It specifies that the output cube resolution should be
corresponding to the 'maxlev' refinement level.

If the maximum level of simulation changes during the evolution, the
parameter allows to read all snapshots with the same output resolution.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-01-05 10:48:00 -03:00
575bcc1d27 PYTHON: Move shrink and interpolation arguments to dataset().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-01-05 10:17:44 -03:00
b3c51569f1 PYTHON: Correct the sign of the electric field components.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-12-30 22:10:07 -03:00
ca1cc9fb31 COORDINATES: Add indices for nb-1, nb+1, ne-1, ne+1.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-12-29 17:42:20 -03:00
4833067bea SCHEMES: Improve detection of numerical degeneracies in HLLD solvers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-12-27 11:37:07 -03:00
052ce7e05a PYTHON: Implement reading subvolume of the datasets in AmunXML.
The subvolume to be read is controlled by additional parameter 'extent'
passed to function dataset(). The syntax of this parameters is

  extent=(xmin, xmax, ymin, ymax, zmin, zmax)

for 3D snapshots or

  extent=(xmin, xmax, ymin, ymax)

for 2D snapshots.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-12-20 22:49:14 -03:00
9ae09e2c46 PYTHON: Fix block counter for 'divv' and 'divb' in AmunXML.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-12-02 17:43:37 -03:00
37fdafd28d EVOLUTION: Fix unphysical cells after applying shapes.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-26 10:50:49 -03:00
f453d86bc4 PYTHON: Add logarithm of density and pressure to AmunXML.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-26 09:52:28 -03:00
2dfcecf789 PYTHON: Fix calculation of temperature in class AmunXML.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-26 09:48:46 -03:00
7fdac2f584 SCHEMES: Make HD ROE solvers numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-17 17:52:22 -03:00
56fb625e59 SCHEMES: Make adiabatic MHD ROE solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-17 17:49:05 -03:00
52747633f8 SCHEMES: Make isothermal MHD ROE solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-17 17:27:37 -03:00
c677e18ecf INTEGRALS: Add subroutine to flush and sync specific IO unit.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-12 07:22:49 -03:00
b2f2b35d0b INTEGRALS: Flush integral buffers.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-12 06:33:56 -03:00
0b1da7d42d INTERPOLATIONS: Enable positivity fix by default.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-11 08:49:27 -03:00
0b6797ede7 INTERPOLATIONS: OCMP7 method requires 6 ghost cells.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 14:54:17 -03:00
16a7842fc0 SCHEMES: Make adiabatic MHD HLLD solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 12:04:03 -03:00
ff7ba2e015 SCHEMES: Make isothermal MHD HLLD solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 12:03:03 -03:00
44b990644d SCHEMES: Make MHD HLLC solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 12:02:16 -03:00
87e748feb7 SCHEMES: Make HD HLLC solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 12:01:28 -03:00
4a4c363660 SCHEMES: Make HLL solver numerically symmetric.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 11:58:17 -03:00
f85c119043 PROBLEMS: Slightly rewrite pressures in tearing problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 10:13:31 -03:00
6150e8bbde BOUNDARIES: Fix bug in specific boundaries along X.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-10 10:05:53 -03:00
277b82ab67 PROBLEMS: Move initialization/finalization of user problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-09 21:41:40 -03:00
18a832f147 PYTHON: Fix detecting magnetic energy in AmunXML.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-11-09 18:58:29 -03:00
fa2b1d05f7 Revert "INTERPOLATION: Improve MP limiting for high-wavenumbers."
This reverts commit 3da83e901d4dd1831d666d3a12e1ee5662f2beb6.
2020-10-25 20:19:20 -03:00
b589ffede9 INTERPOLATIONS: Update coefficients for OCMP methods.
The coefficients were tuned to keep the same maximum disspersion error
while shifting it toward higher frequencies with the scheme order. At
the same time the dissipation decreses at high frequencies with the
scheme order.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-14 08:48:56 -03:00
77b18d0d2b INTERPOLATIONS: Implement 9th order Optimized Compact MP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-11 18:28:56 -03:00
173f13e8b7 INTERPOLATIONS: Remove comments from 5th order OCMP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-09 22:41:10 -03:00
be07cdaaaa INTERPOLATIONS: Implement 7th order Optimized Compact MP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-09 22:27:19 -03:00
7bf0635c6e INTERPOLATIONS: Set central_weight to zero by default.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 22:30:50 -03:00
71290e5731 INTERPOLATIONS: Remove regular compact MP methods.
They are identical to low-dissipation methods with central_weight = 0.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 22:28:55 -03:00
325dc36c44 INTERPOLATIONS: Merge regular explicit MP methods.
They are identical to low-dissipation versions with center_weight = 0.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 22:19:10 -03:00
662f35ff01 INTERPOLATIONS: Implement 9th order low-dissipation Compact MP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 22:08:55 -03:00
20bae9dec1 INTERPOLATIONS: Implement 9th order Compact MP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 19:08:21 -03:00
2d31bb2a4f ALGEBRA: Implement pentadiagona linear equations solver.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 18:31:05 -03:00
352af9a1a6 INTERPOLATIONS: Rewrite compact low-dissipation methods.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 18:19:05 -03:00
a19e8b885f INTERPOLATIONS: Rewrite compact MP methods.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 17:25:19 -03:00
2e2ff48693 INTERPOLATIONS: Rewrite explicit low-dissipation methods.
Introduce a parameter cweight to control the weight toward the central
scheme and so the amount of schemes' dissipation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 16:20:50 -03:00
64cc9794e2 INTERPOLATIONS: Move explicit interpolation coefficients.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 14:17:55 -03:00
a8d85390b4 INTERPOLATIONS: Update coefficients of the MP7LD method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 12:28:17 -03:00
a93bf0daea INTERPOLATIONS: Determine ci5(:) coefficients using di5(:) in OCMP5.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-10-07 08:46:42 -03:00
25dbf151e1 ALGEBRA: Rewrite slightly the tridiagonal solver.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-29 12:04:59 -03:00
168cb03f9f INTERPOLATIONS: Implement 5th order Optimized Compact MP method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-28 11:04:19 -03:00
3da83e901d INTERPOLATION: Improve MP limiting for high-wavenumbers.
This improvement is according to Ahn & Lee, JPC, 2020, 83, 34.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-28 09:58:59 -03:00
098122cf3f INTERPOLATIONS: Add low dissipation version of CRMP7 method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-26 10:05:33 -03:00
ddb2117660 INTERPOLATIONS: Add explicit MP methods with lower dissipation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-25 19:18:13 -03:00
deede9452d EVOLUTION: Determine the relative maximum error.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-09-16 06:40:50 -03:00