print "\nEnter the input file: "; $inputFile = ; chomp $inputFile; unless (open(INPUTFILE, $inputFile)) { print "Cannot read from '$inputFile'.\nProgram closing.\n"; ; exit;} chomp(@dataArray = ); close(INPUTFILE); for ($line = 0; $line <= scalar @dataArray; $line++) { if ($dataArray[$line] =~ m/ATOM\s+(\d+)\s+(\w+)\s+(\w{3})\s+(\w)+\s+(\d+)\s+(\S+\.\S+)\s+(\S+\.\S+)\s+(\S+\.\S+)\s+(.+\S)(.\d\d+\.\d\d.+)/ig) { $m1=$1; $m2=$2; $m3=$3; $m5=$5; $m6=$6; $m7=$7; $m8=$8; $m9=$9; $m10=$10; push(@m3,$m3); push(@m5,$m5); foreach $line (@m3,@m5) {if ($m3[$line] eq $m3[$line+1]) {$m5[$line]=$m5[$line+1];} elsif ($m3[$line] ne $m3[$line+1]){$m5[$line+1]=$m5[$line]+1;}} $~="PDBFORMAT"; format PDBFORMAT = ATOM @|||| @||| @|| @||| @|||||| @|||||| @|||||| @>>>>> @>>>>> $m1, $m2, $m3,$m5, $m6, $m7, $m8, $m9, $m10 . open(PDBFORMAT,">>my2pdb.txt") or die "Can't open anything"; write PDBFORMAT;}} close PDBFORMAT;