Python plot_spectrum_multi Examples

Example #1

def test_constant_current(series_test):
    xdata = 2 * series_test.parameters.get("k3") / series_test.parameters.get("j0")
    p = pylbo.plot_spectrum_multi(
        series_test, xdata=xdata, use_squared_omega=True, markersize=3
    )
    p.ax.set_yscale("symlog", linthresh=1e-8)
    p.ax.set_xlim(1.88, 2.12)
    p.ax.set_ylim(-1e-3, 1e6)
    return p.fig

Example #2

File: test_multi_gravito_MHD.py Project: n-claes/legolas

def test_gravito_mhd(series_test):
    xdata = series_test.get_k0_squared()
    y_scaling = 1 / series_test.get_alfven_speed(which_values="average") ** 2
    p = pylbo.plot_spectrum_multi(
        series_test, xdata=xdata, use_squared_omega=True, markersize=3
    )
    p.set_y_scaling(y_scaling)
    p.ax.set_xlim(0, 550)
    p.ax.set_ylim(0, 550)
    return p.fig

Example #3

File: test_multi_coronal_fluxtube.py Project: n-claes/legolas

def test_coronal_fluxtubes(series_test):
    xdata = (series_test.parameters.get("k3") *
             coronal_tube_setup["config"]["parameters"]["r0"])
    y_scaling = 1 / (xdata * series_test.get_sound_speed("maximum"))
    p = pylbo.plot_spectrum_multi(series_test,
                                  xdata=xdata,
                                  use_squared_omega=False,
                                  markersize=3)
    p.set_y_scaling(y_scaling)
    p.ax.set_ylim(0.84, 5.1)
    return p.fig

Example #4

File: test_multi_interchange.py Project: n-claes/legolas

def test_interchange_flow(series_test):
    xdata = np.linspace(0, np.pi, len(series_test)) / np.pi
    y_scaling = 1 / series_test.get_alfven_speed(which_values="average") ** 2
    p = pylbo.plot_spectrum_multi(
        series_test, xdata=xdata, use_squared_omega=True, markersize=3
    )
    p.set_y_scaling(y_scaling)
    p.ax.set_xlim(-0.01, 1.01)
    p.ax.set_ylim(-4.1, 14.4)
    p.ax.set_xticks(np.arange(0, 1.2, 0.2))
    p.ax.set_yticks(np.arange(-4, 16, 2))
    return p.fig