In [1]:

    

import qcio
import qcutils


    

In [2]:

    

ds_ew=qcio.nc_read_series(qcio.get_filename_dialog())


    

In [3]:

    

ds_pi=qcio.nc_read_series(qcio.get_filename_dialog())


    

In [4]:

    

ds_nc=qcio.nc_read_series(qcio.get_filename_dialog())


    

In [5]:

    

ah_hmp_ew,f=qcutils.GetSeriesasMA(ds_ew,"Ah_HMP_3m")


    

In [6]:

    

ah_hmp_pi,f=qcutils.GetSeriesasMA(ds_pi,"Ah_HMP_3m")


    

In [7]:

    

ah_hmp_nc,f=qcutils.GetSeriesasMA(ds_nc,"Ah_HMP_3m")


    

In [8]:

    

ah_7500_ew,f=qcutils.GetSeriesasMA(ds_ew,"Ah_7500_Av")


    

In [9]:

    

ah_7500_pi,f=qcutils.GetSeriesasMA(ds_pi,"Ah_7500_Av")


    

In [10]:

    

ah_7500_nc,f=qcutils.GetSeriesasMA(ds_nc,"Ah_7500_Av")


    

In [11]:

    

fig_nc=figure()


    

In [12]:

    

plot(ah_7500_nc,ah_hmp_nc,'b.')


    

    
Out[12]:





[<matplotlib.lines.Line2D at 0xbb94fd0>]

In [13]:

    

xlabel("Ah_7500 (g/m3)")


    

    
Out[13]:





<matplotlib.text.Text at 0xbc50bd0>

In [14]:

    

ylabel("Ah_HMP (g/m3)")


    

    
Out[14]:





<matplotlib.text.Text at 0xbc58590>

In [15]:

    

title("Daly Pasture: 2009; No correction")


    

    
Out[15]:





<matplotlib.text.Text at 0xbd59390>

In [16]:

    

fig_ew=figure()


    

In [17]:

    

plot(ah_7500_ew,ah_hmp_ew,"b.")


    

    
Out[17]:





[<matplotlib.lines.Line2D at 0xbbaf3d0>]

In [18]:

    

xlabel("Ah_7500 (g/m3)")


    

    
Out[18]:





<matplotlib.text.Text at 0xbc5ca90>

In [19]:

    

ylabel("Ah_HMP (g/m3)")


    

    
Out[19]:





<matplotlib.text.Text at 0x9c45b10>

In [20]:

    

title("Daly Pasture: 2009; Emma's correction")


    

    
Out[20]:





<matplotlib.text.Text at 0x4c30750>

In [21]:

    

fig_pi=figure()


    

In [22]:

    

plot(ah_7500_pi,ah_hmp_pi,"b.")


    

    
Out[22]:





[<matplotlib.lines.Line2D at 0xbbc6350>]

In [23]:

    

xlabel("Ah_7500 (g/m3)")


    

    
Out[23]:





<matplotlib.text.Text at 0xbbc0650>

In [24]:

    

ylabel("Ah_HMP (g/m3)")


    

    
Out[24]:





<matplotlib.text.Text at 0xbbc6050>

In [25]:

    

title("Daly Pasture: 2009; Peter's correction")


    

    
Out[25]:





<matplotlib.text.Text at 0xc768950>

In [26]:

    

fig_ewpi=figure()


    

In [27]:

    

plot(ah_7500_ew,ah_7500_pi,"b.")


    

    
Out[27]:





[<matplotlib.lines.Line2D at 0xcb20cd0>]

In [28]:

    

xlabel("Ah_7500 (g/m3, Emma's)")


    

    
Out[28]:





<matplotlib.text.Text at 0xc78e150>

In [29]:

    

ylabel("Ah_7500 (g/m3, Peter's)")


    

    
Out[29]:





<matplotlib.text.Text at 0xc78ee90>

In [30]:

    

title("Daly Pasture: 2009; Emma's and Peter's Corrections")


    

    
Out[30]:





<matplotlib.text.Text at 0xc79e890>

In [31]:

    

fe_ew,f=qcutils.GetSeriesasMA(ds_ew,"Fe")


    

In [32]:

    

fe_pi,f=qcutils.GetSeriesasMA(ds_pi,"Fe")


    

In [33]:

    

fe_nc,f=qcutils.GetSeriesasMA(ds_nc,"Fe")


    

In [36]:

    

fig1=figure()


    

In [37]:

    

plot(fe_nc,fe_ew,"b.")


    

    
Out[37]:





[<matplotlib.lines.Line2D at 0xbfcdb90>]

In [38]:

    

xlabel("Fe (W/m2, no correction)")


    

    
Out[38]:





<matplotlib.text.Text at 0x9d47690>

In [39]:

    

ylabel("Fe (W/m2, Peter's)")


    

    
Out[39]:





<matplotlib.text.Text at 0xc673cd0>

In [40]:

    

title("Daly Pasture: 2009; Peter's vs no correction")


    

    
Out[40]:





<matplotlib.text.Text at 0xc680c90>

In [41]:

    

fig2=figure()


    

In [42]:

    

plot(fe_nc,fe_pi,"b.")


    

    
Out[42]:





[<matplotlib.lines.Line2D at 0xc002fd0>]

In [43]:

    

xlabel("Fe (W/m2, no correction)")


    

    
Out[43]:





<matplotlib.text.Text at 0xbfe0c50>

In [44]:

    

ylabel("Fe (W/m2, Peter's)")


    

    
Out[44]:





<matplotlib.text.Text at 0xbfe8690>

In [45]:

    

title("Daly Pasture: 2009; Peter's vs no correction")


    

    
Out[45]:





<matplotlib.text.Text at 0xbff5490>

In [46]:

    

fig3=figure()


    

In [47]:

    

plot(fe_nc,fe_ew,"b.")


    

    
Out[47]:





[<matplotlib.lines.Line2D at 0xc085ad0>]

In [48]:

    

xlabel("Fe (W/m2, no correction)")


    

    
Out[48]:





<matplotlib.text.Text at 0xc058690>

In [49]:

    

ylabel("Fe (W/m2, Emma's)")


    

    
Out[49]:





<matplotlib.text.Text at 0xc061090>

In [50]:

    

title("Daly Pasture: 2009; Emma's vs no correction")


    

    
Out[50]:





<matplotlib.text.Text at 0xc06a910>

In [51]:

    

fig4=figure()


    

In [52]:

    

plot(fe_ew,fe_pi,"b.")


    

    
Out[52]:





[<matplotlib.lines.Line2D at 0xd474590>]

In [53]:

    

xlabel("Fe (W/m2, Emma's)")


    

    
Out[53]:





<matplotlib.text.Text at 0xd44fc50>

In [54]:

    

ylabel("Fe (W/m2, Peter's)")


    

    
Out[54]:





<matplotlib.text.Text at 0xd456650>

In [55]:

    

title("Da;y Pasture: 2009; Emma's vs Peter's")


    

    
Out[55]:





<matplotlib.text.Text at 0xd460e10>

In [ ]: