#!/usr/bin/env perl
# Extract xi.dat and zeta.dat from data.in
# Rok Zitko, rok.zitko@ijs.si, Jun 2010

use warnings;

my $fnin = "data.in";

open(F, "<$fnin") or die "Can't open $fnin: $!";
open(X,">xi.dat") or die "Can't open xi.dat: $!";
open(Z,">zeta.dat") or die "Can't open zeta.dat: $!";

# Find the discretization tables
while(<F>) {
    if (!/^z$/) {
      next; 
    } else {
      last;
    }
}

my $nr1 = <F>;

for (my $i = 0; $i <= $nr1; $i++) {
  my $x = <F>;
  print X $x;
}

my $nr2 = <F>;

for (my $i = 0; $i <= $nr2; $i++) {
  my $x = <F>;
  print Z $x;
}
