#!/bin/sh 
LANG=C

if [ ! -f upr-sel.tsv ]; then
proj +proj=stere +lat_ts=60 +lat_0=90 +lon_0=50 +k_0=1.0 \
	+x_0=-1696000 +y_0=7268100 +ellps=GRS80 \
	< upr-ll.tsv > upr-ps.tsv
ruby -pe 'a = $_.split("\t").map{|s| s.to_f};
next if a[0] < 42_120 or a[0] > 9071_595 \
     or a[1] < 95_076 or a[1] > 11150_302 ' \
	< upr-ps.tsv > upr-sel.tsv
fi

if [ ! -f udic.psv ]; then
sqlite3 /usr/local/share/wmo9/wmo9.sqlite \
  'select idxnum, lat, lon, stnname, obsrems from vola ' \
  > udic.psv
fi

ruby /dev/stdin > upr.tsv <<EOF
ll = {}
sel = {}
File.open('udic.psv') {|fp|
  fp.each {|line|
    stn, lat, lon, name, obsrems = line.split(/\|/)
    rbsn = case obsrems when /RBSN\((\w+)\)/ then \$1 else '-' end
    ll[stn] = [lat, lon, name, rbsn]
  }
}
File.open('upr-sel.tsv') {|fp|
  fp.each {|line|
    x, y, stn, rate = line.split(/\s/)
    sel[stn] = rate
  }
}
for stn in sel.keys.sort
  lon, lat, name, rbsn = ll[stn]
  name = '(not in VolA)' unless name
  rate = sel[stn]
  printf "%s\t%s\t%s\t%4.2f\t%s\t%s\n", lon, lat, stn, rate.to_f, rbsn, name
end
EOF