pro OplotCompareInsitu, t1, f1, density=density, speed=speed, bradial=bradial, btang=btang, _extra=_extra

on_error, 2							; On error, return to caller
COMPILE_OPT hidden

;+
; NAME:
;	OplotCompareInsitu
; PURPOSE:
;	Adds tomographic velocity timeseries to a plot made by href=PlotCompareInsitu=
; CATEGORY:
;	WWW
; CALLING SEQUENCE:
;	OplotCompareInsitu, t1, ff1
; INPUTS:
;	file		scalar or array[*]; type string
;					tomographic files containing velocity data
; OPTIONAL INPUT PARAMETERS:
; OUTPUTS:
;	(none)
; CALLS:
;	ReadEarthFile, PlotCurve
; PROCEDURE:
; MODIFICATION HISTORY:
;	MAR-2000, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
;-

loadct, 12

colors = [30,163,89]
colors = [40,125,190]
colors = [190,125,40]

linestyle=[0,0,0]

x0t = 0.53	; 0.65	;0.4
x0c = 0.60
y0t = 0.90
y0c = 0.15
xw  = 0.12
yh  = 0.075
xs  = 0.30


n = size(f1)
n = n[n[0]+2]/(n[1]*n[2])
k = 0*keyword_set(speed) + 1*keyword_set(density) + 2*keyword_set(bradial) + 3*keyword_set(btang)


; Assume the first series has been plotted with plotcompareinsitu

for i=1,n-1 do begin

	PlotCurve, /oplot, t1, f1[*,k,i], color=colors[i], linestyle=linestyle[i], _extra=_extra

	dx = !x.crange[1]-!x.crange[0]
	dy = !y.crange[1]-!y.crange[0]

	b = !x.crange[0]+x0t*dx
	w = xw*dx
	h = !y.crange[0]+y0t*dy-i*yh*dy
	plots, b+[0.,w],h*[1.,1.], color=colors[i], _extra=_extra
	xyouts, b+(1+xs)*w,h, 'time dep.'+strcompress(i), _extra=_extra

endfor


return  &  end