code: def get_xy_steps(bbox, h_dim):
      r
      x_range, y_range = get_xy_range(bbox)
  
      x_steps = np.ceil(x_range / h_dim)
      y_steps = np.ceil(y_range / h_dim)
  
      return int(x_steps), int(y_steps)
docs: r"""Return meshgrid spacing based on bounding box.
  
      bbox: dictionary
          Dictionary containing coordinates for corners of study area.
      h_dim: integer
          Horizontal resolution in meters.
  
      Returns
      -------
      x_steps, (X, ) ndarray
          Number of grids in x dimension.
      y_steps: (Y, ) ndarray
          Number of grids in y dimension.