rm sunfil.ras 
sioseis << eof
procs diskin plot end
diskin
   ftr 1 ltr 20
   fno 9100 lno 9250 noinc 50
   ipath /home/phenk/data/Tj37/raw_segy.dir/day203-1413z.segy end
end
geom
    gxp 1 -12.5 ggx 12.5 end
end
nmo
    vtp 1500 0 end
end
plot
   ann2 GMT ann sh&tr
   trpin 20 nspace 1 fspace 20 spacei 20
   nibs 2858 srpath sunfil.ras taginc 1
   def .04 nsecs .5 vscale 2.5
   end
end
end
eof
display -rotate 90 sunfil.ras &
convert -rotate 90 sunfil.ras direct_wave.png