示例#1
0
indata = np.loadtxt(infile,
                    dtype='float',
                    delimiter=None,
                    comments='#',
                    usecols=[0, 5])
want = (indata[:, 1] >= 3) & (indata[:, 1] <= 5)
data = indata[want, :]
#%%
span = [55, 100]
for i, v in enumerate(data[span[0]:span[1], 0]):
    epicid = np.int(v)
    acti = i + span[0]
    cfg['campaign'] = int(data[acti, 1])
    print epicid, cfg['campaign']

    output = mS.runOne(epicid, cfg)

    try:
        print output['exception']
        rep = "%u --bad -- %s\n" % (epicid, output['exception'])
        fid.write(rep)
        print ' ! !  EXCEPTION EXCEPTION ! ! ! ! '
        oneout = "/home/sthomp/DAVE/playK2/epic%u.exc" % epicid
        fida = open(oneout, 'w')
        fida.write(output['backtrace'])
        fida.close()

    except KeyError:

        try:
示例#2
0
cfg['stellarFile']='/home/sthomp/DAVE/dave/etc/k2EpicCatalogStellarTable5.txt'

#cfg['modshiftBasename']='/home/stho
#cfg['prfPath']='morejunk/junk';

indata=np.loadtxt(infile,dtype='float',delimiter=',',comments='#',usecols=[0,3])
data=indata
#%%
span=[498,500]
for i,v in enumerate(data[span[0]:span[1],0]):    
    epicid=np.int(v)
    acti=i+span[0]
    cfg['campaign'] = 6
    print epicid,cfg['campaign']

    clip=mS.runOne(epicid, cfg)
     
    try:
        print clip['exception']
        rep="%u --bad -- %s\n" % (epicid, clip['exception'])
        print ' ! !  EXCEPTION EXCEPTION ! ! ! ! '
        oneout="/home/sthomp/daveOutput/c6/epic%s.exc" % epicid
        #fida=open(oneout,'w')
        #fida.write(clip['backtrace'])
        #fida.close()
        
    except KeyError:
        #print "%u   %s" % (epicid,clip.disposition.isCandidate)
        pass
    
    
rundata=rdata[idx,:]

#%%

periods=np.linspace(0.32,2.1,num=25000)
bls=np.zeros((410,25000))
fid=open(outfile,'a')
span=[55,410]
n=0
for i,v in enumerate(rundata[span[0]:span[1],0]):    
    epicid=np.int(v)
    acti=i+span[0];
    cfg['campaign'] = int(rundata[acti,1])
    print epicid,cfg['campaign']

    clip=mS.runOne(epicid, cfg)
    
    #plt.plot(clip.bls.bls_search_periods,clip.bls.convolved_bls,'r')
    
    try:
        f=interp1d(clip.bls.bls_search_periods,clip.bls.convolved_bls)
        ibls=f(periods)
        #bls[n,:]=ibls;
        clip.bls['ibls']=ibls;
        clip.bls['iperiods']=periods
        clip.bls['kepmag']=rundata[acti,2]
        clip=pipe.saveClip(clip)
        n=n+1;
    except AttributeError:
        pass
示例#4
0
cfg['taskList'][-1] = "dpp.saveClip" 
#cfg['modshiftBasename']='/home/stho
#cfg['prfPath']='morejunk/junk';

indata=np.loadtxt(infile,dtype='float',delimiter=None,comments='#',usecols=[0,5])
want=(indata[:,1]>=3) & (indata[:,1] <=5)
data=indata[want,:]
#%%
span=[55,100]
for i,v in enumerate(data[span[0]:span[1],0]):    
    epicid=np.int(v)
    acti=i+span[0]
    cfg['campaign'] = int(data[acti,1])
    print epicid,cfg['campaign']

    output=mS.runOne(epicid, cfg)
     
    try:
        print output['exception']
        rep="%u --bad -- %s\n" % (epicid, output['exception'])
        fid.write(rep)
        print ' ! !  EXCEPTION EXCEPTION ! ! ! ! '
        oneout="/home/sthomp/DAVE/playK2/epic%u.exc" % epicid
        fida=open(oneout,'w')
        fida.write(output['backtrace'])
        fida.close()
        
    except KeyError:
        
        try: