Commit b1e23450 authored by Bert Vandenbroucke's avatar Bert Vandenbroucke
Browse files

Evrard collapse now works with GIZMO.

parent f8546e71
......@@ -9,15 +9,15 @@ InternalUnitSystem:
# Parameters governing the time integration
TimeIntegration:
time_begin: 0. # The starting time of the simulation (in internal units).
time_end: 3. # The end time of the simulation (in internal units).
time_end: 0.8 # The end time of the simulation (in internal units).
dt_min: 1e-7 # The minimal time-step size of the simulation (in internal units).
dt_max: 1e-2 # The maximal time-step size of the simulation (in internal units).
dt_max: 1e-3 # The maximal time-step size of the simulation (in internal units).
# Parameters governing the snapshots
Snapshots:
basename: evrard # Common part of the name of output files
time_first: 0. # Time of the first output (in internal units)
delta_time: 0.03 # Time difference between consecutive outputs (in internal units)
delta_time: 0.1 # Time difference between consecutive outputs (in internal units)
# Parameters governing the conserved quantities statistics
Statistics:
......@@ -39,3 +39,6 @@ Gravity:
# Parameters related to the initial conditions
InitialConditions:
file_name: ./evrard.hdf5 # The file to read
PhysicalConstants:
G: 1.
......@@ -645,18 +645,6 @@ __attribute__((always_inline)) INLINE static void hydro_kick_extra(
/* Make sure the gpart knows the mass has changed. */
p->gpart->mass = p->conserved.mass;
#if !defined(EOS_ISOTHERMAL_GAS)
/* If the energy needs to be updated, we need to do it before the momentum
is updated, as the old value of the momentum enters the equations. */
p->conserved.energy += dt * (p->conserved.momentum[0] * a_grav[0] +
p->conserved.momentum[1] * a_grav[1] +
p->conserved.momentum[2] * a_grav[2]);
p->conserved.energy += dt * (a_grav[0] * p->gravity.mflux[0] +
a_grav[1] * p->gravity.mflux[1] +
a_grav[2] * p->gravity.mflux[2]);
#endif
/* Kick the momentum for half a time step */
/* Note that this also affects the particle movement, as the velocity for
the particles is set after this. */
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment