In [4]:

    

x = range(100)
y = [i ** 2 for i in x]
print y


    

In [4]:

    

1 / 0


    

In [5]:

    

%hist


    

In [6]:

    

%pushd /
!ls
%popd


    

In [7]:

    

!!ls


    

In [8]:

    

for f in _:
    !wc "$f"


    

In [9]:

    

plot(x, y)


    

In [1]:

    

from mpl_toolkits.basemap import Basemap
from netCDF4 import Dataset
with Dataset('/data/cmip5/output1/MOHC/HadGEM2-ES/rcp85/mon/atmos/Amon/r1i1p1/'
             'v20111215/tas/tas_Amon_HadGEM2-ES_rcp85_r1i1p1_200512-203011.nc') as d:
    lat = d.variables['lat'][:]
    lon = d.variables['lon'][:]
    tas = d.variables['tas'][0,...]
m = Basemap(projection = 'merc', llcrnrlat = -80, urcrnrlat = 80,
            llcrnrlon = lon[0], urcrnrlon = lon[-1], lat_ts = 0)
m.drawcoastlines()
m.drawparallels(range(-90, 90, 30))
m.drawmeridians(range(-0, 360, 30))
x, y = m(*np.meshgrid(lon, lat))
m.contourf(x, y, tas, 100)
plt.colorbar()