Skip to content

Instantly share code, notes, and snippets.

@keflavich
Last active June 5, 2020 12:47
Show Gist options
  • Save keflavich/7d28b965fcf46ccc292871dd80a0551c to your computer and use it in GitHub Desktop.
Save keflavich/7d28b965fcf46ccc292871dd80a0551c to your computer and use it in GitHub Desktop.
Moment 0
[########################################] | 100% Completed | 9min 34.0s
[########################################] | 100% Completed | 9min 34.1s
[########################################] | 100% Completed | 9min 34.2s
[########################################] | 100% Completed | 9min 34.2s
[########################################] | 100% Completed | 9min 34.3s
[########################################] | 100% Completed | 9min 34.4s
[########################################] | 100% Completed | 9min 34.4s
/orange/adamginsburg/repos/radio-beam/build/lib.linux-x86_64-3.7/radio_beam/multiple_beams.py:261: UserWarning: Do not use the average beam for convolution! Use the smallest common beam from `Beams.common_beam()`.
warnings.warn("Do not use the average beam for convolution! Use the"
WARNING: BeamAverageWarning: Arithmetic beam averaging is being performed. This is not a mathematically robust operation, but is being permitted because the beams differ by <1 [spectral_cube.base_class]
[ ] | 0% Completed | 4.3s
[ ] | 0% Completed | 4.4s
[ ] | 0% Completed | 4.5s
[ ] | 0% Completed | 4.6s
[ ] | 0% Completed | 4.6s
[ ] | 0% Completed | 4.7s
[ ] | 0% Completed | 4.7s
Traceback (most recent call last):
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 566, in _do_parse
return cls._parse_unit(s, detailed_exception=False)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 487, in _parse_unit
raise ValueError()
ValueError
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 569, in _do_parse
return cls._parser.parse(s, lexer=cls._lexer, debug=debug)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/extern/ply/yacc.py", line 331, in parse
return self.parseopt_notrack(input, lexer, debug, tracking, tokenfunc)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/extern/ply/yacc.py", line 1199, in parseopt_notrack
tok = call_errorfunc(self.errorfunc, errtoken, self)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/extern/ply/yacc.py", line 193, in call_errorfunc
r = errorfunc(token)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 443, in p_error
raise ValueError()
ValueError
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/core.py", line 1854, in __call__
return f.parse(s)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 551, in parse
result = cls._do_parse(s, debug=debug)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/format/generic.py", line 574, in _do_parse
raise ValueError(f"Syntax error parsing unit '{s}'")
ValueError: Syntax error parsing unit 'm/s'
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/orange/adamginsburg/ALMA_IMF/reduction/analysis/line_quicklooks.py", line 93, in <module>
mvcube = mcube.with_spectral_unit(u.km/u.s, velocity_convention='radio')
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/spectral_cube.py", line 1609, in moment0
return self.moment(axis=axis, order=0, how=how)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/base_class.py", line 657, in newfunc
result = function(*args, **kwargs)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/dask_spectral_cube.py", line 906, in moment
out = self._compute(out)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/dask_spectral_cube.py", line 274, in _compute
return array.compute(**self._scheduler_kwargs)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/base.py", line 166, in compute
(result,) = compute(self, traverse=False, **kwargs)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/base.py", line 437, in compute
results = schedule(dsk, keys, **kwargs)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/threaded.py", line 84, in get
**kwargs
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/local.py", line 486, in get_async
raise_exception(exc, tb)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/local.py", line 316, in reraise
raise exc
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/local.py", line 222, in execute_task
result = _execute_task(task, data)
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/core.py", line 121, in _execute_task
return func(*(_execute_task(a, cache) for a in args))
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/dask/array/core.py", line 103, in getter
c = a[b]
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/dask_spectral_cube.py", line 199, in __getitem__
wcs_tolerance=self._cube._wcs_tolerance)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/masks.py", line 238, in _filled
ex = self.exclude(data=data, wcs=wcs, view=view, **kwargs)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/masks.py", line 159, in exclude
self._validate_wcs(data, wcs, **kwargs)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/masks.py", line 399, in _validate_wcs
self._mask1._validate_wcs(new_data=new_data, new_wcs=new_wcs, **kwargs)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/masks.py", line 400, in _validate_wcs
self._mask2._validate_wcs(new_data=new_data, new_wcs=new_wcs, **kwargs)
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/masks.py", line 534, in _validate_wcs
**kwargs):
File "/orange/adamginsburg/repos/spectral-cube/spectral_cube/wcs_utils.py", line 387, in check_equality
u1 = u.Unit(c1[1])
File "/orange/adamginsburg/miniconda3/lib/python3.7/site-packages/astropy/units/core.py", line 1874, in __call__
raise ValueError(msg)
ValueError: 'm/s' did not parse as unit: Syntax error parsing unit 'm/s' If this is meant to be a custom unit, define it with 'u.def_unit'. To have it recognized inside a file reader or other code, enable it with 'u.add_enabled_units'. For details, see http://docs.astropy.org/en/latest/units/combining_and_defining.html
@keflavich
Copy link
Author

This is coming from one of these lines: (between the print statements):

                    pl.clf()
                    dt(); print("Moment 2")
                    mom2 = mcube.with_spectral_unit(u.km/u.s, velocity_convention='radio').linewidth_fwhm()
                    mom2.write('collapse/moment2/{0}'.format(fn.replace(suffix,"_mom2fwhm_kms.fits")),
                               overwrite=True)
                    mom2.quicklook('collapse/moment2/pngs/{0}'.format(fn.replace(suffix,"_mom2fwhm_kms.png")))

                    pl.clf()
                    dt(); print("Spatial max (peak spectrum)")

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment