Josh Borrow · Matthieu Schaller · Matthieu Schaller
--- a/src/hydro/SPHENIX/hydro_io.h
+++ b/src/hydro/SPHENIX/hydro_io.h

  list[10] = io_make_output_field_convert_part(
      "DiffusionParameters", FLOAT, 1, UNIT_CONV_NO_UNITS, 0.f, parts, xparts,
      convert_diffusion, "Diffusion coefficient (alpha_diff) of the particles");
+
+  list[11] = io_make_output_field(
+      "LaplacianInternalEnergies", FLOAT, 1, UNIT_CONV_FREQUENCY_SQUARED,
+      5.f - 3.f * hydro_gamma, parts, diffusion.laplace_u,
+      "Laplacian (del squared) of the Internal Energy per "
+      "unit mass of the particles");
+
+  list[12] = io_make_output_field(
+      "VelocityDivergences", FLOAT, 1, UNIT_CONV_FREQUENCY, 0.f, parts,
+      viscosity.div_v,
+      "Local velocity divergence field around the particles. Provided without "
+      "cosmology, as this includes the Hubble flow.");