In [ ]:

from IPython.core.display import HTML

with open('../../common/creativecommons.html', 'r') as f:

with open('../../common/custom.css', 'r') as f:

HTML(styles)

text = 'Check this post at'
name = get_notebook_name()
link = """<p>%s <a href="%s/%s"><em>nbviewer</em>.</a></p>""" % (text, uri, name)




In [ ]:

import numpy as np
import matplotlib.pyplot as plt
from mpltools import style
style.use('dark_background')




In [ ]:

# Constants.
x, y = 0, 0
n0 = 0.5
g = 9.8
rho = 1025
K = 0.02
l = 0.014
k = 0.014
w = 0.43
H = 100

X = np.arange(1, 21)
t = np.arange(1, 21)



Equations:



In [ ]:

N = n0 * np.cos((k * x) + (l * y) - (w * t))
P0 = (rho * g * n0) / (np.cosh(K * H)) * np.cosh(K * (0 + H)) * np.cos((k * x) + (l * y)-(w * t))
U0 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (0 + H)) * np.cos((k * x) + (l * y)-(w * t))
V0 = ((l * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (0 + H)) * np.cos((k * x) + (l * y)-(w * t))
W0 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.sinh(K * (0 + H)) * np.sin((k * x) + (l * y)-(w * t))

P10 = (rho * g * n0) / (np.cosh(K * H)) * np.cosh(K * (-10 + H)) * np.cos((k * x) + (l * y)-(w * t))
U10 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-10 + H)) * np.cos((k * x) + (l * y)-(w * t))
V10 = ((l * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-10 + H)) * np.cos((k * x) + (l * y)-(w * t))
W10 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.sinh(K * (-10 + H)) * np.sin((k * x) + (l * y)-(w * t))

P50 = (rho * g * n0) / (np.cosh(K * H)) * np.cosh(K * (-50 + H)) * np.cos((k * x) + (l * y)-(w * t))
U50 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-50 + H)) * np.cos((k * x) + (l * y)-(w * t))
V50 = ((l * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-50 + H)) * np.cos((k * x) + (l * y)-(w * t))
W50 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.sinh(K * (-50 + H)) * np.sin((k * x) + (l * y)-(w * t))

P100 = (rho * g * n0) / (np.cosh(K * H)) * np.cosh(K * (-100 + H)) * np.cos((k * x) + (l * y)-(w * t))
U100 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-100 + H)) * np.cos((k * x) + (l * y)-(w * t))
V100 = ((l * g * n0) / w  * (np.cosh(K * H))) * np.cosh(K * (-100 + H)) * np.cos((k * x) + (l * y)-(w * t))
W100 = ((k * g * n0) / w  * (np.cosh(K * H))) * np.sinh(K * (-100 + H)) * np.sin((k * x) + (l * y)-(w * t))




In [ ]:

fig, ax = plt.subplots()
ax.plot(N)
ax.set_ylabel(r'$\eta$ [m]')
ax.set_xlabel('Time [s]')

fig, ax = plt.subplots()
ax.plot(X, U0, X, U10, X, U50, X, U100)
ax.set_ylabel(r'U [m s$^{-1}$]')
ax.set_xlabel('Time [s]')
ax.legend(('0 m', '10 m', '50 m', '100 m'))

fig, ax = plt.subplots()
plot(X, V0, X, V10, X, V50, X, V100)
ax.set_ylabel(r'V [m s$^{-1}$]')
ax.set_xlabel('Time [s]')
ax.legend(('0 m', '10 m', '50 m', '100 m'))

fig, ax = plt.subplots()
ax.plot(X, W0, X, W10, X, W50, X, W100)
ax.set_ylabel(r'W [m s$^{-1}$]')
ax.set_xlabel('Time [s]')
ax.legend(('0 m', '10 m', '50 m', '100 m'))

fig, ax = plt.subplots()
ax.plot(X, P0, X, P10, X, P50, X, P100)
ylabel('Pressure [Pa]')
ax.set_xlabel('Time [s]')
ax.legend(('0 m', '10 m', '50 m', '100 m'))




In [ ]:

HTML(html)