Dear GMT developers and users,
I have XY locations corresponding to seismicity of earthquakes. I want to compute the total amount of seismic moment and earthquake rate by firstly clustering/grouping the events into a certain grid, say every 1deg, as shown by the figure attached. The outputs I expect are a number of grids, with each grid listing the number of earthquakes that fall within the grid, including their XY locations and magnitudes. Is it possible to do that in GMT?