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()
# choose: get_particles for all particles # get_particles_dm for dm particles only gps = "get_particles" i = len(sys.argv) if(i!=2): print("usage: gen_particles_list.py snap") exit(1) facr = 1 # * rvir from AHF for maximal radial distance to include particles snap = int(sys.argv[1]) d = mys.d(snap) pdm = d+"particles_dm" xl,yl,zl,ml,rl = mys.getxyzmr(snap,2) nhalo= mys.get_nhalo(snap) os.system("rm "+d+"particles_dm") os.system("echo '"+str(nhalo)+"' >> "+pdm) for i in range(len(xl)): xc=str(xl[i]); yc=str(yl[i]); zc=str(zl[i]); r=str(facr*rl[i]) tmpdat = d+"tmp_"+str(i+1)+".dat" # all particles (DM/stars mixed) cmd1 = gps+" -inp "+d cmd1 += " -xc " + xc + " -yc " + yc + " -zc " + zc + " -r " + r cmd1 += ">" + tmpdat cmd2 = "true" my.threadif(cmd1,cmd2,run,run,show,run) if(not loop): break