#!/bin/sh
#
# Scans SCSN formatted output and converts it to PDPC input.
# http://www.data.scec.org/ftp/catalogs/SCSN/
# Usage: "create_catalog_anss.sh < input > output"

awk '
{       year  = $1*1.0
        month = $2*1.0
        day   = $3*1.0
        time  = ($4*3600) + ($5*60) + $6
        lat   = substr($0,25,3) + substr($0,29,5) / 60.0
        lon   = substr($0,34,4) - substr($0,39,5) / 60.0
        qual  = substr($0,45,1)
        mag   = $11*1.0
        dep   = $12*1.0
}
(qual == " ") {qual = "Z"}
(lat != 0) && (lon != 0) {print year,month,day,lon,lat,mag}
'