In [ ]:

    

from IPython.core.display import HTML

with open('../../common/creativecommons.html', 'r') as f:
    html = f.read()
    
with open('../../common/custom.css', 'r') as f:
    styles = f.read()
    
HTML(styles)

text = 'Check this post at'
uri = 'http://nbviewer.ipython.org/urls/raw.github.com/ocefpaf/python4oceanographers/master/content/downloads/notebooks'
name = get_notebook_name()
link = """<p>%s <a href="%s/%s"><em>nbviewer</em>.</a></p>""" % (text, uri, name)
html += str(link)


    

In [ ]:

    

import numpy as np
import seawater as sw
import matplotlib.pyplot as plt


    

In [ ]:

    

L = 33.2, 2.5
G = 35.0, 12.5
M = (L[0] + G[0]) * 0.5, (L[1] + G[1]) * 0.5
line = sw.dens0(M[0], M[1]) - 1000


    

In [ ]:

    

t = np.arange(0, 15 + 1, 1)
s = np.arange(32.5, 35 + 0.1, 0.1)

S, T = np.meshgrid(s, t)
sigma_theta = sw.dens0(S, T) - 1000
levels = [24.5, 25, 26, 26.5, 27, 27.5, 28]


    

In [ ]:

    

from matplotlib import rcParams
from matplotlib.ticker import MultipleLocator
rcParams['xtick.direction'] = 'out'
rcParams['ytick.direction'] = 'out'
fig, ax = plt.subplots(figsize=(6, 4))
ax.xaxis.set_major_locator(MultipleLocator(0.5))
ax.xaxis.set_minor_locator(MultipleLocator(0.1))
ax.yaxis.set_major_locator(MultipleLocator(5))
ax.yaxis.set_minor_locator(MultipleLocator(1))
ax.set_ylabel(u"Temperatura \u00b0C")
ax.set_xlabel(r"Salinidade [g kg$^{-1}$]")
ax.axis([32.5, 35, 0, 15])

cs = ax.contour(s, t, sigma_theta, colors='black', levels=levels)
ax.clabel(cs, fontsize=9, inline=1, fmt='%2.1f')
sg = ax.contour(s, t, sigma_theta, linestyle='--', colors='grey', levels=[0, line])
ax.clabel(sg, fontsize=9, inline=1, fmt='%2.1f')
ax.plot(L[0], L[1], 'ko')
ax.plot(G[0], G[1], 'ko')
ax.plot(M[0], M[1], 'ro')
xoffset, yoffset = 0.08, 0.5
ax.text(L[0] + xoffset, L[1] - yoffset, 'L')
ax.text(G[0] - xoffset, G[1] + yoffset, 'G')
ax.text(M[0] + xoffset, M[1] - yoffset, 'M')
fig.savefig("Diagrama-TS.pdf", dpi=150)


    

In [ ]:

    

HTML(html)