;function HHMMSS, axis, index, value ; ;hour = floor(value) ;minute = floor((value - double(hour)) * 60.0) ;sec = floor((value - double(hour) - double(minute) / 60.0) * 3600.0) ;msec = (value - double(hour) - double(minute) / 60.0 - double(sec) / 3600.0) * 3600000.0 ; ;;return, string(hour, minute, sec, msec, $ ;; format = "(i2.2, ':', i2.2, ':', i2.2, '.', d)") ; ;if msec ge 1000.0 then begin ; sec = sec + 1 ; msec = msec - 1000.0 ;endif ; ;if msec ge 500.0 then sec = sec + 1 ; ;if sec ge 60.0 then begin ; minute = minute + 1 ; sec = sec - 60.0 ;endif ;if minute ge 60.0 then begin ; hour = hour + 1 ; minute = minute - 60.0 ;endif ; ;return, string(hour, minute, sec, $ ; format = "(i2.2, ':', i2.2, ':', i2.2)") ;end ; ;;-------------------------------------------------------- ; ;function MMSS, axis, index, value ; ;minute = floor(value * 60.0) ;sec = floor((value - double(minute) / 60.0) * 3600.0) ;msec = (value - double(minute) / 60.0 - double(sec) / 3600.0) * 3600000.0 ; ;if msec ge 1000.0 then begin ; sec = sec + 1 ; msec = msec - 1000.0 ;endif ; ;if msec ge 500.0 then sec = sec + 1 ; ;if sec ge 60.0 then begin ; minute = minute + 1 ; sec = sec - 60.0 ;endif ; ;return, string(minute, sec, $ ; format = "(i.2, ':', i2.2)") ;end ; ;-------------------------------------------------- ; ; Main Program: ; ;pro dwe print, 'Output Medium: ' print, ' ' print, ' output.ps (0)' print, ' output.eps (1)' print, ' Tektronix-Screen (2)' print, ' X-Screen (3)' print, ' Windows (4)' print, ' terminate program (5)' read, medium case medium of 0: begin set_plot, 'PS' device, filename='output.ps', $ ; /landscape, $ ; xoffset=2.5, yoffset=26., $ ; xsize=24., ysize=16. /portrait, $ xoffset=2., yoffset=3., $ xsize=16., ysize=24. end 1: begin set_plot, 'PS' device, filename='output.eps', $ /encapsulated, $ xsize=24., ysize=16. end 2: begin set_plot, 'TEK' end 3: begin set_plot, 'X' ; window, 0, xsize=650.*210/297, ysize=650. ; !p.color=255 ; !p.background=0 end 4: begin window, 0, xsize=600., ysize=400. ; !p.color=255 ; !p.background=0 end else: begin print, 'Programmende' stop end endcase ; t_start and t_end determination openr, 1,'carriera.dat', error = err if err ne 0 then begin print, !err_string stop endif dat = dblarr(2) readf, 1, dat t_start = dat(0) while not eof(1) do readf, 1, dat close, 1 t_end = dat(0) print, 'Time scale:' print, ' ' print, ' t_0 = 0 (1)' print, ' Start of test = 0 (2)' print, ' scale with given time (3)' read, tscale case tscale of 1: begin print, 'Enter moment of t_0 as recorded from HK3 (in hours)' read, t_sub xtit = '!18t - t!D0!N!X [min]' end 2: begin t_sub = t_start xtit = 'TEST ELAPSED TIME [min] (T=0 <=> ' + string(floor(t_start), format = '(I0)') + ':' + string((t_start - floor(t_start)) * 60.0, format = '(D4.1)') + ' SCET)' end else: begin t_sub = 0.0 xtit = 'SCET [min]' end endcase t_start = t_start - t_sub t_end = t_end - t_sub print, 'Which time range (in minutes, remember t0!) ?' print, 'Enter minimum time (-1 for full time range)' read, mintime if mintime ne -1 then begin print, 'Enter maximum time' read, maxtime !x.range(0) = mintime !x.range(1) = maxtime endif else begin !x.range(0) = t_start * 60.0 !x.range(1) = t_end * 60.0 endelse !x.style = 1 dummy = ' ' ;if !x.range(1) - !x.range(0) le 10.0 then begin ; print, 'Connect data points by lines? (y/n)' ; read, dummy ; if dummy eq 'y' then psy = 0 else psy = 3 ;endif else psy = 3 print, 'Connect NCO-data points by lines? (y/n)' read, dummy if dummy eq 'y' then ncopsy = 0 else ncopsy = 3 print, 'Connect data AGC-points by lines? (y/n)' read, dummy if dummy eq 'y' then agcpsy = 0 else agcpsy = 3 xyouts, 0.5, 0.5, /normal, '!5' ; setzt den Font !p.multi = [0, 1, 14, 0, 0] ; legt fest, wieviele plots auf 1 seite kommen ; !p.psym = 0 ; punkte durch linien verbunden ; !p.psym = 3 ; keine linien sondern punkte !p.charsize = 1.0 !p.charthick = 1.5 print, 'Enter test name' read, dummy ptit = 'DWE parameters: ' + dummy print, 'Enter date of test' read, dummy ptit = ptit + ', ' + dummy !x.thick = 3.0 !y.thick = 3.0 ; Abfragen print, 'Which frequency range ?' print, 'Enter minimum frequency (0 for auto detect)' read, minnco if minnco ne 0 then begin print, 'Enter maximum frequency' read, maxnco ystylenco = 1 endif else begin maxnco = 0.0 ystylenco = 3 endelse print, 'Which AGC range ?' print, 'Enter minimum AGC (0 for auto detect)' print, 'recommended: 0' read, minagc if minagc ne 0 then begin print, 'Enter maximum AGC' read, maxagc ystyleagc = 1 endif else begin maxagc = 0.0 ystyleagc = 3 endelse ; NCO-Plot openr, 1,'ncoa.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'ncoa.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 16./18., 1.0, 1.0] plot, dat(0,*), dat(1,*), $ position = [0.1, 16./18., 1.0, 1.0], $ ; psym = psy, $ psym = ncopsy, $ xtitle = ' ', $ yrange = [minnco, maxnco], $ ystyle = ystylenco, $ ytitle = 'f!DR!N [Hz]', $ xtickname = replicate(' ', 30), $ xticklen = 0.03 ; gestrichelte Null-Linie oplot, [!x.crange(0), !x.crange(1)], [0.0, 0.0], $ linestyle = 2, thick = 0.5 ; AGC-Plot openr, 1,'agca.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'agca.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 14./18., 1.0, 16./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 14./18., 1.0, 16./18.], $ ; psym = psy, $ psym = agcpsy, $ xtitle = ' ', $ yrange = [minagc, maxagc], $ ystyle = ystyleagc, $ ytitle = 'AGC [dBm]', $ ; yrange = [-120.0, -100.0], $ xtickname = replicate(' ', 30), $ xticklen = 0.03 ; R_TEMP-Plot openr, 1,'r_temp.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(4,ln) openr, 1,'r_temp.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(4,1) !p.region = [0.0, 12./18., 1.0, 14./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 12./18., 1.0, 14./18.], $ psym = ncopsy, $ xtitle = ' ', $ ytitle = 'R_TEMP [C]', $ ystyle = 1, $ xtickname = replicate(' ', 30), $ yminor = 2, $ yrange = [60., 120.], $ xticklen = 0.03 oplot, dat(0,*), dat(2,*), psym = ncopsy oplot, dat(0,*), dat(3,*), psym = ncopsy ; R_LOCK-Plot openr, 1,'r_lock.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'r_lock.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) ; dat(1,*) = (dat(1,*) gt 2.4) ; setzt alle dat, die groesser als 2.4 (= 0x3C * 10.24 / 255) sind, auf 1 ; seit 22.8.96 nicht mehr noetig, macht schon das Auswerteprogramm !p.region = [0.0, 11./18., 1.0, 12./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 11./18., 1.0, 12./18.], $ xtitle = ' ', $ ytitle = 'R_LOCK', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; R_SELECT-Plot openr, 1,'r_select.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'r_select.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) ; dat(1,*) = (dat(1,*) lt 0.5) ; vertauscht 0 und 1, s. TM/TC ; seit 22.8.96 nicht mehr noetig, macht schon das Programm 'psa.c' !p.region = [0.0, 10./18., 1.0, 11./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 10./18., 1.0, 11./18.], $ xtitle = ' ', $ ytitle = 'R_SELECT', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; R_POWER-Plot openr, 1,'r_power.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'r_power.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 9./18., 1.0, 10./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 9./18., 1.0, 10./18.], $ xtitle = ' ', $ ytitle = 'R_POWER', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 !p.psym = 0 ; T_TEMP-Plot openr, 1,'t_tempa.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(4,ln) openr, 1,'t_tempa.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(4,1) !p.region = [0.0, 7./18., 1.0, 9./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 7./18., 1.0, 9./18.], $ xtitle = ' ', $ ytitle = 'T_TEMP [C]', $ ystyle = 1, $ yrange = [60., 120.], $ yminor = 2, $ xtickname = replicate(' ', 30), $ xticklen = 0.03 oplot, dat(0,*), dat(2,*) oplot, dat(0,*), dat(3,*) openr, 1,'t_tempb.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(4,ln) openr, 1,'t_tempb.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(4,1) oplot, dat(0,*), dat(1,*), linestyle = 2 oplot, dat(0,*), dat(2,*), linestyle = 2 oplot, dat(0,*), dat(3,*), linestyle = 2 ; T_BOXTEMP-Plot openr, 1,'t_box.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_box.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 6./18., 1.0, 7./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 6./18., 1.0, 7./18.], $ xtitle = ' ', $ ytitle = 'T_BOXTEMP [C]', $ ystyle = 3, $ ; yrange = [min(dat(1,*))-1.0, max(dat(1,*))+1.0], $ yrange = [5.0001, 39.9999], $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; T_N_CURRENT-Plot openr, 1,'t_n_curr.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_n_curr.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 5./18., 1.0, 6./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 5./18., 1.0, 6./18.], $ xtitle = ' ', $ ytitle = 'T_N_CURR [A]', $ ; yrange = [0.1, 0.35], $ yrange = [0.0, 0.8], $ ystyle = 3, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; T_N_CL-Plot openr, 1,'t_n_cl.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_n_cl.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 4./18., 1.0, 5./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 4./18., 1.0, 5./18.], $ xtitle = ' ', $ ytitle = 'T_N_CL', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; T_R_CURRENT-Plot openr, 1,'t_r_curr.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_r_curr.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 3./18., 1.0, 4./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 3./18., 1.0, 4./18.], $ xtitle = ' ', $ ytitle = 'T_R_CURR [A]', $ ystyle = 3, $ yrange = [0.0, 0.8], $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; T_R_CL-Plot openr, 1,'t_r_cl.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_r_cl.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 2./18., 1.0, 3./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 2./18., 1.0, 3./18.], $ xtitle = ' ', $ ytitle = 'T_R_CL', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 ; T_LOCK-Plot openr, 1,'t_locka.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_locka.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) !p.region = [0.0, 1./18., 1.0, 2./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 1./18., 1.0, 2./18.], $ xtitle = ' ', $ ytitle = 'T_LOCK', $ ystyle = 1, $ yrange = [-1., 2.], $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ xtickname = replicate(' ', 30), $ xticklen = 0.06 openr, 1,'t_lockb.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_lockb.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) oplot, dat(0,*), dat(1,*), linestyle = 2 ; T_SELECT-Plot openr, 1,'t_select.dat', error = err if err eq 0 then begin ln = long(0.0) while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(2,ln) openr, 1,'t_select.dat' readf, 1, dat close, 1 dat(0,*) = (dat(0,*) - t_sub) * 60.0 ; <=> Minuten endif else dat = dblarr(2,1) ;spawn, /noshell, ['cp', '/user2/duttaroy/IDL/allgemein/hhmmss.pro', '.'] ;form = 'hhmmss' ;form = '/aibn92_0/duttaroy/IDL/allgemein/hhmmss' !p.region = [0.0, 0./18., 1.0, 1./18.] plot, dat(0,*), dat(1,*), $ position = [0.1, 0./18., 1.0, 1./18.], $ ; xtitle = 'TIME AFTER START OF TEST [min]', $ xtitle = xtit, $ ytitle = 'T_SELECT', $ ystyle = 1, $ yticks = 3, $ ytickv = [-1., 0., 1., 2.], $ ytickname = [' ', '0', '1', ' '], $ yminor = -1, $ ; xtickname=[' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' '] ; xtickformat = '/user2/duttaroy/IDL/allgemein/hhmmss', $ ; xtickformat = form, $ ; xtickformat = '/aibn92_0/duttaroy/IDL/allgemein/hhmmss', $ ; xtickformat = strlowcase('/aibn92_0/duttaroy/IDL/allgemein/hhmmss'), $ ; xtickformat = 'HHMMSS', $ yrange = [-1., 2.], $ xticklen = 0.06, $ xcharsize = 1.3 ;spawn, /noshell, ['rm', '-f', 'hhmmss.pro'] xyouts, 0.55, 1.0, '!A' + ptit + '!N', $ /normal, $ alignment = 0.5, $ charsize = 1.3, $ charthick = 2.0 if medium eq 0 then device, /close !x.range = 0 !x.style = 1 !x.thick = 1.0 !y.range = 0 !y.style = 1 !y.thick = 1.0 !p.multi = 0 !p.psym = 0 !p.charsize = 1.0 !p.charthick = 1.0 !p.title = '' !p.region = 0 end