# Get this figure: fig <- get_figure("cbrandon", 7)
# Get this figure's data: data <- get_figure("cbrandon", 7)$data
# Add data to this figure: p <- add_trace(p, x=c(4, 5), y=c(4, 5), kwargs=list(filename="Vermont Lunar CubeSat orbit", fileopt="extend"))
# Get y data of first trace: y1 <- get_figure("cbrandon", 7)$data[[1]]$y

# Get figure documentation: https://plot.ly/r/get-requests/
# Add data documentation: https://plot.ly/r/file-options/

# You can reproduce this figure in R with the following code!

# Learn about API authentication here: https://plot.ly/r/getting-started
# Find your api_key here: https://plot.ly/settings/api

library(plotly)
trace1 <- list(
  x = c(108.0, 114.0, 241.0, 299.0, 308.0, 409.0, 432.0, 462.0, 777.0, 869.0, 922.0, 923.0, 920.0, 920.0, 920.0, 1087.0, 1087.0, 1166.0, 1233.0, 1233.0, 1497.0, 1562.0, 1661.0, 1775.0, 1963.0, 1963.0, 2091.0, 2235.0, 2370.0, 2415.0, 2417.0, 2578.0, 2795.0, 2963.0, 3153.0, 3322.0, 3555.0, 3714.0, 3903.0, 4079.0, 4318.0, 4630.0, 4812.0), 
  y = c(40.5172, 40.5114, 40.519400000000005, 40.5174, 40.5189, 40.519800000000004, 40.5186, 40.518100000000004, 40.519600000000004, 40.5187, 40.5199, 40.5182, 40.516600000000004, 40.516600000000004, 40.516600000000004, 40.5206, 40.5206, 40.519000000000005, 40.5197, 40.5197, 40.5165, 40.5163, 40.514500000000005, 40.519000000000005, 40.517, 40.517, 40.5105, 40.5075, 40.5142, 40.507400000000004, 40.506800000000005, 40.516600000000004, 40.5161, 40.5159, 40.5159, 40.5158, 40.5165, 40.5168, 40.5169, 40.5173, 40.517500000000005, 40.518, 40.5186), 
  name = "Inclination", 
  type = "scatter", 
  yaxis = "y2"
)
trace2 <- list(
  x = c(108.0, 114.0, 241.0, 299.0, 308.0, 409.0, 432.0, 462.0, 777.0, 869.0, 922.0, 923.0, 920.0, 920.0, 920.0, 1087.0, 1087.0, 1166.0, 1233.0, 1233.0, 1497.0, 1562.0, 1661.0, 1775.0, 1963.0, 1963.0, 2091.0, 2235.0, 2370.0, 2415.0, 2417.0, 2578.0, 2795.0, 2963.0, 3153.0, 3322.0, 3555.0, 3714.0, 3903.0, 4079.0, 4318.0, 4630.0, 4812.0), 
  y = c(6035.0, 5162.0, 5877.0, 6397.0, 6112.0, 6051.0, 6195.0, 6073.0, 6056.0, 5349.0, 5267.0, 5612.0, 5323.0, 5323.0, 5323.0, 5237.0, 5237.0, 5517.0, 4593.0, 4593.0, 5226.0, 4327.0, 3506.0, 4833.0, 4527.0, 4527.0, 2598.0, 2462.0, 3438.0, 2817.0, 2487.0, 3644.0, 3218.0, 3007.0, 2934.0, 2985.0, 2967.0, 3004.0, 3037.0, 3130.0, 3366.0, 3707.0, 3924.0), 
  name = "Eccentr.", 
  type = "scatter"
)
trace3 <- list(
  x = c(108.0, 114.0, 241.0, 299.0, 308.0, 409.0, 432.0, 462.0, 777.0, 869.0, 922.0, 923.0, 920.0, 920.0, 920.0, 1087.0, 1087.0, 1166.0, 1233.0, 1233.0, 1497.0, 1562.0, 1661.0, 1775.0, 1963.0, 1963.0, 2091.0, 2235.0, 2370.0, 2415.0, 2417.0, 2578.0, 2795.0, 2963.0, 3153.0, 3322.0, 3555.0, 3714.0, 3903.0, 4079.0, 4318.0, 4630.0, 4812.0), 
  y = c(15.21012737, 15.21036615, 15.21109505, 15.21097009, 15.21113058, 15.21127139, 15.21133141, 15.211344760000001, 15.211971870000001, 15.21212188, 15.2121475, 15.21221282, 15.21217141, 15.21217141, 15.21217141, 15.21236756, 15.21236756, 15.21248494, 15.21248898, 15.21248898, 15.21270573, 15.21279416, 15.212773160000001, 15.21316468, 15.21319895, 15.21319895, 15.21346638, 15.21355108, 15.21366688, 15.21388864, 15.213884890000001, 15.20964887, 15.20980078, 15.209936560000001, 15.21012142, 15.21027329, 15.210460020000001, 15.210612900000001, 15.21081767, 15.211006040000001, 15.21123837, 15.2115199, 15.21167118), 
  name = "M Motion", 
  type = "scatter", 
  yaxis = "y3"
)
data <- list(trace1, trace2, trace3)
layout <- list(
  autosize = TRUE, 
  bargap = 0.2, 
  bargroupgap = 0, 
  barmode = "stack", 
  boxgap = 0.3, 
  boxgroupgap = 0.3, 
  boxmode = "overlay", 
  dragmode = "zoom", 
  font = list(
    color = "#000", 
    family = "'Open sans', verdana, arial, sans-serif", 
    size = 12
  ), 
  height = 565, 
  hovermode = "x", 
  legend = list(
    x = 0.98, 
    y = 0.98, 
    bgcolor = "#fff", 
    bordercolor = "#000", 
    borderwidth = 1, 
    font = list(
      color = "", 
      family = "", 
      size = 0
    ), 
    traceorder = "normal"
  ), 
  margin = list(
    r = 80, 
    t = 100, 
    b = 80, 
    l = 80, 
    pad = 2
  ), 
  paper_bgcolor = "#fff", 
  plot_bgcolor = "#fff", 
  separators = ".,", 
  showlegend = TRUE, 
  title = "Vermont Lunar CubeSat", 
  titlefont = list(
    color = "rgb(0, 0, 255)", 
    family = "", 
    size = 0
  ), 
  width = 1344, 
  xaxis = list(
    anchor = "y", 
    autorange = TRUE, 
    autotick = TRUE, 
    domain = c(0.1, 1), 
    dtick = 500, 
    exponentformat = "e", 
    gridcolor = "#ddd", 
    gridwidth = 1, 
    linecolor = "#000", 
    linewidth = 1, 
    mirror = "all", 
    nticks = 0, 
    overlaying = FALSE, 
    position = 0, 
    range = c(108, 4812), 
    rangemode = "normal", 
    showexponent = "all", 
    showgrid = TRUE, 
    showline = TRUE, 
    showticklabels = TRUE, 
    tick0 = 0, 
    tickangle = "auto", 
    tickcolor = "#000", 
    tickfont = list(
      color = "rgb(0, 0, 255)", 
      family = "", 
      size = 0
    ), 
    ticklen = 5, 
    ticks = "outside", 
    tickwidth = 1, 
    title = "Revolutions since November 19, 2013, 8:15pm EST Launch", 
    titlefont = list(
      color = "rgb(0, 0, 255)", 
      family = "", 
      size = 0
    ), 
    type = "linear", 
    zeroline = TRUE, 
    zerolinecolor = "#000", 
    zerolinewidth = 1
  ), 
  yaxis = list(
    anchor = "x", 
    autorange = TRUE, 
    autotick = TRUE, 
    domain = c(0, 1), 
    dtick = 500, 
    exponentformat = "e", 
    gridcolor = "#ddd", 
    gridwidth = 1, 
    linecolor = "#000", 
    linewidth = 1, 
    mirror = "all", 
    nticks = 0, 
    overlaying = FALSE, 
    position = 0, 
    range = c(2243.3888888888887, 6615.611111111111), 
    rangemode = "normal", 
    showexponent = "all", 
    showgrid = TRUE, 
    showline = TRUE, 
    showticklabels = TRUE, 
    tick0 = 0, 
    tickangle = "auto", 
    tickcolor = "#000", 
    tickfont = list(
      color = "rgb(255, 127, 14)", 
      family = "", 
      size = 0
    ), 
    ticklen = 5, 
    ticks = "outside", 
    tickwidth = 1, 
    title = "Eccentricity x 10^-7", 
    titlefont = list(
      color = "rgb(255, 127, 14)", 
      family = "", 
      size = 0
    ), 
    type = "linear", 
    zeroline = TRUE, 
    zerolinecolor = "#000", 
    zerolinewidth = 1
  ), 
  yaxis2 = list(
    anchor = "x", 
    autorange = TRUE, 
    autotick = TRUE, 
    domain = c(0, 1), 
    dtick = 0.002, 
    exponentformat = "e", 
    gridcolor = "#ddd", 
    gridwidth = 1, 
    linecolor = "#000", 
    linewidth = 1, 
    mirror = "all", 
    nticks = 0, 
    overlaying = "y", 
    position = 0, 
    range = c(40.50603333333334, 40.521366666666665), 
    rangemode = "normal", 
    showexponent = "all", 
    showgrid = TRUE, 
    showline = TRUE, 
    showticklabels = TRUE, 
    side = "right", 
    tick0 = 0, 
    tickangle = "auto", 
    tickcolor = "#000", 
    tickfont = list(
      color = "rgb(31, 119, 180)", 
      family = "", 
      size = 0
    ), 
    ticklen = 5, 
    ticks = "outside", 
    tickwidth = 1, 
    title = "Inclination (degrees)", 
    titlefont = list(
      color = "rgb(31, 119, 180)", 
      family = "", 
      size = 0
    ), 
    type = "linear", 
    zeroline = TRUE, 
    zerolinecolor = "#000", 
    zerolinewidth = 1
  ), 
  yaxis3 = list(
    anchor = "free", 
    autorange = TRUE, 
    autotick = TRUE, 
    domain = c(0, 1), 
    dtick = 0.0005, 
    exponentformat = "e", 
    gridcolor = "#ddd", 
    gridwidth = 1, 
    linecolor = "#000", 
    linewidth = 1, 
    mirror = "all", 
    nticks = 0, 
    overlaying = "y", 
    position = 0, 
    range = c(15.209413327222222, 15.214124182777779), 
    rangemode = "normal", 
    showexponent = "all", 
    showgrid = TRUE, 
    showline = TRUE, 
    showticklabels = TRUE, 
    side = "left", 
    tick0 = 0, 
    tickangle = "auto", 
    tickcolor = "#000", 
    tickfont = list(
      color = "rgb(44, 160, 44)", 
      family = "", 
      size = 0
    ), 
    ticklen = 5, 
    ticks = "outside", 
    tickwidth = 1, 
    title = "Mean Montion (revolutions/day)", 
    titlefont = list(
      color = "rgb(44, 160, 44)", 
      family = "", 
      size = 0
    ), 
    type = "linear", 
    zeroline = TRUE, 
    zerolinecolor = "#000", 
    zerolinewidth = 1
  )
)
p <- plot_ly()
p <- add_trace(p, x=trace1$x, y=trace1$y, name=trace1$name, type=trace1$type, yaxis=trace1$yaxis)
p <- add_trace(p, x=trace2$x, y=trace2$y, name=trace2$name, type=trace2$type)
p <- add_trace(p, x=trace3$x, y=trace3$y, name=trace3$name, type=trace3$type, yaxis=trace3$yaxis)
p <- layout(p, autosize=layout$autosize, bargap=layout$bargap, bargroupgap=layout$bargroupgap, barmode=layout$barmode, boxgap=layout$boxgap, boxgroupgap=layout$boxgroupgap, boxmode=layout$boxmode, dragmode=layout$dragmode, font=layout$font, height=layout$height, hovermode=layout$hovermode, legend=layout$legend, margin=layout$margin, paper_bgcolor=layout$paper_bgcolor, plot_bgcolor=layout$plot_bgcolor, separators=layout$separators, showlegend=layout$showlegend, title=layout$title, titlefont=layout$titlefont, width=layout$width, xaxis=layout$xaxis, yaxis=layout$yaxis, yaxis2=layout$yaxis2, yaxis3=layout$yaxis3)