#!/bin/bash

# Get stations from UTEP_gravdb in box set by N,S,E,W
awk -v NORTH=$1 -v SOUTH=$2 -v EAST=$3 -v WEST=$4 \
  'BEGIN{print "# BBOX: ", NORTH, SOUTH, EAST, WEST}
   /^#/{print $0; next}
   {if($2>=WEST && $2<=EAST && $3>=SOUTH && $3<=NORTH)print}' UTEP_gravdb 
