In [1]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

from hmf import Perturbations


    

In [3]:

    

def get_track(log_abundance=-5.0):
    M = arange(8.0,13.0,0.01)
    z_list = arange(6.0,20.4,0.5)
    n_z = size(z_list)
    mass_track = empty((0))
    delta_track = empty((0))
    for z in z_list:
       pert = Perturbations(M,z=z)
       mass_func_tmp = pert.dndlnm
       delta= abs(log10(mass_func_tmp)-log_abundance)
       i_min = argmin(delta)
       mass_track = append(mass_track, M[i_min])
       delta_track = append(delta_track, delta[i_min])
    return z_list, mass_track, delta_track


    

In [4]:

    

abundances = arange(-6.0,-2.9,0.1)
print abundances
for log_ab in abundances:
    filename= "../data/mass_track_abundance%.1f.dat"%(log_ab)
    z_track, m_track, diff_track  = get_track(log_abundance=log_ab)
    out = open(filename, "w")
    for z, m, d in zip(z_track, m_track, diff_track):
        out.write("%f %f %f\n"%(z,m,d))
    out.close()
    print filename


    

    




[-6.  -5.9 -5.8 -5.7 -5.6 -5.5 -5.4 -5.3 -5.2 -5.1 -5.  -4.9 -4.8 -4.7 -4.6
 -4.5 -4.4 -4.3 -4.2 -4.1 -4.  -3.9 -3.8 -3.7 -3.6 -3.5 -3.4 -3.3 -3.2 -3.1
 -3. ]
../data/mass_track_abundance-6.0.dat
../data/mass_track_abundance-5.9.dat
../data/mass_track_abundance-5.8.dat
../data/mass_track_abundance-5.7.dat
../data/mass_track_abundance-5.6.dat
../data/mass_track_abundance-5.5.dat
../data/mass_track_abundance-5.4.dat
../data/mass_track_abundance-5.3.dat
../data/mass_track_abundance-5.2.dat
../data/mass_track_abundance-5.1.dat
../data/mass_track_abundance-5.0.dat
../data/mass_track_abundance-4.9.dat
../data/mass_track_abundance-4.8.dat
../data/mass_track_abundance-4.7.dat
../data/mass_track_abundance-4.6.dat
../data/mass_track_abundance-4.5.dat
../data/mass_track_abundance-4.4.dat
../data/mass_track_abundance-4.3.dat
../data/mass_track_abundance-4.2.dat
../data/mass_track_abundance-4.1.dat
../data/mass_track_abundance-4.0.dat
../data/mass_track_abundance-3.9.dat
../data/mass_track_abundance-3.8.dat
../data/mass_track_abundance-3.7.dat
../data/mass_track_abundance-3.6.dat
../data/mass_track_abundance-3.5.dat
../data/mass_track_abundance-3.4.dat
../data/mass_track_abundance-3.3.dat
../data/mass_track_abundance-3.2.dat
../data/mass_track_abundance-3.1.dat
../data/mass_track_abundance-3.0.dat

In [ ]: