; load the chi-squared data from filtered lensview log
; grep out lines with "main: Starting minimis" and "main: ImgDoConjGradient returned"
; from the log file. Then load them with this.
; Works for a single value of crit radius. You should chop up the filtered
;	file to unique vel dispersions
;	surf: the big fat array of values
pro		makechiarr_dev_halo, filename, surf, kappa,axratio,dev_oa,dev_scale,halo_b_out,halo_ell_out,halo_oa_out,halo_core_out

	restore,'/home/rwayth/idl/surfacetemp_piep_halo.dat'
	chidata = read_ascii(filename, TEMPLATE = BIGTEMPLATE)
; field 7	"lensoffset"
; field 8	"x:"
; field 9	x offset
; field 10	"offsety:"
; field 11	y offset
; field 12	"Critrad:" or some mass scale
; field 13	mass_scale
; field 14	"shear/Ellipticity" or some ellipticity parameter
; field 15	shear
; field 16	"shear_angle/Orient_Angle"
; field 17	shear angle (degrees)
; field 18	"Devauc_scale:"
; field 19	scale
; field 20  "R_power:"
; field 21	sersic r factor
; field 22	"Critrad:" for halo
; field 23	halo_b
; field 24	"Ellipticity:" for halo
; field 25	ellipticity
; field 26	"Orient_Angle:"
; field 27	angle
; field 28	"Core_rad:"
; field 29	halo_core
; field 30	date (next line)
; field 36-37	"Chi Squ:"
; field 38	chi_squ
; field 39	"Entropy:"
; field 40	entrpoy
; field 40-42	"Reduced Chisqu:"
; field 43	reduced_chisqu
; field 42-44	"Num img pix:"
; field 45	num_img_pix
; field 46-48	"Num src pix:"
; field 49	num_src_pix
	xosize = size([uniq(chidata.FIELD09[sort(chidata.FIELD09)])])
	yosize = size([uniq(chidata.FIELD11[sort(chidata.FIELD11)])])
	devksize = size([uniq(chidata.FIELD13[sort(chidata.FIELD13)])])
	axrsize = size([uniq(chidata.FIELD15[sort(chidata.FIELD15)])])
	devoasize = size([uniq(chidata.FIELD17[sort(chidata.FIELD17)])])
	devscsize = size([uniq(chidata.FIELD19[sort(chidata.FIELD19)])])
	halo_b_sz = size([uniq(chidata.FIELD23[sort(chidata.FIELD23)])])
	halo_ell_sz = size([uniq(chidata.FIELD25[sort(chidata.FIELD25)])])
	halo_oa_sz = size([uniq(chidata.FIELD27[sort(chidata.FIELD27)])])
	halo_core_sz = size([uniq(chidata.FIELD29[sort(chidata.FIELD29)])])
	axrord = fltarr(axrsize[1])
	devoaord = fltarr(devoasize[1])
	devkord = fltarr(devksize[1])
	devscord = fltarr(devscsize[1])
	halo_b = fltarr(halo_b_sz[1])
	halo_ell = fltarr(halo_ell_sz[1])
	halo_oa = fltarr(halo_oa_sz[1])
	halo_core = fltarr(halo_core_sz[1])

print,'Dev kappa size: ',strtrim(devksize[1],2),', Axratio size: ',strtrim(axrsize[1],2),'. Orient ang size: ',strtrim(devoasize[1],2),', Dev scale: ',strtrim(devscsize[1],2),'. Halo b size: ',strtrim(halo_b_sz[1],2),', Halo ellip size: ',strtrim(halo_ell_sz[1],2),', Halo orient ang size: ',strtrim(halo_oa_sz[1],2),', Halo core size: ',strtrim(halo_core_sz[1],2)

	surface = fltarr(devksize[1], axrsize[1],devoasize[1],devscsize[1],halo_b_sz[1],halo_ell_sz[1],halo_oa_sz[1],halo_core_sz[1])
	for k=0L, devksize[1]-1 do begin
		devkord[k] = chidata.FIELD13[k*axrsize[1]*devoasize[1]*devscsize[1]*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1]]
		for l=0L, axrsize[1] -1 do begin
			axrord[l] = chidata.FIELD15[l*devoasize[1]*devscsize[1]*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1]]
			for i=0L,devoasize[1]-1 do begin
				devoaord[i] = chidata.FIELD17[i*devscsize[1]*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1]]
				for j=0L,devscsize[1]-1 do begin
					devscord[j] = chidata.FIELD19[j*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1]]
					for m=0L,halo_b_sz[1]-1 do begin
						halo_b[m] = chidata.FIELD23[m*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1]]
						for n=0L, halo_ell_sz[1]-1 do begin
							halo_ell[n] = chidata.FIELD25[n*halo_oa_sz[1]*halo_core_sz[1]]
							for o=0L,halo_oa_sz[1]-1 do begin
								halo_oa[o] = chidata.FIELD27[o*halo_core_sz[1]]
								for p=0L,halo_core_sz[1]-1 do begin
									halo_core[p] = chidata.FIELD29[p]
									surface[k,l,i,j,m,n,o,p] = chidata.FIELD38[ $
										k*devoasize[1]*axrsize[1]*devscsize[1]*halo_ell_sz[1]*halo_core_sz[1]*halo_oa_sz[1]*halo_b_sz[1] $
										+ l*devoasize[1]*devscsize[1]*halo_b_sz[1]* halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1] $
										+ i*devscsize[1]*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1] $
										+ j*halo_b_sz[1]*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1] $
										+ m*halo_ell_sz[1]*halo_oa_sz[1]*halo_core_sz[1] $
										+ n*halo_oa_sz[1]*halo_core_sz[1] $
										+ o*halo_core_sz[1] $
										+ p $
										]
								endfor
							endfor
						endfor
					endfor
				endfor
			endfor
		endfor
	endfor
	surf = surface
	axratio = axrord
	dev_oa = devoaord
	kappa = devkord
	dev_scale = devscord
	halo_b_out = halo_b
	halo_ell_out = halo_ell
	halo_oa_out = halo_oa
	halo_core_out = halo_core
end