Plot uncertainty ellipse location


maybe this is a more math problem than just GMT’s, here is my issue:
I have a single seismic station, and I know the location with estimate of a source’s location from back-azimuth and distance both in degrees. I have a estimate of the location as well in Long/Lat. On both Back-azimuth and epicentral distance, I have uncertainties (i.e., +/- a value).
What I’d like to plot on the map is the uncertainty ellipse.

The plot -SE|e is not what I need, as this does not solve what I am seeking for. So I can compute outside GMT the ellipse trace, but my question is, do we have in some ways a purely GMT’s solution to such task?