In [1]:
%run basics
In [2]:
l3name=qcio.get_filename_dialog(path="../../Sites")
In [3]:
l4name=qcio.get_filename_dialog(path="../../Sites")
In [4]:
ds3=qcio.nc_read_series(l3name)
ds4=qcio.nc_read_series(l4name)
In [8]:
dt3=ds3.series["DateTime"]["Data"]
dt4=ds4.series["DateTime"]["Data"]
In [5]:
Fsd_l3_data,Fsd_l3_flag,Fsd_l3_attr=qcutils.GetSeriesasMA(ds3,"Fsd")
Fsd_l4_data,Fsd_l4_flag,Fsd_l4_attr=qcutils.GetSeriesasMA(ds4,"Fsd")
In [6]:
data_l3=ds3.series["Fsd"]['Data'].copy()
flag_l3=ds3.series["Fsd"]['Flag'].copy()
data_l4=ds4.series["Fsd_access"]['Data'].copy()
flag_l4=ds4.series["Fsd_access"]['Flag'].copy()
In [9]:
si=qcutils.GetDateIndex(dt3,"2013-02-06 12:00")
ei=qcutils.GetDateIndex(dt3,"2013-02-07 00:00")
In [12]:
index=numpy.where(numpy.mod(flag_l3,10)==0)[0]
flag_l3[index]=0
In [17]:
index=numpy.where(flag_l3!=0)[0]
In [20]:
flag_l3[index]=flag_l4[index]
In [21]:
print flag_l3[si:ei],flag_l4[si:ei]
In [ ]: