Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~rminchin/idl/masdpsall.pro Дата изменения: Tue Oct 2 22:29:01 2012 Дата индексирования: Sun Feb 3 19:01:24 2013 Кодировка: Поисковые слова: earthquakes

pro masdpsall,flist,projid,date,b,srcOn=srcOn,srcOff=srcOff,bpon=bpon,bpoff=bpoff,src_onoff=src_onoff,side=side,flux=flux,dateroll=dateroll,verbose=verbose

; check if flist (if supplied as a list) s a string variable, if not
; assume it is a numeric list and add .fits.

if (n_elements(flist) EQ 4) then begin
if (size(flist,/type) NE 7) then begin
fnums = flist
flist = strarr(4)
for i = 0, 3 do begin
flist[i] = string((fnums[i]),format='(I05)')+'.fits'
endfor
endif
endif

;print,flist

;flist=strarr(4)
;for i = 0, 3 do begin
; flist[i]=string(fnum+(100*i),format='(I05)')+'.fits'
;endfor

if(n_elements(dateroll) EQ 0) then dateroll = intarr(4)
if(n_elements(date) EQ 4) then dodates = 0 else dodates = 1

if(n_elements(flist) EQ 1) then begin
fnum = flist
flist =strarr(4)
masmakeflist,projid,date,fnum,flist,dateroll=dateroll,verbose=verbose
endif


if(n_elements(side) NE 1) then side = 1

for i = 0, 13 do begin

if (dodates) then begin
fbase='/share/pdata'+string(i+1,format='(I0)')+'/pdev/'+projid+'.'+string(date,format='(I8)')+'.b'+string(i-7*fix(i/7),format='(I1)')+'s'+string(fix(side),format='(I1)')+'g'+string(fix(i/7),format='(I1)')+'.'

fbase2='/share/pdata'+string(i+1,format='(I0)')+'/pdev/'+projid+'.'+string(date+1,format='(I8)')+'.b'+string(i-7*fix(i/7),format='(I1)')+'s'+string(fix(side),format='(I1)')+'g'+string(fix(i/7),format='(I1)')+'.'

endif

; masonoff2,flist,bout,fbase=fbase
;print,fbase,flist

files = strarr(4)
for j = 0, 3 do begin
if (dodates EQ 0) then begin
files[j] = '/share/pdata'+string(i+1,format='(I0)')+'/pdev/'+projid+'.'+string(date[j],format='(I8)')+'.b'+string(i-7*fix(i/7),format='(I1)')+'s'+string(fix(side),format='(I1)')+'g'+string(fix(i/7),format='(I1)')+'.'+flist[j]
endif else if (dateroll[j]) then begin
files[j] = fbase2+flist[j]
endif else begin
files[j] = fbase+flist[j]
endelse
endfor

; is=masposonoff(files,bout,sclcal=sclcal,sclJy=sclJy,median=median)

; masdps,flist,bout,fbase=fbase,flux=flux

is=masdpsp(files,bout,srcOn=srcOnout,srcOff=srcOffout,bpon=bponout,bpoff=bpoffout,src_onoff=src_onoffout,flux=flux)

if (i eq 0) then begin
b = replicate(bout,14)
srcOn = replicate(srcOnout[0],14,n_elements(srcOnout))
srcOff = replicate(srcOffout[0],14,n_elements(srcOffout))
bpon = replicate(bponout[0],14,n_elements(bponout))
bpoff = replicate(bpoffout[0],14,n_elements(bpoffout))
src_onoff = replicate(src_onoffout[0],14,n_elements(src_onoffout))
nsum=masfilesum(files[0],fsum,/list)
nrows = fsum.nrows
endif
b[i] = bout
srcOn[i] = srcOnout
srcOff[i] = srcOffout
bpon[i] = bponout
bpoff[i] = bpoffout
src_onoff[i] = src_onoffout

print,'done '+string(i+1,format='(I0)')+'/14'

endfor

restmulti = !p.multi

!p.multi = [0,4,4]

for i = 0, 13 do masplot,b[i]

!p.multi=restmulti

wait,5

masaddwgt,b,nrows
masplotall,b

return
end