code: def _run_crefl(refl, mus, muv, phi, height, sensor_name, *coeffs):
      atm_vars_cls = (
          _VIIRSAtmosphereVariables
          if sensor_name.lower() == "viirs"
          else _MODISAtmosphereVariables
      )
      atm_vars = atm_vars_cls(mus, muv, phi, height, *coeffs)
      sphalb, rhoray, TtotraytH2O, tOG = atm_vars()
      return _correct_refl(refl, tOG, rhoray, TtotraytH2O, sphalb)
  
imports: from __future__ import annotations
  import logging
  from typing import Optional, Type, Union
  import dask.array
  import numpy
  import xarray
  from satpy.dataset.dataid import WavelengthRange
label_desc: Use image node_id 4 for calling a global function with example usage: _run_crefl(refl, mus, muv, phi, height, sensor_name) and returns: _correct_refl