# colocar en el vecto n_magnus los distintos 'angulos que se disponen de coeficientes

[rows,cols] = np.shape(globals.data)

n = cols - 8

n_magnus = globals.data[0:rows,8:8+n]

#Convert deg to rad

magnus_ang = globals.data_raw[rows,8:8+n]*np.pi/180

tot = rows*n

M_0 = globals.data[:,0]

# Mach para doble interpolaci'on

globals.Mpa = np.zeros(tot)

for i in range(rows):

for k in range(n):

globals.Mpa[i*n+k] = M_0[i]

print('globals.Mpa')

print(globals.Mpa)

print(np.shape(globals.Mpa))

# angulos de magnus

globals.ang = []

for ii in range(rows):

globals.ang = np.append(globals.ang,magnus_ang)

print('globals.ang')

print(globals.ang)

print(np.shape(globals.ang))

# Cnpa

globals.Cnpa_proce = np.zeros(tot)

for i in range(rows):

for k in range(n):

globals.Cnpa_proce[i*n+k] = globals.data[i,k+8]

print('globals.Cnpa_proce')

print(globals.Cnpa_proce)

print(np.shape(globals.Cnpa_proce))

#interp2d option

#globals.interp = interpolate.interp2d(globals.Mpa, globals.ang, globals.Cnpa_proce)

#RectBivariateSpline option

n = np.shape(globals.Mpa[0::4])

m = np.shape(globals.ang[0:4])

globals.interp = interpolate.RectBivariateSpline(globals.Mpa[0::4], globals.ang[0:4], np.reshape(np.asarray(globals.Cnpa_proce),(n[0],m[0])))