26171438601

Committed 20 May 2026 02:07AM UTC coverage: 90.571% (+0.01%) from 90.557%

Build # 26171438601

Build Type

push

github

Committed by

jameswilburlewis

Commit Message

Fix crash when a spectrogram is plotted with Y axis interpolation is enabled and Y scaling is linear.
Add regression test for y_interp bug.

Coverage Stats

0 of 4 new or added lines in 2 files covered. (0.0%)

17831 existing lines in 383 files now uncovered.

44379 of 48999 relevant lines covered (90.57%)

1.43 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

78.95

/pyspedas/tplot_tools/tplot_math/divide.py

# Copyright 2020 Regents of the University of Colorado. All Rights Reserved.
# Released under the MIT license.
# This software was developed at the University of Colorado's Laboratory for Atmospheric and Space Physics.
# Verify current version before use at: https://github.com/MAVENSDC/Pytplot

import pyspedas
from pyspedas.tplot_tools import store_data, tinterp
import numpy as np
import copy
import logging

def divide(tvar1,tvar2,newname=None):
    """
    Divides two tplot variables.  Will interpolate if the two are not on the same time cadence.

    Parameters
    ----------
        tvar1 : str
            Name of first tplot variable.
        tvar2 : int/float
            Name of second tplot variable
        newname : str
            Name of new tvar for divided data.  If not set, then the data in tvar1 is replaced.

    Returns
    -------
        None

    Examples
    --------

        >>> pyspedas.store_data('a', data={'x':[0,4,8,12,16], 'y':[1,2,3,4,5]})
        >>> pyspedas.store_data('c', data={'x':[0,4,8,12,16,19,21], 'y':[1,4,1,7,1,9,1]})
        >>> pyspedas.divide('a','c','a_over_c')
        """

    # interpolate tvars
    tv2 = tinterp(tvar1, tvar2)
    # separate and divide data
    data1 = pyspedas.tplot_tools.data_quants[tvar1].values
    data2 = pyspedas.tplot_tools.data_quants[tv2].values
    data = data1 / data2
    # store divided data
    if newname is None:
        pyspedas.tplot_tools.data_quants[tvar1].values = data
        return tvar1
    if 'spec_bins' in pyspedas.tplot_tools.data_quants[tvar1].coords:
        store_data(newname, data={'x': pyspedas.tplot_tools.data_quants[tvar1].coords['time'].values, 'y': data,
                                           'v': pyspedas.tplot_tools.data_quants[tvar1].coords['spec_bins'].values})
        pyspedas.tplot_tools.data_quants[newname].attrs = copy.deepcopy(pyspedas.tplot_tools.data_quants[tvar1].attrs)
    else:
       store_data(newname, data={'x':pyspedas.tplot_tools.data_quants[tvar1].coords['time'].values, 'y': data})
       pyspedas.tplot_tools.data_quants[newname].attrs = copy.deepcopy(pyspedas.tplot_tools.data_quants[tvar1].attrs)

    return newname

1	# Copyright 2020 Regents of the University of Colorado. All Rights Reserved.
2	# Released under the MIT license.
3	# This software was developed at the University of Colorado's Laboratory for Atmospheric and Space Physics.
4	# Verify current version before use at: https://github.com/MAVENSDC/Pytplot
5
6	import pyspedas	2✔
7	from pyspedas.tplot_tools import store_data, tinterp	2✔
8	import numpy as np	2✔
9	import copy	2✔
10	import logging	2✔
11
12	def divide(tvar1,tvar2,newname=None):	2✔
13	"""
14	Divides two tplot variables. Will interpolate if the two are not on the same time cadence.
15
16	Parameters
17	----------
18	tvar1 : str
19	Name of first tplot variable.
20	tvar2 : int/float
21	Name of second tplot variable
22	newname : str
23	Name of new tvar for divided data. If not set, then the data in tvar1 is replaced.
24
25	Returns
26	-------
27	None
28
29	Examples
30	--------
31
32	>>> pyspedas.store_data('a', data={'x':[0,4,8,12,16], 'y':[1,2,3,4,5]})
33	>>> pyspedas.store_data('c', data={'x':[0,4,8,12,16,19,21], 'y':[1,4,1,7,1,9,1]})
34	>>> pyspedas.divide('a','c','a_over_c')
35	"""
36
37	# interpolate tvars
UNCOV 38	tv2 = tinterp(tvar1, tvar2)	1✔
39	# separate and divide data
UNCOV 40	data1 = pyspedas.tplot_tools.data_quants[tvar1].values	1✔
UNCOV 41	data2 = pyspedas.tplot_tools.data_quants[tv2].values	1✔
UNCOV 42	data = data1 / data2	1✔
43	# store divided data
UNCOV 44	if newname is None:	1✔
45	pyspedas.tplot_tools.data_quants[tvar1].values = data	×
46	return tvar1	×
UNCOV 47	if 'spec_bins' in pyspedas.tplot_tools.data_quants[tvar1].coords:	1✔
48	store_data(newname, data={'x': pyspedas.tplot_tools.data_quants[tvar1].coords['time'].values, 'y': data,	×
49	'v': pyspedas.tplot_tools.data_quants[tvar1].coords['spec_bins'].values})
50	pyspedas.tplot_tools.data_quants[newname].attrs = copy.deepcopy(pyspedas.tplot_tools.data_quants[tvar1].attrs)	×
51	else:
UNCOV 52	store_data(newname, data={'x':pyspedas.tplot_tools.data_quants[tvar1].coords['time'].values, 'y': data})	1✔
UNCOV 53	pyspedas.tplot_tools.data_quants[newname].attrs = copy.deepcopy(pyspedas.tplot_tools.data_quants[tvar1].attrs)	1✔
54
UNCOV 55	return newname	1✔

spedas / pyspedas / 26171438601

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous