#!/usr/bin/env perl

use Env;
use warnings;
use strict;

my $Nz = 4;

chomp(my $min =   `getparam broaden_min param`);
chomp(my $max =   `getparam broaden_max param`);
chomp(my $ratio = `getparam broaden_ratio param`);
chomp(my $T =     `getparam T param`);
chomp(my $alpha = `getparam broaden_alpha param`);
chomp(my $gamma = `getparam broaden_gamma param`);

my $pr = "FDM_dens";

avg("spec_${pr}_A_d_u-A_d_u.bin",    "c-imG-u.dat");
avg("spec_${pr}_self_d_u-A_d_u.bin", "c-imF-u.dat");
avg("spec_${pr}_A_d_d-A_d_d.bin",    "c-imG-d.dat");
avg("spec_${pr}_self_d_d-A_d_d.bin", "c-imF-d.dat");

sub avg
{
    my $fn = shift;
    my $out = shift;
    system "$NRGPREFIX broaden -x $gamma -m $min -M $max -r $ratio   $fn $Nz $alpha $T 1e-99";
    system "mv spec.dat $out";
}