gals_txt = (args.slice_gals).split('/')[-1]
fnamesuffix = '_%s%s' % (gals_txt, fnamesuffix)
#integrated
l_gals, _ = C.sort_gals(args.slice_gals)
l_gals = sorted(l_gals)
gals_slice__integr = np.zeros(H.califaIDs_all.shape, dtype = np.bool)
for g in l_gals:
i = H.califaIDs_all.tolist().index(g)
gals_slice__integr[i] = True
mask__Tg = np.zeros((N_T, H.N_zones__g.sum()), dtype = bool)
for iT in xrange(N_T):
mask__Tg[iT] = mask_zones_iT(iT, H, args.bamin, gals_slice__g)
for iGal, califaID in enumerate(H.califaIDs):
K = C.read_one_cube(califaID, EL = True, GP = True, debug = args.debug, v_run = args.v_run)
tipos, tipo, tipo_m, tipo_p = C.get_morfologia(califaID)
my_type = C.my_morf(tipos)
N_zone = K.N_zone
# Setup elliptical-rings geometry
pa, ba = K.getEllipseParams()
K.setGeometry(pa, ba)
minzones = 5
for iT, tSF in enumerate(tSF__T):
mask__z, where_slice = H.get_prop_gal(mask__Tg[iT], califaID, return_slice = True)
N_mask = mask__z.astype(int).sum()
C.debug_var(args.debug, N_mask = N_mask, N_zone = K.N_zone)
if (N_mask < (K.N_zone - minzones)):
mask__yx = K.zoneToYX(mask__z, extensive = False, surface_density = False)
t_init_gal = time.clock()
args = parser_args()
debug = args.debug
iT = args.itSF
H = C.H5SFRData(args.hdf5)
paths = C.CALIFAPaths()
C.debug_var(debug, args = args)
if (len(np.where(H.califaIDs == args.califaID)[0]) == 0):
exit('<<< plot: %s: no data in HDF5 file.' % args.califaID)
K = C.read_one_cube(args.califaID, EL = True, GP = True)
age = H.tSF__T[iT]
ageMyr = age / 1e6
# Setup elliptical-rings geometry
pa, ba = K.getEllipseParams()
K.setGeometry(pa, ba)
tipos, tipo, tipo_m, tipo_p = C.get_morfologia(args.califaID)
galaxyImgFile = paths.get_image_file(args.califaID)
if args.output is None:
output_filename = '%s_dossier_v2.pdf' % args.califaID
else:
