In [1]:

    

%run basics
%matplotlib


    

    




Using matplotlib backend: Qt4Agg

In [8]:

    

l5name = qcio.get_filename_dialog(path="../Sites")
#l3name = qcio.get_filename_dialog(path="../Sites")


    

In [3]:

    

# get the file name
l5name = qcio.get_filename_dialog(path="../Sites")
# read the netCDF file
ds_l5 = qcio.nc_read_series(l5name)
# convert global attributes from Unicode to ASCII
for item in ds_l5.globalattributes.keys():
    #print item,type(ds_l5.globalattributes[item])
    attr = str(ds_l5.globalattributes[item])
    ds_l5.globalattributes.pop(item,None)
    item = item.encode('ascii','ignore')
    ds_l5.globalattributes[item] = attr.encode('ascii','ignore')
# convert variable attributes
for item in ds_l5.series.keys():
    attr_dict = ds_l5.series[item]["Attr"]
    for attr in attr_dict:
        value = str(ds_l5.series[item]["Attr"][attr])
        ds_l5.series[item]["Attr"].pop(attr,None)
        attr = attr.encode('ascii','ignore')
        ds_l5.series[item]["Attr"][attr] = value.encode('ascii','ignore')
# convert variable names from Unicode to ASCII
for item in ds_l5.series.keys():
    tmp = ds_l5.series[item]
    ds_l5.series.pop(item,None)
    ds_l5.series[item.encode('ascii','ignore')] = tmp
outname = l5name.replace(".nc","_ascii.nc")
ncfile = qcio.nc_open_write(outname)
qcio.nc_write_series(ncfile,ds_l5)


    

In [3]:

    

for nl in [0,10,20,30,40,50,60,70]:
#for nl in [0,1,2,3,4,5,6,7,8,9]:
    n = 0
    ds = qcio.nc_read_series(ncname)
    for item in ds.globalattributes.keys():
        n = n + 1
        if n>nl: ds.globalattributes.pop(item,None)
    outname = ncname.replace(".nc","_"+str(nl)+".nc")
    ncfile = qcio.nc_open_write(outname)
    qcio.nc_write_series(ncfile,ds)


    

In [12]:

    

l5name = qcio.get_filename_dialog(path="../Sites")
l3name = qcio.get_filename_dialog(path="../Sites")
ds_l5 = qcio.nc_read_series(l5name)
ds_l3 = qcio.nc_read_series(l3name)
for item in ds_l5.globalattributes.keys():
    ds_l5.globalattributes.pop(item,None)
for item in ds_l3.globalattributes.keys():
    ds_l5.globalattributes[item] = ds_l3.globalattributes[item]
outname = l5name.replace(".nc","_c.nc")
ncfile = qcio.nc_open_write(outname)
qcio.nc_write_series(ncfile,ds_l5)


    

In [ ]: