diff --git a/src/ic.c b/src/ic.c
index c9af122f0e0c39585e9b3caec9fba8d353e867cf..5cb07ff70823f0b5588064274919e447be154ff8 100644
--- a/src/ic.c
+++ b/src/ic.c
@@ -330,9 +330,11 @@ void read_ic ( char* fileName, double dim[3], struct part **parts,  int* N, int*
   readArray(h_grp, "Mass", FLOAT, *N, 1, *parts, mass, COMPULSORY);
   readArray(h_grp, "SmoothingLength", FLOAT, *N, 1, *parts, h, COMPULSORY);
   readArray(h_grp, "InternalEnergy", FLOAT, *N, 1, *parts, u, COMPULSORY);
-  readArray(h_grp, "TimeStep", FLOAT, *N, 1, *parts, dt, OPTIONAL);
   readArray(h_grp, "ParticleIDs", ULONGLONG, *N, 1, *parts, id, COMPULSORY);
+  readArray(h_grp, "TimeStep", FLOAT, *N, 1, *parts, dt, OPTIONAL);
+  readArray(h_grp, "Acceleration", FLOAT, *N, 3, *parts, a, OPTIONAL);
   readArray(h_grp, "TimeStep", FLOAT, *N, 1, *parts, dt, OPTIONAL );
+  readArray(h_grp, "Density", FLOAT, *N, 1, *parts, rho, OPTIONAL );
   readArray(h_grp, "Acceleration", FLOAT, *N, 3, *parts, a, OPTIONAL );
 
   /* Close particle group */