Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~tghosh/software/integrate.pro Дата изменения: Mon Sep 20 22:01:45 2004 Дата индексирования: Sun Apr 10 13:04:49 2016 Кодировка: Поисковые слова: rocket

pro integrate, sfile, sbc, out_ar, integ
; takes a normalized spectrum from structure=sfile, subcorrelator=sbc, and
; integrates between between x-limits defined by 2 cursor clicks.
; Answer in opt_depth*km/sec
;
; To call type, integrate,in_struct,#sbc, varn
sbc=sbc-1
y=sfile.(sbc).d[*]
y=-1.*(y-1.)
y=(1./(1.-y))
y=alog(y)
plot,y
print,'Define the region to integrate over by two cursor clicks.'
cursor,x0,y0,/down
cursor,x1,y1,/down
integ=total(y[x0:x1])
vel=corfrq(sfile.(sbc).h,/retvel)
velperch=(vel(x0)-vel(x1))/(x1-x0)
integ = integ*velperch
print, 'integral = ', integ
out_ar=y
end