In [1]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

import numpy as np


    

In [3]:

    

x=[3280, 3382, 3455]
y=[2331, 1935, 1818]


    

In [4]:

    

dx = np.diff(x)
dy = np.diff(y)
print dx
print dy


    

    




[102  73]
[-396 -117]

In [8]:

    

print (dy * 100) / dx / 1000. # mm/100m * 1m / 1000mm


    

    




[-0.389 -0.161]

In [22]:

    

import matplotlib as mpl
import matplotlib.pyplot as plt


    

In [23]:

    

plt.plot(x,y)
plt.xlabel('Elevation (m)')
plt.ylabel('Ablation (mm)')
plt.title('Ablation measurements at Karabatkak Glacier')


    

    
Out[23]:





<matplotlib.text.Text at 0x110e52ad0>






    

In [24]:

    

dydx = dy.astype(float)/dx
print dydx # mm / m


    

    




[-3.88235294 -1.60273973]

In [25]:

    

plt.savefig('karabatkak_abl_grad.png')


    

    






<matplotlib.figure.Figure at 0x110f0a210>

In [25]: