代码如下

#!/usr/bin/env perl 
use strict;
use threads;
use Thread::Semaphore;

my $trait = $ARGV[0];
my $NumOfTrait = 0;
my $thread;
my $max_cpu = 28;
my $semaphore = new Thread::Semaphore($max_cpu);
my (@line, @line, %traitDict, $i);

open REF, "taxa_".$trait.".txt";
while(<REF>){
  chomp;
  @line = split /\t/;
  $traitDict{$line[0]} = $line[1];
  $NumOfTrait ++;
}

# check num of trait <> 28
for($i=1; $i<=$NumOfTrait; $i++){
  $semaphore->down();
  $thread = threads->new(\&run_gemma, $trait, $i, $traitDict{$i});
  $thread->detach();
}
#
&waitquit();

# functions 
sub waitquit{
  my ($i);
  for($i=0;$i<$max_cpu;$i++){
    $semaphore->down();
  }
  $semaphore->up($max_cpu);
}

sub run_gemma{
  my($phe, $col, $name) = @_;
  system("bash ../02.gemma.sh Bna171_".$phe." ".$col." ".$name);
  print "[INFO] $name is done.\n";
  $semaphore->up();
}