diff --git a/examples/PertubedBox/makeIC.py b/examples/PertubedBox/makeIC.py index 81ae9e5909ec51cb640209ce759d618311ce811f..a5e831eca02463d287ce2c7748eb780ef66aeb33 100644 --- a/examples/PertubedBox/makeIC.py +++ b/examples/PertubedBox/makeIC.py @@ -63,7 +63,7 @@ for i in range(L): v[index,1] = 0. v[index,2] = 0. m[index] = mass - h[index] = 2.251 * boxSize / L + h[index] = 1.1255 * boxSize / L u[index] = internalEnergy ids[index] = index diff --git a/examples/SedovBlast/makeIC.py b/examples/SedovBlast/makeIC.py index 9b4b5443f472edf5bb299ed5e7261d115c96293c..f3ed3288306a1ff7bb0387896a7e4e0f803a1f35 100644 --- a/examples/SedovBlast/makeIC.py +++ b/examples/SedovBlast/makeIC.py @@ -67,7 +67,7 @@ for i in range(L): v[index,1] = 0. v[index,2] = 0. m[index] = mass - h[index] = 2.251 / 2 * boxSize / L + h[index] = 1.1255 * boxSize / L u[index] = internalEnergy ids[index] = index if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L: diff --git a/examples/SedovBlast/makeIC_fcc.py b/examples/SedovBlast/makeIC_fcc.py index 88cbaf8042323ea91ed7dd09b1bd63418aff3e3f..8a5d50031d06247afdd9b51ab9fe43bcca87a963 100644 --- a/examples/SedovBlast/makeIC_fcc.py +++ b/examples/SedovBlast/makeIC_fcc.py @@ -70,7 +70,7 @@ for i in range(L): v[index,1] = 0. v[index,2] = 0. m[index] = mass - h[index] = 2.251 / 2 * hbox + h[index] = 1.1255 * hbox u[index] = internalEnergy ids[index] = index if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 1.2 * hbox: diff --git a/examples/UniformBox/makeIC.py b/examples/UniformBox/makeIC.py index 9f740398164a3411453f5a55a011428ececaeb19..9f5e3fd7eeafcc6db824a668fd7fdb11c7b472f5 100644 --- a/examples/UniformBox/makeIC.py +++ b/examples/UniformBox/makeIC.py @@ -74,7 +74,7 @@ ds = grp.create_dataset('Masses', (numPart,1), 'f') ds[()] = m m = zeros(1) -h = full((numPart, 1), 2.251 * boxSize / L) +h = full((numPart, 1), 1.1255 * boxSize / L) ds = grp.create_dataset('SmoothingLength', (numPart,1), 'f') ds[()] = h h = zeros(1)