pro shiftplot
gauss_pfad='/data/etnas1/ulysses/ket/rate/shifttab/gauss-10.dat'
shift_pfad='/data/etnas1/ulysses/ket/rate/shifttab/shiftplot.dat'
gauss=read_ascii(gauss_pfad)
shift=read_ascii(shift_pfad)

days=gauss.field01[0,*]
a4000d1=gauss.field01[1,*]
a4000d2=gauss.field01[3,*]
a4000c2=gauss.field01[5,*]
a4000s2=gauss.field01[7,*]
p4000d1=gauss.field01[9,*]
p4000d2=gauss.field01[11,*]
p4000c2=gauss.field01[13,*]
p4000s2=gauss.field01[15,*]

sdays=shift.field1[0,*]
c2shift=shift.field1[1,*]
s2shift=shift.field1[2,*]


liniendicke = 1
schriftgroesse = 1
schriftdicke = 1

Set_Plot, 'PS'

;----------C2:---------------
Device, filename='shiftplot_c2_inv.eps', $
set_font=schriftart, $
          /encapsulated, $
          /color, $
          bits_per_pixel=24, $ 
          /landscape, $
          xsize=29.7, $
          ysize=21.0, $
          yoffset=29

plot, days-10, p4000c2*1.023, $
        /nodata, $
        xstyle=1, ystyle=1, $
        yrange=[100,125], $
        xrange=[290, 6500], $
        ;xticks=25, $
        xminor=50, $
        ytitle='C2 channel', $
        charsize=schriftgroesse, charthick=schriftdicke, $
        font=schrifttyp, $
        position=[0.05, 0.05, 0.95, 0.95], $
        psym=10, $
        /noerase, $
        xgridstyle=5, $
        ygridstyle=5, $
        color=fsc_color('black')

oplot, days-10, p4000c2*1.023, $
        color=fsc_color('darkgreen')

oplot, sdays, (-1*c2shift)+120.5, $
       color=fsc_color('red'), $
       psym=1


device, /close_file


;--------------S2:----------------
Device, filename='shiftplot_s2_inv.eps', $
set_font=schriftart, $
          /encapsulated, $
          /color, $
          bits_per_pixel=24, $ 
          /landscape, $
          xsize=29.7, $
          yoffset=29, $
          ysize=21.0

plot, days-10, p4000s2*1.0224, $
        /nodata, $
        xstyle=1, ystyle=1, $
        yrange=[120,140], $
        xrange=[290, 6500], $
        ;xticks=25, $
        xminor=50,$
        ytitle='S2 channel', $
        charsize=schriftgroesse, charthick=schriftdicke, $
        font=schrifttyp, $
        position=[0.05, 0.05, 0.95, 0.95], $
        psym=10, $
        xgridstyle=5, $
        ygridstyle=5, $
        /noerase, $
        color=fsc_color('black')

oplot, days-10, p4000s2*1.0224, $
        color=fsc_color('darkgreen')

oplot, sdays, (-1*s2shift)+115, $
       color=fsc_color('red'), $
       psym=1





device, /close_file
set_plot, 'X'
end