In [2]:

    

%run basics


    

In [3]:

    

tower_name=qcio.get_filename_dialog(path="../../Sites/")
aws_name=qcio.get_filename_dialog(path="../../Sites/")


    

In [4]:

    

ds_tower=qcio.nc_read_series(tower_name)


    

In [5]:

    

ds_aws=qcio.nc_read_series(aws_name)


    

In [6]:

    

dt_tower=ds_tower.series["DateTime"]["Data"]
dt_aws=ds_aws.series["DateTime"]["Data"]


    

In [18]:

    

startdate = "2011-01-01 00:30"
enddate = "2012-01-01 00:00"
si_tower=qcutils.GetDateIndex(dt_tower,startdate)
ei_tower=qcutils.GetDateIndex(dt_tower,enddate)
si_aws=qcutils.GetDateIndex(dt_aws,startdate)
ei_aws=qcutils.GetDateIndex(dt_aws,enddate)
ldt_towere=dt_tower[si_tower:ei_tower+1]
ldt_aws=dt_aws[si_aws:ei_aws+1]


    

In [19]:

    

Ts_tower,f,a=qcutils.GetSeriesasMA(ds_tower,"Ts")
Flu_tower,f,a=qcutils.GetSeriesasMA(ds_tower,"Flu")
Ts_bios,f,a=qcutils.GetSeriesasMA(ds_bios,"Ts",si=si_bios,ei=ei_bios)
Flu_bios,f,a=qcutils.GetSeriesasMA(ds_bios,"Flu",si=si_bios,ei=ei_bios)


    

In [22]:

    

fig=plt.figure(1)
ax1=plt.subplot(411)
ax1.plot(dt_tower,Ts_tower,'b.')
plt.ylabel("Ts_tower")
ax2=plt.subplot(412)
ax2.plot(dt_bios,Ts_bios,'b.')
plt.ylabel("Ts_BIOS2")
ax3=plt.subplot(413)
ax3.plot(dt_tower,Flu_tower,'r+')
plt.ylabel("Flu_tower")
ax4=plt.subplot(414)
ax4.plot(dt_bios,Flu_bios,'r+')
plt.ylabel("Flu_BIOS2")
plt.show()


    

In [ ]: