コード例 #1
0
ファイル: velhi.py プロジェクト: sbird/DLA_script
 def __init__(self,snap_dir,snapnum,minpart,reload_file=False,savefile=None):
     if savefile==None:
         savefile_s=path.join(snap_dir,"snapdir_"+str(snapnum).rjust(3,'0'),"velocity_grid.hdf5")
     else:
         savefile_s = savefile
     HaloHI.__init__(self,snap_dir,snapnum,minpart=minpart,reload_file=reload_file,savefile=savefile_s)
     return
コード例 #2
0
ファイル: boxhi.py プロジェクト: sbird/DLA_script
 def save_file(self, save_grid=False, LLS_cut = 17., DLA_cut = 20.3):
     """Save the file, by default without the grid"""
     HaloHI.save_file(self,save_grid)
     #Save a list of DLA positions instead
     f=h5py.File(self.savefile,'r+')
     ind = np.where(self.sub_nHI_grid > DLA_cut)
     ind_LLS = np.where((self.sub_nHI_grid > LLS_cut)*(self.sub_nHI_grid < DLA_cut))
     grp = f.create_group("abslists")
     grp.create_dataset("DLA",data=ind)
     grp.create_dataset("DLA_val",data=self.sub_nHI_grid[ind])
     grp.create_dataset("LLS",data=ind_LLS)
     grp.create_dataset("LLS_val",data=self.sub_nHI_grid[ind_LLS])
     f.close()
コード例 #3
0
ファイル: boxhi.py プロジェクト: sbird/DLA_script
 def save_file(self, save_grid=False, LLS_cut=17., DLA_cut=20.3):
     """Save the file, by default without the grid"""
     HaloHI.save_file(self, save_grid)
     #Save a list of DLA positions instead
     f = h5py.File(self.savefile, 'r+')
     ind = np.where(self.sub_nHI_grid > DLA_cut)
     ind_LLS = np.where(
         (self.sub_nHI_grid > LLS_cut) * (self.sub_nHI_grid < DLA_cut))
     grp = f.create_group("abslists")
     grp.create_dataset("DLA", data=ind)
     grp.create_dataset("DLA_val", data=self.sub_nHI_grid[ind])
     grp.create_dataset("LLS", data=ind_LLS)
     grp.create_dataset("LLS_val", data=self.sub_nHI_grid[ind_LLS])
     f.close()