Vertical derivative of Bouguer anomalies

I would like to compute the vertical derivative of Bouguer anomalies using data in both geographical coordinates and projected coordinates (UTM). How should I proceed?

grdfft -D?