Also write the physical constants when writing files in parallel

src/parallel_io.c

@@ -1153,6 +1153,9 @@ void prepare_file(struct engine* e, const char* baseName, long long N_total[6],
   /* Print the run's policy */
   io_write_engine_policy(h_file, e);
 
+  /* Print the physical constants */
+  phys_const_print_snapshot(h_file, e->physical_constants);
+
   /* Print the SPH parameters */
   if (e->policy & engine_policy_hydro) {
     h_grp = H5Gcreate(h_file, "/HydroScheme", H5P_DEFAULT, H5P_DEFAULT,

src/serial_io.c

@@ -1018,6 +1018,9 @@ void write_output_serial(struct engine* e, const char* baseName,
     /* Print the run's policy */
     io_write_engine_policy(h_file, e);
 
+    /* Print the physical constants */
+    phys_const_print_snapshot(h_file, e->physical_constants);
+
     /* Print the SPH parameters */
     if (e->policy & engine_policy_hydro) {
       h_grp = H5Gcreate(h_file, "/HydroScheme", H5P_DEFAULT, H5P_DEFAULT,