gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity
1,913 Commits 5 Branches 0 Tags
Commit Graph

214 Commits

Author SHA1 Message Date
Grzegorz Kowal
87216b1b73 EQUATIONS: Add the vector of variable indices for update_flux(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:56 -03:00
Grzegorz Kowal
60539e0558 SCHEMES: Move subroutines update_flux() up. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:52 -03:00
Grzegorz Kowal
be548fde57 SCHEMES: Separate HLLC Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:45 -03:00
Grzegorz Kowal
dc44a3ba41 SCHEMES: Separate HLLC Riemann solver for adiabatic HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:40 -03:00
Grzegorz Kowal
f3aa475f63 SCHEMES: Unify HLL Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:33 -03:00
Grzegorz Kowal
bf9ef52fb9 SCHEMES: Rename riemann pointer to numerical_flux. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:28 -03:00
Grzegorz Kowal
c54505917a SCHEMES: Use rather subroutines instead of pointers from equations. 
Use compatible (for the same set of equations and equation of state) in
Riemann solvers, not the procedure pointers.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:23 -03:00
Grzegorz Kowal
9801ac6d8b SCHEMES: Fix description of riemann_hll(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:19 -03:00
Grzegorz Kowal
fed3900810 SCHEMES: Separate universal HLL solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:14 -03:00
Grzegorz Kowal
1154cb5ec0 SCHEMES: Use higher_order_flux_correction(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:10 -03:00
Grzegorz Kowal
1acb7cca04 SCHEMES: Add a separate subroutine for higher order flux corrections. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:04 -03:00
Grzegorz Kowal
ad37de3a99 EQUATIONS, SCHEMES: Merge speeds cm and cp into one variable. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:43:55 -03:00
Grzegorz Kowal
f3670ddf37 EQUATIONS: Reorganize variable indices. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:43:45 -03:00
Grzegorz Kowal
ad0303d3ba Update copyright. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-02 13:18:53 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
7578f1e9f2 Revert "SCHEMES: Let user control the high order flux correction." 
This reverts commit 0dcc3cc95c384c8ed63d6b62aa1cdbcaa58e6672.
 2019-10-04 11:11:05 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
8e9e8f196c EVOLUTION: Update the increment for passive scalars. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-02 21:16:29 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
5847a4f086 SCHEMES: Correct numerical fluxes for higher-orders. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-30 22:50:08 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
4b7615e72f INTEGRALS: Use 'magnetized' instead of 'ibx > 0'. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-08 07:03:39 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
7b241c90a8 MESH: Add status argument to prolong_block(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-07 10:51:00 -03:00
Grzegorz Kowal
81294d958f INTERPOLATIONS: Implement 9th order MP9 reconstruction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-07 10:32:37 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
24f2098163 PROBLEMS: Rewrite the current sheet test problem. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-26 09:25:14 -03:00
Grzegorz Kowal
e68a068d56 SOURCES: Implement anomalous resistivity. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-25 09:53:05 -03:00
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
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
Grzegorz Kowal
a66cd2e27f DRIVER: Fix ETA and EXECUTION time calculation. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-18 16:57:35 -03:00
Grzegorz Kowal
5086bf6bd4 DRIVER: Print problem name before initializations. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-18 12:07:02 -03:00
Grzegorz Kowal
00967c153e EVOLUTION: Small changes to evolve_rk43(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-13 15:11:47 -02:00