In [2]:
import sys
sys.path.append('../scripts')
import constants as c
import glob
import qcio
import qcutils
import netCDF4
import matplotlib.pyplot as plt
In [3]:
bios_name=qcio.get_filename_dialog()
In [47]:
bios_ncfile=netCDF4.Dataset(bios_name)
time=bios_ncfile.variables["time"][:]
time_units=getattr(bios_ncfile.variables["time"],"units")
dt_bios=netCDF4.num2date(time,time_units)
Sws_bios=bios_ncfile.variables["SoilMoist"][:,0,0]
Precip_bios=1800*bios_ncfile.variables["Rainf"][:,0]
Fsd_bios=bios_ncfile.variables["SWdown"][:,0]
In [35]:
si_bios=qcutils.GetDateIndex(list(dt),"2012-01-01 00:00")
ei_bios=qcutils.GetDateIndex(list(dt),"2013-01-01 00:00")
In [36]:
access_name=qcio.get_filename_dialog()
In [49]:
ds_access=qcio.nc_read_series(access_name)
dt_access=ds_access.series["DateTime"]["Data"]
Sws_access,f,a=qcutils.GetSeriesasMA(ds_access,"Sws_11")
Precip_access,f,a=qcutils.GetSeriesasMA(ds_access,"Precip_11")
Fsd_access,f,a=qcutils.GetSeriesasMA(ds_access,"Fsd_11")
In [38]:
tower_name=qcio.get_filename_dialog()
In [54]:
ds_tower=qcio.nc_read_series(tower_name)
dt_tower=ds_tower.series["DateTime"]["Data"]
Sws_tower,f,a=qcutils.GetSeriesasMA(ds_tower,"Sws")
Precip_tower,f,a=qcutils.GetSeriesasMA(ds_tower,"Precip")
Fsd_tower,f,a=qcutils.GetSeriesasMA(ds_tower,"Fsd")
In [51]:
fig=plt.figure()
plt.plot(list(dt_bios[si_bios:ei_bios]),Sws_bios[si_bios:ei_bios],'b.')
plt.plot(dt_access,Sws_access,'r+')
plt.plot(dt_tower,Sws_tower,'go')
Out[51]:
In [52]:
fig=plt.figure()
plt.plot(list(dt_bios[si_bios:ei_bios]),Precip_bios[si_bios:ei_bios],'b.')
plt.plot(dt_access,Precip_access,'r+')
plt.plot(dt_tower,Precip_tower,'go')
Out[52]:
In [55]:
fig=plt.figure()
plt.plot(list(dt_bios[si_bios:ei_bios]),Fsd_bios[si_bios:ei_bios],'b.')
plt.plot(dt_access,Fsd_access,'r+')
plt.plot(dt_tower,Fsd_tower,'go')
Out[55]:
In [ ]: