# see also README.pdf for instruction
import h5py
import numpy as np
import matplotlib.pyplot as plt

f = h5py.File("R2D2_difrot_ana.h5",'r')

# choose case from 'low', 'middle', and 'high'
datacase = 'low'

# define 2D geometry
RR, TH = np.meshgrid(f[datacase+'/rr'][...],f[datacase+'/th'][...],indexing='ij')
X, Y = RR*cos(TH), RR*sin(TH)

plt.pcolormesh(Y,X,f[datacase+'/vph_m'][:,:,200])

f.close()