return surface def make_arrays(resolution='coarse', surface=False, case=None, period=False, homedir=None): wpfiles = get_filenames(case, homedir=homedir + '/WINDPROF') wp = [] ncols = 0 # number of timestamps for f in wpfiles: if resolution == 'fine': wp.append(mf.parse_windprof(f, 'fine')) elif resolution == 'coarse': wp.append(mf.parse_windprof(f, 'coarse')) else: print 'Error: resolution has to be "fine" or "coarse"' ncols += 1 ''' creates 2D arrays with spd and dir ''' nrows = len( wp[0].HT.values) # number of altitude gates (fine same as coarse) hgt = wp[0].HT.values # print len(hgt) wspd = np.empty([nrows, ncols]) wdir = np.empty([nrows, ncols]) timestamp = [] for i, p in enumerate(wp): timestamp.append(p.timestamp)