In [1]:

    

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


    

In [4]:

    

nc_name = qcio.get_filename_dialog(path="../../Sites")


    

In [15]:

    

ds = qcio.nc_read_series(nc_name)


    

In [6]:

    

for var in ds.series.keys():
    ds.series[var]["Attr"]["missing_value"] = -9999


    

In [16]:

    

for var in ds.series.keys():
    if "valid_range" in ds.series[var]["Attr"].keys():
        valid_range = ds.series[var]["Attr"]["valid_range"]
        if '"' in valid_range: valid_range = valid_range.replace('"','')


    

In [19]:

    

for attr in ds.series["Fsd"]["Attr"].keys():
    print attr,ds.series["Fsd"]["Attr"][attr]


    

    




instrument Kipp and Zonen CNR1(not CNR4!)
ancillary_variables Fsd QC flag
height 70 m
coverage_L2 84
long_name Down-welling short wave
standard_name surface_downwelling_shortwave_flux_in_air
serial_number 
units W/m2
coverage_L3 84
valid_range -10,1500
rangecheck_upper [1500]*12
rangecheck_lower [-10]*12

In [20]:

    

out_name = nc_name.replace(".nc","_vr.nc")


    

In [21]:

    

ncfile = qcio.nc_open_write(out_name)


    

In [22]:

    

qcio.nc_write_series(ncfile,ds)


    

In [ ]: