The read_ic routine now reads the densities as well if present in the IC file.

Former-commit-id: a10e6828e22416893ff0528e7fe6b8f88031995f

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 */