The R Markdown code is available on GitHub

Dependencies

library(lubridate)
library(tidyverse)
library(hrbrthemes)
library(viridis)

Read Data

rnhions <- read.csv('../../../tmp/nasa/spdf/new-horizons/swap/ions/new_horizons_swap_pickup-ions_20081116180800_v1.0.1.csv')

Convert date time

nhions <- rnhions %>% mutate(ts = as_datetime(time / 1000))

Range

Time range

summary(nhions$ts)
##                  Min.               1st Qu.                Median 
## "2008-11-16 18:08:26" "2013-10-30 18:08:27" "2015-06-07 18:08:28" 
##                  Mean               3rd Qu.                  Max. 
## "2015-01-30 05:54:24" "2016-06-09 06:08:28" "2017-03-31 18:08:28"

Distance range

summary(nhions$distance)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   11.63   27.70   32.61   31.48   35.65   38.07

Plots

Distance from sun

ggplot(data = nhions) +
  geom_point(aes(x = ts, y = distance)) +
  labs(title = "New Horizons SWAP ions data", x = "Time", y = "Distance from Sun (AU)") +
  theme_light()

Density

ggplot(data = nhions) +
  geom_point(aes(x = ts, y = density)) +
  labs(title = "New Horizons SWAP ions density", x = "Time", y = "cm**-3") +
  theme_light()

ggplot(data = nhions) +
  geom_point(aes(x = distance, y = density)) +
  labs(title = "New Horizons SWAP ions density", x = "Distance from Sun (AU)", y = "cm**-3") +
  theme_light()

Temperature

ggplot(data = nhions) +
  geom_point(aes(x = ts, y = temperature)) +
  labs(title = "New Horizons SWAP ions Temperature", x = "Time", y = "K") +
  theme_light()

ggplot(data = nhions) +
  geom_point(aes(x = distance, y = temperature)) +
  labs(title = "New Horizons SWAP ions Temperature", x = "Distance from Sun (AU)", y = "K") +
  theme_light()

Density and Pressure

ggplot(data = nhions) +
  geom_point(aes(x = density, y = temperature, color = pressure)) +
  labs(title = "New Horizons SWAP ions", x = "Density (cm**-3)", y = "Temperature (K)") +
  theme_light()

ggplot(data = nhions) +
  geom_point(aes(x = pressure, y = temperature, color = density)) +
  labs(title = "New Horizons SWAP ions", x = "Pressure (pPa)", y = "Temperature (K)") +
  theme_light()

ggplot(data = nhions) +
  geom_point(aes(x = density, y = pressure, color = temperature)) +
  labs(title = "New Horizons SWAP ions", x = "Density (cm**-3)", y = "Pressure (pPa)") +
  theme_light()