how to do terrain correction to gravity data in Matlab

I want to do terrain correction to a gravity data grid. I have grid files of topography and free air gravity. I got a Matlab code, named nagyprism.m to calculate the topographic correction due to a prism. But I am not clear about how to apply this to a grid data. The only inputs to the code are coordinates of a prism. Please explain this.

Answer 1

[Pierre845]

Not sure what you're after here .. you got a matlab file and don't know how to use it? well post the code for broader interest and potential answers.

My help is that you need two surfaces or grids of data: the topographic one and the gravity one; then you just need to do a wise substraction of the two, making sure you are using the same units in the substraction (either mass or gravity acceleration).

Hence, a 3D matrix will do the magic - one for the topography and one for the gravity; the two matrix need to have the same x,y coordinates if not you will need to do an interpolation to align coordinates.