fig = pygmt.Figure()
pygmt.makecpt(cmap="roma", series=[-30,60], reverse=True)
fig.grdimage(grid=vel, cmap=True)
fig.basemap(frame="a")
fig.colorbar(frame="a5")
fig.show()