Paste #326

Welcome To LodgeIt

Welcome to the LodgeIt pastebin. In order to use the notification feature a 31 day cookie with an unique ID was created for you. The lodgeit database does not store any information about you, it's just used for an advanced pastebin experience :-). Read more on the about lodgeit page. Have fun :-)

hide this notification

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
import caesar
import numpy as np

caesarfile = '/orange/narayanan/desika.narayanan/gizmo_runs/simba/m100n1024/Groups/m100n1024_078.hdf5'

obj = caesar.load(caesarfile)

ngal = obj.ngalaxies
sfr = [obj.galaxies[i].sfr for i in range(ngal)]
mstar = [obj.galaxies[i].masses['stellar'] for i in range(ngal)]
mmol = [obj.galaxies[i].masses['H2'] for i in range(ngal)]
metallicity = [obj.galaxies[i].metallicity for i in range(ngal)]

np.savez('popping_shivaei.npz',ngal=ngal,sfr=sfr,mstar=mstar,mmol=mmol,metallicity=metallicity)