#!/usr/local/bin/perl -w use strict; use Getopt::Long; my $pdbfile; my $odatfile; GetOptions( "pdb=s" => \$pdbfile, "odat=s" => \$odatfile, ) or die; my $pdb=`grep "^ATOM" $pdbfile`; my $odat=`grep -v " -" $odatfile |grep -v "INDO"`; my @pdblines=split "\n",$pdb; my @odatlines=split "\n",$odat; my $prev_res=-1; my $res_cnt=-1; my $entropy=0.0; foreach my $line (@pdblines){ my $res=substr($line,22,4); $res=$res+0; if($res>$prev_res){ $res_cnt++; $entropy=$odatlines[$res_cnt]; $entropy=substr($entropy,0,7); $prev_res=$res; } print substr($line,0,60)." $entropy\n"; } 1;