For 2-dimensional data you can use
IMOdel IMOdel [# # [#]] [# # [#]]
where 2nd triplet of numbers specify the lower and upper limits, and the number of steps to use in the Y direction.
The method only uses a simple rectangle rule to calculate the integral.