Dynamic power spectra calculations
using CDAWeb
using SPEDAS
using DimensionalData
using CairoMakie, SpacePhysicsMakie
t0 = "2015-10-16"
t1 = "2015-10-16T03"
ds = CDAWeb.get_dataset("MMS4_SCM_SRVY_L2_SCSRVY", t0, t1; clip=true)
da = DimArray(ds["mms4_scm_acb_gse_scsrvy_srvy_l2"])
tplot(da)
using DSP
metadata = Dict("DISPLAY_TYPE" => "spectrogram", "SCALETYP" => log10, :ylabel => "Frequency (Hz)", :yscale => log10)
pvar = setmeta(pspectrum(da; nfft = 512), metadata)
f = tplot(pvar[:,:,1])
ylims!(5e-2, 2e1)
f