import gpplot as gplt
import numpy as np
import matplotlib.pyplot as plt

# Set fontsizes
plt.rcParams.update({'axes.titlesize':16,'axes.labelsize':14,'xtick.labelsize':12,'ytick.labelsize':12})

# Import data and extract fields into dictionary
data = gplt.dataFromFile('RayleighTaylor0128_3_002975.h5')

plt.figure()

# 2 side-by-side subplots
plt.subplot(1,2,1)
plt.plot(data['dens'][:,64],data['y'],'-b',linewidth=2)
plt.xlabel(r'$\rho$')
plt.ylabel(r'$y$')
plt.title(r'Density at $y=1/8$')
plt.grid('both')

plt.subplot(1,2,2)
plt.plot(data['x'],data['dens'][300,:],'-b',linewidth=2)
plt.xlabel(r'$x$')
plt.ylabel(r'$\rho$')
plt.title(r'Density at $x=0.586$')
plt.grid('both')

# Show plot on screen
plt.show()