In [1]:

    

%run basics
%matplotlib


    

    




Using matplotlib backend: Qt4Agg

In [2]:

    

cf = qcio.load_controlfile(path='controlfiles')
infilename = qcio.get_infilenamefromcf(cf)
ds4 = qcio.nc_read_series(infilename)


    

In [3]:

    

ds5 = qcio.copy_datastructure(cf,ds4)
qcutils.UpdateGlobalAttributes(cf,ds5,"L5")
ds5.cf = cf
ds_alt = {}
qcgf.ImportSeries(cf,ds5)
qcck.do_qcchecks(cf,ds5)
for ThisOne in cf["Fluxes"].keys():
    qcgf.GapFillParseControlFile(cf,ds5,ThisOne,ds_alt)


    

In [4]:

    

Fc_data_nf,Fc_flag_nf,Fc_attr_nf = qcutils.GetSeriesasMA(ds5,"Fc")


    

In [5]:

    

qcck.ApplyTurbulenceFilter(cf,ds5)


    

    




WARNING:qc.rp: CPD results filename not in control file

In [6]:

    

Fc_data_f,Fc_flag_f,Fc_attr_f = qcutils.GetSeriesasMA(ds5,"Fc")


    

In [8]:

    

print numpy.ma.count(Fc_data_nf),numpy.ma.count(Fc_data_f)


    

    




44630 22569

In [11]:

    

print len(numpy.where(Fc_flag_nf==0)[0]),len(numpy.where(Fc_flag_f==0)[0])


    

    




41144 21272

In [ ]: