ezCon240108a Errors with ezConAstroMath = 2

[nathanielbutts]

ezCon240108a Errors with ezConAstroMath = 2

Howdy Ted,

Had a strange error pop up trying to use ezConAstroMath = 2. I didn't do any digging into the code yet as I found this on my lunch break, but wanted to make you aware.

Setup

System

ezcon is aliased to ~/bin/ezRA/ezRA/ezCon.py
using programRevision = ezCon240108a.py

OS: Linux Mint 21.2 x86_64
Host: OptiPlex 3050
Kernel: 5.15.0-91-generic
Shell: bash 5.1.16
DE: Cinnamon 5.8.4

Commands

Failed when using commands:

ezcon . -ezConAstroMath 2
ezcon . -ezConAstroMath 2 -ezConGalCrossingGLonCenterL -154 -149 6 -ezConGalCrossingGLonCenterL 35 40 6
ezcon sokyrad-1m_240126_00.txt -ezConAstroMath 2

Passes when using commands:

ezcon . -ezConGalCrossingGLonCenterL -154 -149 6 -ezConGalCrossingGLonCenterL 35 40 6
ezcon sokyrad-1m_240126_00.txt -ezConAstroMath 1

Output on Failure

sokyrad-1m_240130_00.txt  43 plotting ezCon382antRBTByFreqBinAvgRfi.png ================================

   createEzConOutEzb ===============
                         ezConAstroMath = 2
Traceback (most recent call last):
  File "/home/nathan/bin/ezRA/ezRA/ezCon.py", line 10345, in <module>
    main()
  File "/home/nathan/bin/ezRA/ezRA/ezCon.py", line 10226, in main
    createEzConOutEzb()                 # using ezConAstroMath, creates partial ezConOut,
  File "/home/nathan/bin/ezRA/ezRA/ezCon.py", line 4363, in createEzConOutEzb
    vStarEarthAroundSun = -float(pointingTelescope.radial_velocity_correction(kind='barycentric') / (1000. * u.m / u.s))   # into km/s
  File "/usr/lib/python3/dist-packages/astropy/coordinates/sky_coordinate.py", line 1942, in radial_velocity_correction
    gcrs_p, gcrs_v = location.get_gcrs_posvel(obstime)
  File "/usr/lib/python3/dist-packages/astropy/coordinates/earth.py", line 734, in get_gcrs_posvel
    return self._get_gcrs_posvel(obstime,
  File "/usr/lib/python3/dist-packages/astropy/coordinates/earth.py", line 708, in _get_gcrs_posvel
    pos = itrs_cart.transform(itrs_to_gcrs)
  File "/usr/lib/python3/dist-packages/astropy/coordinates/representation.py", line 1402, in transform
    p = erfa_ufunc.rxp(matrix, self.get_xyz(xyz_axis=-1))
  File "/usr/lib/python3/dist-packages/astropy/coordinates/representation.py", line 1354, in get_xyz
    return np.stack([self._x, self._y, self._z], axis=xyz_axis)
  File "/usr/lib/python3/dist-packages/astropy/units/quantity.py", line 1696, in __array_function__
    return super().__array_function__(function, types, args, kwargs)
  File "/home/nathan/.local/lib/python3.10/site-packages/numpy/core/shape_base.py", line 456, in stack
    return _nx.concatenate(expanded_arrays, axis=axis, out=out,
  File "/usr/lib/python3/dist-packages/astropy/units/quantity.py", line 1701, in __array_function__
    args, kwargs, unit, out = function_helper(*args, **kwargs)
TypeError: concatenate() got an unexpected keyword argument 'dtype'