PRO ADIlines

;Use the screen as an output
SET_PLOT,'ps'
DEVICE, FILE='Errors.ps', /COLOR
;Setup screen and plot positions
!P.POSITION=[0.1,0.3,0.9,0.9]
!P.COLOR=1   
!P.THICK=3.0
!P.FONT=1
!X.THICK=3.0
!Y.THICK=3.0
r = [255,0,0,0,0,0,66,184,245,255,255,255,230,191,148,105,161,184]
g = [255,0,84,199,255,255,255,255,255,209,135,28,0,0,0,0,0,0]
b = [255,0,255,255,255,178,0,0,0,0,0,0,0,0,0,20,135,186]
TVLCT, r,g,b




;Set some constants
loops  = 5
trials = 7

;Make floating point array
err = FLTARR(loops)

;Open the input data file
OPENR, unit, 'Data', /GET_LUN

;Set some variables
mu        = 0.
temp      = 0.
blankline = ''

FOR trial = 1, trials DO BEGIN
  ;Read the acceleration parameter
  READF, unit, mu
  PRINT, 'mu = ', mu
  ;Read the data
  FOR loop = 0, loops-1 DO BEGIN
    READF, unit, temp
    err(loop) = temp
    PRINT, 'error = ', err(loop), temp
  ENDFOR
  ;draw the graph
  IF (trial EQ 1) THEN BEGIN
    PLOT, INDGEN(loops), err, XSTYLE=5, YSTYLE=5, COLOR=13,$
      YRANGE=[0., 3.0e-6], THICK=0.1
  ENDIF ELSE BEGIN
    OPLOT, INDGEN(loops), err, COLOR=13, THICK=0.1
  ENDELSE
  IF (trial LT trials) THEN READF, unit, blankline
ENDFOR




;Close and free the file
CLOSE, unit
FREE_LUN, unit



; add the axes
myticks=['0', '1', '2', '3', '4']
myticklen=-0.03
mytickvals=[0,1,2,3,4]
AXIS, XAXIS=0, XTITLE='ADI loop', TICKLEN=myticklen,$
         XTICKS=5, XTICKNAME=myticks, XTICKV=mytickvals
myticks(*)=' '
AXIS, XAXIS=1, TICKLEN=myticklen,$
         XTICKS=5, XTICKNAME=myticks, XTICKV=mytickvals

myticks=['0', '1e-6', '2e-6','3e-6']
myticklen=-0.015
mytickvals=[0.0, 1.0e-6, 2.0e-6, 3.0e-6]
AXIS, YAXIS=0, YTITLE='Error', TICKLEN=myticklen,$
         YTICKS=4, YTICKNAME=myticks, YTICKV=mytickvals
myticks(*)=' '
AXIS, YAXIS=1, TICKLEN=myticklen,$
         YTICKS=4, YTICKNAME=myticks, YTICKV=mytickvals



;Close the postscript file
DEVICE, /CLOSE

END