In [1]:

    

import os
import sys
# check the scripts directory is present
if not os.path.exists("../scripts/"):
    print "compare_access: the scripts directory is missing"
    sys.exit()
# since the scripts directory is there, try importing the modules
sys.path.append('../scripts')
import qcio
import qcutils
import qcplot


    

In [9]:

    

drname=qcio.get_filename_dialog(path="../../Sites")
ds_dr=qcio.nc_read_series(drname)


    

    




No handlers could be found for logger "qc.utils"

In [4]:

    

dsname=qcio.get_filename_dialog(path="../../Sites")
ds_ds=qcio.nc_read_series(dsname)


    

In [6]:

    

ldt_dr=ds_dr.series["DateTime"]["Data"]
ldt_ds=ds_ds.series["DateTime"]["Data"]
print ldt_dr[0],ldt_ds[0]
print ldt_dr[-1],ldt_ds[-1]


    

    




2009-01-01 00:00:00 2009-01-01 00:00:00
2009-12-31 23:30:00 2009-12-31 23:30:00

In [7]:

    

Fsd_dr,f,a=qcutils.GetSeriesasMA(ds_dr,"Fsd")
Fsd_ds,f,a=qcutils.GetSeriesasMA(ds_ds,"Fsd")


    

In [8]:

    

fig=figure()
qcplot.xyplot(Fsd_dr,Fsd_ds,regr=1,
              title='2009',
              xlabel='Fsd Regrowth (W/m2)',
              ylabel='Fsd Savanna (W/m2)')


    

In [ ]: