a=rtable("itgood.txt",3)
itgood=fix(a[0,*])
tgood=reform(a[1,*])
nj=n_elements(itgood)
nz=n_elements(xya.z)
z=xya.z
;
abmz=fltarr(nz,nj)
jbmz=fltarr(nz,nj)
b2mz=fltarr(nz,nj)
epsMmz=fltarr(nz,nj)
;
for j=0,nj-1 do begin
  print,j,itgood[j]
  abmz[*,j]=xya.abmz[*,itgood[j]]
  jbmz[*,j]=xya.jbmz[*,itgood[j]]
  b2mz[*,j]=xya.b2mz[*,itgood[j]]
  epsMmz[*,j]=xya.epsMmz[*,itgood[j]]
endfor
;
;pc_read_xyaver,obj=xya
help,xya.epsMmz,epsMmz
help,xya.b2mz,b2mz
;contour,transpose(xya.b2mz),xya.t,xya.z,/fill,nlev=30
;contour,transpose(xya.epsMmz),xya.t,xya.z,/fill,nlev=30,xtit='t',ytit='z'
contour,alog(transpose(epsMmz)),tgood,z,nlev=30,/fil
;
save,file='xyaver.sav',z,tgood,b2mz,epsMmz,abmz,jbmz
spawn,"cvs add -kb xyaver.sav"
END