I was trying to plot a field with a clip file so that it only plots where the data exists which is an oblique rectangular shape. However, I could not find the equivalent of psclip in pygmt. What would be the simplest way to plot a grd file inside a polygon?