function grid2d,xmin,xmax,xnum,ymin,ymax,ynum
  x=linterp(0,xmin,xnum-1,xmax,dindgen(xnum))
  x=rebin(x,xnum,ynum)
  y=linterp(0,ymin,ynum-1,ymax,dindgen(ynum))
  y=rebin(transpose(y),xnum,ynum)
  return,[[x[*]],[y[*]]]
end