from pylab import * from matplotlib import rc rc('text', usetex=True) # use latex in figures # check i = len(sys.argv) if i != 3: print "Usage: plot_Mstar_rhalf.py snap outfile" print "Example: plot_Mstar_rhalf.py 270 mr2_00270.png" exit(1) snap = int(sys.argv[1]) outfile = sys.argv[2] # read in Mstar, rhalf mstar = array(mys.get_M_star(snap)) mdm = array(mys.get_M_dm(snap)) issub = array(mys.get_is_sub(snap)) rhalf = array(mys.get_rhalf_star(snap)) * 1000 mtot = mdm + mstar cc = mstar / mtot print cc rh1 = [] ms1 = [] rh2 = [] ms2 = [] pylab.figure()
if i != 3: print "Usage: plot_sfr_halo.py snap_start snap_stop" print "Example: plot_sfr_halo.py 200 270" exit(1) snapstart = int(sys.argv[1]) snapstop = int(sys.argv[2]) dnz = [] z = [] i = 0 for snap in range(snapstart, snapstop): if (not mys.exists_snap(snap)): continue i = i + 1 mstar = mys.get_M_star(snap)[0] zsnap = mys.get_z(snap) print mstar, zsnap if (i == 1): zold = zsnap nold = mstar continue znew = zsnap nnew = mstar z.append(znew) dnz.append((nnew - nold) / (zold - znew)) #looks clearer without log, for whole simulation zold = znew nold = nold #sort by z so lines between points in plot are useful
i=len(sys.argv) if i!=3: print "Usage: plot_sfr_halo.py snap_start snap_stop" print "Example: plot_sfr_halo.py 200 270" exit(1) snapstart = int(sys.argv[1]) snapstop = int(sys.argv[2]) dnz=[]; z=[] i=0; for snap in range(snapstart,snapstop): if(not mys.exists_snap(snap)): continue i=i+1 mstar=mys.get_M_star(snap)[0] zsnap=mys.get_z(snap) print mstar,zsnap if(i==1): zold=zsnap nold=mstar continue znew=zsnap nnew=mstar z.append(znew) dnz.append((nnew-nold)/(zold-znew)) #looks clearer without log, for whole simulation zold=znew nold=nold #sort by z so lines between points in plot are useful z=array(z); dnz=array(dnz)
from pylab import * from matplotlib import rc rc('text',usetex=True)# use latex in figures # check i=len(sys.argv) if i!=3: print "Usage: plot_Mstar_rhalf.py snap outfile" print "Example: plot_Mstar_rhalf.py 270 mr2_00270.png" exit(1) snap=int(sys.argv[1]) outfile=sys.argv[2] # read in Mstar, rhalf mstar= array(mys.get_M_star(snap)) mdm = array(mys.get_M_dm(snap)) issub= array(mys.get_is_sub(snap)) rhalf= array(mys.get_rhalf_star(snap))*1000 mtot=mdm+mstar cc=mstar/mtot print cc rh1=[];ms1=[]; rh2=[]; ms2=[] pylab.figure() for i in range(len(cc)): if(cc[i]>0.5):