function mk_ideal_shd, njmx, rawx, rawy
;
;  Make shadowgram in idealized grid
;
;  2004-08-31/NJW
;

if n_params() lt 1 then begin
   print,'Syntax: mk_ideal_shd, njmx, rawx, rawy'
   return
endif

sj = strtrim(njmx,2)

; Open the pixel redistribution matrix
   
fxbopen, unit, '/right_place/JMX'+sj+'-PRED-MOD_001.fits', 1, access='R'


ideal_shd = fltarr(256,256)
      
for j = 0L, 255L do begin
for i = 0L, 255L do begin
   index = j * 256L + i
   fxbread, unit, num, 1, index+1
   fxbread, unit, idx, 2, index+1
   fxbread, unit, wei, 3, index+1
   if num ne n_elements(idx) then begin
      print,' Problem for index: ', index
   endif
   for k = 0, num-1 do begin
      l = idx(k) mod 256L
      m = idx(k) / 256L
      ideal_shd(l,m) = ideal_shd(l,m) + rshd(i,j)*wei(k)
   endfor
endfor
endfor
      
fxbclose, unit

return, ideal_shd

end