PRO PlottOverviewMRT1
sunfile ='H:\MyDocs\Callisto\IHY\Korea\Essay\PlotOV\OVS_MRT_20090427_083402.prn'
sundata = read_ascii(sunfile, DATA_START=1, DELIMITER=";")
freq = reform(sundata.field1[0,*])
ysun = reform(sundata.field1[1,*])
reffile = 'H:\MyDocs\Callisto\IHY\Korea\Essay\PlotOV\OVS_MRT_20090427_083957.prn';T0
refdata = read_ascii(reffile, DATA_START=1, DELIMITER=";")
freq = reform(refdata.field1[0,*])
yref = reform(refdata.field1[1,*])
title = 'Callisto spectral overview Mauritius Radio Telescope (vertical log-per)'
xtitle = 'Frequency [MHz]'
Ytitle = 'Intensity [mV]'
Filename = 'H:\MyDocs\Callisto\IHY\Korea\Essay\PlotOV\MRTovplot100pwm.eps'
freqrange = [0,870];frequency range to plot
freqstep = 50; (50) stepsize frequency-axis
dbrange = [-5,45];y-axis range
;set_plot,'PS'
;DEVICE, /ENCAPSULATED, BITS_PER_PIXEL=8, /COLOR, FILENAME=filename,xsize=20, ysize=10
plot,freq,ysun,$
xrange=freqrange,xtickinterval=freqstep,$
yrange=[1200,2500],ytickinterval=100,title=title,xtitle=xtitle,ytitle=ytitle
loadct,5
oplot,freq,ysun,color=133
oplot,freq,yref,color=33
xyouts,100,2300,'Sky (red)' ,color=133,charthick=3, charsize=1
xyouts,400,2300,'50 Ohms (blue)',color=33 ,charthick=3, charsize=1
;device,/close
;set_plot,'Win'
image3d = TVRD(TRUE=1)
WRITE_JPEG, 'H:\MyDocs\Callisto\IHY\Korea\Essay\PlotOV\MRTovplot100pwm.jpg', image3d, TRUE=1, QUALITY=100
end
In the above program you need to edit pathnames and filenames suitable for local situation. Similar examples can be downloaded form here:
http://www.astro.phys.ethz.ch/instrument/callisto/ecallisto/applidocs.htm