CATS GPR example, plot 2

    Assign ranges to every trace according to distance from
the southern survey boundary.  SEG-Y only allows integers,
so multiply the range by 1000.  Thus meter 26 becomes 26000.
    This GPR dataset was collected with a scan recorded by
time rather by distance, so some the distance between scans
is not uniform.  For instance, there are 75 scans in the
first meter and 57 scans in the second meter.  We could use
plot parameter HSCALE 1000 to plot by distance. 


sioseis << eof
procs diskin geom header diskoa prout end
diskin
   fno 26 lno 33
   ipath file59.segy end
end
diskoa
    opath data end
end
prout
   fno 1 lno 51 indices l3 l4 l11 l12 l19 l21 end
end
header
    fno 0 lno 9999 ftr 0 ltr 9999 l19 = l10 l21 = l10 end
end
plot
   srpath sunfil ann sh&tr taginc 10 wiggle 0 trpin 300
    nibs 2859 def .005  vscale 113 nsecs .03 tlines .005 .01 end
end
geom
  dfls 0 lprint 7
  fs 1 ls 1 gxp 1 1000 75 1987 end
  fs 2 ls 2 gxp 1 2000 57 2982 end
  fs 3 ls 3 gxp 1 3000 59 3983 end
  fs 4 ls 4 gxp 1 4000 67 4985 end
  fs 5 ls 5 gxp 1 5000 58 5983 end
  fs 6 ls 6 gxp 1 6000 59 6983 end
  fs 7 ls 7 gxp 1 7000 57 7982 end
  fs 8 ls 8 gxp 1 8000 59 8983 end
  fs 9 ls 9 gxp 1 9000 63 9984 end
  fs 10 ls 10 gxp 1 10000 59 10983 end
  fs 11 ls 11 gxp 1 11000 62 11984 end
  fs 12 ls 12 gxp 1 12000 62 12984 end
  fs 13 ls 13 gxp 1 13000 56 13982 end
  fs 14 ls 14 gxp 1 14000 58 14983 end
  fs 15 ls 15 gxp 1 15000 61 15984 end
  fs 16 ls 16 gxp 1 16000 60 16983 end
  fs 17 ls 17 gxp 1 17000 56 17982 end
  fs 18 ls 18 gxp 1 18000 62 18984 end
  fs 19 ls 19 gxp 1 19000 54 19981 end
  fs 20 ls 20 gxp 1 20000 53 20981 end
  fs 21 ls 21 gxp 1 21000 66 21985 end
  fs 22 ls 22 gxp 1 22000 60 22983 end
  fs 23 ls 23 gxp 1 23000 59 23983 end
  fs 24 ls 24 gxp 1 24000 58 24983 end
  fs 25 ls 25 gxp 1 25000 56 25982 end
  fs 26 ls 26 gxp 1 26000 53 26981 end
  fs 27 ls 27 gxp 1 27000 52 27981 end
  fs 28 ls 28 gxp 1 28000 61 28984 end
  fs 29 ls 29 gxp 1 29000 60 29983 end
  fs 30 ls 30 gxp 1 30000 61 30984 end
  fs 31 ls 31 gxp 1 31000 65 31985 end
  fs 32 ls 32 gxp 1 32000 68 32985 end
  fs 33 ls 33 gxp 1 33000 52 33981 end
  fs 34 ls 34 gxp 1 34000 54 34981 end
  fs 35 ls 35 gxp 1 35000 52 35981 end
  fs 36 ls 36 gxp 1 36000 54 36981 end
  fs 37 ls 37 gxp 1 37000 54 37981 end
  fs 38 ls 38 gxp 1 38000 61 38984 end
  fs 39 ls 39 gxp 1 39000 54 39981 end
  fs 40 ls 40 gxp 1 40000 52 40981 end
  fs 41 ls 41 gxp 1 41000 59 41983 end
  fs 42 ls 42 gxp 1 42000 54 42981 end
  fs 43 ls 43 gxp 1 43000 60 43983 end
  fs 44 ls 44 gxp 1 44000 73 44986 end
  fs 45 ls 45 gxp 1 45000 59 45983 end
  fs 46 ls 46 gxp 1 46000 53 46981 end
  fs 47 ls 47 gxp 1 47000 53 47981 end
  fs 48 ls 48 gxp 1 48000 48 48979 end
  fs 49 ls 49 gxp 1 49000 52 49980 end
  fs 50 ls 50 gxp 1 50000 57 50982 end
  fs 51 ls 51 gxp 1 51000 end
end
end
eof
##################################################
# script to plot so that north is on the right.
# This doesn't use plot parameter DIR LTR since it
# drops the time line annotation.  We just sort the
# data in reverse shot and trace order.
sioseis << eof1
procs sort diskin plot end
sort
   rev1 yes rev2 yes
   ipath data.geom lkey1 3 lkey2 4 opath sfile end
end
diskin
   spath sfile
   ipath data.geom  end
end
shift
    datume 600 datumv 200000 end # elevation is in mm
end
plot 
    ann fanno fanno -4 taginc 60
    srpath sunfil wiggle 0 trpin 20 
    clip .05 tlines .005 .01 
     nibs 75 def .08  vscale 150 nsecs .03 end 
 end 
end
eof1
xloadimage -r 90 sunfil &