code: def sgli_pol_file(tmp_path_factory):
      """Create a POL stub file."""
      filename = tmp_path_factory.mktemp("data") / "test_pol_file.h5"
      with h5py.File(filename, "w") as h5f:
          global_attributes = h5f.create_group("Global_attributes")
          global_attributes.attrs["Scene_start_time"] = np.array(
              [START_TIME.strftime("%Y%m%d %H:%M:%S.%f")[:-3]],
              dtype="|S21",
          )
          global_attributes.attrs["Scene_end_time"] = np.array(
              [END_TIME.strftime("%Y%m%d %H:%M:%S.%f")[:-3]],
              dtype="|S21",
          )
  
          image_data = h5f.create_group("Image_data")
          image_data.attrs["Number_of_lines"] = 1854
          image_data.attrs["Number_of_pixels"] = 1250
  
          p1_0 = image_data.create_dataset(
              "Lt_P1_0", data=FULL_KM_ARRAY, chunks=(116, 157)
          )
          p1_0.attrs["Slope_reflectance"] = np.array(
              [5e-05], dtype=np.float32
          )
          p1_0.attrs["Offset_reflectance"] = np.array(
              [0.0], dtype=np.float32
          )
          p1_0.attrs["Slope"] = np.array([0.02], dtype=np.float32)
          p1_0.attrs["Offset"] = np.array([-25], dtype=np.float32)
          p1_0.attrs["Mask"] = np.array([16383], dtype=np.uint16)
          p1_0.attrs["Bit00(LSB)-13"] = np.array(
              [
                  b"Digital Number\n16383 : Missing value\n16382 : Saturation value"
              ],
              dtype="|S61",
          )
  
          p1_m60 = image_data.create_dataset(
              "Lt_P1_m60", data=FULL_KM_ARRAY, chunks=(116, 157)
          )
          p1_m60.attrs["Slope_reflectance"] = np.array(
              [5e-05], dtype=np.float32
          )
          p1_m60.attrs["Offset_reflectance"] = np.array(
              [-60.0], dtype=np.float32
          )
          p1_m60.attrs["Slope"] = np.array([0.0012], dtype=np.float32)
          p1_m60.attrs["Offset"] = np.array([-1.65], dtype=np.float32)
          p1_m60.attrs["Mask"] = np.array([16383], dtype=np.uint16)
          p1_m60.attrs["Bit00(LSB)-13"] = np.array(
              [
                  b"Digital Number\n16383 : Missing value\n16382 : Saturation value"
              ],
              dtype="|S61",
          )
  
          p1_60 = image_data.create_dataset(
              "Lt_P1_60", data=FULL_KM_ARRAY, chunks=(116, 157)
          )
          p1_60.attrs["Slope_reflectance"] = np.array(
              [5e-05], dtype=np.float32
          )
          p1_60.attrs["Offset_reflectance"] = np.array(
              [60.0], dtype=np.float32
          )
          p1_60.attrs["Slope"] = np.array([0.0012], dtype=np.float32)
          p1_60.attrs["Offset"] = np.array([-1.65], dtype=np.float32)
          p1_60.attrs["Mask"] = np.array([16383], dtype=np.uint16)
          p1_60.attrs["Bit00(LSB)-13"] = np.array(
              [
                  b"Digital Number\n16383 : Missing value\n16382 : Saturation value"
              ],
              dtype="|S61",
          )
  
          geometry_data = h5f.create_group("Geometry_data")
          longitude = geometry_data.create_dataset(
              "Longitude",
              data=FULL_KM_ARRAY.astype(np.float32),
              chunks=(47, 63),
          )
          longitude.attrs["Resampling_interval"] = 1
          latitude = geometry_data.create_dataset(
              "Latitude",
              data=FULL_KM_ARRAY.astype(np.float32),
              chunks=(47, 63),
          )
          latitude.attrs["Resampling_interval"] = 1
  
          return filename
  
imports: import sys
  from datetime import datetime, timedelta
  import dask
  import h5py
  import numpy
  import pytest
  from satpy.readers.sgli_l1b import HDF5SGLI
label_desc: Use image node_id 3 for calling a global function with example usage: sgli_pol_file(tmp_path_factory) and returns: filename