; TASK Identify the orders

datestr = '201130'
echo=1
NCP = nsxng_load_ncp()
redo = 0

c = nsxng_constants()
ncp = nsxng_load_ncp()

; Eris
fn = fnsets( datestr, 45, 2, ok)
nsxng_read_image, IA, echo, file=fn[1]
nxnyn, *(IA.image), nc, nr
plotsub, *(IA.image), bin=2
color = [x11rgb('yellow'), x11rgb('red')]
for nso = 3, 7 do begin
   if nso eq 7 then cols = indgen(nc/2) else cols = indgen(nc)
   for mode = 1, 2 do begin
      rows = nsxng_find_real_image_row(mode, cols, nso)
      oplot, cols, rows, color=color[mode-1]
      if mode eq 1 then xyouts, cols[10], rows[10], string(nso, format='(I1)'), color=color[1]
   endfor
end
tv2png, 'fig/1_id_orders_'+IA.root+'.png', /tr, /co

; HD19061
fn = fnsets( datestr, 93, 1, ok)
nsxng_read_image, IA, echo, file=fn[0]
nxnyn, *(IA.image), nc, nr
plotsub, *(IA.image), bin=2
color = [x11rgb('yellow'), x11rgb('red')]
for nso = 3, 7 do begin
   if nso eq 7 then cols = indgen(nc/2) else cols = indgen(nc)
   for mode = 1, 2 do begin
      rows = nsxng_find_real_image_row(mode, cols, nso)
      oplot, cols, rows, color=color[mode-1]
      if mode eq 1 then xyouts, cols[10], rows[10], string(nso, format='(I1)'), color=color[1]
   endfor
end
tv2png, 'fig/1_id_orders_'+IA.root+'.png', /tr, /co

; HD7983
fn = fnsets( datestr, 33, 1, ok)
nsxng_read_image, IA, echo, file=fn[0]
nxnyn, *(IA.image), nc, nr
plotsub, *(IA.image), 200, 900, 200, 200, bin=2, /mag, zmax=80, zmin = -10
color = [x11rgb('yellow'), x11rgb('red')]
for nso = 3, 7 do begin
   if nso eq 7 then cols = indgen(nc/2) else cols = indgen(nc)
   for mode = 1, 2 do begin
      rows = nsxng_find_real_image_row(mode, cols, nso)
      oplot, cols, rows, color=color[mode-1]
      if mode eq 1 then xyouts, cols[10], rows[10], string(nso, format='(I1)'), color=color[1]
   endfor
end
nsxng_wide_profile_mash, IA, ncp
nso = 3
mode = 1
rows = nsxng_find_real_image_row(mode, cols, nso, ncp, IA.SlitOffset)
oplot, cols, rows, color=color[mode-1], line=2
mode = 2
rows = nsxng_find_real_image_row(mode, cols, nso, ncp, IA.SlitOffset)
oplot, cols, rows, color=color[mode-1], line=2

tv2png, 'fig/1_id_orders_'+IA.root+'.png', /tr, /co
;stop

; Look at shifts
seqnostart = 31
seqnoend = 120
nd =  1+seqnoend-seqnostart
if n_elements(d) ne nd or keyword_set(redo) then begin
   d = rdsets(datestr, seqnostart, nd, /img)
   for id = 0, nd-1 do begin
      img = d[id]
      nsxng_wide_profile_mash, img, ncp, /silent
      d[id] = img
   endfor
end
indxseq = where(d.object ne 'HD7983' and d.object ne 'HD6038')
SlitOffset_median = median(d[indxseq].slitoffset)
str_SlitOffset_median =  datestr+' SlitOffset_median ' +string(SlitOffset_median)
print, datestr, ' SlitOffset_median ', SlitOffset_median

a=d.object & s = a[sort(a)] & object_u = s[uniq(s)]
color_u = [255, x11rgb('DarkGreen'), x11rgb('DarkBlue'), x11rgb('DarkViolet'), x11rgb('orange')]

plot, d.airmass, d.slitoffset, ps=4, xtit='Airmass', yr=[2.2, 3.2], ytit='SlitOffset'
xyouts, !x.crange[0], mean(!y.crange),  str_SlitOffset_median
for iu = 0, n_elements(object_u)-1 do begin
   indx = where(d.object eq object_u[iu])
   oplot, d[indx].airmass, d[indx].slitoffset, ps=4, color=color_u[iu]
   xyouts, 1.2, 2.3 + iu/10.,  object_u[iu], color=color_u[iu]
endfor
oplot, !x.crange, [0,0]+SlitOffset_median
tv2png, 'fig/1_id_orders_SlitOffset_v_airmass.png', /tr, /co

seqno = indgen(nd) + seqnostart
plot, seqno, d.slitoffset, ps=4, xtit='SeqNo', yr=[2.2, 3.2], ytit='SlitOffset'
xyouts, !x.crange[0], mean(!y.crange),  str_SlitOffset_median
for iu = 0, n_elements(object_u)-1 do begin
   indx = where(d.object eq object_u[iu])
   oplot, seqno[indx], d[indx].slitoffset, ps=4, color=color_u[iu]
   xyouts, !x.crange[1], 2.3 + iu/10.,  object_u[iu], color=color_u[iu], al=1
endfor
oplot, !x.crange, [0,0]+SlitOffset_median
tv2png, 'fig/1_id_orders_SlitOffset_v_seqno.png', /tr, /co

end


;201129 SlitOffset_median        3.0408224