function qs_norh, radius, center = center

if n_elements(center) ne 2 then center = [256.5, 256.5]

widget_control, /hour

xt=(yt=fltarr(300))

;openr, lun, '/home/grechnev/idl_lib/grlib/rdmodel1.dat', /get
openr, lun, 'c:\rsi\idl52\lib\istp\rdmodel1.dat', /get
readu, lun, xt, yt
free_lun,lun

if !version.OS ne 'windows' and !version.OS ne 'Win32' then byteorder, xt, yt, /lsw

M=300
dM=15

pix=4.911

yt=[yt, fltarr(dM)]

x=fltarr(512, 512)

solr=radius/pix
solrad=solr*1.2

r=(findgen(M+dM)+1)*solrad/M

xt=r/solr

        for j=0, M+dM-1 do begin

N=fix(4*!pi*r(j) > 72*4)

t=findgen(N)/(N-1)*2*!pi

xx=cos(t)*r(j)+Center(0)
yy=sin(t)*r(j)+Center(1)

x(xx, yy)=yt(j)

        endfor

x(256-5:256+5, 256-5:256+5) = median(x(256-5:256+5, 256-5:256+5), 3)

return, smooth(x,3)


	end