Tidy hydro debug

src/hydro/Minimal/hydro_debug.h

@@ -36,15 +36,18 @@
 __attribute__((always_inline)) INLINE static void hydro_debug_particle(
     const struct part* p, const struct xpart* xp) {
   printf(
-      "x=[%.3e,%.3e,%.3e], "
-      "v=[%.3e,%.3e,%.3e],v_full=[%.3e,%.3e,%.3e] \n a=[%.3e,%.3e,%.3e], "
-      "u=%.3e, du/dt=%.3e v_sig=%.3e, P=%.3e\n"
-      "h=%.3e, dh/dt=%.3e wcount=%d, m=%.3e, dh_drho=%.3e, rho=%.3e, "
-      "time_bin=%d\n",
-      p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2], xp->v_full[0],
-      xp->v_full[1], xp->v_full[2], p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
-      p->u, p->u_dt, p->force.v_sig, hydro_get_comoving_pressure(p), p->h,
-      p->force.h_dt, (int)p->density.wcount, p->mass, p->density.rho_dh, p->rho,
+      "\n "
+      "x=[%.6g, %.6g, %.6g], v=[%.3g, %.3g, %.3g], \n "
+      "v_full=[%.3g, %.3g, %.3g], a=[%.3g, %.3g, %.3g], \n "
+      "m=%.3g, u=%.3g, du/dt=%.3g, P=%.3g, c_s=%.3g, \n "
+      "v_sig=%.3g, h=%.3g, dh/dt=%.3g, wcount=%.3g, rho=%.3g, \n "
+      "dh_drho=%.3g, time_bin=%d \n",
+      p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2],
+      xp->v_full[0], xp->v_full[1], xp->v_full[2],
+      p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
+      p->mass, p->u, p->u_dt, hydro_get_comoving_pressure(p),
+      p->force.soundspeed, p->force.v_sig, p->h, p->force.h_dt,
+      p->density.wcount, p->rho, p->density.rho_dh,
       p->time_bin);
 }
 

src/hydro/MinimalMultiMat/hydro_debug.h

@@ -38,15 +38,18 @@
 __attribute__((always_inline)) INLINE static void hydro_debug_particle(
     const struct part* p, const struct xpart* xp) {
   printf(
-      "x=[%.3e,%.3e,%.3e], "
-      "v=[%.3e,%.3e,%.3e],v_full=[%.3e,%.3e,%.3e] \n a=[%.3e,%.3e,%.3e], "
-      "u=%.3e, du/dt=%.3e v_sig=%.3e, P=%.3e\n"
-      "h=%.3e, dh/dt=%.3e wcount=%d, m=%.3e, dh_drho=%.3e, rho=%.3e, "
-      "time_bin=%d, mat_id=%d\n",
-      p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2], xp->v_full[0],
-      xp->v_full[1], xp->v_full[2], p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
-      p->u, p->u_dt, p->force.v_sig, hydro_get_comoving_pressure(p), p->h,
-      p->force.h_dt, (int)p->density.wcount, p->mass, p->density.rho_dh, p->rho,
+      "\n "
+      "x=[%.6g, %.6g, %.6g], v=[%.3g, %.3g, %.3g], \n "
+      "v_full=[%.3g, %.3g, %.3g], a=[%.3g, %.3g, %.3g], \n "
+      "m=%.3g, u=%.3g, du/dt=%.3g, P=%.3g, c_s=%.3g, \n "
+      "v_sig=%.3g, h=%.3g, dh/dt=%.3g, wcount=%.3g, rho=%.3g, \n "
+      "dh_drho=%.3g, time_bin=%d, mat_id=%d \n",
+      p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2],
+      xp->v_full[0], xp->v_full[1], xp->v_full[2],
+      p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
+      p->mass, p->u, p->u_dt, hydro_get_comoving_pressure(p),
+      p->force.soundspeed, p->force.v_sig, p->h, p->force.h_dt,
+      p->density.wcount, p->rho, p->density.rho_dh,
       p->time_bin, p->mat_id);
 }