In [1]:

    

%run basics
%matplotlib


    

    




Using matplotlib backend: Qt4Agg

In [2]:

    

nc_name = qcio.get_filename_dialog(path='../Sites/')
ds = qcio.nc_read_series(nc_name)


    

In [3]:

    

dt = ds.series["DateTime"]["Data"]
ustar,f,a = qcutils.GetSeriesasMA(ds,"ustar")
ustar_CR3,f,a = qcutils.GetSeriesasMA(ds,"ustar_CR3")
Fh,f,a = qcutils.GetSeriesasMA(ds,"Fh")
Fh_CR3,f,a = qcutils.GetSeriesasMA(ds,"Fh_CR3")
Fe,f,a = qcutils.GetSeriesasMA(ds,"Fe")
Fe_CR3,f,a = qcutils.GetSeriesasMA(ds,"Fe_CR3")
Fc,f,a = qcutils.GetSeriesasMA(ds,"Fc")
Fc_CR3,f,a = qcutils.GetSeriesasMA(ds,"Fc_CR3")


    

In [4]:

    

xlabel = 'Fh_CR3 (W/m2)'; ylabel = 'Fh (W/m2)'
qcplot.xyplot(Fh_CR3,Fh,sub=[2,2,1],regr=1,xlabel=xlabel,ylabel=ylabel)
xlabel = 'Fe_CR3 (W/m2)'; ylabel = 'Fe (W/m2)'
qcplot.xyplot(Fe_CR3,Fe,sub=[2,2,2],regr=1,xlabel=xlabel,ylabel=ylabel)
xlabel = 'Fc_CR3 (g/m2/s)'; ylabel = 'Fc (g/m2/s)'
qcplot.xyplot(Fc_CR3,Fc,sub=[2,2,3],regr=1,xlabel=xlabel,ylabel=ylabel)
xlabel = 'u*_CR3 (m/s)'; ylabel = 'u* (m/s)'
qcplot.xyplot(ustar_CR3,ustar,sub=[2,2,4],regr=1,xlabel=xlabel,ylabel=ylabel)


    

In [4]:

    

fig=plt.figure()
plt.plot(dt,Fh,'b.')
plt.plot(dt,Fh_CR3,'r+')
plt.show()


    

In [ ]: