import Meteoframes as mf
import matplotlib.pyplot as plt
from matplotlib import dates
import datetime
import seaborn as sns
import pandas as pd
f='/home/rvalenzuela/GPS_IWV/case03/bby_gpsiwv_010123-24.dat'
df = mf.parse_gps_iwv(f)
st = datetime.datetime(2001, 1, 23, 0, 0)
en = datetime.datetime(2001, 1, 25, 0, 0)
xticks = pd.date_range(start=st,end=en, freq='3H')
fig,ax=plt.subplots()
ax.plot(df.index, df['IPW'],'-o')
ax.set_xticks(xticks)
ax.invert_xaxis()
datefmt = dates.DateFormatter('%H\n%d')
ax.xaxis.set_major_formatter(datefmt)
ax.set_xlabel(r'$\Leftarrow$'+' Time [UTC]')
ax.set_ylabel('IWV [cm]')
t1='GPS-IWV at Bodega Bay'
t2='\nDate: ' + st.strftime('%b-%Y')
plt.title(t1+t2)
plt.show(block=False)
