docs: def headgrid2(self, x1, x2, nx, y1, y2, ny, layers=None, printrow=False):
          """Grid of heads
          
          Parameters
          ----------
          x1, x2, nx : 
              x values are generated as linspace(x1, x2, nx)
          y1, y2, ny : 
              y values are generated as linspace(y1, y2, ny)
          layers : integer, list or array, optional
              layers for which grid is returned
          printrow : boolean, optional
              prints dot to screen for each row of grid if set to `True`
          
          Returns
          -------
          h : array size `nlayers, ny, nx`
          
          See also
          --------
          
          :func:`~timml.model.Model.headgrid`
          
          """
          
          xg, yg = np.linspace(x1, x2, nx), np.linspace(y1, y2, ny)
          return self.headgrid(xg, yg, layers=layers, printrow=printrow)