Fig11.ipynb

Kernel: Python 3 (ipykernel)

import numpy as np
import h5py as h5
import matplotlib.pyplot as plt
from matplotlib import rcParams
from numpy import matrix

rcParams['font.family'] = 'serif'
rcParams['font.serif'] = ['Computer Modern Roman']
rcParams['text.usetex'] = True
rcParams['xtick.major.size']=7
rcParams['ytick.major.size']=7
rcParams['xtick.minor.size']=3.5
rcParams['ytick.minor.size']=3.5
rcParams['xtick.labelsize']=10
rcParams['ytick.labelsize']=10
rcParams['contour.negative_linestyle'] = 'solid'

F = h5.File('../data/C0093.h5','r')
C0093_D_uu    =F['D_uu'][()]
C0093_E_uu    =F['E_uu'][()]
C0093_P_uu    =F['P_uu'][()]
C0093_PI1_uu  =F['PI1_uu'][()]
C0093_T_uu    =F['T_uu'][()]
C0093_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/C0220.h5','r')
C0220_D_uu    =F['D_uu'][()]
C0220_E_uu    =F['E_uu'][()]
C0220_P_uu    =F['P_uu'][()]
C0220_PI1_uu  =F['PI1_uu'][()]
C0220_T_uu    =F['T_uu'][()]
C0220_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/C0500.h5','r')
C0500_D_uu    =F['D_uu'][()]
C0500_E_uu    =F['E_uu'][()]
C0500_P_uu    =F['P_uu'][()]
C0500_PI1_uu  =F['PI1_uu'][()]
C0500_T_uu    =F['T_uu'][()]
C0500_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/P0550.h5','r')
P0550_D_uu    =F['D_uu'][()]
P0550_E_uu    =F['E_uu'][()]
P0550_P_uu    =F['P_uu'][()]
P0550_PI1_uu  =F['PI1_uu'][()]
P0550_T_uu    =F['T_uu'][()]
P0550_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/P1000.h5','r')
P1000_D_uu    =F['D_uu'][()]
P1000_E_uu    =F['E_uu'][()]
P1000_P_uu    =F['P_uu'][()]
P1000_PI1_uu  =F['PI1_uu'][()]
P1000_T_uu    =F['T_uu'][()]
P1000_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/P2000.h5','r')
P2000_D_uu    =F['D_uu'][()]
P2000_E_uu    =F['E_uu'][()]
P2000_P_uu    =F['P_uu'][()]
P2000_PI1_uu  =F['PI1_uu'][()]
P2000_T_uu    =F['T_uu'][()]
P2000_yp_loc  =F['yp_loc'][()]
F.close()

F = h5.File('../data/P5200.h5','r')
P5200_D_uu    =F['D_uu'][()]
P5200_E_uu    =F['E_uu'][()]
P5200_P_uu    =F['P_uu'][()]
P5200_PI1_uu  =F['PI1_uu'][()]
P5200_T_uu    =F['T_uu'][()]
P5200_yp_loc  =F['yp_loc'][()]
F.close()

plt.plot(P0550_yp_loc, P0550_P_uu, 'C0-x')
plt.plot(P1000_yp_loc, P1000_P_uu, 'C1-x')
plt.plot(P2000_yp_loc, P2000_P_uu, 'C2-x')
plt.plot(P5200_yp_loc, P5200_P_uu, 'C3-x')
plt.plot(C0093_yp_loc, C0093_P_uu, 'C4-x')
plt.plot(C0220_yp_loc, C0220_P_uu, 'C5-x')
plt.plot(C0500_yp_loc, C0500_P_uu, 'C6-x')
plt.plot([-1,-2],[0,0],'k-x',label=r'$P^+$')

plt.plot(P0550_yp_loc, P0550_T_uu, 'C0-.')
plt.plot(P1000_yp_loc, P1000_T_uu, 'C1-.')
plt.plot(P2000_yp_loc, P2000_T_uu, 'C2-.')
plt.plot(P5200_yp_loc, P5200_T_uu, 'C3-.')
plt.plot(C0093_yp_loc, C0093_T_uu, 'C4-.')
plt.plot(C0220_yp_loc, C0220_T_uu, 'C5-.')
plt.plot(C0500_yp_loc, C0500_T_uu, 'C6-.')
plt.plot([-1,-2],[0,0],'k-.',label=r'$T^+$')

plt.plot(P0550_yp_loc, P0550_PI1_uu, 'C0--')
plt.plot(P1000_yp_loc, P1000_PI1_uu, 'C1--')
plt.plot(P2000_yp_loc, P2000_PI1_uu, 'C2--')
plt.plot(P5200_yp_loc, P5200_PI1_uu, 'C3--')
plt.plot(C0093_yp_loc, C0093_PI1_uu, 'C4--')
plt.plot(C0220_yp_loc, C0220_PI1_uu, 'C5--')
plt.plot(C0500_yp_loc, C0500_PI1_uu, 'C6--')
plt.plot([-1,-2],[0,0],'k--',label=r'$\Pi^+$')

plt.plot(P0550_yp_loc, P0550_D_uu, 'C0:')
plt.plot(P1000_yp_loc, P1000_D_uu, 'C1:')
plt.plot(P2000_yp_loc, P2000_D_uu, 'C2:')
plt.plot(P5200_yp_loc, P5200_D_uu, 'C3:')
plt.plot(C0093_yp_loc, C0093_D_uu, 'C4:')
plt.plot(C0220_yp_loc, C0220_D_uu, 'C5:')
plt.plot(C0500_yp_loc, C0500_D_uu, 'C6:')
plt.plot([-1,-2],[0,0],'k:',label=r'$D^+$')

plt.plot([-1,-2],[0,0],'k-',label=r'$-\epsilon^+$')
plt.plot(P0550_yp_loc, -P0550_E_uu, 'C0-', label=r'$Re_\tau = 550$ (P)')
plt.plot(P1000_yp_loc, -P1000_E_uu, 'C1-', label=r'$Re_\tau = 1000$ (P)')
plt.plot(P2000_yp_loc, -P2000_E_uu, 'C2-', label=r'$Re_\tau = 2000$ (P)')
plt.plot(P5200_yp_loc, -P5200_E_uu, 'C3-', label=r'$Re_\tau = 5200$ (P)')
plt.plot(C0093_yp_loc, -C0093_E_uu, 'C4-', label=r'$Re_\tau = 93$ (C)')
plt.plot(C0220_yp_loc, -C0220_E_uu, 'C5-', label=r'$Re_\tau = 220$ (C)')
plt.plot(C0500_yp_loc, -C0500_E_uu, 'C6-', label=r'$Re_\tau = 500$ (C)')
plt.plot([0,100],[0,0],'k-',linewidth=0.5)
plt.ylabel(r'$B^+$')
plt.xlabel(r'$y^+$')
plt.legend(frameon=False,ncol=1,loc='upper right', bbox_to_anchor=(1.35, 1.0))
plt.xlim([0,20])
plt.ylim([-0.6,0.6])
#plt.text(-4,0.1,'(e)')


plt.savefig('figure11.eps',dpi=300,bbox_inches='tight')
plt.show()
plt.close()