Skip to content
Snippets Groups Projects

Stefan add readme files

Merged Matthieu Schaller requested to merge Stefan_add_README_files into stats_include_e_pot_ext
+ 425
21
Compare changes
  • Side-by-side
  • Inline
Files
+ 25
0
;
;To make the initial conditions we distribute gas particles randomly in a cube
;with a side length twice that of the virial radius. The density profile
;of the gas is proportional to r^(-2) where r is the distance from the centre
;of the cube.
;
;The parameter v_rot (in makeIC.py and cooling.yml) sets the circular velocity
;of the halo, and by extension, the viral radius, viral mass, and the
;internal energy of the gas such that hydrostatic equilibrium is achieved.
;
;While the system is initially in hydrostatic equilibrium, the cooling of the
;gas means that the halo will collapse.
;
;To run this example, make such that the code is compiled with either the
;isothermal potential or softened isothermal potential, and 'const_lambda'
;cooling, set in src/const.h. In
;the latter case, a (small) value of epsilon needs to be set in cooling.yml.
;0.1 kpc should work well.
;
;The plotting scripts produce a plot of the density, internal energy and radial
;velocity profile for each snapshot. test_energy_conservation.py shows the
;evolution of energy with time. These can be used to check if the example
;has run properly.
;
;
\ No newline at end of file
Loading