diff --git a/examples/SodShock/makeIC.py b/examples/SodShock/makeIC.py index 027e727d38cf18647a7247ef79b37423dc16aa59..105acce134ae21069e9ecd1278e541516943827b 100644 --- a/examples/SodShock/makeIC.py +++ b/examples/SodShock/makeIC.py @@ -37,13 +37,13 @@ fileName = "sodShock.hdf5" #--------------------------------------------------- #Read in high density glass -glass1 = h5py.File("../Glass/glass_12500.hdf5") -pos1 = glass1["/PartType0/Coordinates"][:,:] +glass1 = h5py.File("../Glass/glass_200000.hdf5") +pos1 = glass1["/PartType1/Coordinates"][:,:] pos1 = pos1 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5] #Read in high density glass -glass2 = h5py.File("../Glass/glass_3125.hdf5") -pos2 = glass2["/PartType0/Coordinates"][:,:] +glass2 = h5py.File("../Glass/glass_50000.hdf5") +pos2 = glass2["/PartType1/Coordinates"][:,:] pos2 = pos2 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5] #Generate high density region