from matplotlib import gridspec
import math
import matplotlib
import nanopores as nano
import nanopores.geometries.pughpore as pughpore
import matplotlib.pyplot as plt
from matplotlib.ticker import FormatStrFormatter
import numpy as np
import os
import sys
import nanopores.tools.fields as f
HOME = os.path.expanduser("~")
PAPERDIR = os.path.join(HOME, "papers", "paper-howorka")
FIGDIR = os.path.join(PAPERDIR, "figures", "")
DATADIR = os.path.join(HOME, "Dropbox", "nanopores", "fields")
f.set_dir_mega()
number = False
geop = nano.Params(pughpore.params)
hpore = geop.hpore
fieldsname = 'events3_onlyone_5'
params = dict(avgbind1=2e7,
avgbind2=3e4,
P_bind1=8.e-2,
P_bind2=0 * 3e-1,
z0=hpore / 2. + 0.)
#cmap=matplotlib.cm.get_cmap('viridis')
data = f.get_fields(fieldsname, **params)
figname = fieldsname + '_%.1e_%.1e_%.1e_%.1e' % (
# (c) 2016 Gregor Mitscha-Baude
import os
HOME = os.path.expanduser("~")
#PAPERDIR = os.path.join(HOME, "papers", "paper-howorka")
#FIGDIR = os.path.join(PAPERDIR, "figures", "")
FIGDIR = os.path.join(HOME, "Dropbox", "nanopores", "figures")
FIGDIR_HOWORKA = os.path.join(HOME, "Dropbox", "Paper Howorka", "figures")
#DATADIR = os.path.join(HOME, "Dropbox", "nanopores", "fields")
import nanopores.tools.fields as fields
fields.set_dir_mega()
#fields.set_dir(DATADIR)