1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34
| #
# Calcul des ratios EJMA pour
#
vx=1.82*(q/4)/(Dm*t*(Db/Dm)**0.5)
vy=2*(q/4)/w
print("Valeur vx: \n", vx)
print("Valeur vy: \n", vy)
#
#
# chargement du CSV
#
path="cp_ejma.csv"
file=open(path,"r")
test=csv.reader(file,quoting=csv.QUOTE_NONNUMERIC)
file.close()
#
#Conversion du CSV en tableau
data2 = [data for data in test]
data_array = np.asarray(data2, dtype = float)
#
# attribution d'une colonne du tableau à une variable (x,y et z)
x=data_array[:,0]
y=data_array[:,1]
z=data_array[:,2]
################################################
# Définition du maillage pour l'interpolation
xi = np.linspace(0.2,4.1,200)
yi = np.linspace(0.00,1.1,200)
# grid the data.
#En fonction des entrées utilisateurs, la fonction renvoi les coeficients
# Cp, Cf ou Cd
C = griddata((x,y),z,(vx,vy), method='cubic')
print ("Valeur C: \n", C) |
Partager