Skip to content
Snippets Groups Projects
Commit 5c296e9a authored by Folkert Nobels's avatar Folkert Nobels
Browse files

Update starformation law plotter

parent a1f4b098
Branches
Tags
1 merge request!705Star formation following Schaye08
...@@ -4,10 +4,11 @@ import h5py ...@@ -4,10 +4,11 @@ import h5py
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
# for the moment only use a single snapshot # for the moment only use a single snapshot
i = 110 i = 29
# Load the snapshot # Load the snapshot
f = h5py.File('output_%04dcopy.hdf5'%i, 'r') #f = h5py.File('output_%04dcopy.hdf5'%i, 'r')
f = h5py.File('output_%04d.hdf5'%i, 'r')
# Load the stars # Load the stars
starparticle = f['PartType4'] starparticle = f['PartType4']
...@@ -60,9 +61,10 @@ for i in range(1,len(r_array)): ...@@ -60,9 +61,10 @@ for i in range(1,len(r_array)):
Sigmagas[i-1] = massGAS/area_loop Sigmagas[i-1] = massGAS/area_loop
# Put SF law in correct units # Put SF law in correct units
SigmaSFR *= 1e10/1.10e8 SigmaSFR *= 1e10/0.29e8
Sigmagas *= 1e10/(1e3**2) Sigmagas *= 1e10/(1e3**2)
SigmaSFR[SigmaSFR==0] = 1e-6 SigmaSFR[SigmaSFR==0] = 1e-6
print(Sigmagas)
Sgas_array =np.logspace(-1,3,40) Sgas_array =np.logspace(-1,3,40)
Delta_Sgas = np.log10(Sgas_array[1]) - np.log10(Sgas_array[0]) Delta_Sgas = np.log10(Sgas_array[1]) - np.log10(Sgas_array[0])
...@@ -90,7 +92,7 @@ print(SigmaSFR_avg) ...@@ -90,7 +92,7 @@ print(SigmaSFR_avg)
print(Sigmagas_avg) print(Sigmagas_avg)
# Put SF law in correct units # Put SF law in correct units
SigmaSFR_avg *= 1e10/1.10e8 SigmaSFR_avg *= 1e10/0.29e8
Sigmagas_avg *= 1e10/(1e6) Sigmagas_avg *= 1e10/(1e6)
SigmaSFR_avg[SigmaSFR_avg==0]=1e-6 SigmaSFR_avg[SigmaSFR_avg==0]=1e-6
...@@ -106,7 +108,7 @@ ncrit = 10 # M_sun /pc^2 ...@@ -106,7 +108,7 @@ ncrit = 10 # M_sun /pc^2
# plot SF law # plot SF law
plt.loglog(Sigmagas,SigmaSFR,'o',label='Data') plt.loglog(Sigmagas,SigmaSFR,'o',label='Data')
plt.loglog(Sigmagas_avg,SigmaSFR_avg,'or',label='binned Data') #plt.loglog(Sigmagas_avg,SigmaSFR_avg,'o',label='binned Data')
plt.loglog(Sigmagaslaw, SigmaSFRlaw,'r-.',label='KS law') plt.loglog(Sigmagaslaw, SigmaSFRlaw,'r-.',label='KS law')
#plt.loglog(np.logspace(-1,3,40),1e-5*np.ones(40),'o') #plt.loglog(np.logspace(-1,3,40),1e-5*np.ones(40),'o')
plt.axvline(x=10,linestyle='--', color='r',label='$\Sigma_c = 10$') plt.axvline(x=10,linestyle='--', color='r',label='$\Sigma_c = 10$')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment