open(IN, "/path/outModified.pl") or die "$!"; while (my $line = ) { chomp($line); my @array = (split (/\s+/, $line))[6, 7, 8]; # print "@array\n"; push @points, [ @array ]; } close(IN); print '@points : ', Dumper \@points; open my $out_file, '>', "/path/dist_modified.pl" or die "cannot open: $!\n"; for my $i1 ( 0 .. $#points -1 ){ my ( $x1, $y1, $z1 ) = @{ $points[$i1] }; for my $i2 ( 1 .. $#points){ my ( $x2, $y2, $z2 ) = @{ $points[$i2] }; my $dist = sqrt(($x2 - $x1)**2 + ($y2 - $y1)**2 + ($z2 - $z1)**2); print $out_file "153L.pdb: the distance between CA $i1 and CA $i2 = $dist\n"; } } close $out_file or die "cannot close file: $!\n"; #### ATOM 2 CA ARG A 1 6.324 32.707 50.379 ATOM 13 CA THR A 2 5.197 32.618 46.826 ATOM 20 CA ASP A 3 4.020 36.132 46.259 ATOM 28 CA CYS A 4 7.131 38.210 45.919 ATOM 34 CA TYR A 5 6.719 38.935 42.270 ATOM 46 CA GLY A 6 2.986 39.221 41.892 ATOM 50 CA ASN A 7 -0.269 37.184 41.565 ATOM 58 CA VAL A 8 -1.140 35.549 38.341 ATOM 65 CA ASN A 9 -4.817 35.710 39.211 ATOM 73 CA ARG A 10 -4.704 39.489 39.013 ATOM 84 CA ILE A 11 -3.243 39.663 35.555 #### Index Distance 1 2.69 2 2.97032 3 3.25064 4 3.53096 5 3.81128 6 4.0916 7 4.37192 8 4.65224 9 4.93256 10 5.21288 11 5.4932 12 5.77352 13 6.05384 14 6.33416 15 6.61448 16 6.8948 17 7.17512 18 7.45544 19 7.73576 20 8.01608 21 8.2964 22 8.57672 23 8.85704 24 9.13736 25 9.41768 26 9.698 27 9.97832 28 10.25864 29 10.53896 30 10.81928 31 11.0996 32 11.37992 33 11.66024 34 11.94056 35 12.22088 36 12.5012 37 12.78152 38 13.06184 39 13.34216 40 13.62248 41 13.9028 42 14.18312 43 14.46344 44 14.74376 45 15.02408 46 15.3044 47 15.58472 48 15.86504 49 16.14536 50 16.42568 51 16.706 52 16.98632 53 17.26664 54 17.54696 55 17.82728 56 18.1076 57 18.38792 58 18.66824 59 18.94856 60 19.22888 61 19.5092 62 19.78952 63 20.06984 64 20.35016 65 20.63048 66 20.9108 67 21.19112 68 21.47144 69 21.75176 70 22.03208 71 22.3124 72 22.59272 73 22.87304 74 23.15336 75 23.43368 76 23.714 77 23.99432 78 24.27464 79 24.55496 80 24.83528 81 25.1156 82 25.39592 83 25.67624 84 25.95656 85 26.23688 86 26.5172 87 26.79752 88 27.07784 89 27.35816 90 27.63848 91 27.9188 92 28.19912 93 28.47944 94 28.75976 95 29.04008 96 29.3204 97 29.60072 98 29.88104 99 30.16136 100 30.44168 101 30.722 102 31.00232 103 31.28264 104 31.56296 105 31.84328 106 32.1236 107 32.40392 108 32.68424 109 32.96456 110 33.24488 111 33.5252 112 33.80552 113 34.08584 114 34.36616 115 34.64648 116 34.9268 117 35.20712 118 35.48744 119 35.76776 120 36.04808 121 36.3284 122 36.60872 123 36.88904 124 37.16936 125 37.44968 126 37.73 127 38.01032 128 38.29064 129 38.57096 130 38.85128 131 39.1316 132 39.41192 133 39.69224 134 39.97256 135 40.25288 136 40.5332 137 40.81352 138 41.09384 139 41.37416 140 41.65448 141 41.9348 142 42.21512 143 42.49544 144 42.77576 145 43.05608 146 43.3364 147 43.61672 148 43.89704 149 44.17736 150 44.45768 151 44.738 152 45.01832 153 45.29864 154 45.57896 155 45.85928 156 46.1396 157 46.41992 158 46.70024 159 46.98056 160 47.26088 161 47.5412 162 47.82152 163 48.10184 164 48.38216 165 48.66248 166 48.9428 167 49.22312 168 49.50344 169 49.78376 170 50.06408 171 50.3444 172 50.62472 173 50.90504 174 51.18536 175 51.46568 176 51.746 177 52.02632 178 52.30664 179 52.58696 180 52.86728 181 53.1476 182 53.42792 183 53.70824