ENH: Create survey plot of simultaneous THEMIS data (use magfig() as template for inputs)

[KCollins]

Concurrent plots of THEMIS data (velocity, dynamic pressure, e.g.) - Copy panels C and D (or M/N) in Figure 1 of doi.org/10.1029/2022GL099768.

Optionally, may call this from magall() function.

[KCollins]

Here's the figure, for reference:

Panels (a–k) show (a) magnetic field in GSE coordinates, (b) ion density, (c) ion velocity, (d) X component of dynamic pressure in GSE coordinates, (e) ion energy flux, (f) detrended ground magnetic field at Vieze Island (VIZ), (g) Pc5 band-pass filtered ground magnetic field data at VIZ, (h) wavelet power spectrum of urn:x-wiley:00948276:media:grl65024:grl65024-math-0015 at VIZ, (i) urn:x-wiley:00948276:media:grl65024:grl65024-math-0016 at VIZ, (j) wavelet power spectrum of urn:x-wiley:00948276:media:grl65024:grl65024-math-0017 at VIZ, and (k) urn:x-wiley:00948276:media:grl65024:grl65024-math-0018 at VIZ, respectively, measured on 28 June 2016. Panels (l–v) are in the same format as panels (a–k), except that the ground magnetic field in panels (q–v) was measured at Ittoqqortoormiit (SCO) and all of the data shown in these panels were measured on 16 June 2008.