Added a missing division by E(z) for the line of sight velocity

This commit is contained in:
Guilhem Lavaux 2013-07-18 17:34:52 +02:00
parent cd2bef15b3
commit 46b9d3ec40

View file

@ -141,7 +141,7 @@ void metricTransform(SimuData *data, int axis, bool reshift, bool pecvel, double
expfact[i] = z / z_old; expfact[i] = z / z_old;
// Add peculiar velocity // Add peculiar velocity
if (pecvel) if (pecvel)
z += v/100; z += v/(100*e_computer(z0));
} }
catch(const InvalidRangeException& e) { catch(const InvalidRangeException& e) {
cout << "Trying to interpolate out of the tabulated range." << endl; cout << "Trying to interpolate out of the tabulated range." << endl;