Skip to content
Snippets Groups Projects
Commit 2c057c42 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Updated the glass pre-ICs and the sod shock script. 

Former-commit-id: 367c2893c4514c30f1570637f950a57f9e184a5c
parent dc623367
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
examples/SodShock/makeIC.py
+ 4
− 4
View file @ 2c057c42
...@@ -37,13 +37,13 @@ fileName = "sodShock.hdf5" ...@@ -37,13 +37,13 @@ fileName = "sodShock.hdf5"
#--------------------------------------------------- #---------------------------------------------------
#Read in high density glass #Read in high density glass
glass1 = h5py.File("../Glass/glass_12500.hdf5") glass1 = h5py.File("../Glass/glass_200000.hdf5")
pos1 = glass1["/PartType0/Coordinates"][:,:] pos1 = glass1["/PartType1/Coordinates"][:,:]
pos1 = pos1 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5] pos1 = pos1 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5]
#Read in high density glass #Read in high density glass
glass2 = h5py.File("../Glass/glass_3125.hdf5") glass2 = h5py.File("../Glass/glass_50000.hdf5")
pos2 = glass2["/PartType0/Coordinates"][:,:] pos2 = glass2["/PartType1/Coordinates"][:,:]
pos2 = pos2 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5] pos2 = pos2 / 20. # Particles are in [0:0.5, 0:0.5, 0:0.5]
#Generate high density region #Generate high density region
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment