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 medium = 0 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 set_plot, 'WIN' ; window, 0, xsize=600., ysize=400. ; !p.color=255 ; !p.background=0 end else: begin print, 'Program terminated' stop end endcase !x.thick = 3.0 !y.thick = 3.0 !p.charsize = 1.5 !p.charthick = 2.0 xyouts, 0.5, 0.5, /normal, '!5' ; set font to font-no.5 ; set font to font-no.5 dummy = ' ' fn = 'recovtest_rw0.out' ln = long(0.) col = 16 openr, 1, fn while not eof(1) do begin readf, 1, dummy ln = ln + 1 endwhile close, 1 dat = dblarr(col, ln) openr, 1, fn readf, 1, dat close, 1 fn = 'recovtest_rw-50.out' datoff = dblarr(col, ln) openr, 1, fn readf, 1, datoff close, 1 ;for i=1, 14 do begin ; case i of ; 1: ytit = 'Residual gamma_DWC' ; 2: ytit = 'Residual cos(gamma_DWC)' ; 3: ytit = 'Residual gamma_C' ; 4: ytit = 'Residual cos(gamma_C)' ; 5: ytit = 'Residual v_C * cos(gamma_C)' ; 6: ytit = 'Residual Gamma_des' ; 7: ytit = 'Residual cos(gamma_des)' ; 8: ytit = 'Residual v_des * cos(gamma_des)' ; 9: ytit = 'Residual Gamma_NS' ; 10: ytit = 'Residual cos(gamma_NS)' ; 11: ytit = 'Residual v_NS * cos(gamma_NS)' ; 12: ytit = 'Residual v_EW * cos(gamma_DWC)' ; 13: ytit = 'Residual v_EW' ; 14: ytit = 'Residual lon_H' ; endcase ; plot, dat[0,*], datoff[i,*]-dat[i,*], $ ; title = fn, $ ; xtitle = '!18t - t!D0!N!X [min]', $ ;; ytitle = 'Longitudinal Offset [deg]', $ ; ytitle = ytit, $ ; ystyle = 3, $ ; xstyle = 1 ;endfor plot, dat[0,*], datoff[5,*]/datoff[2,*]-dat[5,*]/dat[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Residual: v_C * cos(gamma_C) / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 plot, dat[0,*], datoff[15,*]/datoff[2,*]-dat[15,*]/dat[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Residual: v_LS / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 plot, dat[0,*], datoff[5,*]/datoff[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Retrieval: v_C * cos(gamma_C) / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 plot, dat[0,*], datoff[15,*]/datoff[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Retrieval: v_LS / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 plot, dat[0,*], (datoff[15,*] + datoff[5,*]) / datoff[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Retrieval: (v_LS + v_C * cos(gamma_C)) / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 plot, dat[0,*], (datoff[15,*] + datoff[5,*]) / datoff[2,*] - (dat[15,*] + dat[5,*]) / dat[2,*], $ title = fn, $ xtitle = '!18t - t!D0!N!X [min]', $ ytitle = 'Residual: (v_LS + v_C * cos(gamma_C)) / cos(gamma_DWC)', $ ystyle = 3, $ xstyle = 1 if (medium eq 0 or medium eq 1) then begin device, /close ; close 'output.ps' print, 'Rename "output.ps" or "output.eps" before you call this program again !!' endif end