code: def _expandGLQ(self, zeros, lmax, lmax_calc):
          
          if self.normalization == '4pi':
              norm = 1
          elif self.normalization == 'schmidt':
              norm = 2
          elif self.normalization == 'unnorm':
              norm = 3
          elif self.normalization == 'ortho':
              norm = 4
          else:
              raise ValueError(
                  "Normalization must be '4pi', 'ortho', 'schmidt', or " +
                  "'unnorm'. Input value was {:s}"
                  .format(repr(self.normalization)))
  
          if zeros is None:
              zeros, weights = _shtools.SHGLQ(self.lmax)
  
          data = _shtools.MakeGridGLQC(self.coeffs, zeros, norm=norm,
                                       csphase=self.csphase, lmax=lmax,
                                       lmax_calc=lmax_calc)
          gridout = SHGrid.from_array(data, grid='GLQ', copy=False)
          return gridout
docs: Evaluate the coefficients on a Gauss-Legendre quadrature grid.