From fd335842ee6d61d2209faa3de11437577669b5b7 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Tue, 22 Jan 2013 22:51:21 +0000
Subject: [PATCH] The read_ic routine now reads the densities as well if
 present in the IC file.

Former-commit-id: a10e6828e22416893ff0528e7fe6b8f88031995f
---
 src/ic.c | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/ic.c b/src/ic.c
index c9af122f0e..5cb07ff708 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 */
-- 
GitLab