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]