Example #1
0
from bgpe.util.stellarpop import n_component
from pystarlight.util.redenninglaws import calc_redlaw
from bgpe.io.hdf5util import inithdf5, read_filterhdf5
from numpy.lib.recfunctions import drop_fields, merge_arrays, rename_fields
from bgpe.io.readfilterset import readfilterset
import cosmocalc
from pystarlight.util.constants import L_sun
from bgpe.util.cosmo import zcor
from bgpe.photometry.syntphot import spec2filterset
import os

#logging
import logging
from bgpe.core.log import setConsoleLevel
setConsoleLevel(logging.CRITICAL)

def argparser():
    
    ''' Defines the input and help options of the program... '''
    
    parser = argparse.ArgumentParser(description='Convert STARLIGHT-SDSS input + output + table output to a hdf5 library database.')
    
    parser.add_argument('-B', metavar='BaseFile.txt', type=str, nargs=1,
                        help='Input base file given on STARLIGHT standards. See STARLIGHT manual for details.', required=True)    
    parser.add_argument('-Bd', metavar='/path/to/base_directory', type=str, nargs=1,
                        help='The directory where Base files are stored', required=True)
#    parser.add_argument('-si', metavar='/path/to/sdss_directory', type=str, nargs=1,
#                        help='The directory where SDSS Input .txt files are stored', required=True)
#    parser.add_argument('-st', metavar='/path/to/tables_directory', type=str, nargs=1,
#                        help='The directory where output Tables are stored', required=True)
Example #2
0
    parser.add_argument('-dz', metavar='0.01', type=float, nargs=1, help='Delta redshift', required=True)
    parser.add_argument('-v', '-v', action='count')
    parser.add_argument('--version', action='version', version='%s version %s' % (_bgpe_name_, _bgpe_version_))
    
    args = parser.parse_args()
    return args

if __name__ == '__main__':
    
    t0 = time.time() # To measure execution times.
    
    args = argparser()
    
    #logging
    from bgpe.core.log import setConsoleLevel
    setConsoleLevel(logging.DEBUG)
    
    log = logging.getLogger('bgpe.starlight2photdb')
    
    # 0 - Main definitions
    in_dir = args.si[0]
    out_dir = args.so[0]
    tables_dir = args.st[0]
    aux = np.loadtxt(args.i[0], dtype = np.str).T
    infiles = aux[0]
    outfiles = aux[1]
    el_file = '%s/sample.F%s.%s.f.lines.dat.BS.bz2' % (tables_dir, args.d[0], args.sa[0]) 
    syn01_file = '%s/sample.F%s.%s.f.Starlight.SYN01.tab.BS.bz2' % (tables_dir, args.d[0], args.sa[0]) 
    syn02_file = '%s/sample.F%s.%s.f.Starlight.SYN02.tab.BS.bz2' % (tables_dir, args.d[0], args.sa[0]) 
    syn03_file = '%s/sample.F%s.%s.f.Starlight.SYN03.tab.BS.bz2' % (tables_dir, args.d[0], args.sa[0]) 
    syn04_file = '%s/sample.F%s.%s.f.Starlight.SYN04.tab.BS.bz2' % (tables_dir, args.d[0], args.sa[0])