open(TABLEA, "mcp_list");
@tablea = <TABLEA>;
# Specify the correspoding output file
open(OUT,">output_inversion");


for ($i = 0; $i < @tablea; $i++) {
    chomp ($tablea[$i]);
    ($mcpFile) = (split /\s+/,$tablea[$i])[0];
    
    system("wc $mcpFile > crap");
    open(TABLEB,'crap');
    @tableb = <TABLEB>;
    chomp ($tableb[0]);
    ($count) = (split /\s+/,$tableb[0])[1];
    $numObs = $count - 9;
    close(TABLEB);
    unlink('crap');
    #print $mcpFile," ",$numObs,"\n";
    
    $numLines = int($numObs/6);
    $remainder = $numObs - ($numLines*6);
    if ($numLines eq 0) {
        $numLines = $numLines + 1;
    }
    #print $numLines," ",$remainder,"\n";
    # Now begin with the output file
    open(TABLEB, $mcpFile);
    @tableb = <TABLEB>;
    for ($j = 0; $j < @tableb; $j++) {
        chomp ($tableb[$j]);
        ($PDE,$year,$month,$day,$hour,$minute,$second,$eqlat,$eqlong,$eqdepth,$mag) = (split /\s+/,$tableb[$j])[0,1,2,3,4,5,6,7,8,9,11];
        if ($PDE eq "PDE") {
            printf OUT "%2d%2d%2d%2d%2d %s %s%7.3f %s%8.3f %s%6.2f %s %s %s %s \n", $year%100,$month,$day,$hour,$minute,$second,"0.00",$eqlat,"0.00",$eqlong,"0.00",$eqdepth,"0.00",$numObs,$mag, "0.00", "\n";
        }
for ($k = 0; $k < @tableb; $k++) {
            chomp ($tableb[$k]);
            ($netsta, $delay_time) = (split /\s+/,$tableb[$j])[1,9];
            ($net, $sta) = (split /\./, $netsta)[0,1];
            print $net, " ", $sta, "\n";