#!/bin/bash
#
# Run UTEP station blocks through USGS terrain_correct

OUTPUT=UTEP.usgs

function run_tc() {
  NORTH=$1
  SOUTH=$2
  EAST=$3
  WEST=$4
  echo "# Running block $NORTH to $SOUTH, $WEST to $EAST" >> $OUTPUT
  ~/devel/gravity/GravDB/getstations ./UTEP_gravdb.pf \
    $NORTH $WEST $SOUTH $EAST | grep -v "^#" > test_stns.pf
  ~/devel/gravity/terrain/terrain_correct <<EOI
test_stns.pf

0
$SOUTH 0 $EAST 0
$NORTH 0 $WEST 0
0
15
0
0


h
EOI
  cat test_stns.new >> $OUTPUT
}

echo "# Automated processing from script by blocks" > $OUTPUT
#	North South East West
run_tc  43	35	-112	-115
run_tc  47	42	-112	-120
run_tc  47	44	-101	-110
run_tc  35	32	-109	-112
run_tc  45	42	-109	-112
