pro gencompeit,nxx,nyy,rad

;  routine to generate a multiplicative array rad
;  of size nxx,nyy intended to compensate the radial
;  fall-off of the EIT data.  Assumes 
;  an r^14 multiplicative function

;  coordinates for sun center derived from least-squares
;  fit to the limb
xc = 154.462 & yc = 490.518 & rsun = 381.728

;  radius array
rad = fltarr(nxx,nyy)

for j = 0,nyy-1 do for i = 0,nxx-1 do begin
	rad(i,j) = sqrt( (float(i)-xc)^2 + (float(j)-yc)^2)/rsun
endfor

;  get multiplicative factor
whr = where (rad gt 1.)
rad(whr) = rad(whr)^14.

; set remainder to 1
whr = where(rad le 1.)
rad(whr) = 1.0
rad(whr) = 2.5            ;;;;  Alice's mod


end