code: def reduce_mda(mda, max_size=100):
      """Recursively remove arrays with more than `max_size` elements from the given metadata dictionary."""
      reduced = {}
      for key, val in mda.items():
          if isinstance(val, dict):
              reduced[key] = reduce_mda(val, max_size)
          elif not (
              isinstance(val, np.ndarray) and val.size > max_size
          ):
              reduced[key] = val
      return reduced
  
imports: from __future__ import annotations
  import bz2
  import logging
  import os
  import shutil
  import tempfile
  import warnings
  from contextlib import closing, contextmanager
  from io import BytesIO
  from shutil import which
  from subprocess import PIPE, Popen
  import numpy
  import pyproj
  import xarray
  from pyresample.geometry import AreaDefinition
  from satpy import config
  from satpy.readers import FSFile
  from satpy.utils import get_legacy_chunk_size
label_desc: Use image node_id 20 for calling a global function with example usage: reduce_mda(mda, max_size) and returns: reduced