#! /usr/local/bin/perl 

use XML::Simple;

$XMLfname=$ARGV[0];
$slip_limit=$ARGV[1];

my $fault_model = XMLin($XMLfname, searchpath=>'.');

my $nfault = @{$fault_model->{fault}};

for ($i=0;$i<$nfault;$i++) {
  if ( $fault_model->{fault}->[$i]->{rectangular_fault}->{slip_rate}->{value} >= $slip_limit) {
    $f1="fault.".$i;
    open (F,">$f1") or die "cant open $f1\n";
    print F "$fault_model->{fault}->[$i]->{rectangular_fault}->{location}->{end1}->{c1}->{value} $fault_model->{fault}->[$i]->{rectangular_fault}->{location}->{end1}->{c2}->{value} $fault_model->{fault}->[$i]->{fname}\n";
    print F "$fault_model->{fault}->[$i]->{rectangular_fault}->{location}->{end2}->{c1}->{value} $fault_model->{fault}->[$i]->{rectangular_fault}->{location}->{end2}->{c2}->{value} $fault_model->{fault}->[$i]->{rectangular_fault}->{length}->{value} $i\n";
    close (F);
  }
}