diff --git a/csds b/csds
index 62c0b35675e88a76615154dcbdf0122755be7e56..30524e2c425e276289218025e4afc68425a65176 160000
--- a/csds
+++ b/csds
@@ -1 +1 @@
-Subproject commit 62c0b35675e88a76615154dcbdf0122755be7e56
+Subproject commit 30524e2c425e276289218025e4afc68425a65176
diff --git a/src/csds_io.c b/src/csds_io.c
index d77812abe23bb760cdb4b4289ace127238a5c695..1c0145a9e5672af4b86d31b5d55712a4a66a2bd6 100644
--- a/src/csds_io.c
+++ b/src/csds_io.c
@@ -102,6 +102,19 @@ void csds_write_description(struct csds_writer* log, struct engine* e) {
   fprintf(f, "  NumberGParts: %li\n", e->s->nr_gparts);
   fprintf(f, "\n");
 
+  /* Write the cosmology */
+  fprintf(f, "Cosmology:\n");
+  fprintf(f, "  Omega_cdm: %g\n", e->cosmology->Omega_cdm);
+  fprintf(f, "  Omega_lambda: %g\n", e->cosmology->Omega_lambda);
+  fprintf(f, "  Omega_b: %g\n", e->cosmology->Omega_b);
+  fprintf(f, "  Omega_r: %g\n", e->cosmology->Omega_r);
+  fprintf(f, "  Omega_k: %g\n", e->cosmology->Omega_k);
+  fprintf(f, "  Omega_nu_0: %g\n", e->cosmology->Omega_nu_0);
+  fprintf(f, "  w_0: %g\n", e->cosmology->w_0);
+  fprintf(f, "  w_a: %g\n", e->cosmology->w_a);
+  fprintf(f, "  Hubble0: %g\n", e->cosmology->H0);
+  fprintf(f, "\n");
+
   /* Write unit system */
   const struct unit_system* us = e->internal_units;
   fprintf(f, "InternalUnitSystem:\n");