header, data_dict = ltools.from_HDF(filepath[0], datatypes) try: BRdat = data_dict[datatypes[0]] except KeyError: sys.exit("No BR532 data available!") try: PRdat = data_dict[datatypes[1]] except KeyError: sys.exit("No PR532 data available!") altrange = np.arange(150,10000,10) BRdat = ltools.alt_resample(BRdat,altrange) PRdat = ltools.alt_resample(PRdat,altrange) altrange = BRdat.columns.values copolvals = np.hstack(BRdat.values).astype('float32') depolvals = np.hstack(PRdat.values).astype('float32') numbins = 100 depolmin = 0.0 depolmax = 0.5 copolmin = 0.0 copolmax = 5.0 copolhist=h2d.fullhist(copolvals,numbins,copolmin,copolmax,-9999.,-8888.)
#double check to make sure the mask files match up with the data files for d,m in zip(datafiles, maskfiles): [d_stat,d_date,d_type] = d.split('_') [m_stat,m_date,m_type] = m.split('_') print 'Checking mask/data match for %s'%(d_date) if d_date == m_date and d_stat == m_stat: print 'Check!' continue else: sys.exit("Error: Mask files don't match data files") #open, altitude resample, and concatenate data and mask files for d,m in zip(datafiles, maskfiles): d_temp, data_prod = LNC.lnc_reader(d) d_realt = LNC.alt_resample(d_temp,altrange) try: d_event = pan.concat([d_event,d_realt]) except NameError: d_event = d_realt m_temp, data_prod = LNC.lnc_reader(m) m_realt = LNC.alt_resample(m_temp,altrange) try: m_event = pan.concat([m_event,m_realt]) except NameError: m_event = m_realt
BRfiles.append(f) #make sure the files are in a common order of ascending date (assuming they're all #from the same station BRfiles.sort() counter = 0 header = {} #open, altitude resample, and concatenate data and mask files dataprods = [] counter = 0 for br in BRfiles: br_temp, tempprod = LNC.lnc_reader(br) br_realt = LNC.alt_resample(br_temp,altrange) if counter == 0: br_event = br_realt dataprods.append(BR_type) else: br_event = pan.concat([br_event,br_realt]) counter +=1 altrange = br_realt.columns.values #sort by index to make certain data is in order then set date ranges to match br_event = br_event.sort_index()