pro write_reg,fname,x,y,r,h=h

	if n_params(0) eq 0 then begin
	  print,'write_reg,fname,x,y,r'
	  return
	endif

	openw,1,fname
	for i=0,n_elements(h)-1 do printf,1,h(i)
	if n_elements(r) gt 1 then s = reverse(sort(r)) else s=0
	for i = 0,n_elements(x)-1 do begin
	  n = 1
	  if r(s(i)) gt 0 then a = ' ' else begin
	    a = '-'
	    n = -1
	  endelse
	  printf,1,a,'CIRCLE(',strn(x(s(i))),',',strn(y(s(i))),',',$
	    strn(n*r(s(i))),')'
	endfor
	close,1
	return
	end