import sys
import os
#===========================================
"""
***********************************
--> Lectura Subhalos
***********************************
"""
print("-----------------------------------\n")
print("Inicia asignacion de los subfind\n")
sub = gadget.Subfind('/home/dmontenegro/Data/Sims512/cosmobh01/groups_015/fof_subhalo_tab_015.0.hdf5' ,combineFiles=True)
Mass_halo= sub.SubhaloMass*1e10/0.7
vel_halo = sub.SubhaloVel ##km/s
vel_dis_halo = sub.SubhaloVelDisp ##km/s
vel_max_halo = sub.SubhaloVmax ##km/s
spin_halo=sub.SubhaloSpin ##ckpc/h
r_halo = sub.SubhaloPos ## posicion en x,y,z ckp/h
#print("size data subhalo",len(Mass_halo))
print("-----------------------------------\n")
print("Termina asignacion de los subfind\n")
from struct import *
import sys
import os
#===========================================
"""
***********************************
--> Lectura Subhalos
***********************************
"""
print("-----------------------------------\n")
print("Inicia asignacion de los subfind\n")
sub = gadget.Subfind('../../Data/groups_015/fof_subhalo_tab_015.0.hdf5',
combineFiles=True)
Mass_halo = sub.SubhaloMass * 10e10 / 0.7
vel_halo = sub.SubhaloVel ##km/s
vel_dis_halo = sub.SubhaloVelDisp ##km/s
vel_max_halo = sub.SubhaloVmax ##km/s
spin_halo = sub.SubhaloSpin ##ckpc/h
r_halo = sub.SubhaloPos ## posicion en x,y,z ckp/h
print("-----------------------------------\n")
print("Termina asignacion de los subfind\n")
print(":::::::::::::::::::::::::::::::::::::\n")
print("Inicia graficas de control halos\n")
"""
********************************