From 890c0b0a9b54e47bebc65b6b93860ffa1950e72b Mon Sep 17 00:00:00 2001
From: Jan Benda <jan.benda@uni-tuebingen.de>
Date: Fri, 2 Apr 2021 09:30:28 +0200
Subject: [PATCH] added efishtitle with two fish

---
 efishtitle.pdf | Bin 254317 -> 254317 bytes
 efishtitle.py  |  22 ++++++++++++++++++++++
 2 files changed, 22 insertions(+)

diff --git a/efishtitle.pdf b/efishtitle.pdf
index 60413563638623e9712a6ad157ceae69325f72bb..65473762df0099025524b5d100357ec5f9413507 100644
GIT binary patch
delta 26
icmaERi2v;&{)QID7N#xCx9+hRS(;dGzju$hp9uh^I}2I>

delta 26
icmaERi2v;&{)QID7N#xCx9+hRnwT1Hzju$hp9uh@)C)}j

diff --git a/efishtitle.py b/efishtitle.py
index 38b7774..24c4c12 100644
--- a/efishtitle.py
+++ b/efishtitle.py
@@ -154,6 +154,27 @@ def plot_threewavefish(ax):
     ax.set_aspect('equal')
 
     
+def plot_twowavefishside(ax):
+    maxx = 30.0
+    maxy = maxx * aspect_ratio(ax)
+    x = np.linspace(-maxx, maxx, 600)
+    y = np.linspace(-maxy, maxy, 200)
+    xx, yy = np.meshgrid(x, y)
+    fish1 = (('Alepto_male', 'side'), (-12, 0), (1, 0.2), 23.0, 15)
+    fish2 = (('Alepto', 'top'), (14, -1.5), (1.0, 0.45), 23.0, -20)
+    poles1 = efish_monopoles(*fish1[1:])
+    poles2 = efish_monopoles(*fish2[1:])
+    pot = epotential_meshgrid(xx, yy, None, poles1, poles2)
+    mz = 0.65
+    zz = squareroot_transform(pot/200, mz)
+    levels = np.linspace(-mz, mz, 16)
+    ax.contourf(x, y, -zz, levels)
+    ax.contour(x, y, -zz, levels, **csLine)
+    plot_fish(ax, *fish1, bodykwargs=fishBody, finkwargs=fishFins)
+    plot_fish(ax, *fish2, bodykwargs=fishBody, finkwargs=fishFins)
+    ax.set_aspect('equal')
+
+    
 if __name__ == "__main__":
     # hide spines:
     plt.rcParams['axes.spines.left'] = False
@@ -172,4 +193,5 @@ if __name__ == "__main__":
     fig, ax = plt.subplots(figsize=(15.2/2.54, 4.0/2.54))
     fig.subplots_adjust(left=0, right=1, top=1, bottom=0)
     plot_threewavefish(ax)
+    #plot_twowavefishside(ax)
     fig.savefig('efishtitle.pdf')