code: def _get_cos_sza(utc_time, lons, lats):
      cos_sza = da.map_blocks(
          _cos_zen_ndarray,
          lons,
          lats,
          utc_time,
          meta=np.array((), dtype=lons.dtype),
          dtype=lons.dtype,
          chunks=lons.chunks,
      )
      return cos_sza
  
imports: from __future__ import annotations
  import hashlib
  import os
  import shutil
  import warnings
  from datetime import datetime
  from functools import update_wrapper
  from glob import glob
  from typing import Any, Callable, Optional, Union
  import dask
  import numpy
  import xarray
  from dask import array
  from pyorbital.astronomy import cos_zen
  from pyorbital.astronomy import get_alt_az
  from pyorbital.orbital import get_observer_look
  from pyresample.geometry import AreaDefinition, StackedAreaDefinition, SwathDefinition
  import satpy
  from satpy.utils import PerformanceWarning, get_satpos, ignore_invalid_float_warnings
label_desc: Use image node_id 16 for calling a global function with example usage: _get_cos_sza(utc_time, lons, lats) and returns: cos_sza