Re-factored the whole external potential mechanism to match what is done for cooling. Renamed 'disk' to 'disc'.

Give the Gizmo volume the correct dimension in the outputs.

Is this correct ? 

UNIT_CONV_VOLUME is U_L^3
UNIT_CONV_INV_VOLUME is U_L^-3


See merge request !248

Radiated energy



See merge request !247

Align critical structures on 128 bytes and make all their sizes multiples of 32 bytes.

Here is another one for you.

Following the discussion we had earlier about the task structure I have gone through the rest of the code and made sure everything important is aligned on 128 bytes. 

See merge request !240

Collect the radiated energy of all particles when computing statistics and output the result in the log file.

Moved the structure related to cooling into a separate file. Added a structure to the xpart to contain the accumulated radiative information.

Use the cooling time-step condition when computing dt. Better documentation. More uniform interface accross cooling functions.

Drift all particles before repartitioning

Drifting them afterwards can result in some cells having inconsistent particle lists, so this must be done first. This fixes the crashes I was seeing in MPI runs of the SodShock with 4x12 cores.


See merge request !241

Fixed bug in GIZMO_SPH hydro_get_internal_energy.

The energy stored in `conserved.energy` is not the thermal energy (per unit mass) `u`, but the total thermal energy `u*conserved.mass`. Updated `hydro_get_internal_energy` and `hydro_set_entropy` to reflect this.

See merge request !244

Set the default value of the subsize to 64M

Fix to #212. 

See merge request !243

Add the cooling infrastructure and the const_du and const_lambda

@dc-arri1 I have opened this merge request so that we can push the main changes into the master branch to make everything available to everyone.

I'll have a few questions for you here. I'd be great if you could answer them so that I can merge things.

See merge request !242

This shouldn't be necessary when we are drifting all, we have just repartitioned or the particles are just being initialised (before the first step)

Added the new cooling parameters to the parameter documentation. Renamed the parameters to differentiate them.