论坛 › 【分子模拟编程】 › pl程序修改

625557018 发表于 2017-2-21 16:23:24

pl程序修改

程序如下：
#!/usr/bin/perl
##############################################################################
### SECTION:3: splits multi-MOLECULE input mol2 into individual mol2 files ###
##############################################################################

print("Splitting the multiple Ligand mol2 file into single mol2 files.\n");

#print("The single mol2 files will have the prefix: $prefix \n");
my($mol2_file) = "";
my($global_count) = 0;
my($local_count) = 0;

$mol2_file = $ARGV;
open (MOL2, $mol2_file)
    || die "\nError, cannot open ligand mol2 file: $mol2_file\n";

# assign basename for output files based on input MOL2 filename:

while (<MOL2>)
{
    chomp($_);
    if ($_ =Name)
    {
      $global_count++;
      open (OUT, ">$global_count.mol2")
            || die "\nError, cannot open output mol2 file: $global_count.mol2\n";
      print OUT "$name\n";
      $local_count++;
    }

    if ($_ =~ /^(@<TRIPOS>MOLECULE)/ && $local_count > 1)
    {
      $local_count = 1;
      next;
    }

    if ($_ !~ /^(@<TRIPOS>MOLECULE)/ && $local_count == 1 )
    {
      print OUT "$_\n" ;
      next;
    }
}

print("The number of single Ligand mol2 files generated is $global_count.\n");
close OUT;

问题：这个程序是拆分多分子mol2文件成单一分子的，但是输出文件名不是按原mol2中分子名称输出，而是按1,2,3...顺序输出，要怎样修改程序才能按原mol2中分子名称输出？？？

查看完整版本: pl程序修改