In [1]:
%run basics
%matplotlib
In [26]:
l3_name = qcio.get_filename_dialog()
ds_l3 = qcio.nc_read_series(l3_name)
access_name = qcio.get_filename_dialog()
ds_access = qcio.nc_read_series(access_name)
In [29]:
ldt_l3 = ds_l3.series["DateTime"]["Data"]
ldt_access = ds_access.series["DateTime"]["Data"]
start_date = max([ldt_l3[0],ldt_access[0]])
end_date = min([ldt_l3[-1],ldt_access[-1]])
print start_date,end_date
si_access = qcutils.GetDateIndex(ldt_access,str(start_date))
ei_access = qcutils.GetDateIndex(ldt_access,str(end_date))
si_l3 = qcutils.GetDateIndex(ldt_l3,str(start_date))
ei_l3 = qcutils.GetDateIndex(ldt_l3,str(end_date))
print si_access,ei_access,si_l3,ei_l3
In [30]:
Wd_16m,f,a = qcutils.GetSeriesasMA(ds_l3,"Wd_16m",si=si_l3,ei=ei_l3)
Ws_16m,f,a = qcutils.GetSeriesasMA(ds_l3,"Ws_16m",si=si_l3,ei=ei_l3)
Wd_8m,f,a = qcutils.GetSeriesasMA(ds_l3,"Wd_8m",si=si_l3,ei=ei_l3)
Wd_4m,f,a = qcutils.GetSeriesasMA(ds_l3,"Wd_4m",si=si_l3,ei=ei_l3)
Wd_access,f,a = qcutils.GetSeriesasMA(ds_access,"Wd_11",si=si_access,ei=ei_access)
Wd_diff_16m = Wd_16m-Wd_access
Wd_diff_8m = Wd_8m-Wd_access
Wd_diff_4m = Wd_4m-Wd_access
Wd_diff_16m = numpy.ma.masked_where(Ws_16m<5,Wd_diff_16m)
Wd_diff_16m = numpy.ma.masked_where(abs(Wd_diff_16m)>180,Wd_diff_16m)
Wd_diff_8m = numpy.ma.masked_where(Ws_16m<5,Wd_diff_8m)
Wd_diff_4m = numpy.ma.masked_where(Ws_16m<5,Wd_diff_4m)
In [7]:
fig=plt.figure()
plt.plot(ldt_l3,Wd_16m,'b.')
plt.plot(ldt_l3,Wd_8m,'g^')
plt.plot(ldt_l3,Wd_4m,'yo')
plt.plot(ldt_l3,Wd_access,'r+')
plt.show()
In [15]:
fig=plt.figure()
plt.plot(ldt_l3,Wd_diff_16m,'b.')
plt.plot(ldt_l3,Wd_diff_8m,'g^')
plt.plot(ldt_l3,Wd_diff_4m,'yo')
plt.show()
In [32]:
bins = range(-180,180,10)
fig = plt.figure()
n, bins, patches = plt.hist(Wd_diff_16m, bins, normed=1, facecolor='green', alpha=0.75)
plt.xlabel("Wd Difference (16m - ACCESS)")
plt.ylabel("Frequency")
plt.title("Yanco: 2015")
plt.show()
In [11]:
si1_l3 = qcutils.GetDateIndex(ldt_l3,"2012-07-02 00:30")
ei1_l3 = qcutils.GetDateIndex(ldt_l3,"2012-07-10 12:00")
si1_access = qcutils.GetDateIndex(ldt_access,"2012-07-02 00:30")
ei1_access = qcutils.GetDateIndex(ldt_access,"2012-07-10 12:00")
In [13]:
ldt_l3_1 = ldt_l3[si1_l3:ei1_l3+1]
Wd_16m_1,f,a = qcutils.GetSeriesasMA(ds_l3,"Wd_16m",si=si1_l3,ei=ei1_l3)
Wd_8m_1,f,a = qcutils.GetSeriesasMA(ds_l3,"Wd_8m",si=si1_l3,ei=ei1_l3)
Wd_access_1,f,a = qcutils.GetSeriesasMA(ds_access,"Wd_11",si=si1_access,ei=ei1_access)
In [14]:
fig=plt.figure()
plt.plot(ldt_l3_1,Wd_16m_1,'b.')
plt.plot(ldt_l3_1,Wd_8m_1,'g^')
plt.plot(ldt_l3_1,Wd_access_1,'r+')
plt.show()
In [16]:
avgdiff_16 = numpy.ma.mean(Wd_16m_1-Wd_access_1)
print avgdiff_16
In [ ]: