code: def test_load_surf_data_file_glob(tmp_path):
      data2D = np.ones((20, 3))
      fnames = []
      for f in range(3):
          fd, filename = tempfile.mkstemp(
              prefix=f"glob_{f}_", suffix=".gii", dir=str(tmp_path)
          )
          os.close(fd)
          fnames.append(filename)
          data2D[:, f] *= f
          darray = gifti.GiftiDataArray(
              data=data2D[:, f], datatype="NIFTI_TYPE_FLOAT32"
          )
          gii = gifti.GiftiImage(darrays=[darray])
          nb.save(gii, fnames[f])
  
      assert_array_equal(
          load_surf_data(
              os.path.join(os.path.dirname(fnames[0]), "glob*.gii")
          ),
          data2D,
      )
  
      # make one more gii file that has more than one dimension
      fd, filename = tempfile.mkstemp(
          prefix="glob_3_", suffix=".gii", dir=str(tmp_path)
      )
      os.close(fd)
      fnames.append(filename)
      darray1 = gifti.GiftiDataArray(
          data=np.ones((20,)), datatype="NIFTI_TYPE_FLOAT32"
      )
      gii = gifti.GiftiImage(darrays=[darray1, darray1, darray1])
      nb.save(gii, fnames[-1])
  
      data2D = np.concatenate((data2D, np.ones((20, 3))), axis=1)
      assert_array_equal(
          load_surf_data(
              os.path.join(os.path.dirname(fnames[0]), "glob*.gii")
          ),
          data2D,
      )
  
      # make one more gii file that has a different shape in axis=0
      fd, filename = tempfile.mkstemp(
          prefix="glob_4_", suffix=".gii", dir=str(tmp_path)
      )
      os.close(fd)
      fnames.append(filename)
      darray = gifti.GiftiDataArray(
          data=np.ones((15, 1)), datatype="NIFTI_TYPE_FLOAT32"
      )
      gii = gifti.GiftiImage(darrays=[darray])
      nb.save(gii, fnames[-1])
  
      with pytest.raises(
          ValueError,
          match="files must contain data with the same shape",
      ):
          load_surf_data(
              os.path.join(os.path.dirname(fnames[0]), "*.gii")
          )
      for f in fnames:
          os.remove(f)
  
imports: import os
  import tempfile
  import warnings
  import nibabel
  import numpy
  import pytest
  from nibabel import gifti
  from numpy.testing import assert_array_almost_equal, assert_array_equal
  from scipy.spatial import Delaunay
  from scipy.stats import pearsonr
  from nilearn import datasets, image
  from nilearn._utils import data_gen
  from nilearn.image import resampling
  from nilearn.surface import Mesh, Surface, load_surf_data, load_surf_mesh, surface
  from nilearn.surface.surface import _gifti_img_to_mesh, _load_surf_files_gifti_gzip
  from nilearn.surface.tests._testing import flat_mesh, generate_surf, z_const_img
label_desc: Use image node_id 16 for calling a global function with example usage: test_load_surf_data_file_glob(tmp_path) without return types