#!/bin/sh
#
# Scans SCEDC formatted output and converts it to PDPC input.
# http://www.data.scec.org/catalog_search/date_mag_loc.php
# Usage: "create_catalog_anss.sh < input > output"

awk '
{	year  = substr($0, 1,4) * 1.0
	month = substr($0, 6,2) * 1.0
	day   = substr($0, 9,2) * 1.0
	hour  = substr($0,12,2) * 1.0
	min   = substr($0,15,2) * 1.0
	sec   = substr($0,19,5) * 1.0
	et    = substr($0,24,2)
	mag   = substr($0,27,4) * 1.0
	mt    = substr($0,32,1)
	lat   = substr($0,36,7) * 1.0
	lon   = substr($0,44,8) * 1.0
	dep   = substr($0,53,5) * 1.0
	qual  = substr($0,59,1)
	evid  = substr($0,61,8) * 1.0
	nph   = substr($0,70,3) * 1.0
        ngrm  = substr($0,75,4) * 1.0
}
(et   == " ") {et   = "uk"}
(mt   == " ") {mt   = "n"}
(qual == " ") {qual = "Z"}

(lat != 0) && (lon != 0) {print year,month,day,lon,lat,mag}
'