Ejemplo n.º 1
0
        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)
Ejemplo n.º 2
0
    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: