3ecfc9bf92
Update .gitlab-ci.yml
2019-10-04 20:15:49 +00:00
090c4d7f48
Update .gitlab-ci.yml
2019-10-04 20:12:31 +00:00
c8bfe39ad5
Add .gitlab-ci.yml
2019-10-04 20:05:13 +00:00
f36979d74a
SCHEMES: Let user control the high order flux correction - corrected.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-04 11:12:05 -03:00
7578f1e9f2
Revert "SCHEMES: Let user control the high order flux correction."
...
This reverts commit 0dcc3cc95c384c8ed63d6b62aa1cdbcaa58e6672.
2019-10-04 11:11:05 -03:00
0dcc3cc95c
SCHEMES: Let user control the high order flux correction.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-04 11:02:05 -03:00
0806bfc1e0
EQUATIONS: Do not specify variable range if not necessary.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 22:17:09 -03:00
92b905cecd
SCHEMES: Pass only flux variables to Riemann solvers.
...
There is no need to interpolate or pass all variables to the Riemann
solvers. Just pass those related to flux calculation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 21:48:04 -03:00
e26d9b23e4
PROBLEMS: Implement passive scalar for KH test problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 21:14:05 -03:00
6506083785
EQUATIONS: Update scalars in correct_unphysical_states() too.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:52:34 -03:00
c1f8233613
BOUNDARIES: Update conservative scalars after boundary update too.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:47:33 -03:00
fe4a712c76
PROBLEMS: Calculate conservative representation of scalars in RT.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:42:22 -03:00
1037377414
PROBLEMS: Add passive scalar initialization for RT problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:39:16 -03:00
31ac085d70
EQUATIONS: Take into account passive scalars in prim2cons().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:38:27 -03:00
f22cc55b58
EQUATIONS: Take into account passive scalars in cons2prim().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-03 13:27:47 -03:00
8e9e8f196c
EVOLUTION: Update the increment for passive scalars.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-02 21:16:29 -03:00
a80b8cf41f
BOUNDARIES, DRIVER, EVOLUTION: Use nf for the number of block fluxes.
...
Now, since the number of variables and fluxes may be different, it is
important to keep the size of the block flux array using nf instead of
nv.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-02 15:39:25 -03:00
a5094d1e6b
EQUATIONS: Print number of fluxes and passive scalars.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-02 14:03:12 -03:00
5760b6c40c
EQUATIONS: Add parameters for the number of fluxes and passive scalars.
...
Initialize these parameters and generate the names of passive scalars.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-02 12:28:33 -03:00
0133adfc53
MPITOOLS, BOUNDARIES: Simplify arguments of exchange_real_arrays().
...
Some arguments, like sproc or stag, are redundant, so remove them.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-10-01 11:15:31 -03:00
05ec954e16
COORDINATES: Print zmn & zmx, not ymn & ymx if zmn > zmx.
...
Reported-by: Lorenzo Nativi <lorenzo.nativi@astro.su.se>
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-09-23 22:32:45 -03:00
5847a4f086
SCHEMES: Correct numerical fluxes for higher-orders.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-30 22:50:08 -03:00
368894560b
INTERPOLATIONS: Add order module parameter.
...
This parameter keeps the order of the used reconstruction.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-30 21:12:28 -03:00
f0ca778eb6
PYTHON: Fix subroutines to read integrals.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-11 11:53:06 -03:00
4b7615e72f
INTEGRALS: Use 'magnetized' instead of 'ibx > 0'.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-08 07:03:39 -03:00
055f67c874
MESH: Replace allocatable array with module one in prolong_block().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-07 11:12:49 -03:00
7b241c90a8
MESH: Add status argument to prolong_block().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-07 10:51:00 -03:00
81294d958f
INTERPOLATIONS: Implement 9th order MP9 reconstruction.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-07 10:32:37 -03:00
cc29921715
MESH: Skip higher levels, if no need for refinement in generate_mesh().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-04 22:49:49 -03:00
55ad33889c
INTEGRALS: Introduce 'integrals_append' and 'statistics_append'.
...
These parameters allow to control if the integrals and statistics files
should be appended or new files should be created for the restarted
jobs.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-03-04 22:42:24 -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
a982743d4c
PROBLEMS: Change the magnetic field transition profile in current sheet.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-26 20:56:21 -03:00
24f2098163
PROBLEMS: Rewrite the current sheet test problem.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-26 09:25:14 -03:00
e68a068d56
SOURCES: Implement anomalous resistivity.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-25 09:53:05 -03:00
9a53103f14
PYTHON: Determine the maximum level from the blocks.
...
In some situations, the block refinement didn't reach the maximum level
available. There is no need to rescale all block to this level then.
Also the memory consumption in such case is significantly reduced.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-23 08:51:07 -03:00
b13e728819
MESH: Improve printing of the currently generated mesh level.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-22 20:17:24 -03:00
76d7cf8826
IO, DRIVER, PYTHON: Fix domain block dimensions with job restart.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-21 15:58:38 -03:00
e1ee6f39bb
DRIVER: There is no need to limit tm_exec.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 17:07:52 -03:00
a66cd2e27f
DRIVER: Fix ETA and EXECUTION time calculation.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 16:57:35 -03:00
5086bf6bd4
DRIVER: Print problem name before initializations.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 12:07:02 -03:00
0ceca17a95
PYTHON: Add refinement 'level' to be read as dataset.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 11:27:47 -03:00
81b94a8009
PYTHON: Remove wrong condition in amun_compatible().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 11:11:16 -03:00
001748ff92
PYTHON: Support older snapshots with 'rdims'.
...
The new snapshots store 'domain_base_dims' instead.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 11:04:55 -03:00
eccc933c20
PYTHON: Check if the shrink parameter is correct.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-18 10:53:47 -03:00
00967c153e
EVOLUTION: Small changes to evolve_rk43().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-13 15:11:47 -02:00
abb1f43528
Rename nb => bcells to nn => bcells.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-13 11:44:51 -02:00
230b8eb9bf
EVOLUTION: Fix SSPRK(10,4) method.
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-12 17:36:27 -02:00
da9e895406
DRIVER: Use check_status() to verify subroutine status flag.
...
In addition, print error in case of any initialization or finalization
subroutine indicates problems.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-12 13:59:17 -02:00
cbc1a8335f
MPITOOLS: Add function check_status().
...
This function takes one logical input, and if MPI is used reduces it
among all MPI processes with logical OR. Otherwise, it returns the input
value.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-12 12:15:14 -02:00
bf6847e5fd
MESH: Add status argument to refine/derefine_selected_blocks().
...
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2019-02-12 11:06:34 -02:00
