def onShrink(self, event):
nstart,nstop=mys.get_range()
# have xs,xs_star set in the end
for nc in range(nstop-nstart+1):
nsnap = nc + nstart
if(not mys.snap_exists(nsnap)):
continue
print(mys.get_nhalo(nsnap))
for i in range(3):#range(mys.get_nhalo(nsnap)):
if(not mys.exists_snap(nsnap)):
print("snapshot "+str(nsnap)+" missing")
continue
cmd = "shrink_sphere.py "+str(nsnap)+" "+str(i+1)
my.run(cmd)
if(not mys.is_dmonly()):
cmd = "shrink_sphere_stars.py "+str(nsnap)+" "+str(i+1)
my.run(cmd)
# TODO: if empty xs (shrinking sphere algorithm did not succeed: use xc)
my.done()
def onShrink(self, event):
nstart, nstop = mys.get_range()
# have xs,xs_star set in the end
for nc in range(nstop - nstart + 1):
nsnap = nc + nstart
if (not mys.snap_exists(nsnap)):
continue
print(mys.get_nhalo(nsnap))
for i in range(3): #range(mys.get_nhalo(nsnap)):
if (not mys.exists_snap(nsnap)):
print("snapshot " + str(nsnap) + " missing")
continue
cmd = "shrink_sphere.py " + str(nsnap) + " " + str(i + 1)
my.run(cmd)
if (not mys.is_dmonly()):
cmd = "shrink_sphere_stars.py " + str(nsnap) + " " + str(
i + 1)
my.run(cmd)
# TODO: if empty xs (shrinking sphere algorithm did not succeed: use xc)
my.done()
from numpy import array, inner, sqrt
import lib.mysql as mys
import lib.initialize as my
i = len(sys.argv)
if (i != 3):
print("shrink_sphere.py snap hid")
exit(0)
snap = int(sys.argv[1])
hid = int(sys.argv[2])
eps = 1e-5
frad = 0.90 # shrinked sphere has radius frad*maxr, i.e. frad = 0.90 means 10% smaller
if (not mys.exists_snap(snap)):
print("snapshot " + str(snap) + " missing")
exit(0)
xc, yc, zc, mvir, rvir = mys.getxyzmr(snap, 1)
halo = my.open_file(mys.d(snap) + "dm/dm_" + str(hid) + ".dat", "r")
x = []
y = []
z = []
m = []
for line in halo:
val = line.split()
m.append(float(val[0]))
x.append(float(val[1]))
y.append(float(val[2]))
z.append(float(val[3]))
import lib.mysql as mys
import lib.initialize as my
i=len(sys.argv)
if(i!=3):
print("shrink_sphere.py snap hid")
exit(0)
snap=int(sys.argv[1])
hid=int(sys.argv[2])
eps=1e-5
frad = 0.9 # shrinked sphere has radius frad*maxr, i.e. frad = 0.90 means 10% smaller
#print("missing file")
if(not mys.exists_snap(snap)):
print("snapshot "+str(snap)+" missing")
exit(0)
xc,yc,zc,mvir,rvir=mys.getxyzmr(snap,1)
halo = my.open_file(mys.d(snap)+"stars/stars_"+str(hid)+".dat","r")
x=[];y=[];z=[];m=[]
for line in halo:
val=line.split()
m.append(float(val[0]))
x.append(float(val[1]))
y.append(float(val[2]))
z.append(float(val[3]))
x = array(x); y = array(y); z = array(z); m = array(m)
#print(x)
