a45a380d0c
SOURCES: Consider -B.div(B) terms for isothermal KEPES too.
...
The -B.div(B) source term in the momentum equation was dropped for
the isothermal MHD case. Apparently, its lack causes some numerical
instabilities related to the accumulation of the divergence of B.
Therefore, take it into account for the isothermal case too.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-07-30 16:46:29 -03:00
e76e875004
Update the copyright year to 2024.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-03-07 09:34:43 -03:00
81de98d9e2
Update the copyright year to 2023.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2023-02-01 18:36:37 -03:00
6ba8d5d34f
SOURCES: Reduce the number of arguments to update_sources().
...
Pointer pdata already has a reference to the increment array dU(), so
remove the redundant argument du.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-02-09 16:25:07 -03:00
9829505650
Update copyright year to 2022.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-02-02 09:51:41 -03:00
5c96f3fd5a
SOURCES: Add KEPES source terms for isothermal GLM-MHD.
...
These source terms guarantee the KEPES solver to conserve the entropy.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-01-22 11:57:01 -03:00
a634d9bed5
VARIOUS: Make thread number thread private.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-01-08 11:52:16 -03:00
d8e959b219
SOURCES: Make variable initialization OpenMP conformant.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-01-08 10:52:05 -03:00
226998e80e
SOURCES: Add source terms related to the MHD KEPES solver.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-01-06 16:46:36 -03:00
c6fc8d5194
OPENMP: Declare some local variables as threadprivate.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-12-07 19:55:30 -03:00
ea0b9b83dc
WORKSPACE: Make workspace thread safe.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-12-07 10:46:18 -03:00
67835ffdff
SOURCES: Remove empty if statement.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-12-03 19:55:53 -03:00
25ad944f44
SOURCES: Add procedure interface for extra source terms.
...
A user defined source terms can be pointer to the procedure pointer
'update_extra_sources'.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-26 13:37:25 -03:00
d961a080ce
SOURCES: Use print_message().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-19 12:39:36 -03:00
f659000e8c
PROFILE: Remove all custom profiling.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-16 15:22:15 -03:00
1b1e3ab034
SOURCES: Use resize_workspace() in update_sources().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-13 20:38:27 -03:00
a381fc0cce
MESH: Remove workspace from this module.
...
Use everywhere the workspace provided by module WORKSPACE.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-12 22:35:05 -03:00
d4e442ead3
SOURCES: Use workspace in update_sources().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-11-12 11:53:43 -03:00
94a621a3c5
USER_PROBLEM: Rename use problem subroutines.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-10-26 08:52:03 -03:00
3d16d5c13d
Update copyrights.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2021-02-04 17:35:04 -03:00
9e507b283c
SOURCES: Remove unused variable for 2D case.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-08-15 01:12:41 -03:00
78d520a87c
SOURCES: Remove unused variables.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-08-06 14:33:42 -03:00
ad0303d3ba
Update copyright.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2020-03-02 13:18:53 -03:00
e263e4a4d8
SOURCES: Print anomalous resistivity parameters only if enabled.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-27 10:07:46 -03:00
e68a068d56
SOURCES: Implement anomalous resistivity.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-25 09:53:05 -03:00
abb1f43528
Rename nb => bcells to nn => bcells.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-13 11:44:51 -02:00
b2009f03be
SOURCES: Make initialize_sources() resistant to wrong parameters.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-08 15:37:06 -02:00
969afab57f
SOURCES: Use bcells instead of im, jm, km.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-05 11:17:59 -02:00
9a43c5b19d
SOURCES, HELPERS: Add print_parameter_double() to HELPERS.
...
Use helper subroutines in print_sources().
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-01-30 18:55:41 -02:00
cf3fbf6902
DRIVER, SOURCES, EQUATIONS: Add print_sources().
...
Also add magnetized flag in EQUATIONS to indicate if the equations
system is magnetized. This helps to reduce printing some unnecessary
info.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-01-30 12:20:18 -02:00
b0953016a5
DRIVER, SOURCES: Initialize sources after equations.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-01-29 10:05:36 -02:00
93342fb333
SOURCES: Replace explicit calls to get_parameter_*() with unified one.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-01-28 21:32:37 -02:00
f8b65842af
Move source files from src/ to sources/.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-01-28 09:14:20 -02:00
