

In [1]:

    
%run basics
%matplotlib
import qcrp
import os









    



Using matplotlib backend: Qt4Agg



In [43]:

    
nc_name = qcio.get_filename_dialog(path="../../Sites/")
path_name = os.path.split(nc_name)
data_path = path_name[0]
plot_path = data_path.replace("Data","Plots")
if not os.path.exists(plot_path): os.makedirs(plot_path)
ds = qcio.nc_read_series(nc_name)
qcrp.CalculateET(ds)
series_dict = {"daily":{},"annual":{},"cumulative":{},"lists":{}}
sdl = series_dict["lists"]
sdl["nee"] = ["NEE"]
sdl["gpp"] = ["GPP"]
sdl["fre"] = ["Fre"]
sdl["co2"] = sdl["nee"]+sdl["gpp"]+sdl["fre"]
for item in sdl["co2"]:
    series_dict["daily"][item] = {}
    series_dict["cumulative"][item] = {}
    series_dict["daily"][item]["operator"] = "sum"
    series_dict["daily"][item]["format"] = "0.00"
    series_dict["cumulative"][item]["operator"] = "sum"
    series_dict["cumulative"][item]["format"] = "0.00"
series_dict["daily"]["Fe"] = {"operator":"average","format":"0.0"}
series_dict["daily"]["Fh"] = {"operator":"average","format":"0.0"}
series_dict["daily"]["Fg"] = {"operator":"average","format":"0.0"}
series_dict["daily"]["Fn"] = {"operator":"average","format":"0.0"}
series_dict["daily"]["ET"] = {"operator":"sum","format":"0.0"}
series_dict["daily"]["Precip"] = {"operator":"sum","format":"0.0"}
series_dict["cumulative"]["ET"] = series_dict["daily"]["ET"]
series_dict["cumulative"]["Precip"] = series_dict["daily"]["Precip"]
series_dict["annual"] = series_dict["daily"]
series_dict["monthly"] = series_dict["daily"]
xl_name = nc_name.replace(".nc","_Summary.xls")
xl_file = qcio.xl_open_write(xl_name)
daily_dict = qcrp.L6_summary_daily(xl_file,ds,series_dict)
monthly_dict = qcrp.L6_summary_monthly(xl_file,ds,series_dict)
annual_dict = qcrp.L6_summary_annual(xl_file,ds,series_dict)
cumulative_dict = qcrp.L6_summary_cumulative(xl_file,ds,series_dict)
xl_file.save(xl_name)
cf = {"Options":{"call_mode":"interactive"},"Files":{"plot_path":plot_path}}
qcrp.L6_summary_plotdaily(cf,ds,daily_dict)
qcrp.L6_summary_plotcumulative(cf,ds,cumulative_dict)