mvn_kp_altplot.pro

;+
;
; :Name: mvn_kp_altplot
; 
; Copyright 2017 Regents of the University of Colorado. All Rights Reserved.
; Released under the MIT license.
; This software was developed at the University of Colorado's Laboratory for Atmospheric and Space Physics.
; Verify current version before use at: https://lasp.colorado.edu/maven/sdc/public/pages/software.html
; :Author: Kristopher Larsen
;   
; :Description:
;   This simple routine plots one or more altitude profiles from the insitu 
;   KP data structure.  Any data fields may be plotted together, on 
;   individual or single plots, using both direct and function graphics. 
;
; :Params:
;    kp_data: in, required, type=structure
;       the INSITU KP data structure from which to plot data
;    parameter: in, required, type=strarr,intarr
;       the INSITU kp data fields to plot, maybe an integer or string array 
;       for multiple choices
;    time: in, optional, can be a scalar or a two item array of type:
;         long(s)        orbit number
;         string(s)      format:  YYYY-MM-DD/hh:mm:ss       
;       A start or start & stop time (or orbit #) range for reading kp data. 
; 
; :Keywords:
;    list: in, optional, type=boolean
;       Used to print out the contents of the input data structure.
;           If set as a keyword, /list, this is printed to the screen.
;           If set as a variable, list=list, a string array is returned 
;           containing the structure index and tag names.
;    range: in, optional, type=boolean
;       if selected, will list the beginning and end times of kp_data
;    plot_object: out, optional
;       if provided, the name of the variable to which the plot created
;       will be assigned.  Ignored if /directgraphic is also set.
;    directgraphic: in, optional, type=boolean
;       if selected, will override teh default Graphics plot procedure and 
;       use direct graphics instead
;    davin: in, optional, type=boolean
;       As requested by Davin Larson, this keyword will flip the X and Y 
;       axis of each plot.
;
; :Obsolete:
;    The following are kept in case we wish to return this functionality
;
;    xrange: in, optional, type=fltarr
;       Minimum and maximum range for the x-axis. If multiple plots are
;       included, the number of xrange arrays must match.
;    yrange: in, optional, type=fltarr
;       Minimum and maximum range for the y-axis. If multiple plots are
;       included, the number of yrange arrays must match.
;
; :Version:   1.1 (2015-Apr 28) McGouldrick  
;
; :History:
;   v1.0 (July 8, 2014)
;-
@mvn_kp_tag_parser
@mvn_kp_tag_list
@mvn_kp_range
@mvn_kp_range_select
@mvn_kp_tag_verify

pro MVN_KP_ALTPLOT, kp_data, parameter, time=time, list=list, range=range, $
                    directgraphic=directgraphic, $
                    davin=davin, y_labels=y_labels, $
                    plot_object=plot_object, $
                    error=error, _extra = e, help=help


  ;provide help for those who don't have IDLDOC installed
  if keyword_set(help) then begin
    mvn_kp_get_help,'mvn_kp_altplot'
    return
  endif

  ;CHECK THAT THE INPUT PARAMETERS ARE VALID
  ;DETERMINE ALL THE PARAMETER NAMES THAT MAY BE USED LATER
  
  MVN_KP_TAG_PARSER, kp_data, base_tag_count, first_level_count, $
    second_level_count, base_tags,  first_level_tags, second_level_tags

  ;LIST OF ALL POSSIBLE PLOTABLE PARAMETERS IF /LIST IS SET

  if keyword_set(list) then begin
    if arg_present(list) then begin
      ; return the KP list to the provided variable
      mvn_kp_get_list, kp_data, list=list
    endif else begin
      ; print the KP list to screen
      mvn_kp_get_list, kp_data, /list
    endelse
    return
  endif

  ;PROVIDE THE TEMPORAL RANGE OF THE DATA SET IN BOTH DATE/TIME AND ORBITS 
  ;IF REQUESTED.
  if keyword_set(range) then begin
    MVN_KP_RANGE, kp_data
    return
  endif
  
  ;SET THE VARIOUS PLOT OPTIONS, SHOULD THEY REQUIRE IT
  if keyword_set(title) eq 0 then begin
   if n_elements(parameter) eq 1 then title=''
   if n_elements(parameter) ne 1 then title=strarr(n_elements(parameter))
  endif
  
  if keyword_set(directgraphic) eq 0 then begin
   if Float(!Version.Release) lt 8.0 THEN $
      directgraphic = 1  ;FORCE DIRECT GRAPHICS IF USER HAS OLD VERSION OF IDL
  endif
  
  ;CHECK THAT, IF INCLUDED, XRANGE AND YRANGE CONTAIN SUFFICIENT DATA
;  if keyword_set(yrange) then begin
;    if (n_elements(yrange)/n_elements(parameter)) ne 2 then begin
;      print,'When using the YRANGE keyword, ranges must be set for each plot'
;      return    
;    endif
;  endif
;  if keyword_set(xrange) then begin
;    if (n_elements(xrange)/n_elements(parameter)) ne 2 then begin
;      print,'When using the XRANGE keyword, ranges must be set for each plot'
;      return    
;    endif
;  endif  
  
  ;IF THE USER SUPPLIES A TIME RANGE, SET THE BEGINNING AND END INDICES
  
  if keyword_set(time) then begin ;determine the start and end indices to plot
    MVN_KP_RANGE_SELECT, kp_data, time, kp_start_index, kp_end_index
  endif else begin                ;otherwise plot all data within structure
   kp_start_index = 0
   kp_end_index = n_elements(kp_data.orbit)-1
  endelse
  ;
  ;  Check for invalid time ranges
  ;
  if kp_start_index eq -1 or kp_end_index eq -1 then begin
    print,'Sorry, the times you requested are not contained within the data.'
    print,'Check your time range and try again.'
    return
  endif
  if kp_start_index eq kp_end_index then begin
    print,'Sorry, start and end times are the same. Nothing to plot!'
    return
  endif
  if( kp_start_index gt kp_end_index )then begin
    print,'WARNING: Start time provided is later than end time.'
    print,'         Start time = ',time[0]
    print,'           End time = ',time[1]
    do_swap = ''
    while( do_swap ne 's' and do_swap ne 'S' and $
      do_swap ne 'q' and do_swap ne 'Q' )do begin
      read,prompt="Please press 's' to swap, or 'q' to exit and try again: ",$
           do_swap
    endwhile
    if( do_swap eq 's' or do_swap eq 'q' )then begin
      temp = kp_start_index
      kp_start_index = kp_end_index
      kp_end_index = temp
    endif
    if( do_swap eq 'q' or do_swap eq 'Q' )then $
      message,'Exiting as per user request to re-enter time range'
  endif
;
;  Copy the overplot check from mvn_kp_plot rather than the one used here
;
  overplot_check = 0
  for i = 0,n_elements(parameter)-1 do begin
    pos = strpos(parameter[i],',')
    if pos ne -1 then overplot_check = 1
  endfor

  ;CREATE THE PLOT VECTORS
  
  if overplot_check ne 1 then begin
;
; Insert code form plot here
;
    if n_elements(parameter) eq 1 then begin
    ;only going to plot a single altitude plot

      MVN_KP_TAG_VERIFY, kp_data, parameter,base_tag_count, $
        first_level_count, base_tags, first_level_tags, $
        check, level0_index, level1_index, tag_array

      if check eq 0 then begin  ;CHECK THAT THE REQUESTED PARAMETER EXISTS

        y = kp_data[kp_start_index:kp_end_index].spacecraft.altitude
        x = kp_data[kp_start_index:kp_end_index].(level0_index).(level1_index)

        if keyword_set(error) then begin
          mvn_kp_tag_verify, kp_data, error, base_tag_count, $
            first_level_count, base_tags, first_level_tags, $
            err_check, err_level0, err_level1, temp_tag
          if err_check eq 0 then begin

            mvn_kp_define_error_bars, kp_data[kp_start_index:kp_end_index], $
              level0_index, level1_index, base_tags, $
              err_level0, err_level1, temp_tag, $
              x_error

          endif else begin

            print,'Requested error parameter is not included in the data.'
            print,'Try /LIST to check for it.'
            print,'Creating requested plot WITHOUT error bars'

          endelse

        endif else err_check = 1 ; if not plotting error bars...

        ;SET THE YLABELS
        if keyword_set(y_labels) ne 1 then begin
          y_labels = strupcase(string(tag_array[0]+'.'+tag_array[1]))
        endif

      endif else begin
        print,'Requested plot parameter is not included in the data.'
        print,'Try /LIST to confirm your parameter choice.'
        return
      endelse
    endif ;end of single altitude variable creation loop

  ;CREATE SINGLE ALTITUDE PLOT
  
  if keyword_set(directgraphic) eq 0 then begin 
    ; PLOT USING THE NEW IDL GRAPHICS PLOT FUNCTION
    if n_elements(parameter) eq 1 then begin
      w = window(window_title='Maven KP Data Altitude Plots')
      if keyword_set(davin) then begin
        if err_check eq 0 then begin
          plot1 = errorplot( y,x,x_error, $
                     xtitle='Spacecraft Altitude [km]', $
                     ytitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
                     /current, xstyle=1, ystyle=1, _extra = e)                             
        endif else begin
          plot1 = plot(y,x,xtitle='Spacecraft Altitude, km',$
               ytitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
               /current, xstyle=1, ystyle=1, _extra = e)
        endelse
      endif else begin
        if err_check eq 0 then begin
          null_error = fltarr(2,n_elements(y))
          plot1 = errorplot( x,y,x_error,null_error, $
                  ytitle='Spacecraft Altitude [km]', $
                  xtitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
                  /current, xstyle=1, ystyle=1, _extra = e)
        endif else begin
          plot1 = plot(x,y,ytitle='Spacecraft Altitude, km',$
                xtitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
                /current, xstyle=1, ystyle=1, _extra = e)
        endelse
      endelse
    endif
  endif
  if keyword_set(directgraphic) ne 0 then begin
    ;USE THE OLD DIRECT GRAPHICS PLOT PROCEDURES
    if n_elements(parameter) eq 1 then begin
      device,decomposed=0
      loadct,0,/silent
      !P.MULTI = [0, n_elements(parameter), 1]
      if keyword_set(davin) then begin
        plot,y,x,xtitle='Spacecraft Altitude, km', $
             ytitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
             background=255, color=0,$
             xrange=yrange,yrange=xrange,xstyle=1,ystyle=1, _extra = e
        if err_check eq 0 then begin
          errplot,y,x-x_error[0,*],x+x_error[1,*],color=0,_extra=e
        endif
      endif else begin
        plot,x,y,ytitle='Spacecraft Altitude, km', $
             xtitle=strupcase(string(tag_array[0]+'.'+tag_array[1])),$
             background=255, color=0,xstyle=1,ystyle=1, _extra = e
        if err_check eq 0 then begin
          print,'***WARNING***'
          print,'mvn_kp_altplot'
          print,'Old direct graphics errplot cannot make horizontal error bars'
          print,'Suggest using either /davin keyword to flip x and y.'
          print,'  or do not use the /directgraphics keyword.'
          print,'Will continue plotting but will ignore error bars.'
        endif
      endelse
    endif
  endif
  
  ;CREATE MULTIPLE ALITIUDE PLOT VECTORS

  if n_elements(parameter) gt 1 then begin
;
;  Now, insert the create multi vectors code
;
    mvn_kp_create_multi_vectors, kp_data[kp_start_index:kp_end_index], $
                                 parameter, x, x_error, x_axis_title, $
                                 error=error, y_labels=y_labels, $
                                 err_check=err_check
    y = kp_data[kp_start_index:kp_end_index].spacecraft.altitude
  endif
  
  ;CREATE DUMMY YRANGE IF NOT DEFINED
  temp_yrange = dblarr(2,n_elements(parameter))
  if n_elements(yrange) ne 0 then begin
    temp_yrange = yrange
  endif else begin
    for i=0,n_elements(parameter)-1 do begin
      temp_yrange[0,i] = min(y[*],/NaN)
      temp_yrange[1,i] = max(y[*],/NaN)
    endfor
  endelse
  ;CREATE DUMMY XRANGE IF NOT DEFINED
  temp_xrange = dblarr(2,n_elements(parameter))
  if n_elements(xrange) ne 0 then begin
    temp_xrange = xrange
  endif else begin
    for i=0,n_elements(parameter)-1 do begin
      temp_xrange[0,i] = min(x[i,*],/NaN)
      temp_xrange[1,i] = max(x[i,*],/NaN)
    endfor
  endelse     
          
  ;CREATE THE MULTPLE ALTITUDE PLOT  
  if keyword_set(directgraphic) eq 0 then begin  
    ;PLOT USING THE NEW IDL GRAPHICS PLOT FUNCTION
    if n_elements(parameter) gt 1 then begin
      w = window(window_title='Maven KP Data Altitude Plots')
      plot1 = []; define empty vector (IDL *+ only)
      for i = 0, n_elements(parameter) -1 do begin
        if keyword_set(davin) then begin
          if err_check[i] ne 0 then begin
            plot1 = [plot1, $
                     plot(y,x[i,*], xtitle='Spacecraft Altitude, km', $
                         ytitle=x_axis_title[i], $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],$
                         xstyle=1,ystyle=1,xrange=temp_yrange[*,i],$
                         yrange=temp_xrange[*,i], _extra = e) ]
          endif else begin
            plot1 = [plot1, $
                     errorplot(y,x[i,*], reform(x_error[*,i,*]), $
                              xtitle='Spacecraft Altitude, km', $
                              ytitle=x_axis_title[i], $
                              layout=[n_elements(parameter),1,i+1],/current,$
                              title=title[i],$
                              xstyle=1,ystyle=1,xrange=temp_yrange[*,i],$
                              yrange=temp_xrange[*,i], _extra = e) ]
          endelse
        endif else begin ; do not swap the axis
          if err_check[i] ne 0 then begin
            plot1 = [plot1, $
                     plot(x[i,*], y, ytitle='Spacecraft Altitude, km', $
                         xtitle=x_axis_title[i], $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],$
                         xstyle=1,ystyle=1,yrange=temp_yrange[*,i],$
                         xrange=temp_xrange[*,i], _extra = e) ]
          endif else begin
            null_error = fltarr(2,n_elements(y))
            plot1 = [plot1, $
                     errorplot(x[i,*], y, reform(x_error[*,i,*]), null_error, $
                              ytitle='Spacecraft Altitude, km', $
                              xtitle=x_axis_title[i], $
                              layout=[n_elements(parameter),1,i+1],/current,$
                              title=title[i],$
                              xstyle=1,ystyle=1,yrange=temp_yrange[*,i],$
                              xrange=temp_xrange[*,i], _extra = e) ]
          endelse ; err_check
        endelse   ; davin check
      endfor ; loop over all parameters
    endif    ; if there is more than on parameter
  endif      ; if object oriented graphics
  
  if keyword_set(directgraphic) ne 0 then begin 
    ;PLOT USING THE OLD IDL DIRECT GRAPHICS
          device,decomposed=0
          loadct,0,/silent
    if n_elements(parameter) gt 1 then begin
      !P.MULTI = [0, n_elements(parameter), 1]
      if keyword_set(davin) then begin
        plot,y,x[0,*],xtitle='Spacecraft Altitude, km', $
             ytitle=x_axis_title[0],xstyle=1,ystyle=1,$
             title=title[0],background=255,color=0,$
             charsize=2,font=-1,$
             xrange=temp_yrange[*,0],yrange=temp_xrange[*,0], _extra = e
        ;
        ;  If requested, plot the error bars
        ;
        if err_check[0] eq 0 then begin
          errplot,y,x[0,*]-x_error[0,0,*],x[0,*]+x_error[1,0,*], $
                  color=0, _extra=e
        endif
        for i=1,n_elements(parameter)-1 do begin
          plot,y,x[i,*],xtitle='Spacecraft Altitude, km', $
               ytitle=x_axis_title[i],xstyle=1,ystyle=1,$
               title=title[i],color=0,charsize=2,font=-1,$
               xrange=temp_yrange[*,i],yrange=temp_xrange[*,i], _extra = e
          ;
          ;  If requested, plot the error bars
          ;
          if err_check[i] eq 0 then begin
            errplot,y,x[i,*]-x_error[0,i,*],x[i,*]+x_error[1,i,*], $
                    color=0, _extra=e
          endif
        endfor 
      endif else begin
        if( where(err_check eq 1) lt 0 )then begin
          ;
          ;  Then at least one error bars has been requested
          ;
          print,'***WARNING***'
          print,'mvn_kp_altplot'
          print,'Old direct graphics errplot cannot make horizontal error bars'
          print,'Suggest using either /davin keyword to flip x and y.'
          print,'  or do not use the /directgraphics keyword.'
          print,'Will continue plotting but will ignore error bars.'
        endif
        ;
        ;  Now, plot the relevant parameters
        ;
        plot,x[0,*],y,ytitle='Spacecraft Altitude, km', $
             xtitle=x_axis_title[0],xstyle=1,ystyle=1,$
             title=title[0],background=255,color=0,charsize=2,font=-1,$
             yrange=temp_yrange[*,0],xrange=temp_xrange[*,0], _extra = e
        for i=1,n_elements(parameter)-1 do begin
          plot,x[i,*],y,ytitle='Spacecraft Altitude, km', $
               xtitle=x_axis_title[i],xstyle=1,ystyle=1,$
               title=title[i], color=0,charsize=2,font=-1,$
               yrange=temp_yrange[*,i],xrange=temp_xrange[*,i], _extra = e
        endfor 
      endelse ; if davin then else
    endif     ; if there is more tha one parameter
  endif       ; if using direct graphics
 endif else begin ; Finished with no over plots; now do overplots

    ;BEGIN SEPARATE ROUTINES IF ANY OVERPLOTTING IS REQUIRED.

  ;ANALYZE THE INPUT STRINGS TO DETERMINE PARAMETERS AND SIZES
    
    plot_count =intarr(n_elements(parameter))
    total_lines = 0
    true_index = intarr(50)    
    
    for i=0, n_elements(parameter)-1 do begin
      check = strmatch(parameter[i],'*,*')
      if check eq 1 then begin                      ;over plots 
        extract = strmid(strsplit(parameter[i],',',/extract) ,0,1)
        new_param = strsplit(parameter[i],',',/extract)
        for j=0,n_elements(extract)-1 do begin
          if strmatch(extract[j],'[0123456789]') eq 1 then begin 
            ;structure index call
            true_index[total_lines] = fix(new_param[j])
          endif else begin
            ;call based on structure names
            mvn_kp_structure_index, kp_data, new_param[j], new_index, $
                                    first_level_tags
            true_index[total_lines] = new_index            
          endelse
          total_lines = total_lines+1
          plot_count[i] = plot_count[i]+1
        endfor    
      endif else begin                              ;single plots
        extract = strmid(parameter[i],0,1)
        new_param = strsplit(parameter[i],',',/extract)
        if strmatch(extract,'[0123456789]') eq 1 then begin
          ;structure index call
          true_index[total_lines] = fix(parameter[i])
        endif else begin
          ;structure name call
          mvn_kp_structure_index, kp_data, new_param, new_index, $
                                  first_level_tags
          true_index[total_lines] = new_index
        endelse
        total_lines = total_lines + 1
        plot_count[i] = 1
      endelse
    endfor
 
    true_index = true_index[0:total_lines-1]

  ;CHECK PARAMETER VALIDITY AND EXTRACT DATA
  
  ;
  ; Now, use the create_multi_vector code with true_index in place
  ;  of parameter to fill the data, error, and label vectors
  ;
  y = kp_data[kp_start_index:kp_end_index].spacecraft.altitude
  mvn_kp_create_multi_vectors, kp_data[kp_start_index:kp_end_index], $
                               true_index, x, x_error, x_axis_title, $
                               error=error, y_labels=y_labels, $
                               err_check=err_check

  ;CREATE DUMMY YRANGE IF NOT DEFINED
  temp_yrange = dblarr(2,n_elements(true_index))
  if n_elements(yrange) ne 0 then begin
    temp_yrange = yrange
  endif else begin
    for i=0,n_elements(true_index)-1 do begin
      temp_yrange[0,i] = min(y,/NaN)
      temp_yrange[1,i] = max(y,/NaN)
    endfor
  endelse
  ;CREATE DUMMY XRANGE IF NOT DEFINED
  temp_xrange = dblarr(2,n_elements(true_index))
  if n_elements(xrange) ne 0 then begin
    temp_xrange = xrange
  endif else begin
    for i=0,n_elements(true_index)-1 do begin
      temp_xrange[0,i] = min(x[i,*],/NaN)
      temp_xrange[1,i] = max(x[i,*],/NaN)
    endfor
  endelse  

  ;CREATE THE PLOTS
  
  if keyword_set(directgraphic) eq 0 then begin
  ;   if n_elements(parameter) gt 1 then begin
      oplot_index = 0
       w = window(window_title='Maven KP Data Altitude Plots')
      plot1 = [] ; define empty array for plots (IDL 8+ only)
      for i = 0, n_elements(parameter) -1 do begin
        if keyword_set(davin) then begin
          if plot_count[i] eq 1 then begin
            plot1 = [plot1, $
                     plot(y,x[oplot_index,*], $
                         xtitle='Spacecraft Altitude [km]', $
                         ytitle=x_axis_title[oplot_index], $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],$
                         xstyle=1,ystyle=1,yrange=temp_xrange[*,oplot_index],$
                         xrange=temp_yrange[*,oplot_index], _extra = e) ]
            oplot_index= oplot_index+1
          endif else begin
            plot1 = [plot1, $
                     plot(y,x[oplot_index,*], $
                         xtitle='Spacecraft Altitude [km]', $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],linestyle=0,$
                         yrange=temp_xrange[*,oplot_index],xstyle=1,ystyle=1,$
                         xrange=temp_yrange[*,oplot_index],$
                         name=x_axis_title[oplot_index], _extra = e) ]
            l = legend(target=plot1,$
                       position=[(i*(1./n_elements(parameter)))$
                                 +(.5/(n_elements(parameter))), .85],$
                       /normal,linestyle=0,font_size=8)
            oplot_index = oplot_index+1
            hold_index = oplot_index-1
            for j=1,plot_count[i]-1 do begin      
              plot1 = [plot1, $
                       plot(y,x[oplot_index,*], $
                           xtitle='Spacecraft Altitude[ km]', $
                           layout=[n_elements(parameter),1,i+1],/current,$
                           title=title[i],linestyle=j,$
                           xstyle=1,ystyle=1,$
                           yrange=temp_xrange[*,hold_index],$
                           name=x_axis_title[oplot_index], $
                           xrange=temp_yrange[*,hold_index], _extra = e) ]
              l = legend(target=plot1,$
                         position=[(i*(1./n_elements(parameter))) $
                                   +(.5/(n_elements(parameter))), $
                                   .85+(j*0.05)],$
                         /normal,linestyle=j,font_size=8)
              oplot_index=oplot_index+1
            endfor    
          endelse
        endif else begin
          if plot_count[i] eq 1 then begin
            plot1 = [plot1, $
                     plot(x[oplot_index,*], y, $
                         ytitle='Spacecraft Altitude [km]', $
                         xtitle=x_axis_title[oplot_index], $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],$
                         xstyle=1,ystyle=1,xrange=temp_xrange[*,oplot_index],$
                         yrange=temp_yrange[*,oplot_index], _extra = e) ]
            oplot_index= oplot_index+1
          endif else begin
            plot1 = [plot1, $
                     plot(x[oplot_index,*], y, $
                         ytitle='Spacecraft Altitude [km]', $
                         layout=[n_elements(parameter),1,i+1],/current,$
                         title=title[i],linestyle=0,$
                         xrange=temp_xrange[*,oplot_index],xstyle=1,ystyle=1,$
                         yrange=temp_yrange[*,oplot_index],$
                         name=x_axis_title[oplot_index], _extra = e) ]
            l = legend(target=plot1,$
                       position=[(i*(1./n_elements(parameter))) $
                                 +(.5/(n_elements(parameter))),.85],$
                       /normal,linestyle=0,font_size=8)
            oplot_index = oplot_index+1
            hold_index = oplot_index-1
            for j=1,plot_count[i]-1 do begin      
              plot1 = [plot1, $
                       plot(x[oplot_index,*], y, $
                           ytitle='Spacecraft Altitude [km]', $
                           layout=[n_elements(parameter),1,i+1],/current,$
                           title=title[i],linestyle=j,$
                           xstyle=1,ystyle=1,$
                           xrange=temp_xrange[*,hold_index],$
                           name=x_axis_title[oplot_index], $
                           yrange=temp_yrange[*,hold_index], _extra = e) ]
               l = legend(target=plot1,$
                          position=[(i*(1./n_elements(parameter)))$
                                    +(.5/(n_elements(parameter))),$
                                    .85+(j*0.05)],$
                          /normal,linestyle=j,font_size=8)
              oplot_index=oplot_index+1
            endfor    
          endelse
        endelse 
      endfor
  endif 
  if keyword_set(directgraphic) eq 1 then begin
    device,decomposed=1
  ;  if n_elements(parameter) gt 1 then begin
      !P.MULTI = [0, n_elements(parameter), 1]
      oplot_index = 0 
      if keyword_set(davin) then begin
        for i = 0, n_elements(parameter) -1 do begin
          if plot_count[i] eq 1 then begin
            plot,y,x[oplot_index,*],xtitle='Spacecraft Altitude [km]', $
                 xtitle=x_axis_title[oplot_index],title=title[i],$
                 background='FFFFFF'x,color=0,$
                 font=-1,xstyle=1,ystyle=1,$
                 yrange=temp_xrange[*,oplot_index],$
                 xrange=temp_yrange[*,oplot_index], _extra = e
            oplot_index = oplot_index+1
          endif else begin 
            plot,y,x[oplot_index,*],xtitle='Spacecraft Altitude [km]',$
                 title=title[oplot_index],thick=thick,linestyle=linestyle,$
                 xlog=xaxis_log,ylog=yaxis_log, background='FFFFFF'x,$
                 color=0,xstyle=1,ystyle=1,$
                 yrange=temp_xrange[*,oplot_index],$
                 xrange=temp_yrange[*,oplot_index], _extra = e
            plots,[(i*(1./n_elements(parameter)))$
                     +(.25/(n_elements(parameter))),$
                   (i*(1./n_elements(parameter)))$
                     +(.48/(n_elements(parameter)))],$
                   [.81,.81],linestyle=0,color=0,/normal
            xyouts,(i*(1./n_elements(parameter)))$
                    +(.5/(n_elements(parameter))),$
                    .8,x_axis_title[oplot_index],color=0,/normal
            oplot_index = oplot_index+1
            for j=1,plot_count[i]-1 do begin      
              oplot,y,x[oplot_index,*],linestyle=j,thick=thick,color=0
              plots,[(i*(1./n_elements(parameter)))$
                     +(.25/(n_elements(parameter))),$
                     (i*(1./n_elements(parameter)))$
                     +(.48/(n_elements(parameter)))],$
                     [.81+(j*0.03),.81+(j*0.03)],linestyle=j,color=0,/normal
              xyouts,(i*(1./n_elements(parameter)))$
                     +(.5/(n_elements(parameter))),$
                     .8+(j*.03),x_axis_title[oplot_index],color=0,/normal
              oplot_index=oplot_index+1
            endfor        
          endelse 
         endfor         
      endif else begin
        for i = 0, n_elements(parameter) -1 do begin
          if plot_count[i] eq 1 then begin
            plot,x[oplot_index,*],y,ytitle='Spacecraft Altitude [km]', $
                 xtitle=x_axis_title[oplot_index],title=title[i],$
                 background='FFFFFF'x,color=0,$
                 font=-1,xstyle=1,ystyle=1,$
                 xrange=temp_xrange[*,oplot_index],$
                 yrange=temp_yrange[*,oplot_index], _extra = e
            oplot_index = oplot_index+1
          endif else begin 
            plot,x[oplot_index,*],y,ytitle='Spacecraft Altitude [km]',$
                 title=title[oplot_index],$
                 background='FFFFFF'x,color=0,xstyle=1,ystyle=1,$
                 xrange=temp_xrange[*,oplot_index],$
                 yrange=temp_yrange[*,oplot_index], _extra = e
            plots,[(i*(1./n_elements(parameter)))$
                   +(.25/(n_elements(parameter))),$
                   (i*(1./n_elements(parameter)))$
                   +(.48/(n_elements(parameter)))],$
                   [.81,.81],linestyle=0,color=0,/normal
            xyouts,(i*(1./n_elements(parameter)))$
                    +(.5/(n_elements(parameter))),$
                    .8,x_axis_title[oplot_index],color=0,/normal
            oplot_index = oplot_index+1
            for j=1,plot_count[i]-1 do begin      
              oplot,x[oplot_index,*],y,linestyle=j,color=0, _extra=e
              plots,[(i*(1./n_elements(parameter)))$
                    +(.25/(n_elements(parameter))),$
                    (i*(1./n_elements(parameter)))$
                    +(.48/(n_elements(parameter)))],$
                    [.81+(j*0.03),.81+(j*0.03)],linestyle=j,color=0,/normal
              xyouts,(i*(1./n_elements(parameter)))$
                     +(.5/(n_elements(parameter))),$
                     .8+(j*.03),x_axis_title[oplot_index],color=0,/normal
              oplot_index=oplot_index+1
            endfor        
          endelse 
         endfor 
      endelse
    endif
  endelse ; end overplots if..then..else
;
;  If using oo graphics and a variable was provided, return plot to that var
;
  if( ~keyword_set(directgraphic) and $
      arg_present(plot_object) )then plot_object=plot1

end