pro bp_sh3,instr,p

; v1 2/9/4 implements SH3 bad pixel algorithm on MRI (176 by 254), 
;    SLI (128 by 254) and HRI (160 by 254)
;    image array indices are based on Tables 5.2-2,3,4 in the image calibration
;    document.  The row indices in that table refer to physical rows on the 
;    CCD, which are displaced by -1 from the row numbers in an MRI, SLI or HRI image
;    because they add a row below (as well as above).  Columns preserve
;    their indices.

if instr eq 'MRI' then begin
  p(0,0:253)=p(13,0:253)
  p(1,0:112)=p(13,0:112)
  p(27,162:163)=p(25,162:163)
  p(28,162:163)=p(25,162:163)
  p(172,193)=p(170,193)
endif
if instr eq 'SLI' then begin
  p(87,73:253)=p(86,73:253)
  p(0,0:7)=p(122,0:7)
  p(0,34:253)=p(3,34:253)
  p(1,93:253)=p(3,93:253)
  p(2,151:253)=p(3,151:253)
  p(76,62:63)=p(74,62:63)
  p(77,62:63)=p(74,62:63)
  p(1,0:6)=p(122,0:6)
  p(2,0:5)=p(122,0:5)
  p(3,0:4)=p(122,0:4)
  p(4,0:4)=p(122,0:4)
  p(5,0:4)=p(122,0:4)
  p(6,0:3)=p(122,0:3)
  p(7,0:3)=p(122,0:3)
  p(8,0:3)=p(122,0:3)
  p(9,0:2)=p(122,0:2)
  for i=10,86 do begin
    p(i,0:1)=p(122,0:1)
  endfor
  p(87:121,0)=p(122,0)
endif
if instr eq 'HRI' then begin
  p(0,34:78)=p(2,34:78)
  p(0,131:196)=p(3,131:196)
  p(0,197:253)=p(4,197:253)
  p(1,44:74)=p(2,44:74)
  p(1,183:253)=p(4,183:253)
  p(2,183:217)=p(4,183:217)
  p(3,202:216)=p(4,202:216)
endif

return
end

; old MRI
; p(1,0:111)=p(13,0:111)
; p(27,161:162)=p(25,161:162)
; p(28,161:162)=p(25,161:162)
; p(172,192)=p(170,192)

; old SLI
; p(87,72:253)=p(86,72:253)
; p(0,0:6)=p(122,0:6)
; p(0,33:252)=p(3,33:252)
; p(1,92:252)=p(3,92:252)
; p(2,150:252)=p(3,150:252)
; p(76,61:62)=p(74,61:62)
; p(77,61:62)=p(74,61:62)
; p(1,0:5)=p(122,0:5)
; p(2,0:4)=p(122,0:4)
; p(3,0:3)=p(122,0:3)
; p(4,0:3)=p(122,0:3)
; p(5,0:3)=p(122,0:3)
; p(6,0:2)=p(122,0:2)
; p(7,0:2)=p(122,0:2)
; p(8,0:2)=p(122,0:2)
; p(9,0:1)=p(122,0:1)
; for i=10,86 do begin
;   p(i,0)=p(122,0)
; endfor

; old HRI
; p(0,33:77)=p(2,33:77)
; p(0,130:195)=p(3,130:195)
; p(0,196:253)=p(4,196:253)
; p(1,43:74)=p(2,43:74)
; p(1,182:253)=p(4,182:253)
; p(2,182:216)=p(4,182:216)
; p(3,201:215)=p(4,201:215)

; old SLI