format='csv',
devices=[device],
measure="hourly",
order="time",
end=end,
start=sta,
limit=10000)
df = StringIO(df)
df = pd.read_csv(df, sep=',', error_bad_lines=False)
# GPS doesn't read initially
df = df[3:]
df['time'] = pd.to_datetime(df['time'])
df.index = df['time']
df['doy'] = doy_(df.index)
td = pd.to_timedelta(df.long.iloc[0] / 360. * 24., unit='h').to_pytimedelta()
df['solartime'] = df.time + td
df['solartime'] = pd.to_datetime(df['solartime'])
df.index = df['solartime']
# Some useful variables
SWP = SWpot_(df.time.dt, df.lat, df.long)
df['k'] = df.SWdw / SWP
psi = solar_psi_(df.time.dt, df.lat, df.long) # Solar Zenith Angle
