In [1]:

    

%run basics
%matplotlib


    

    




Using matplotlib backend: Qt4Agg

In [2]:

    

name_dingo = "../Sites/Whroo/Data/DINGO/V12a/Whroo_L6.nc"
ds_dingo = qcio.nc_read_series(name_dingo)
ldt_dingo = ds_dingo.series["DateTime"]["Data"]
fsd_dingo,f,a = qcutils.GetSeriesasMA(ds_dingo,"Fsd")
gpp_dingo,f,a = qcutils.GetSeriesasMA(ds_dingo,"GPP")


    

In [3]:

    

name_ofqc = "../Sites/Whroo/Data/Portal/Whroo_L6.nc"
ds_ofqc = qcio.nc_read_series(name_ofqc)
ldt_ofqc = numpy.array(ds_ofqc.series["DateTime"]["Data"])
fsd_ofqc,f,a = qcutils.GetSeriesasMA(ds_ofqc,"Fsd")
fsdsyn_ofqc,f,a = qcutils.GetSeriesasMA(ds_ofqc,"Fsd_syn")
sa,f,a = qcutils.GetSeriesasMA(ds_ofqc,"solar_altitude")
gpp_ofqc,f,a = qcutils.GetSeriesasMA(ds_ofqc,"GPP_LT")
nee_ofqc,f,a = qcutils.GetSeriesasMA(ds_ofqc,"NEE_LT")
er_ofqc,f,a = qcutils.GetSeriesasMA(ds_ofqc,"ER_LT")


    

In [4]:

    

name_rep = "../Sites/Whroo/Data/REddyProc/Whroo_L6.nc"
ds_rep = qcio.nc_read_series(name_rep)
ldt_rep = ds_rep.series["DateTime"]["Data"]
fsd_rep,f,a = qcutils.GetSeriesasMA(ds_rep,"Fsd")
gpp_rep,f,a = qcutils.GetSeriesasMA(ds_rep,"GPP")


    

In [5]:

    

fig=plt.figure()
plt.plot(fsd_ofqc,gpp_ofqc,'b.',label="OzFluxQC")
plt.plot(fsd_rep,gpp_rep,'r+',label="REddyProc")
plt.xlim([-20,50])
plt.ylim([-20,20])
plt.xlabel("Fsd (W/m2)")
plt.ylabel("GPP (umol/m2/s)")
plt.title("Whroo: 2011 to 2015")
plt.legend()
plt.show()


    

In [48]:

    

fig=plt.figure()
plt.plot(ldt_ofqc,fsd_ofqc,'b.')
plt.plot(ldt_ofqc,fsdsyn_ofqc,'r+')
plt.plot(ldt_ofqc,sa,'go')
plt.show()


    

In [68]:

    

fig=plt.figure()
plt.plot(fsd_rep,gpp_rep,'b.')
plt.xlim([-20,50])
plt.ylim([-20,20])
plt.xlabel("Fsd (W/m2)")
plt.ylabel("GPP (umol/m2/s)")
plt.title("Whroo: REddyProc")
plt.show()


    

In [29]:

    

fig=plt.figure()
plt.plot(fsd_dingo,gpp_dingo,'b.')
plt.show()


    

In [56]:

    

idx = numpy.ma.where((gpp_ofqc==0)&(fsd_ofqc>10))[0]
print len(idx)


    

    




575

In [57]:

    

for i in range(100):
    print ldt_ofqc[idx[i]],fsd_ofqc[idx[i]],fsdsyn_ofqc[idx[i]],gpp_ofqc[idx[i]],nee_ofqc[idx[i]],er_ofqc[idx[i]]


    

    




2011-12-02 05:30:00 33.46975 0.0 0.0 1.06671805683 1.06671805683
2011-12-02 06:00:00 115.734 0.000190508251318 0.0 -0.292528000038 -0.292528000038
2011-12-03 06:00:00 110.6781 0.000402314924586 0.0 0.326941931202 0.326941931202
2011-12-04 05:30:00 32.11527 0.0 0.0 1.3082102054 1.3082102054
2011-12-04 06:00:00 116.8991 0.00069930294491 0.0 -1.915910045 -1.915910045
2011-12-09 05:30:00 31.66507 0.0 0.0 -1.30396509679 -1.30396509679
2011-12-09 06:00:00 107.4047 0.00137346329411 0.0 -5.30347370867 -5.30347370867
2011-12-12 05:30:00 33.86893 0.0 0.0 -0.201591104165 -0.201591104165
2011-12-12 06:00:00 116.2462 0.000441409012986 0.0 -1.80538903502 -1.80538903502
2011-12-13 05:30:00 14.00895 0.0 0.0 -10.2460708626 -10.2460708626
2011-12-13 06:00:00 58.29335 0.000214344806385 0.0 -1.21422403655 -1.21422403655
2011-12-14 05:30:00 31.98307 0.0 0.0 0.616399174078 0.616399174078
2011-12-14 06:00:00 113.3059 8.15839334686e-05 0.0 -0.208930064837 -0.208930064837
2011-12-17 06:00:00 87.24711 3.55919749401e-07 0.0 3.63946477321 3.63946477321
2011-12-18 06:00:00 33.70848 1.11591234051e-08 0.0 1.11559123503 1.11559123503
2011-12-20 05:30:00 35.73178 0.0 0.0 2.85765367183 2.85765367183
2011-12-20 06:00:00 91.50012 6.84399333617e-15 0.0 0.105309915866 0.105309915866
2011-12-21 06:00:00 87.5268 7.68873193287e-23 0.0 0.359683575209 0.359683575209
2011-12-25 05:30:00 11.21276 0.0 0.0 2.10706296021 2.10706296021
2011-12-25 06:00:00 49.31652 0.0 0.0 1.54687185488 1.54687185488
2011-12-26 06:00:00 66.51482 0.0 0.0 9.18729072984 9.18729072984
2011-12-27 05:30:00 20.52332 0.0 0.0 1.95402309848 1.95402309848
2011-12-27 06:00:00 93.17953 0.0 0.0 0.0391420005066 0.0391420005066
2011-12-28 05:30:00 17.48715 0.0 0.0 2.12359388986 2.12359388986
2011-12-28 06:00:00 98.85627 0.0 0.0 1.02517237472 1.02517237472
2011-12-29 05:30:00 20.5801 0.0 0.0 1.73855261531 1.73855261531
2011-12-29 06:00:00 95.69814 0.0 0.0 0.190625257201 0.190625257201
2011-12-30 05:30:00 17.37488 0.0 0.0 0.754916585345 0.754916585345
2011-12-30 06:00:00 85.27795 0.0 0.0 1.78188376149 1.78188376149
2012-01-05 05:30:00 11.43538 0.0 0.0 2.16664815238 2.16664815238
2012-01-05 06:00:00 43.99102 0.0 0.0 1.21831609838 1.21831609838
2012-01-06 05:30:00 12.00184 0.0 0.0 1.81648065021 1.81648065021
2012-01-11 06:00:00 40.69904 0.0 0.0 -0.751963066762 -0.751963066762
2012-01-13 06:00:00 31.46369 0.0 0.0 2.96944975683 2.96944975683
2012-01-14 06:00:00 56.40198 0.0 0.0 0.110827271811 0.110827271811
2012-01-15 06:00:00 48.7319 0.0 0.0 1.00229113836 1.00229113836
2012-01-16 06:00:00 43.90872 0.0 0.0 2.65752730606 2.65752730606
2012-01-17 06:00:00 44.13012 0.0 0.0 2.32716231143 2.32716231143
2012-01-17 06:30:00 127.9109 6.5106237054 0.0 0.0042765724439 0.0042765724439
2012-01-19 06:00:00 39.79247 0.0 0.0 1.64946404935 1.64946404935
2012-01-19 06:30:00 124.6088 1.47051344247 0.0 0.61325559965 0.61325559965
2012-01-20 06:00:00 37.48366 0.0 0.0 2.78578142399 2.78578142399
2012-01-20 06:30:00 119.1935 0.480766734307 0.0 1.56429529671 1.56429529671
2012-01-21 06:00:00 31.28969 0.0 0.0 2.84671958735 2.84671958735
2012-01-21 06:30:00 113.9841 0.100637852785 0.0 -0.0609903806067 -0.0609903806067
2012-01-24 06:00:00 27.7181 0.0 0.0 0.413154771981 0.413154771981
2012-01-24 06:30:00 106.3347 9.94110943244e-08 0.0 -0.95748747171 -0.95748747171
2012-01-26 06:00:00 24.33282 0.0 0.0 0.926780084836 0.926780084836
2012-01-26 06:30:00 99.41067 0.0 0.0 -0.427456008191 -0.427456008191
2012-01-28 06:00:00 12.04299 0.0 0.0 2.53721299149 2.53721299149
2012-01-28 06:30:00 36.48751 0.0 0.0 2.72132011793 2.72132011793
2012-01-31 06:00:00 13.93341 0.0 0.0 3.95144342136 3.95144342136
2012-01-31 06:30:00 78.74403 0.0 0.0 2.01908927714 2.01908927714
2012-02-01 06:30:00 93.17162 0.0 0.0 0.839746707702 0.839746707702
2012-02-05 06:30:00 61.85131 0.0 0.0 3.34540667112 3.34540667112
2012-02-08 06:30:00 70.20588 0.0 0.0 1.75480057427 1.75480057427
2012-02-11 06:30:00 57.03128 0.0 0.0 3.74712810473 3.74712810473
2012-02-13 06:30:00 50.1085 0.0 0.0 1.3040535174 1.3040535174
2012-02-15 07:00:00 131.7922 3.34307878677 0.0 1.67556865255 1.67556865255
2012-02-16 07:00:00 39.00459 1.37973313591 0.0 1.74736683241 1.74736683241
2012-02-21 06:30:00 23.75285 0.0 0.0 3.7983171999 3.7983171999
2012-02-21 07:00:00 83.45065 5.2869620351e-08 0.0 2.40043274946 2.40043274946
2012-02-26 06:30:00 20.26641 0.0 0.0 2.56571379644 2.56571379644
2012-02-26 07:00:00 97.48089 0.0 0.0 0.817812993163 0.817812993163
2012-02-29 07:00:00 39.46045 0.0 0.0 1.57071533828 1.57071533828
2012-03-02 07:00:00 26.35533 0.0 0.0 -1.92331414447 -1.92331414447
2012-03-03 07:00:00 16.39093 0.0 0.0 3.91662642989 3.91662642989
2012-03-06 06:30:00 10.02979 0.0 0.0 3.52380293201 3.52380293201
2012-03-06 07:00:00 75.71791 0.0 0.0 2.49808807223 2.49808807223
2012-03-07 07:00:00 78.27703 0.0 0.0 2.01135633074 2.01135633074
2012-03-08 07:00:00 72.49073 0.0 0.0 1.37317787022 1.37317787022
2012-03-10 07:00:00 62.99122 0.0 0.0 3.34475589103 3.34475589103
2012-03-12 07:00:00 37.54271 0.0 0.0 1.77171247136 1.77171247136
2012-03-15 07:00:00 14.77107 0.0 0.0 4.46057838584 4.46057838584
2012-03-15 07:30:00 62.89638 7.56325890346 0.0 2.59585038993 2.59585038993
2012-03-16 07:30:00 14.32858 4.44311358953 0.0 1.47158651589 1.47158651589
2012-03-17 07:00:00 46.25202 0.0 0.0 3.17973385094 3.17973385094
2012-03-17 07:30:00 140.0092 2.33015041632 0.0 2.99871495815 2.99871495815
2012-03-18 07:00:00 47.59256 0.0 0.0 3.11116385799 3.11116385799
2012-03-20 07:00:00 31.79199 0.0 0.0 2.37540118208 2.37540118208
2012-03-20 07:30:00 90.98168 0.100892202151 0.0 0.0528838122437 0.0528838122437
2012-03-21 07:00:00 31.23541 0.0 0.0 3.60388130454 3.60388130454
2012-03-21 07:30:00 61.20091 0.0164921350966 0.0 1.68399888768 1.68399888768
2012-03-22 07:00:00 37.24337 0.0 0.0 2.44901088456 2.44901088456
2012-03-22 07:30:00 127.5248 0.00120438145128 0.0 2.46983156917 2.46983156917
2012-03-24 07:00:00 16.18912 0.0 0.0 2.12389792556 2.12389792556
2012-03-24 07:30:00 81.80565 1.21206109431e-08 0.0 1.26522904221 1.26522904221
2012-04-06 18:00:00 34.382 0.012282503669 0.0 1.8031500352 1.8031500352
2012-04-09 18:00:00 35.58956 3.2121747839e-30 0.0 0.525523944361 0.525523944361
2012-04-10 18:00:00 20.49554 0.0 0.0 0.131139958001 0.131139958001
2012-04-11 18:00:00 20.74329 0.0 0.0 0.580038328259 0.580038328259
2012-04-23 17:30:00 28.62116 8.11313597773 0.0 1.55596646343 1.55596646343
2012-04-24 17:30:00 24.64674 4.45596657816 0.0 0.475240776776 0.475240776776
2012-04-25 17:30:00 24.86754 2.12791689302 0.0 2.99846091356 2.99846091356
2012-04-29 17:30:00 32.44363 0.00482614565141 0.0 -0.346475693942 -0.346475693942
2012-05-01 17:30:00 26.2957 3.45106443863e-07 0.0 -0.216411890995 -0.216411890995
2012-05-02 17:30:00 23.52782 1.19376060319e-12 0.0 1.17422992024 1.17422992024
2012-05-05 17:30:00 10.1285 0.0 0.0 0.250207497433 0.250207497433
2012-05-07 17:30:00 17.47918 0.0 0.0 -0.699395862589 -0.699395862589
2012-05-23 17:00:00 45.44139 8.03247199896 0.0 -0.375011133123 -0.375011133123

In [ ]: