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)
