;+
;PROCEDURE: swe_a0_snap
;PURPOSE:
; Plots 3D snapshots in a separate window for times selected with the cursor in
; a tplot window. Hold down the left mouse button and slide for a movie effect.
;
; If housekeeping data exist (almost always the case), then they are displayed
; as text in a small separate window.
;
;USAGE:
; swe_3d_snap
;
;INPUTS:
;
;KEYWORDS:
; LAYOUT: A named variable to specify window layouts.
;
; 0 --> Default. No fixed window positions.
; 1 --> Macbook Air with Dell 1920x1200 external screen.
; 2 --> HP Z220 with twin Dell 1920x1200 screens.
; 3 --> Dell 2009W 1680x1050 with Macbook Air
;
; This puts up snapshot windows in convenient, non-overlapping
; locations, depending on display hardware.
;
; MODEL: Plot a model of the 3D data product with the test pulser on in a
; separate window. (An analytic approximation to the test pulser
; signal is used. See 'swe_testpulser_model.pro' for details.)
;
; KEEPWINS: If set, then don't close the snapshot window(s) on exit.
;
; ZRANGE: Sets color scale range. Default = [1,3000].
;
; ZLOG: Sets log color scaling. Default = 1.
;
; ARCHIVE: If set, show snapshots of archive data (A1).
;
;CREATED BY: David L. Mitchell 07-24-12
;FILE: swe_a0_snap.pro
;VERSION: 1.0
;LAST MODIFICATION: 07/24/12
;-
pro swe_a0_snap, layout=layout, model=model, keepwins=keepwins, zrange=zrange, zlog=zlog, $
archive=archive
@mvn_swe_com
if keyword_set(archive) then aflg = 1 else aflg = 0
if (aflg) then begin
if (size(swe_3d_arc,/type) ne 8) then begin
print,"No valid 3D archive data."
return
endif
endif else begin
if (size(swe_3d,/type) ne 8) then begin
print,"No valid 3D survey data."
return
endif
endelse
if keyword_set(model) then mflg = 1 else mflg = 0
if (size(swe_hsk,/type) ne 8) then hflg = 0 else hflg = 1
if keyword_set(keepwins) then kflg = 0 else kflg = 1
if (size(zrange,/type) eq 0) then zrange = [1.,3000.]
if (size(zlog,/type) eq 0) then zlog = 1
; gudsum = ['CC'X, '1E'X] ; up to early cruise
gudsum = ['C0'X, 'DE'X] ; once the new sweep tables are uploaded
nchan = [64, 32, 16] ; number of energy channels, indexed by the GROUP tag in A0
period = 2D^(dindgen(6) + 1D) ; sampling interval (sec), indexed by the PERIOD tag in A0
; Put up snapshot window(s)
Twin = !d.window
print,'Use button 1 to select time; button 3 to quit.'
if not keyword_set(layout) then layout = 0
hsk_ysize = 545
case layout of
0 : begin
window, /free, xsize=800, ysize=500
Swin = !d.window
if (hflg) then begin
window, /free, xsize=225, ysize=hsk_ysize
Hwin = !d.window
endif
if (mflg) then begin
window, /free, xsize=800, ysize=500
Mwin = !d.window
endif
end
1 : begin
window, /free, xpos=225, ypos=-600
Swin = !d.window
if (hflg) then begin
window, /free, xsize=225, ysize=hsk_ysize, xpos=1400, ypos=-550
Hwin = !d.window
endif
end
2 : begin
window, /free, xsize=800, ysize=500, xpos=1120, ypos=640
Swin = !d.window
if (hflg) then begin
window, /free, xsize=200, ysize=hsk_ysize, xpos=880, ypos=500
Hwin = !d.window
endif
if (mflg) then begin
window, /free, xsize=800, ysize=500, xpos=1120, ypos=100
Mwin = !d.window
endif
end
3 : begin
window, /free, xsize=800, ysize=500, xpos=1680, ypos=640
Swin = !d.window
if (hflg) then begin
window, /free, xsize=200, ysize=hsk_ysize, xpos=2600, ypos=500
Hwin = !d.window
endif
if (mflg) then begin
window, /free, xsize=800, ysize=500, xpos=1680, ypos=100
Mwin = !d.window
endif
end
else : begin
window, /free
Swin = !d.window
if (hflg) then begin
window, /free, xsize=225, ysize=hsk_ysize
Hwin = !d.window
endif
end
endcase
; Select the first time, then get the 3D spectrum closest that time
wset,Twin
ctime2,trange,npoints=1,/silent,button=button
if (size(trange,/type) eq 2) then begin ; Abort before first time select.
wdelete,Swin ; Don't keep empty windows.
if (hflg) then wdelete,Hwin
if (mflg) then wdelete,Mwin
wset,Twin
return
endif
if (aflg) then dt = min(abs(swe_3d_arc.time - trange[0]), iref) $
else dt = min(abs(swe_3d.time - trange[0]), iref)
if (hflg) then dt = min(abs(swe_hsk.time - trange[0]), jref)
ok = 1
while (ok) do begin
if (aflg) then begin
ddd = swe_3d_arc[iref]
pmsg = "A1"
endif else begin
ddd = swe_3d[iref]
pmsg = "A0"
endelse
n_e = ddd.n_e
x = findgen(80) + 0.5
y = findgen(n_e) + 0.5
z = ddd.data[*,0:(n_e-1)]
wset, Swin
; Put up a 3D spectrogram
title = string(pmsg, time_string(ddd.time), ddd.group, ddd.period, ddd.lut, ddd.npkt, $
format='(a2,4x,a19,4x,"G ",i1,2x,"P ",i1,4x,"LUT: ",i1,4x,"NPKT: ",i3)')
lim = {x_no_interp:1, y_no_interp:1, xrange:[0,80], yrange:[0,n_e], zrange:zrange, $
xmargin:[10,10], charsize:1.2, xtitle:'Angle Bin', ytitle:'Energy Bin', $
ztitle:'Counts', title:title, xstyle:1, ystyle:1, zlog:zlog, xticks:5, xminor:4, $
yticks:4, yminor:4}
indx = where(z eq 0.,count)
if (count gt 0L) then z[indx] = !values.f_nan
specplot,x,y,z,limits=lim
; Put up test pulser model
if (mflg) then begin
wset, Mwin
swe_testpulser_model,group=ddd.group,result=tpmod
title = 'Test Pulser Model: ' + pmsg
lim = {x_no_interp:1, y_no_interp:1, xrange:[0,80], yrange:[0,n_e], zrange:zrange, $
xmargin:[10,10], charsize:1.2, xtitle:'Angle Bin', ytitle:'Energy Bin', $
ztitle:'Counts', title:title, xstyle:1, ystyle:1, zlog:zlog, xticks:5, xminor:4, $
yticks:4, yminor:4}
specplot,x,y,tpmod.a0,limits=lim
endif
; Print out housekeeping in another window
if (hflg) then begin
wset, Hwin
csize = 1.4
nlines = 28
x1 = 0.05
x2 = 0.75
x3 = x2 - 0.12
y1 = 0.95 - (0.9/(float(nlines-1)))*findgen(nlines)
fmt1 = '(f7.2," V")'
fmt2 = '(f7.2," C")'
fmt3 = '(Z2.2)'
j = jref
if (swe_hsk[j].CHKSUM[0] eq gudsum[0]) then col0 = 4 else col0 = 6
if (swe_hsk[j].CHKSUM[1] eq gudsum[1]) then col1 = 4 else col1 = 6
k = 0
erase
xyouts,x1,y1[k++],/normal,"SWEA Housekeeping",charsize=csize
xyouts,x1,y1[k++],/normal,time_string(swe_hsk[j].time),charsize=csize
xyouts,x1,y1[k],/normal,"NPKT",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].NPKT,format='(i3)'),charsize=csize,align=1.0
k++
xyouts,x1,y1[k],/normal,"P28V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P28V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"MCP28V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].MCP28V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"NR28V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].NR28V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"MCPHV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].MCPHV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"NRV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].NRV,format=fmt1),charsize=csize,align=1.0
k++
xyouts,x1,y1[k],/normal,"P12V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P12V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"N12V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].N12V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"P5AV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P5AV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"N5AV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].N5AV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"P5DV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P5DV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"P3P3DV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P3P3DV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"P2P5DV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].P2P5DV,format=fmt1),charsize=csize,align=1.0
k++
xyouts,x1,y1[k],/normal,"ANALV",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].ANALV,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"DEF1V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].DEF1V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"DEF2V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].DEF2V,format=fmt1),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"V0V",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].V0V,format=fmt1),charsize=csize,align=1.0
k++
xyouts,x1,y1[k],/normal,"ANALT",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].ANALT,format=fmt2),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"LVPST",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].LVPST,format=fmt2),charsize=csize,align=1.0
xyouts,x1,y1[k],/normal,"DIGT",charsize=csize
xyouts,x2,y1[k++],/normal,string(swe_hsk[j].DIGT,format=fmt2),charsize=csize,align=1.0
k++
xyouts,x1,y1[k],/normal,"CHKSUM",charsize=csize
xyouts,x2,y1[k],/normal,string(swe_hsk[j].CHKSUM[1],format=fmt3),charsize=csize,align=1.0,$
color=col1
xyouts,x3,y1[k],/normal,string(swe_hsk[j].CHKSUM[0],format=fmt3),charsize=csize,align=1.0,$
color=col0
endif
; Get the next button press
wset,Twin
ctime2,trange,npoints=1,/silent,button=button
if (size(trange,/type) eq 5) then begin
if (aflg) then dt = min(abs(swe_3d_arc.time - trange[0]), iref) $
else dt = min(abs(swe_3d.time - trange[0]), iref)
if (hflg) then dt = min(abs(swe_hsk.time - trange[0]), jref)
ok = 1
endif else ok = 0
endwhile
if (kflg) then begin
wdelete, Swin
if (hflg) then wdelete, Hwin
if (mflg) then wdelete, Mwin
endif
wset, Twin
return
end