code: def plotradtime(TR: xarray.Dataset, c1: Dict[str, Any], log: bool = False):
      
  
      for t in TR.time:  # for each time
          plotirrad(TR.sel(time=t), c1, log)
docs: make one plot per time for now.
  
      TR: 3-D array: time, wavelength, [transmittance, radiance]
  
      radiance is currently single-scatter solar