This file describes how to insert shotpoint coordinates into a segy file
with ministreamer data. Note that the inserted coordinates are for the 
shotpoints, not the reflections point (CMP). Also note that if the water 
depth is not very large compared to the offset, it might be necessary to
make NMO corrections before interpretation.

First we check the content of the segy file:

segyread tape=4392.sgy endian=0 | surange
1+0 records in
6+1 records out
3200 bytes (3.2 kB) copied, 7.8e-05 seconds, 41.0 MB/s
2898 traces:
fldr     4392 5840 (4392 - 5840)
tracf    1 2 (1 - 2)
trid     1
nvs      1
nhs      1
scalco   1
sx       -49
gx       0 10 (0 - 10)
counit   1
ns       6000
dt       1000
hcf      250
year     9
day      148 150 (148 - 150)
hour     0 23 (12 - 1)
minute   0 59 (46 - 2)
sec      0 59 (19 - 15)
timbas   1

This shows that the first and last record numbers are 4392 and 5840.
There are 2 channels (tracf) in each record.

From the UKOOA file we find that the first and last shotpoints are
47 and 1451.

Comparing with the log file we find that shotpoint 47 is record 4436
and shotpoint 1451 is record 5840.
The difference is 4389 = 4436-47 = 5840-1451.

We now pick out the shot numbers and UTM coordinates from the UKOOA file:
cut -c21-25,47-54,56-62 < TOPO-SCANDIA-DEEP-2009-line1-UKOOA.txt | tr "." " " > SP_E_N.dat

Edit the file SP_E_N.dat: remove lines before the first shot (47) and
after the last shot (1451).

Convert the shot coordinates to binary form:
awk '{print $2, $3}' SP_E_N.dat | a2b n1=2 > geom.bin

Generate segy binary header record:
segyhdrs ns=6000 dt=1000 lino=1 ntrpr=1 format=1 mfeet=1 hfile=/dev/null

We are now read to make a new segy file with shotpoint coordinates:
segyread tape=4392.sgy endian=0 hfile=/dev/null bfile=/dev/null | \
suwind key=tracf min=2 | \                 # pick channel 2
suchw key1=fldr key2=fldr a=-4389 b=1 | \  # change record number to shot number
suwind key=fldr min=47 max=1451 | \        # pick shot 47 to 1451
sushw key=sx,sy infile=geom.bin | \        # insert shotpoint coordinates
segywrite tape=Line1.segy endian=0         # write new segy file

This script is the syntese of the steps above

#!/bin/bash -x
#This script inserts shotpoint coordinates into a segy file
segyread tape=5000.sgy endian=0 hfile=/dev/null bfile=/dev/null | \
# pick channel 2
suwind key=tracf min=2 | \
# change record number to shot number
suchw key1=fldr key2=fldr a=-1836 b=1 | \
# pick shot 3170 to 4201
suwind key=fldr min=3170 max=4201 | \
# insert shotpoint coordinates
sushw key=sx,sy infile=geom.bin | \
# write new segy file
segywrite tape=Line3.segy endian=0