#!/bin/sh
#
# Scans NCEDC formatted output and converts it to PDPC input.
# http://www.ncedc.org/ncedc/catalog-search.html
# 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
        lat   = $3*1.0
        lon   = $4*1.0
        mag   = $6*1.0
}
(lat != 0) && (lon != 0) {print year,month,day,lon,lat,mag}
'