#!/usr/bin/perl

@files = ('10_1_2','10_1_3','10_1_4','10_1_5','10_1_6','10_1_7','10_1_8','10_2_2','10_2_3','10_2_4','10_2_5','10_2_6','10_2_7','10_2_8','10_3_1','10_3_2','10_3_3','10_3_4','10_3_5','10_3_6','10_3_7','10_3_8','10_4_1','10_4_2','10_4_3','10_4_4','10_4_5','10_4_6','10_4_7','10_4_8','10_5_1','10_5_2','10_5_3','10_5_4','10_5_5','10_5_6','10_5_7','10_5_8','10_6_1','10_6_2','10_6_3','10_6_4','10_6_5','10_6_6','10_6_7','10_6_8','10_7_1','10_7_2','10_7_3','10_7_4','10_7_5','10_7_6','10_7_7','10_7_8','11_1_2','11_1_3','11_1_4','11_1_5','11_1_6','11_1_7','11_1_8','11_2_2','11_2_3','11_2_4','11_2_5','11_2_6','11_2_7','11_2_8','11_3_1','11_3_2','11_3_3','11_3_4','11_3_5','11_3_6','11_3_7','11_3_8','11_4_1','11_4_2','11_4_3','11_4_4','11_4_5','11_4_6','11_4_7','11_4_8','11_5_1','11_5_2','11_5_3','11_5_4','11_5_5','11_5_6','11_5_7','11_5_8','11_6_1','11_6_2','11_6_3','11_6_4','11_6_5','11_6_6','11_6_7','11_6_8','11_7_1','11_7_2','11_7_3','11_7_4','11_7_5','11_7_6','11_7_7','11_7_8','12_1_2','12_1_3','12_1_4','12_1_5','12_1_6','12_1_7','12_1_8','12_2_2','12_2_3','12_2_4','12_2_5','12_2_6','12_2_7','12_2_8','12_3_1','12_3_2','12_3_3','12_3_4','12_3_5','12_3_6','12_3_7','12_3_8','12_4_1','12_4_2','12_4_3','12_4_4','12_4_5','12_4_6','12_4_7','12_4_8','12_5_1','12_5_2','12_5_3','12_5_4','12_5_5','12_5_6','12_5_7','12_5_8','12_6_1','12_6_2','12_6_3','12_6_4','12_6_5','12_6_6','12_6_7','12_6_8','12_7_1','12_7_2','12_7_3','12_7_4','12_7_5','12_7_6','12_7_7','12_7_8','13_1_2','13_1_3','13_1_4','13_1_5','13_1_6','13_1_7','13_1_8','13_2_2','13_2_3','13_2_4','13_2_5','13_2_6','13_2_7','13_2_8','13_3_1','13_3_2','13_3_3','13_3_4','13_3_5','13_3_6','13_3_7','13_3_8','13_4_1','13_4_2','13_4_3','13_4_4','13_4_5','13_4_6','13_4_7','13_4_8','13_5_1','13_5_2','13_5_3','13_5_4','13_5_5','13_5_6','13_5_7','13_5_8','13_6_1','13_6_2','13_6_3','13_6_4','13_6_5','13_6_6','13_6_7','13_6_8','13_7_1','13_7_2','13_7_3','13_7_4','13_7_5','13_7_6','13_7_7','13_7_8','14_1_2','14_1_3','14_1_4','14_1_5','14_1_6','14_1_7','14_1_8','14_2_2','14_2_3','14_2_4','14_2_5','14_2_6','14_2_7','14_2_8','14_3_1','14_3_2','14_3_3','14_3_4','14_3_5','14_3_6','14_3_7','14_3_8','14_4_1','14_4_2','14_4_3','14_4_4','14_4_5','14_4_6','14_4_7','14_4_8','14_5_1','14_5_2','14_5_3','14_5_4','14_5_5','14_5_6','14_5_7','14_5_8','14_6_1','14_6_2','14_6_3','14_6_4','14_6_5','14_6_6','14_6_7','14_6_8','14_7_1','14_7_2','14_7_3','14_7_4','14_7_5','14_7_6','14_7_7','14_7_8','15_1_2','15_1_3','15_1_4','15_1_5','15_1_6','15_1_7','15_1_8','15_2_2','15_2_3','15_2_4','15_2_5','15_2_6','15_2_7','15_2_8','15_3_1','15_3_2','15_3_3','15_3_4','15_3_5','15_3_6','15_3_7','15_3_8','15_4_1','15_4_2','15_4_3','15_4_4','15_4_5','15_4_6','15_4_7','15_4_8','15_5_1','15_5_2','15_5_3','15_5_4','15_5_5','15_5_6','15_5_7','15_5_8','15_6_1','15_6_2','15_6_3','15_6_4','15_6_5','15_6_6','15_6_7','15_6_8','15_7_1','15_7_2','15_7_3','15_7_4','15_7_5','15_7_6','15_7_7','15_7_8','16_1_2','16_1_3','16_1_4','16_1_5','16_1_6','16_1_7','16_1_8','16_2_2','16_2_3','16_2_4','16_2_5','16_2_6','16_2_7','16_2_8','16_3_1','16_3_2','16_3_3','16_3_4','16_3_5','16_3_6','16_3_7','16_3_8','16_4_1','16_4_2','16_4_3','16_4_4','16_4_5','16_4_6','16_4_7','16_4_8','16_5_1','16_5_2','16_5_3','16_5_4','16_5_5','16_5_6','16_5_7','16_5_8','16_6_1','16_6_2','16_6_3','16_6_4','16_6_5','16_6_6','16_6_7','16_6_8','16_7_1','16_7_2','16_7_3','16_7_4','16_7_5','16_7_6','16_7_7','16_7_8','17_1_2','17_1_3','17_1_4','17_1_5','17_1_6','17_1_7','17_1_8','17_2_2','17_2_3','17_2_4','17_2_5','17_2_6','17_2_7','17_2_8','17_3_1','17_3_2','17_3_3','17_3_4','17_3_5','17_3_6','17_3_7','17_3_8','17_4_1','17_4_2','17_4_3','17_4_4','17_4_5','17_4_6','17_4_7','17_4_8','17_5_1','17_5_2','17_5_3','17_5_4','17_5_5','17_5_6','17_5_7','17_5_8','17_6_1','17_6_2','17_6_3','17_6_4','17_6_5','17_6_6','17_6_7','17_6_8','17_7_1','17_7_2','17_7_3','17_7_4','17_7_5','17_7_6','17_7_7','17_7_8','18_1_2','18_1_3','18_1_4','18_1_5','18_1_6','18_1_7','18_1_8','18_2_2','18_2_3','18_2_4','18_2_5','18_2_6','18_2_7','18_2_8','18_3_1','18_3_2','18_3_3','18_3_4','18_3_5','18_3_6','18_3_7','18_3_8','18_4_1','18_4_2','18_4_3','18_4_4','18_4_5','18_4_6','18_4_7','18_4_8','18_5_1','18_5_2','18_5_3','18_5_4','18_5_5','18_5_6','18_5_7','18_5_8','18_6_1','18_6_2','18_6_3','18_6_4','18_6_5','8_1_2','8_1_3','8_1_4','8_1_5','8_1_6','8_1_7','8_1_8','8_2_2','8_2_3','8_2_4','8_2_5','8_2_6','8_2_7','8_2_8','8_3_1','8_3_2','8_3_3','8_3_4','8_3_5','8_3_6','8_3_7','8_3_8','8_4_1','8_4_2','8_4_3','8_4_4','8_4_5','8_4_6','8_4_7','8_4_8','8_5_1','8_5_2','8_5_3','8_5_4','8_5_5','8_5_6','8_5_7','8_5_8','8_6_1','8_6_2','8_6_3','8_6_4','8_6_5','8_6_6','8_6_7','8_6_8','8_7_1','8_7_2','8_7_3','8_7_4','8_7_5','8_7_6','8_7_7','8_7_8','9_1_2','9_1_3','9_1_4','9_1_5','9_1_6','9_1_7','9_1_8','9_2_2','9_2_3','9_2_4','9_2_5','9_2_6','9_2_7','9_2_8','9_3_1','9_3_2','9_3_3','9_3_4','9_3_5','9_3_6','9_3_7','9_3_8','9_4_1','9_4_2','9_4_3','9_4_4','9_4_5','9_4_6','9_4_7','9_4_8','9_5_1','9_5_2','9_5_3','9_5_4','9_5_5','9_5_6','9_5_7','9_5_8','9_6_1','9_6_2','9_6_3','9_6_4','9_6_5','9_6_6','9_6_7','9_6_8','9_7_1','9_7_2','9_7_3','9_7_4','9_7_5','9_7_6','9_7_7','9_7_8','18_6_6','18_6_7','18_6_8','18_7_1','18_7_2','18_7_3','18_7_4','18_7_5','18_7_6','18_7_7','18_7_8','19_1_2','19_1_3','19_1_4','19_1_5','19_1_6','19_1_7','19_1_8','19_2_2','19_2_3','19_2_4','19_2_5','19_2_6','19_2_7','19_2_8','19_3_1','19_3_2','19_3_3','19_3_4','19_3_5','19_3_6','19_3_7','19_3_8','19_4_1','19_4_2','19_4_3','19_4_4','19_4_5','19_4_6','19_4_7','19_4_8','19_5_1','19_5_2','19_5_3','19_5_4','19_5_5','19_5_6','19_5_7','19_5_8','19_6_1','19_6_2','19_6_3','19_6_4','19_6_5','19_6_6','19_6_7','19_6_8','19_7_1','19_7_2','19_7_3','19_7_4','19_7_5','19_7_6','19_7_7','19_7_8','20_1_2','20_1_3','20_1_4','20_1_5','20_1_6','20_1_7','20_1_8','20_2_2','20_2_3','20_2_4','20_2_5','20_2_6','20_2_7','20_2_8','20_3_1','20_3_2','20_3_3','20_3_4','20_3_5','20_3_6','20_3_7','20_3_8','20_4_1','20_4_2','20_4_3','20_4_4','20_4_5','20_4_6','20_4_7','20_4_8','20_5_1','20_5_2','20_5_3','20_5_4','20_5_5','20_5_6','20_5_7','20_5_8','20_6_1','20_6_2','20_6_3','20_6_4','20_6_5','20_6_6','20_6_7','20_6_8','20_7_1','20_7_2','20_7_3','20_7_4','20_7_5','20_7_6','20_7_7','20_7_8','21_1_2','21_1_3','21_1_4','21_1_5','21_1_6','21_1_7','21_1_8','21_2_2','21_2_3','21_2_4','21_2_5','21_2_6','21_2_7','21_2_8','21_3_1','21_3_2','21_3_3','21_3_4','21_3_5','21_3_6','21_3_7','21_3_8','21_4_1','21_4_2','21_4_3','21_4_4','21_4_5','21_4_6','21_4_7','21_4_8','21_5_1','21_5_2','21_5_3','21_5_4','21_5_5','21_5_6','21_5_7','21_5_8','21_6_1','21_6_2','21_6_3','21_6_4','21_6_5','21_6_6','21_6_7','21_6_8','21_7_1','21_7_2','21_7_3','21_7_4','21_7_5','21_7_6','21_7_7');

print "#Filename Average I, Average U, Average P, dlogP/dU\n";
print "#f\tAverI\tAverU\tAverP\tdlnP/dU\n";


foreach $f(@files) {

	$filename1 = $f . "_ivp_der.dat";

	open (WAVE1, "<", "$filename1");
	
	$idx = 0;
	while ($l1=<WAVE1>) {

		$l1 =~ s/\n//;

		if ($l1 =~ s/#/#/) { } else {
			@lines[$idx] = $l1;
			$idx++;
		}
	}

	close WAVE1;

	$maxline = '';
	$lastMax = 0;

	for ($i=0; $i<$idx-1; $i++) {
	
		$thisline = @lines[$i];
		
#		print $thisline;
	
		($f,$aI,$aU,$aP,$di,$du,$dp,$dlnpdu)=split ("\t", $thisline);
	
		if ($dlnpdu > $lastMax) {
			$lastMax = $dlnpdu;
			$maxline = "$f\t$aI\t$aU\t$aP\t$dlnpdu\n";
		}
	}
	print $maxline;
}