; this originally was /data/eric/OLD/kingkong-eric/fernando-dr/setup

sc=readtext('data/szcirc.txt')


ncma=readtext('data/circles/ms-A5p5_1.txt',/quiet)
ncmb=readtext('data/circles/ms-B5p5_1.txt',/quiet)
ncmc=readtext('data/circles/ms-C5p5_1.txt',/quiet)
ncmd=readtext('data/circles/ms-D5p5_1.txt',/quiet)
ncme=readtext('data/circles/ms-E5p5_1.txt',/quiet)
ncmf=readtext('data/circles/ms-F5p5_1.txt',/quiet)
ncmg=readtext('data/circles/ms-G5p5_1.txt',/quiet)
ncmh=readtext('data/circles/ms-H5p5_1.txt',/quiet)
ncmi=readtext('data/circles/ms-I5p5_1.txt',/quiet)
ncmj=readtext('data/circles/ms-J5p5_1.txt',/quiet)

nema=readtext('data/ellipse/ms-A5p5_1.txt',/quiet)
nemb=readtext('data/ellipse/ms-B5p5_1.txt',/quiet)
nemc=readtext('data/ellipse/ms-C5p5_1.txt',/quiet)
nemd=readtext('data/ellipse/ms-D5p5_1.txt',/quiet)
neme=readtext('data/ellipse/ms-E5p5_1.txt',/quiet)
nemf=readtext('data/ellipse/ms-F5p5_1.txt',/quiet)
nemg=readtext('data/ellipse/ms-G5p5_1.txt',/quiet)
nemh=readtext('data/ellipse/ms-H5p5_1.txt',/quiet)
nemi=readtext('data/ellipse/ms-I5p5_1.txt',/quiet)
nemj=readtext('data/ellipse/ms-J5p5_1.txt',/quiet)

w=where(ncma(0,*) ge 0.2*60.0 and ncma(0,*) le 10.0*60.0)

i=0 & ncd=fltarr(6,10) & e=0
m=ncma & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmb & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmc & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmd & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncme & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmf & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmg & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmh & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmi & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=ncmj & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
i=0 & e=1
m=ncma & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmb & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmc & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmd & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncme & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmf & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmg & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmh & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmi & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=ncmj & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
i=0 & e=2
m=nema & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemb & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemc & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemd & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=neme & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemf & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemg & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemh & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemi & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
m=nemj & ncd[e,i]=exp(mean(alog(m(1,w)/m(0,w)))) & i++
i=0 & e=3
m=nema & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemb & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemc & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemd & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=neme & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemf & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemg & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemh & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemi & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++
m=nemj & ncd[e,i]=exp(mean(alog(m(2,w)/m(0,w)))) & i++

ncd(4,*)=reverse([1.5,1.7,2.6,3.1,3.6,4.7,6.2,7.0,8.7,10.2]); circle masses
ncd(5,*)=reverse([1.4,1.8,2.9,3.8,5.0,6.5,8.7,10.6,12.7,14.6]); ellipse mass

