test_dewi_tauris.pl

#!/usr/bin/env perl
use strict;
use rob_misc;
use binary_grid2;
use binary_grid::C;

# 
# script to test the new polytrope comenv lambda 
# code vs Dewi & Tauris 2000
#

my $outdir = '/tmp/dewi';
mkdirhier($outdir);
my $datafp;

my $population = binary_grid2->new(defaults());



$population->set(
    parse_bse_function_pointer=>\&parse_bse,
    );

sub parse_bse
{

    while(1)
    {

        my $la = $population->tbse_line();
        my $header = shift @$la;
        last if($header eq 'fin');
	if($header eq 'COMENVLOG')
	{
            #if($la->[1] == 5 || $la->[1] == 6)
            #if($la->[1] != 4)
            {
                print {$datafp} "@$la\n";
            }
        }
    }
}

# reconstruct their Fig.1 which is lamdba as a 
# function of R/Rsun for various masses
my $gpfile = {
    dewi=>"$outdir/Dewi_Tauris_Fig1.plt",
    vds=>"$outdir/VDSluys_Fig4.plt"
};
my $gp = {};

foreach my $k (keys %$gpfile)
{    
    open($gp->{$k},'>',$gpfile->{$k})||die;

    my $pdf = $gpfile->{$k};
    $pdf =~ s/plt$/pdf/;

    print {$gp->{$k}} "set terminal pdfcairo
set output \"$pdf\"
set ylabel \"Lambda\"
set xtics
set ytics nomirror
set pointsize 0.2
";
}

print {$gp->{dewi}} "
set xrange[0:500]
set yrange[0.1:100]