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

638 Commits

Author SHA1 Message Date
Grzegorz Kowal
0177811f55 Check subroutine check_mesh() to debug refinement. 
- subroutine check_mesh() verifies, before and after the mesh update,
   if the meta blocks pointer fields point to correct blocks;
 2011-04-21 08:39:31 -03:00
Grzegorz Kowal
c24e09fedc Move coordinate variables generation to another subroutine. 2011-04-15 00:46:19 +02:00
Grzegorz Kowal
aa80869645 Move generation if the level resolutions to init_blocks(). 2011-04-15 00:25:22 +02:00
Grzegorz Kowal
3df814ab8a Store evolution parameters in the restart files. 2011-04-15 00:09:50 +02:00
Grzegorz Kowal
2d513508db Write restart files at the end of program execution. 2011-04-14 16:16:58 +02:00
Grzegorz Kowal
1e7dec177c Restore metablock structure during job restart. 2011-04-14 16:14:03 +02:00
Grzegorz Kowal
8348f9ca64 Restore all datablocks for job restart. 2011-04-14 15:01:17 +02:00
Grzegorz Kowal
75e31649fc Add block_array array to store pointers to restored blocks. 2011-04-14 14:34:51 +02:00
Grzegorz Kowal
98c8471123 Fix appending new meta and data blocks. 2011-04-14 13:40:04 +02:00
Grzegorz Kowal
2ea08c6205 Store the block dimensions as a global attribute. 2011-04-14 12:46:38 +02:00
Grzegorz Kowal
342858f12b Allocate all datablocks in read_attributes_h5(). 2011-04-14 10:26:43 +02:00
Grzegorz Kowal
a617de78ba Allocate all metablocks in read_attribute_h5(). 2011-04-14 00:37:41 +02:00
Grzegorz Kowal
9b46eb5c61 Move init_blocks() to the driver. 2011-04-14 00:19:45 +02:00
Grzegorz Kowal
9b8f30b222 Read global attributes in read_data_h5(). 2011-04-13 18:55:55 +02:00
Grzegorz Kowal
2bf24510d2 Implement opening and closing files in read_data_hdf5(). 2011-04-12 18:30:39 +02:00
Grzegorz Kowal
33881a02bc Add checking if restart file exists and is in HDF5 format. 2011-04-11 17:03:06 +02:00
Grzegorz Kowal
5d8efdcdb3 Add and call the initial subroutine to restore data from HDF5 file. 2011-04-11 16:27:08 +02:00
Grzegorz Kowal
944e680553 Add new parameter 'nres' to determine the restart file number. 
- this parameter determines which restart dump should be read in order
   to restart current job;
 2011-04-11 16:08:00 +02:00
Grzegorz Kowal
23ae07b2f6 Add viscous and resistive terms for adiabatic EOS. 
- also slightly optimize the dissipation term calculations;
 2011-04-05 00:42:24 +02:00
Grzegorz Kowal
eb2f8dae50 Add missing parentheses in reconstruct() subroutine. 2011-03-28 11:10:30 -03:00
Grzegorz Kowal
d9ad9c4232 Improve calculation of the refinement criterion. 
- now the refinement criterion uses the second derivative and is
   normalized properly to lay between 0.0 and 1.0; it is calculated in a
   similar way as in the Flash code;
 2011-03-26 19:33:37 -03:00
Grzegorz Kowal
a0f7d5118f Introduce the compression of the HDF5 output files. 2011-03-26 13:31:52 -03:00
Grzegorz Kowal
0f863de622 Optimize LimO3 reconstruction. 2011-03-25 01:05:58 -03:00
Grzegorz Kowal
f9b00e79af Add the Orszag-Tang MHD test problem. 2011-03-24 19:26:57 -03:00
Grzegorz Kowal
2e4273326f Add minmod4 function and use it in MP interpolation. 2011-03-24 16:15:51 -03:00
Grzegorz Kowal
ba92cd4f7a Change the default compiler to GNU Fortran. 2011-03-23 23:34:48 -03:00
Grzegorz Kowal
d05f1d4116 Remove linking to HDF4 libraries. 2011-03-23 22:22:57 -03:00
Grzegorz Kowal
1841ae56f9 Optimize MP reconstruction methods. 2011-03-23 18:01:16 -03:00
Grzegorz Kowal
b3b4b42ad8 Change flag RESIS to RESISTIVITY. 2011-03-22 17:24:28 -03:00
Grzegorz Kowal
1a26c3e204 Optimize hydro Roe solvers and prepare MHD Roe solvers. 2011-03-22 15:36:59 -03:00
Grzegorz Kowal
2320bafe93 Print some info about solved set of equations. 2011-03-21 18:24:55 -03:00
Grzegorz Kowal
ca30b0674c Implement Roe Riemann solver for isothermal hydrodynamics. 2011-03-21 14:58:18 -03:00
Grzegorz Kowal
cacce96f52 Implement Roe Riemann solver for adiabatic hydrodynamics. 
- the method is based on the Athena implementation;
 2011-03-21 14:29:40 -03:00
Grzegorz Kowal
2e77f327b1 Implement viscous terms. 2011-03-18 15:33:24 -03:00
Grzegorz Kowal
5a8cd7d184 Implement 'multi current sheet' problem. 2011-03-18 14:00:27 -03:00
Grzegorz Kowal
4e1e535ec6 Reorganize compilation flags and move them to make.config. 
- in this way all flags are controlled from make.config and host
   specific files choose the right compilers, options, and libraries,
   depending on the chosen flags;
 2011-03-11 19:55:31 -03:00
Grzegorz Kowal
2c3e4eb509 Slightly optimize subroutine fourier_transform(). 2011-03-10 15:56:21 -03:00
Grzegorz Kowal
f07457e084 Reorganize printed informations, parameters and timings. 2011-03-10 15:22:10 -03:00
Grzegorz Kowal
4176930759 Add more timers to forcing subroutines. 
- now each step of forcing (initialization, evolution, Fourier
   transforms) has separate time; this helps to optimize the forcing
   algorithm;
 2011-03-10 14:57:25 -03:00
Grzegorz Kowal
cd91d350ab Introduce a compilation flag REFINE. 
- the flag REFINE determines if the code should be compiler with the
   adaptive mesh refining/derefining; otherwise, the static mesh is
   generated;
 2011-03-10 01:08:00 -03:00
Grzegorz Kowal
1169454cde Correct calculation of the forcing component amplitudes. 
- correct calculating the amplitudes of the forcing components;
   micro optimizations of the force initialization;

 - correct the evolution of the forcing components so the power input
   corresponds exactly to the defined power fpow;
 2011-03-10 01:04:06 -03:00
Grzegorz Kowal
aa5ecd5d38 FORCE: reduce number of components and fix constants. 
- we don't need all components from a sphere in the Fourier space, only
   half components are independent, so reduce it's number;
 - correct the amplitude calculation so the proper amount of energy is
   injected in the system;
 2011-03-08 23:47:18 -03:00
Grzegorz Kowal
01a1f671bb FORCE: adjust the power coefficient for forcing. 2011-03-08 22:39:38 -03:00
Grzegorz Kowal
0638a74d7b INTEGRALS: add forcing correlations to the integrals file. 
- calculate the force-force correlation in the real and Fourier spaces
   and store them in the integrals file;
 2011-03-08 21:56:47 -03:00
Grzegorz Kowal
dfe54acf48 FORCE: fix Fourier transform for multiblock domains. 2011-03-07 17:23:49 -03:00
Grzegorz Kowal
c597a3ef43 FORCE: correct sign in the inverse Fourier transform. 2011-03-07 01:06:00 -03:00
Grzegorz Kowal
28c525bcfc FORCE: implement condition for minimizing velocity-force correlation. 2011-03-07 00:54:46 -03:00
Grzegorz Kowal
cfb35e7cd9 MPI: implement reduction of a complex array. 
- reduce the velocity Fourier components from all processors;
 2011-03-07 00:28:58 -03:00
Grzegorz Kowal
de3d89bb1b FORCE: implement and use Fourier transform of velocity. 
- we need to calculate Fourier transform of velocity for all forcing
   components in order to minimize the velocity-force correlation;
 2011-03-07 00:08:31 -03:00
Grzegorz Kowal
605cb1e370 FORCE: keep the forcing amplitudes in a separate array. 2011-03-06 23:53:36 -03:00