In [1]:

    

# import numpy as np

# # !/usr/bin/env python3
# # -*- coding: utf-8 -*-
# """
# Created on 20181219

# @author: zhangji

# Trajection of a ellipse, Jeffery equation. 
# """

# %pylab inline
# pylab.rcParams['figure.figsize'] = (25, 11)
# fontsize = 40

# import numpy as np
# import scipy as sp
# from scipy.optimize import leastsq, curve_fit
# from scipy import interpolate
# from scipy.interpolate import interp1d
# from scipy.io import loadmat, savemat
# # import scipy.misc

# import matplotlib
# from matplotlib import pyplot as plt
# from matplotlib import animation, rc
# import matplotlib.ticker as mtick
# from mpl_toolkits.axes_grid1.inset_locator import inset_axes, zoomed_inset_axes
# from mpl_toolkits.mplot3d import Axes3D, axes3d

# from sympy import symbols, simplify, series, exp
# from sympy.matrices import Matrix
# from sympy.solvers import solve

# from IPython.display import display, HTML
# from tqdm import tqdm_notebook as tqdm
# import pandas as pd
# import re
# from scanf import scanf
# import os
# import glob

# from codeStore import support_fun as spf
# from src.support_class import *
# from src import stokes_flow as sf

# rc('animation', html='html5')
# PWD = os.getcwd()
# font = {'size': 20}
# matplotlib.rc('font', **font)
# np.set_printoptions(linewidth=90, precision=5)

%load_ext autoreload
%autoreload 2

import os
import glob
import natsort 
import numpy as np
import scipy as sp
from scipy.optimize import leastsq, curve_fit
from scipy import interpolate, integrate
from scipy import spatial, signal
# from scipy.interpolate import interp1d
from scipy.io import loadmat, savemat
# import scipy.misc
# import importlib
from IPython.display import display, HTML
import pandas as pd
import pickle
import re
from scanf import scanf

import matplotlib
# matplotlib.use('agg')
from matplotlib import pyplot as plt
import matplotlib.colors as colors
from matplotlib import animation, rc
import matplotlib.ticker as mtick
from mpl_toolkits.axes_grid1.inset_locator import inset_axes, zoomed_inset_axes
from mpl_toolkits.mplot3d import Axes3D, axes3d
from mpl_toolkits.axes_grid1.axes_divider import make_axes_locatable
from mpl_toolkits.mplot3d.art3d import Line3DCollection
from matplotlib import cm

from tqdm.notebook import tqdm as tqdm_notebook
from tqdm import tqdm
from time import time
from src.support_class import *
from src import jeffery_model as jm
from codeStore import support_fun as spf
from codeStore import support_fun_table as spf_tb

# %matplotlib notebook
%matplotlib inline
rc('animation', html='html5')
fontsize = 40
PWD = os.getcwd()


    

In [2]:

    

fig = plt.figure(figsize=(2, 2))
fig.patch.set_facecolor('white')
ax0 = fig.add_subplot(1, 1, 1)


    

In [3]:

    

job_dir = 'ecoC01B05_wt0.7_psi-0c'
table_name = 'ecoC01B05_tao1_wm0.7'
t_headle = '(.*?).pickle'


    

In [4]:

    

t_dir = os.path.join(PWD, job_dir)

data = spf_tb.load_table_data_pickle_dir(t_dir, t_headle)
lst_eta = data.lst_eta
theta_max_fre = data.theta_max_fre
phi_max_fre = data.phi_max_fre
psi_max_fre = data.psi_max_fre
eta_max_fre = data.eta_max_fre
data_idx = data.data_idx.fillna(-1).astype(int)


    

In [5]:

    

# put images with same frequence into a subdirect
check_fre_list = [0.0140, 0.0150]
atol_fre_list =  [0.0002, 0.0002]
Table_t_range = np.array((9500, np.inf))

def sub_seperate_a(i1, type_fre, iidx, Table_t_range):
    t1 = np.zeros_like(iidx[0])
    return t1
    
sub_seperate = {0: sub_seperate_a, 
                1: sub_seperate_a, 
                'other': sub_seperate_a, }

tfre = theta_max_fre.copy()
type_fre = tfre.copy()
type_fre.iloc[:, :] = -1
i0 = 0
process_total = 0
for i1, (check_fre, atol_fre) in enumerate(zip(check_fre_list[:], atol_fre_list)):
    use_idx = np.isclose(tfre, check_fre, rtol=0, atol=atol_fre)
    iidx = np.where(use_idx)
    process_total = process_total + np.sum(use_idx)
    t1 = sub_seperate[i1](i1, type_fre, iidx, Table_t_range)
    type_fre.iloc[use_idx] = i0 + t1
    i0 = i0 + 1 + np.max(t1)
    print(i0, np.sum(use_idx))

# assert 1 == 2
# plot one of the remaind cases
if np.any(type_fre.values == -1):
    iidx = np.where(type_fre.values == -1)
    t1 = sub_seperate['other'](-1, type_fre, iidx, Table_t_range)
    for (i3, i4), t11 in zip(np.vstack(iidx).T, t1):
        type_fre.iloc[i3].iloc[i4] = i0 + t11

assert np.all(type_fre.values >= 0)
type_fre2 = type_fre.copy()
# type_fre2.iloc[type_fre.values == 0] = 0
# type_fre2.iloc[type_fre.values == 1] = 4
# type_fre2.iloc[type_fre.values == 2] = 4
# type_fre2.iloc[type_fre.values == 3] = 2
# type_fre2.iloc[type_fre.values == 4] = 3
# type_fre2.iloc[type_fre.values == 5] = 2
# type_fre2.iloc[type_fre.values == 6] = 1
spf_tb.show_traj_phase_map_type(type_fre2)
# spf_tb.save_separate_angleList_fft(job_dir, tfre, check_fre_list, atol_fre_list)


    

    




1 101
2 931






    
Out[5]:





True






    

In [6]:

    

t1 = []
for i0 in np.arange(type_fre2.values.max() + 1):
    t1.append(np.isclose(type_fre2.values, i0).sum())

print(np.sum(t1))
print(t1)
print(t1 / np.sum(t1))


    

    




1035
[101, 931, 3]
[0.09758454 0.89951691 0.00289855]

In [8]:

    

fast_mode1 = 0
figsize = np.array((16, 9)) * 0.5
dpi = 200

tipical_th_ph_list = []
for i0 in np.arange(type_fre2.values.max() + 1)[:]:
    iidx = np.where(np.isclose(type_fre2.values, i0))
    spf_tb.phase_map_show_idx(type_fre2, tipical_th_ph_list, iidx, job_dir, table_name, fast_mode=fast_mode1, )


    

    




-ini_theta 0.143000 -ini_phi 1.142000
-ini_theta 0.142800 -ini_phi 1.142397 -ini_psi 1.570796






    




/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/figure.py:2369: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
  warnings.warn("This figure includes Axes that are not compatible "






    




-ini_theta 0.000000 -ini_phi 0.000000
-ini_theta 0.000000 -ini_phi 0.000000 -ini_psi 1.570796






    




/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/figure.py:2369: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
  warnings.warn("This figure includes Axes that are not compatible "






    




-ini_theta 1.571000 -ini_phi 1.428000
-ini_theta 1.570796 -ini_phi 1.427997 -ini_psi 1.570796
Error in callback <function install_repl_displayhook.<locals>.post_execute at 0x7f5fad59a268> (for post_execute):






    




---------------------------------------------------------------------------
KeyboardInterrupt                         Traceback (most recent call last)
/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/pyplot.py in post_execute()
    107             def post_execute():
    108                 if matplotlib.is_interactive():
--> 109                     draw_all()
    110 
    111             # IPython >= 2

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/_pylab_helpers.py in draw_all(cls, force)
    130         for f_mgr in cls.get_all_fig_managers():
    131             if force or f_mgr.canvas.figure.stale:
--> 132                 f_mgr.canvas.draw_idle()
    133 
    134 atexit.register(Gcf.destroy_all)

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/backend_bases.py in draw_idle(self, *args, **kwargs)
   1897         if not self._is_idle_drawing:
   1898             with self._idle_draw_cntx():
-> 1899                 self.draw(*args, **kwargs)
   1900 
   1901     def draw_cursor(self, event):

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/backends/backend_agg.py in draw(self)
    400         toolbar = self.toolbar
    401         try:
--> 402             self.figure.draw(self.renderer)
    403             # A GUI class may be need to update a window using this draw, so
    404             # don't forget to call the superclass.

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/artist.py in draw_wrapper(artist, renderer, *args, **kwargs)
     48                 renderer.start_filter()
     49 
---> 50             return draw(artist, renderer, *args, **kwargs)
     51         finally:
     52             if artist.get_agg_filter() is not None:

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/figure.py in draw(self, renderer)
   1647 
   1648             mimage._draw_list_compositing_images(
-> 1649                 renderer, self, artists, self.suppressComposite)
   1650 
   1651             renderer.close_group('figure')

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/image.py in _draw_list_compositing_images(renderer, parent, artists, suppress_composite)
    136     if not_composite or not has_images:
    137         for a in artists:
--> 138             a.draw(renderer)
    139     else:
    140         # Composite any adjacent images together

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/artist.py in draw_wrapper(artist, renderer, *args, **kwargs)
     48                 renderer.start_filter()
     49 
---> 50             return draw(artist, renderer, *args, **kwargs)
     51         finally:
     52             if artist.get_agg_filter() is not None:

/anaconda3/envs/py35/lib/python3.5/site-packages/mpl_toolkits/mplot3d/axes3d.py in draw(self, renderer)
    285                 sorted(self.collections,
    286                        key=lambda col: col.do_3d_projection(renderer),
--> 287                        reverse=True)):
    288             col.zorder = zorder_offset + i
    289         for i, patch in enumerate(

/anaconda3/envs/py35/lib/python3.5/site-packages/mpl_toolkits/mplot3d/axes3d.py in <lambda>(col)
    284         for i, col in enumerate(
    285                 sorted(self.collections,
--> 286                        key=lambda col: col.do_3d_projection(renderer),
    287                        reverse=True)):
    288             col.zorder = zorder_offset + i

/anaconda3/envs/py35/lib/python3.5/site-packages/mpl_toolkits/mplot3d/art3d.py in do_3d_projection(self, renderer)
    245             self._segments3d]
    246         segments_2d = [np.column_stack([xs, ys]) for xs, ys, zs in xyslist]
--> 247         LineCollection.set_segments(self, segments_2d)
    248 
    249         # FIXME

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/collections.py in set_segments(self, segments)
   1239             _segments = self._add_offsets(_segments)
   1240 
-> 1241         self._paths = [mpath.Path(_seg) for _seg in _segments]
   1242         self.stale = True
   1243 

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/collections.py in <listcomp>(.0)
   1239             _segments = self._add_offsets(_segments)
   1240 
-> 1241         self._paths = [mpath.Path(_seg) for _seg in _segments]
   1242         self.stale = True
   1243 

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/path.py in __init__(self, vertices, codes, _interpolation_steps, closed, readonly)
    147         self._codes = codes
    148         self._interpolation_steps = _interpolation_steps
--> 149         self._update_values()
    150 
    151         if readonly:

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/path.py in _update_values(self)
    191         return pth
    192 
--> 193     def _update_values(self):
    194         self._simplify_threshold = rcParams['path.simplify_threshold']
    195         self._should_simplify = (

KeyboardInterrupt: 





    













    













    




Error in callback <function flush_figures at 0x7f5fad5a3bf8> (for post_execute):






    




---------------------------------------------------------------------------
KeyboardInterrupt                         Traceback (most recent call last)
/anaconda3/envs/py35/lib/python3.5/site-packages/ipykernel/pylab/backend_inline.py in flush_figures()
    119         # ignore the tracking, just draw and close all figures
    120         try:
--> 121             return show(True)
    122         except Exception as e:
    123             # safely show traceback if in IPython, else raise

/anaconda3/envs/py35/lib/python3.5/site-packages/ipykernel/pylab/backend_inline.py in show(close, block)
     41             display(
     42                 figure_manager.canvas.figure,
---> 43                 metadata=_fetch_figure_metadata(figure_manager.canvas.figure)
     44             )
     45     finally:

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/display.py in display(include, exclude, metadata, transient, display_id, *objs, **kwargs)
    311             publish_display_data(data=obj, metadata=metadata, **kwargs)
    312         else:
--> 313             format_dict, md_dict = format(obj, include=include, exclude=exclude)
    314             if not format_dict:
    315                 # nothing to display (e.g. _ipython_display_ took over)

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/formatters.py in format(self, obj, include, exclude)
    178             md = None
    179             try:
--> 180                 data = formatter(obj)
    181             except:
    182                 # FIXME: log the exception

</anaconda3/envs/py35/lib/python3.5/site-packages/decorator.py:decorator-gen-9> in __call__(self, obj)

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/formatters.py in catch_format_error(method, self, *args, **kwargs)
    222     """show traceback on failed format call"""
    223     try:
--> 224         r = method(self, *args, **kwargs)
    225     except NotImplementedError:
    226         # don't warn on NotImplementedErrors

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/formatters.py in __call__(self, obj)
    339                 pass
    340             else:
--> 341                 return printer(obj)
    342             # Finally look for special method names
    343             method = get_real_method(obj, self.print_method)

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/pylabtools.py in <lambda>(fig)
    242 
    243     if 'png' in formats:
--> 244         png_formatter.for_type(Figure, lambda fig: print_figure(fig, 'png', **kwargs))
    245     if 'retina' in formats or 'png2x' in formats:
    246         png_formatter.for_type(Figure, lambda fig: retina_figure(fig, **kwargs))

/anaconda3/envs/py35/lib/python3.5/site-packages/IPython/core/pylabtools.py in print_figure(fig, fmt, bbox_inches, **kwargs)
    126 
    127     bytes_io = BytesIO()
--> 128     fig.canvas.print_figure(bytes_io, **kw)
    129     data = bytes_io.getvalue()
    130     if fmt == 'svg':

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/backend_bases.py in print_figure(self, filename, dpi, facecolor, edgecolor, orientation, format, bbox_inches, **kwargs)
   2047                         orientation=orientation,
   2048                         dryrun=True,
-> 2049                         **kwargs)
   2050                     renderer = self.figure._cachedRenderer
   2051                     bbox_artists = kwargs.pop("bbox_extra_artists", None)

/anaconda3/envs/py35/lib/python3.5/site-packages/matplotlib/backends/backend_agg.py in print_png(self, filename_or_obj, *args, **kwargs)
    521                 cbook.open_file_cm(filename_or_obj, "wb") as fh:
    522             _png.write_png(renderer._renderer, fh,
--> 523                             self.figure.dpi, metadata=metadata)
    524 
    525     def print_to_buffer(self):

KeyboardInterrupt: 

In [7]:

    

Table_t_range1 = np.array((9500, np.inf))
Table_t_range2 = np.array((9500, np.inf))
nshow = np.inf
# nshow = 2
fast_mode2 = True
figsize = np.array((16, 9)) * 0.5
dpi = 100

for i0 in np.arange(type_fre2.values.max() + 1)[:]:
    iidx = np.where(np.isclose(type_fre2.values, i0))
    spf_tb.phase_map_show_idx_list(type_fre2, iidx, job_dir, nshow=nshow, figsize=figsize, dpi=dpi, 
                                   Table_t_range1=Table_t_range1, Table_t_range2=Table_t_range2, fast_mode=fast_mode2, )


    

In [43]:

    

# for i0 in np.arange(type_fre2.values.max() + 1)[:]:
i0 = 1
Table_t_range = (9500, np.inf)

iidx = np.where(np.isclose(type_fre2.values, i0))
theta = type_fre.index.values[iidx[0]]
phi = type_fre.columns.values[iidx[1]]
theta_phi_list = np.vstack((theta, phi)).T
# print(theta_phi_list)

t1 = []
for theta, phi in theta_phi_list:
    tpick, _ = spf_tb.load_table_date_pickle(job_dir, theta, phi)
    Table_t = tpick['Table_t']
    idx = np.logical_and(Table_t >= Table_t_range[0], Table_t <= Table_t_range[1])
    Table_theta = tpick['Table_theta'][idx]
    Table_phi = tpick['Table_phi'][idx]
    Table_psi = tpick['Table_psi'][idx]
    t2 = (np.max(Table_theta), np.min(Table_theta), np.mean(Table_theta), 
          np.max(Table_phi), np.min(Table_phi), Table_psi[np.argmax(Table_phi)], 
          np.max(Table_psi), np.min(Table_psi), np.mean(Table_psi), )
    t1.append(t2)


    

In [44]:

    

figsize = np.array((16, 9))
dpi = 50

t1 = np.vstack(t1)
# t1[t1[:, 3] > 5.1950, :] = np.nan 
fig, axs = plt.subplots(3, 3, figsize=figsize, dpi=dpi)
for i0, axi in enumerate(axs.ravel()):
#     print(t1[:, i0])
    axi.plot(t1[:, i0], '.')


    

In [22]:

    

tipical_th_ph_list = []
for i0 in np.arange(type_fre2.values.max() + 1)[:]:
    iidx = np.where(np.isclose(type_fre2.values, i0))
    theta = type_fre2.index.values[iidx[0][0]]
    phi = type_fre2.columns.values[iidx[1][0]]
    tipical_th_ph_list.append((theta, phi))
tipical_th_ph_list


    

    
Out[22]:





[(0.143, 1.142), (0.0, 0.0), (1.571, 1.428)]

In [23]:

    

tipical_th_ph_list = [[0.143, 1.142],
                      [0.000, 0.000],
                      [1.571, 1.428]]

tipical_th_ph_use_idx = [0, 1, 2]
Table_t_range = np.array((9000, np.inf))
figsize = np.array((16, 9)) * 0.7
dpi = 200

tipical_th_ph_list_use = np.array(tipical_th_ph_list)[tipical_th_ph_use_idx]
label_list = np.arange(len(tipical_th_ph_list))[tipical_th_ph_use_idx]
spf_tb.show_table_theta_phi_psi_fft_list(tipical_th_ph_list_use, job_dir, label_list=label_list, figsize=figsize, dpi=dpi, 
                                         use_welch=False)
spf_tb.show_table_theta_phi_psi_fft_list(tipical_th_ph_list_use, job_dir, label_list=label_list, figsize=figsize, dpi=dpi, 
                                         use_welch=True)
spf_tb.show_table_result_list(tipical_th_ph_list_use, job_dir, label_list=label_list, Table_t_range=Table_t_range,
                              figsize=figsize, dpi=dpi)


    

    
Out[23]:





True






    













    













    

In [8]:

    

psi_lim_fct = 20
resampling_fct = 10
figsize = np.array((16, 9)) * 0.7
dpi = 200

# case_path_list = separate_fre_path_v2(type_fre2)
# color_list =     [np.array(plt.get_cmap('tab20', len(check_fre_list))(i0)[:3]).reshape((1, 3)) 
#                   for i0 in np.arange(type_fre2.values.max() + 1)]
for idx, psi_lim1 in enumerate(np.linspace(0, 2 * np.pi, psi_lim_fct * 16, endpoint=False)[::psi_lim_fct]):
    psi_lim = (psi_lim1, psi_lim1 + 2 * np.pi / (psi_lim_fct * 16))
    # create fig
    desc = '$\\psi \\in [%.3f \\pi, %.3f \\pi)$' % ((psi_lim[0] / np.pi), (psi_lim[1] / np.pi))
    fig, ax1 = plt.subplots(1, 1, figsize=figsize, dpi=dpi, subplot_kw=dict(polar=True))
    n_xticks = 32
    xticks = np.arange(n_xticks)
    ax1.set_xticks(xticks / n_xticks * 2 * np.pi)
    ax1.set_xticklabels(['$\dfrac{%d}{%d}2\pi$' % (i0, n_xticks) for i0 in xticks])
    ax1.set_yticklabels([])
    ax1.set_ylim(0, np.pi)
    fig.suptitle(desc)
    fig.tight_layout(rect=[0, 0, 1, 0.95])
    
    for i0 in np.arange(type_fre2.values.max() + 1):
        # get case_path
        iidx = np.where(np.isclose(type_fre2.values, i0))
        theta = type_fre2.index.values[iidx[0]]
        phi = type_fre2.columns.values[iidx[1]]
        theta_phi_list = np.vstack((theta, phi)).T
        case_path = []
        for theta, phi in theta_phi_list[:]:
            t_headle = 'th%5.3f_ph%5.3f_(.*?).pickle' % (theta, phi)
            filename = [filename for filename in os.listdir(job_dir)
                        if re.match(t_headle, filename) is not None][0]
            case_path.append(os.path.join(PWD, job_dir, filename))
            
        thandle = ''
        color = np.array(plt.get_cmap('tab20', int(type_fre2.values.max() + 1))(int(i0))[:3]).reshape((1, 3)) 
#         print(int(type_fre2.values.max()), i0, color)
        spf_tb.draw_phase_map_theta(case_path, color, psi_lim, axs=[ax1, ], thandle=thandle, 
                                    resampling=True, resampling_fct=resampling_fct)
    tdir = os.path.join(PWD, job_dir, 'phase_mape_fre')
    if not os.path.exists(tdir):
        os.makedirs(tdir)
    figname = os.path.join(tdir, '%04d.png' % (idx))
    fig.savefig(os.path.join(tdir, figname))
    print('save to %s' % figname)
    plt.close(fig)


    

    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0000.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0001.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0002.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0003.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0004.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0005.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0006.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0007.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0008.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0009.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0010.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0011.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0012.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0013.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0014.png






    






 
 










    










    






 
 










    










    






 
 










    




save to /home/zhangji/stokes_flow_master/head_Force/do_calculate_table/ecoC01B05_wt0.7_psi-0c/phase_mape_fre/0015.png

In [ ]: