else:
q = (solidity-soln1)/(soln2-soln1)
param = np.array([men1,sdv1,rat1,wdt1,spr1,scl1,men2,sdv2,rat2,wdt2,spr2,scl2,men3,sdv3,rat3,wdt3,spr3,scl3,men4,sdv4,rat4,wdt4,spr4,scl4,p,q])
elif cfd_data == 'vort2':
loc1,loc2,loc3,spr1,spr2,skw1,skw2,scl1,scl2,scl3 = vorticity2(tsr,solidity)
param = np.array([loc1,loc2,loc3,spr1,spr2,skw1,skw2,scl1,scl2,scl3])
# print param
# import time
#
# time.sleep(10)
elif cfd_data == 'velo2':
spr1,pow1,pow2,spr2,skw,scl1,scl2,scl3 = velocity2(tsr,solidity)
param = np.array([spr1,pow1,pow2,spr2,skw,scl1,scl2,scl3])
vel = velocity_field(xt,yt,xt + x0,yt + y0,velf,dia,tsr,solidity,cfd_data,param)
print '\nNormalized velocity at (',x0,',',y0,') from the turbine =',vel,'\n' # output velocity (normalized by free stream wind speed)
## Plotting
fs = 25 # 18# font size for plots
# Option to plot velocity profiles
vel_slice = True
vel_slice = False # comment this out if desired on
sdv1 = np.zeros_like(tsr)
sdv2 = np.zeros_like(tsr)
sdv3 = np.zeros_like(tsr)
sdv4 = np.zeros_like(tsr)
rat1 = np.zeros_like(tsr)
wdt1 = np.zeros_like(tsr)
spr1 = np.zeros_like(tsr)
spr2 = np.zeros_like(tsr)
spr3 = np.zeros_like(tsr)
spr4 = np.zeros_like(tsr)
scl1 = np.zeros_like(tsr)
scl2 = np.zeros_like(tsr)
scl3 = np.zeros_like(tsr)
for i in range(np.size(tsr)):
men,sdv,rat,wdt,spr,scl = velocity2(tsr[i],solidity[k])
men1[i] = men[0]
sdv1[i] = sdv[0]
sdv2[i] = sdv[1]
sdv3[i] = sdv[2]
sdv4[i] = sdv[3]
rat1[i] = rat[0]
wdt1[i] = wdt[0]
spr1[i] = spr[0]
spr2[i] = spr[1]
spr3[i] = spr[2]
spr4[i] = spr[3]
scl1[i] = scl[0]
scl2[i] = scl[1]
scl3[i] = scl[2]