cfd_data = 'vort'
cfd_data = 'vort2'
# cfd_data = 'velo'
# cfd_data = 'velo2'
if cfd_data == 'vort':
loc,spr,skw,scl = vorticity(tsr,solidity)
param = np.array([loc,spr,skw,scl])
print param
import time
time.sleep(10)
elif cfd_data == 'velo':
men1,sdv1,rat1,wdt1,spr1,scl1,tsrn1,_ = velocity(tsr-0.1249,solidity)
men2,sdv2,rat2,wdt2,spr2,scl2,tsrn2,_ = velocity(tsr+0.1249,solidity)
if solidity >= 0.35:
men3,sdv3,rat3,wdt3,spr3,scl3,_,soln1 = velocity(tsr,solidity-0.1249)
men4,sdv4,rat4,wdt4,spr4,scl4,_,soln2 = velocity(tsr,solidity+0.1249)
elif solidity >=0.25:
men3,sdv3,rat3,wdt3,spr3,scl3,_,soln1 = velocity(tsr,solidity-0.049)
men4,sdv4,rat4,wdt4,spr4,scl4,_,soln2 = velocity(tsr,solidity+0.1249)
else:
men3,sdv3,rat3,wdt3,spr3,scl3,_,soln1 = velocity(tsr,solidity-0.049)
men4,sdv4,rat4,wdt4,spr4,scl4,_,soln2 = velocity(tsr,solidity+0.049)
if tsrn1 == tsrn2:
p = 0.
else:
p = (tsr-tsrn1)/(tsrn2-tsrn1)
if soln1 == soln2:
sdv3 = np.zeros((n,m))
sdv4 = np.zeros((n,m))
rat1 = np.zeros((n,m))
wdt1 = np.zeros((n,m))
spr1 = np.zeros((n,m))
spr2 = np.zeros((n,m))
spr3 = np.zeros((n,m))
spr4 = np.zeros((n,m))
scl1 = np.zeros((n,m))
scl2 = np.zeros((n,m))
scl3 = np.zeros((n,m))
for i in range(np.size(solidity)):
for j in range(np.size(tsr)):
men,sdv,rat,wdt,spr,scl,_,_ = velocity(X[i,j],Y[i,j])
men1[i,j] = men[0]
sdv1[i,j] = sdv[0]
sdv2[i,j] = sdv[1]
sdv3[i,j] = sdv[2]
sdv4[i,j] = sdv[3]
rat1[i,j] = rat[0]
wdt1[i,j] = wdt[0]
spr1[i,j] = spr[0]
spr2[i,j] = spr[1]
spr3[i,j] = spr[2]
spr4[i,j] = spr[3]
scl1[i,j] = scl[0]
scl2[i,j] = scl[1]
scl3[i,j] = scl[2]
