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: