Back to grdout example 1

#----------------------------------------------------------
# generates grd output
#----------------------------------------------------------

#! /bin/csh

set plt = sioseis       #select which plot to use, pltsegy or sioseis
set infile = pregrd.out
set outfile = togrd.grd
/bin/rm siofil sunfil stack.ps $outfile

@ rp1 = $1
@ rp2 = $2
@ nrp = $rp2 - $rp1
@ ninc = 1

sioseis << eof
procs diskin prout plot diskoa grdout end
#debug procs diskin prout end
diskin
   fno $rp1 lno $rp2
   noinc $ninc
   ipath $infile
   set 2.54 6.5
end 
end

prout
 fno 1 lno 9999999 noinc 100
end 
end

diskoa
 opath togrd.out
end
end

grdout
opath $outfile
#xmin $rp1 xmax $rp2 xinc $ninc xunits rps
xmin 0 xmax 63.125 xinc 0.0125 xunits km
ymin 2.54 ymax 6.5 yinc 0.012
 
end 
end
  

plot
  stime 2.5 nsecs 4
  tlines 0.5 1 
  nibs 7224  vscale 2.5 
  wiggle 0 
!  anninc 15 anntyp 5   !gmt time
   ann rpno  taginc 100 !rp number
!  wiggle 30 scalar 1.0e-2 def 1.5 trpin 4 
  scalar 8.0e-4  clip 0.03 trpin 200 !gains type 5 alpha 1.5 with wbt
!  colors -.01 red -.005 magenta 0 white .003 cyan .007 blue    !plots colors
  pctfil 100
  opath siofil
!  dir ltr
  srpath sunfil end
end

end

eof