From cf0d43482f98d5a83e544a6e642f8e5142677029 Mon Sep 17 00:00:00 2001
From: elserj <elserj@localhost>
Date: Fri, 4 Sep 2015 18:49:01 +0000
Subject: [PATCH] Added script that takes a species list and copies only
 relevant sqltable files (from inparanoid) to a new directory

svn path=/; revision=632
---
 interactome_scripts/copy_sqltables_subset.pl | 36 ++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100755 interactome_scripts/copy_sqltables_subset.pl

diff --git a/interactome_scripts/copy_sqltables_subset.pl b/interactome_scripts/copy_sqltables_subset.pl
new file mode 100755
index 0000000..4f28600
--- /dev/null
+++ b/interactome_scripts/copy_sqltables_subset.pl
@@ -0,0 +1,36 @@
+#!/usr/bin/perl
+
+use strict;
+use warnings;
+
+# check for arguments and explain usage
+if ($#ARGV !=1) {
+	print "usage: copy_sqltables_subset.pl species_list_file output_dir\n";
+	exit;
+}
+
+my $speciesfile = $ARGV[0];
+my $outdir = $ARGV[1];
+
+my @speciesarray;
+open (SPECIESFILE, $speciesfile);
+while(<SPECIESFILE>) {
+		my $line = $_;
+		chomp $line;
+		
+		push (@speciesarray, $line);
+}
+close(SPECIESFILE);
+
+# Get list of files in input dir
+my @file_list = glob("sqltable*");
+
+for (my $i = 0; $i < scalar @speciesarray; $i++) {
+		for (my $j = $i + 1; $j < scalar @speciesarray; $j++) {
+				foreach my $file (@file_list) {
+						if($file =~ m/$speciesarray[$i]/ && $file =~ m/$speciesarray[$j]/ ) {
+								system("cp $file $outdir/$file");
+						}
+				}
+		}
+}
-- 
2.34.1