Python get_M_star Examples

Example #1

File: plot_Mstar_rhalf.py Project: hasanmoudud/darcoda

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()

Example #2

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

Example #3

File: plot_sfr_halo.py Project: PascalSteger/darcoda

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)

Example #4

File: plot_Mstar_rhalf.py Project: PascalSteger/darcoda

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):