#REDVOL = ['70']
#REDVOL = ['65']
#fs = 24
#plt.figure(figsize=(14,12))
fs = 10
plt.figure(figsize=(6,4))
plt.rc('text',usetex=True)
plt.rc('font', family='serif')
#for redvol in REDVOL :
for redvol in REDVOL :
color = iter(cm.brg(np.linspace(0.0,0.8, len(CONFIN))))
for confin in CONFIN :
c = next(color)
#(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA) = readMetrics(confin)
(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA, CURV) = readMetrics(redvol,confin)
v = round(REDVOL[0],2)
lam = LAMBDA[0]
# calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0
lamc = 1.0
if (abs(v - 1.00) < 1e-12) :
lamc = 1.000000000000000000000000
if (abs(v - 0.99) < 1e-12) :
lamc = 1.090275107548953613318768
if (abs(v - 0.98) < 1e-12) :
lamc = 1.133537870150389282270677
if (abs(v - 0.97) < 1e-12) :
lamc = 1.169545928031930431602916
if (abs(v - 0.96) < 1e-12) :
CAPNUM = ['100']
CAPNUM = ['1']
#fs = 24
#plt.figure(figsize=(14,12))
fs = 10
plt.figure(figsize=(6,4))
plt.rc('text',usetex=True)
plt.rc('font', family='serif')
#for redvol in REDVOL :
for capnum in CAPNUM :
color = iter(cm.brg(np.linspace(0.0,0.8, len(REDVOL))))
for redvol in REDVOL :
c = next(color)
#(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA) = readMetrics(confin)
(STR, EPS, VOLRAT, LAMBDA, AREA, VLME, DP) = readMetrics(redvol,capnum)
v = round(float(redvol)/100.0,2)
# calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0
lamc = 1.0
if (abs(v - 0.99) < 1e-12) :
lamc = 1.090275107548953613318768
if (abs(v - 0.98) < 1e-12) :
lamc = 1.133537870150389282270677
if (abs(v - 0.97) < 1e-12) :
lamc = 1.169545928031930431602916
if (abs(v - 0.96) < 1e-12) :
lamc = 1.202032009320693326333173
if (abs(v - 0.95) < 1e-12) :
lamc = 1.232435708492447912665580
CAPNUM = ['1']
#fs = 24
#plt.figure(figsize=(14,12))
fs = 10
plt.figure(figsize=(6,4))
plt.rc('text',usetex=True)
plt.rc('font', family='serif')
#for redvol in REDVOL :
for capnum in CAPNUM :
color = iter(cm.brg(np.linspace(0.0,0.8, len(CONFIN))))
for confin in CONFIN :
c = next(color)
#(STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, DP, ETA) = readMetrics(confin)
(STR, EPS, REDVOL, LAMBDA, AREA, VLME, DP) = readMetrics(confin,capnum)
# CE = CA
# DPE = DP
# for i in range(len(CE)) :
# CE[i] = CA[i]/(EPS[i]*EPS[i])
# #CE[i] = CA[i]/(delta*delta)
# DPE[i] = DP[i]*EPS[i]
# #DPE[i] = DP[i]*delta
# CA = CE
# DP = DPE
REDVOL2 = REDVOL
for i in range(len(REDVOL)) :
v = REDVOL[i]
v = float(v)/100
REDVOL2[i] = v
from matplotlib.pyplot import cm from scipy.interpolate import splrep, splev, interp1d from math import log10, floor MAXRAD = ['90', '92', '94', '96'] MAXRAD = ['90', '99'] #MAXRAD = ['80','82','84','86','88','90', '92', '94', '96', '98'] MAXRAD = ['90','91','92','93','94','95','96','97','98','99'] #MAXRAD = ['99'] color = iter(cm.brg(np.linspace(0.0,0.8, len(MAXRAD)))) fs = 24 plt.figure(figsize=(14,12)) for maxrad in MAXRAD : c = next(color) (STR, CA, EPS, REDVOL, LAMBDA, AREA, VLME, LENGTH, DP, ETA) = readMetrics(maxrad) v = round(REDVOL[0],2) lam = LAMBDA[0] # calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0 if (abs(v - 0.95) < 1e-12) : lamc = 1.2324 if (abs(v - 0.90) < 1e-12) : lamc = 1.3712 if (abs(v - 0.85) < 1e-12) : lamc = 1.5050 if (abs(v - 0.80) < 1e-12) : lamc = 1.6437 if (abs(v - 0.75) < 1e-12) : lamc = 1.7925
maxL = 0 nspline = 1000 vecRedVol = [] vecExponent = [] vecCoeff = [] vecRSquared = [] for redvol in REDVOL : c = next(color) x = [] y = [] L = [] for maxrad in MAXRAD : (STR, Ca, eps, area, vlme, Dp, l) = readMetrics(redvol, maxrad, capnum) lam = sqrt(area/(4.0*pi)) v = round(float(redvol)/100,2) # calculate lamc from the cubic equuation: 2 v x^3 - 3 x + 1 = 0 if (abs(v - 0.95) <= 0.01) : lamc = 1.2324 if (abs(v - 0.90) <= 0.01) : lamc = 1.3712 if (abs(v - 0.85) <= 0.01) : lamc = 1.5050 if (abs(v - 0.80) <= 0.01) : lamc = 1.6437 if (abs(v - 0.75) <= 0.01) : lamc = 1.7925
LAVG = [] for redvol in REDVOL : if (redvol == '80') : MAXRAD = ['80','82','83','84','85','86','87','88','89','90','91','92','93','94','95','96','97','98','99','100','101','102','103'] if (redvol == '75') : MAXRAD = ['80','82','83','84','85','86','87','88','89','90','91','92','93','94','95','96','97','98','99','100','101','102','103','104','105','106','107'] if (redvol == '70') : MAXRAD = ['80','82','83','84','85','86','87','88','89','90','91','92','93','94','95','96','97','98','99','100','101','102','103','104','105','106','107','108','109','110','111'] c = next(color) x = [] y = [] L = [] for maxrad in MAXRAD : (STR, Ca, eps, area, vlme, Dp, eta, l, gammax) = readMetrics(redvol, maxrad, capnum) lam = sqrt(area/(4.0*pi)) v = round(float(redvol)/100,2) # calculate lamc from the cubic equuation: 2 v x^3 - 3 x^2 + 1 = 0 if (abs(v - 0.95) <= 0.01) : lamc = 1.2324 if (abs(v - 0.90) <= 0.01) : lamc = 1.3712 if (abs(v - 0.85) <= 0.01) : lamc = 1.5050 if (abs(v - 0.80) <= 0.01) : lamc = 1.6437 if (abs(v - 0.75) <= 0.01) : lamc = 1.7925