Ejemplos de extract_galaxy en Python

Ejemplo n.º 1

Archivo: find_disks.py Proyecto: jgsuresh/Arepo-Disks

def calc_profile(base_id, base, snap_num, sub_id , interact=0):
	import os
	from read_halo import read_halo
	from center_halo import center_halo
	#from J_align_halo import J_align_halo
	from extract_galaxy import extract_galaxy
	from rotate_halo import rotate_halo
	#from visualize_halo import visualize_halo
	#from J_profile_halo import J_profile_halo
	from J_profile_halo2 import J_profile_halo
	from gal_halo_comparison import gal_halo_comparison
	from bdd_halo import bdd_halo
	from save_halo import save_halo


	save_dir = base_id+"/snap"+str(snap_num).zfill(3)+"/"

	halo = Halo()
	halo.base_id = base_id
	halo.base = base
	halo.snap_num = snap_num
	halo.sub_id = sub_id
	halo.save_dir = save_dir
	halo.interact = interact

	# Minimum threshold for # of particles
	type_thresh = 25
	halo.type_thresh = type_thresh

	# Error flag
	halo.err_flag = 0

	if False:
		read_halo(halo)
		center_halo(halo)
		extract_galaxy(halo)
		rotate_halo(halo)
		gal_halo_comparison(halo)
		#visualize_halo(halo)
		J_profile_halo(halo)
		bdd_halo(halo)
		save_halo(halo)

	if True:
		read_halo(halo)
		center_halo(halo)
		extract_galaxy(halo)
		rotate_halo(halo)
		J_profile_halo(halo)
	
	
	#print "Did rotation fail? ",halo.rotate_flag

	print "J_profile test done"
	return halo

Ejemplo n.º 2

Archivo: J_profile_test.py Proyecto: jgsuresh/Arepo-Disks

def calc_profile(base_id, base, snap_num, sub_id , interact=0):
	import os
	from read_halo import read_halo
	from center_halo import center_halo
	#from J_align_halo import J_align_halo
	from extract_galaxy import extract_galaxy
	#from rotate_halo import rotate_halo
	#from visualize_halo import visualize_halo
	from J_profile_halo2 import J_profile_halo


	save_dir = base_id+"/"

	halo = Halo()
	halo.base_id = base_id
	halo.base = base
	halo.snap_num = snap_num
	halo.sub_id = sub_id
	halo.save_dir = save_dir
	halo.interact = interact

	# Minimum threshold for # of particles
	type_thresh = 25
	halo.type_thresh = type_thresh

	# Error flag
	halo.err_flag = 0


	read_halo(halo)
	center_halo(halo)
	extract_galaxy(halo)
	rotate_halo(halo)
	J_profile_halo2(halo)
	
	
	#print "Did rotation fail? ",halo.rotate_flag

	print "J_profile test done"
	return halo