import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import welch, hanning
plt.rcParams.update({'font.size': 22})

# ***** read u

data = np.genfromtxt("u_v_time_4nodes_hybrid.dat",comments="%")

u1=data[:,0]   
u2=data[:,1]   
u3=data[:,2]   
u4=data[:,3]   

u5=data[:,4]   
u6=data[:,4]   
u7=data[:,4]   
u8=data[:,4]  

# time step
dt= 0.0049
t_tot=dt*len(u1)
t = np.linspace(0,t_tot,len(u1))


#%%%%%%%%%%%%%%% plotting section %%%%%%%%%%%%%%%%%%%%%%%%%%
# plot v
fig1 = plt.figure("Figure 1")
plt.clf() #clear the figure
plt.plot(t,u2,'r--')
plt.xlabel('t')
plt.ylabel('u')
plt.savefig('u_time_hybrid_python.eps')