Hi,
I still need some help on the use of mdtoolchest. The dialsandwindows
stopped at a particular frame when processing one of my trajectories,
giving the following messages.
>fmt: read unexpected character
>apparent state: unit 2 named dna.dna
>last format: (A4,3F10.5,2I2,F8.4,1X,A4,1X,I3,1X,I4)
>Unit 2 is a sequential formatted external file
>*** Execution Terminated (115) ***
I also attached this particular frame and the pdb file. Any suggestion
will be appreciated.
regards,
--
Guanglei Cui
State University of New York
Department of Chemistry
Stony Brook, 11794-3400rdparm transformed trajectory, parm was ../vic.top
-15.265 4.369 -3.649 -15.943 5.048 -3.651 -17.182 4.394 -3.963 -17.176
4.249 -5.043 -18.020 5.039 -3.699 -17.357 2.995 -3.415 -18.439 2.904
-3.328 -16.828 2.015 -4.280 -16.414 0.914 -3.488 -17.258 0.280 -3.215
-15.428 0.040 -4.251 -14.798 0.587 -5.398 -15.047 1.564 -5.786 -13.847
-0.209 -6.055 -13.199 0.107 -6.860 -13.504 -1.441 -5.426 -12.579 -2.204
-5.961 -12.817 -3.173 -5.803 -12.246 -1.902 -6.865 -14.191 -1.973 -4.412
-15.124 -1.241 -3.831 -15.736 -1.909 -2.979 -16.827 2.806 -2.002 -16.325
3.701 -1.634 -15.808 1.670 -2.287 -14.862 2.167 -2.503 -15.637 1.013
-1.434 -17.890 2.268 -1.124 -17.802 2.353 0.536 -19.130 2.191 1.104
-16.989 3.555 0.830 -17.024 1.036 1.009 -17.612 -0.263 0.678 -18.510
-0.324 1.293 -17.864 -0.285 -0.382 -16.745 -1.436 1.040 -17.350 -2.335
1.159 -15.749 -1.698 0.017 -14.793 -2.496 0.590 -15.287 -3.447 0.790
-13.593 -2.650 -0.314 -13.020 -1.568 -0.828 -13.448 -0.581 -0.733 -11.865
-1.576 -1.533 -11.504 -0.671 -2.000 -11.341 -2.887 -1.742 -10.408 -3.043
-2.594 -9.994 -3.954 -2.728 -10.102 -2.263 -3.158 -11.845 -3.985 -1.190
-13.012 -3.883 -0.475 -13.456 -4.861 0.075 -15.899 -1.211 2.394 -15.670
-0.187 2.688 -14.589 -1.944 1.976 -13.801 -1.191 2.012 -14.285 -2.753
2.641 -16.543 -1.881 3.457 -16.095 -1.807 4.959 -17.183 -2.418 5.641
-15.672 -0.422 5.270 -14.806 -2.769 5.036 -15.083 -4.124 5.293 -15.532
-4.167 6.285 -15.803 -4.501 4.566 -13.743 -4.911 5.376 -14.040 -5.959
5.436 -13.071 -4.774 4.128 -11.729 -4.797 4.334 -11.296 -5.787 4.478
-11.008 -4.286 3.178 -11.162 -3.035 2.650 -11.831 -2.273 3.022 -10.365
-2.781 1.641 -9.726 -4.020 1.416 -8.831 -4.592 0.456 -8.355 -3.990
-0.589 -7.749 -4.566 -1.156 -8.576 -3.036 -0.835 -8.524 -5.930 0.407
-8.790 -6.607 1.540 -8.406 -7.611 1.643 -9.627 -6.224 2.499 -10.065
-4.930 2.371 -12.686 -4.559 6.511 -13.034 -3.760 7.165 -11.514 -4.020
5.670 -11.647 -2.961 5.447 -10.569 -4.280 6.147 -12.391 -5.733 7.294
-11.404 -5.699 8.549 -11.681 -6.900 9.324 -11.415 -4.413 9.201 -10.067
-5.883 7.802 -9.627 -7.146 7.338 -9.547 -7.866 8.152 -10.400 -7.505
6.658 -8.215 -7.111 6.640 -7.888 -8.138 6.475 -8.366 -6.294 5.466
-7.080 -5.688 5.329 -6.361 -6.472 5.091 -6.936 -4.612 4.248 -7.302
-3.286 4.235 -7.818 -2.780 5.038 -6.981 -2.664 3.155 -6.291 -3.592
2.402 -5.635 -3.536 1.148 -5.542 -2.496 0.353 -5.114 -2.609 -0.555
-6.143 -1.712 0.565 -4.931 -4.618 0.746 -4.948 -5.714 1.470 -4.231
-6.441 1.120 -5.591 -5.946 2.613 -6.290 -4.830 3.025 -7.198 -6.395
7.571 -7.645 -6.084 8.516 -6.752 -5.188 6.774 -7.451 -4.364 6.922
-5.688 -4.972 6.870 -6.205 -7.441 7.765 -4.858 -7.343 8.726 -4.431
-8.617 9.282 -4.955 -6.217 9.615 -3.778 -6.900 7.617 -3.360 -7.796
6.649 -2.754 -8.576 7.110 -4.216 -8.329 6.235 -2.578 -7.137 5.491
-2.353 -7.969 4.825 -3.172 -6.083 4.789 -2.219 -5.305 4.160 -2.135
-5.610 3.117 -2.665 -3.849 4.296 -3.309 -3.365 5.464 -3.628 -4.077
6.211 -3.299 -2.004 5.629 -3.633 -1.616 6.580 -2.805 -1.226 4.576
-2.648 0.045 4.691 -2.156 0.600 4.005 -2.952 0.531 5.522 -2.265
-1.658 3.463 -2.176 -3.014 3.274 -1.690 -3.468 2.262 -1.237 -6.649
6.026 -1.341 -6.155 6.992 -0.907 -5.649 4.931 -0.317 -4.804 5.288
-0.202 -6.199 4.307 -0.295 -7.646 6.157 1.252 -7.510 6.724 1.812
-8.879 6.996 1.093 -6.578 7.884 2.115 -6.782 5.546 2.485 -7.464
4.391 2.874 -8.452 4.638 1.651 -7.783 3.765 3.433 -6.650 3.511
3.778 -7.209 2.641 2.745 -5.449 3.074 3.662 -4.422 3.035 4.234
-4.510 2.111 2.919 -3.076 2.942 1.896 -2.514 3.649 1.468 -2.992
4.519 1.342 -1.474 3.078 2.118 -1.291 1.930 2.048 -0.259 0.930
1.305 0.661 0.932 3.019 -0.396 -0.144 3.103 0.303 -0.869 4.039
-1.371 -0.085 4.829 -1.581 -1.062 4.771 -1.049 -1.919 5.325 -2.458
-1.003 4.112 -2.315 0.855 3.173 -2.182 1.866 4.702 -6.151 4.137
4.895 -6.449 5.168 4.715 -4.628 4.126 4.326 -4.290 5.086 5.669
-4.184 3.840 5.819 -6.558 3.335 7.379 -6.342 3.809 7.978 -7.741
3.788 7.382 -5.712 5.144 8.081 -5.424 2.786 8.214 -5.718 1.379
8.885 -6.570 1.274 7.233 -6.011 1.004 8.726 -4.528 0.534 8.851
-4.773 -0.521 7.714 -3.450 0.515 8.410 -2.146 0.645 8.712 -1.898
-0.373 7.492 -1.158 1.294 7.019 -1.351 2.537 7.251 -2.267 3.060
6.260 -0.394 3.153 5.712 -0.628 4.493 6.479 -1.240 4.969 5.707
0.277 5.100 4.684 -0.990 4.478 5.843 0.843 2.433 5.064 1.701
2.809 6.309 0.898 1.121 6.091 1.672 0.509 7.189 0.011 0.535
7.727 0.241 -0.561 10.094 -3.978 1.040 10.425 -4.537 1.915 9.630
-2.533 1.488 9.303 -2.591 2.527 10.524 -1.937 1.307 11.101 -4.006
0.008 12.608 -3.677 0.288 13.486 -4.517 -0.578 12.942 -3.648 1.786
12.573 -2.178 -0.326 12.565 -1.895 -1.746 13.425 -2.418 -2.163 11.652
-2.229 -2.240 12.603 -0.404 -2.189 12.533 -0.444 -3.276 11.455 0.173
-1.601 11.774 1.386 -0.966 11.918 2.219 -1.653 10.782 1.884 0.003
10.787 1.477 1.333 11.499 0.755 1.704 9.961 2.062 2.268 9.979
1.508 3.702 10.381 2.218 4.425 8.951 1.235 3.939 10.604 0.616
3.669 9.062 3.169 1.889 8.153 3.647 2.545 9.343 3.630 0.639
8.794 4.454 0.438 10.137 3.060 -0.344 10.093 3.529 -1.492 13.833
0.338 -1.654 14.514 -0.409 -1.246 13.233 1.079 -0.468 13.148 0.420
0.395 13.748 2.021 -0.277 14.369 1.261 -2.636 15.735 2.051 -2.537
16.325 2.103 -3.920 16.518 1.511 -1.445 15.308 3.559 -2.183 14.419
4.315 -3.014 14.819 4.542 -4.002 13.506 3.753 -3.212 13.918 5.618
-2.417 13.438 6.194 -3.208 12.928 5.453 -1.427 12.947 6.603 -0.643
12.467 7.437 -1.155 12.213 6.504 0.607 12.203 5.494 1.518 12.764
4.575 1.430 11.266 5.617 2.421 10.763 6.879 2.167 9.777 7.650
2.867 9.037 7.319 3.795 9.469 8.885 2.315 8.852 9.474 2.855
10.322 9.375 1.396 10.102 10.558 0.871 9.360 11.146 1.221 10.790
10.837 0.186 11.188 8.733 0.667 11.361 7.451 1.113 14.996 6.492
-1.833 15.900 5.909 -1.659 14.382 6.850 -0.480 14.833 6.212 0.281
14.511 7.909 -0.257 15.204 7.540 -2.733 16.253 8.735 -2.481 16.522
9.168 -3.907 17.342 8.193 -1.649 15.535 10.001 -1.750 14.337 10.498
-2.272 14.602 10.920 -3.242 13.669 9.644 -2.378 13.607 11.588 -1.430
12.858 12.055 -2.069 12.982 11.110 -0.227 13.254 11.981 0.788 12.644
12.880 0.696 12.986 11.370 2.100 13.388 10.097 2.439 13.983 9.533
1.736 12.870 9.611 3.553 12.144 10.746 4.060 11.239 10.929 5.159
10.988 10.233 6.149 10.500 12.074 5.120 9.912 12.457 5.846 10.699
13.037 4.289 10.002 14.086 4.501 9.417 14.153 5.321 10.181 14.889
3.915 11.552 13.027 3.228 12.144 11.795 3.125 14.533 12.684 -1.007
15.458 12.415 -1.518 14.689 12.437 0.465 15.289 11.549 0.666 15.105
13.225 1.094 14.045 13.942 -1.278 14.341 14.293 -2.121 3.657 14.979
12.072 4.021 15.105 11.192 4.275 16.479 11.023 5.042 16.807 11.724
3.343 17.028 11.155 4.818 16.813 9.685 5.027 17.882 9.698 6.127
16.146 9.691 6.270 15.505 8.403 6.691 16.275 7.757 7.211 14.389
8.518 7.088 13.461 9.549 6.570 13.734 10.457 7.746 12.294 9.380
7.819 11.597 10.201 8.596 12.075 8.282 9.462 11.080 8.331 10.043
10.833 7.543 9.613 10.601 9.207 8.809 12.990 7.353 8.175 14.176
7.477 8.396 14.999 6.659 3.993 16.348 8.440 2.989 16.027 8.718
4.789 15.131 7.984 4.478 14.185 8.427 4.722 15.049 6.900 4.114
17.428 7.556 3.019 17.593 6.376 2.995 18.949 5.822 1.686 17.142
6.861 3.564 16.554 5.171 4.763 16.893 4.434 4.512 17.672 3.714
5.565 17.181 5.114 5.381 15.747 3.600 6.156 16.354 3.132 5.842
14.677 4.319 6.116 13.665 3.457 6.939 13.922 2.790 6.433 12.397
4.176 5.684 12.053 5.261 4.983 12.790 5.623 5.902 10.838 5.837
5.205 10.609 6.630 7.006 9.984 5.358 7.335 8.822 5.790 7.810
8.267 5.093 6.687 8.407 6.443 7.669 10.299 4.314 7.399 11.495
3.672 8.001 11.660 2.651 4.515 15.130 2.480 3.463 15.273 2.724
4.840 13.626 2.648 3.984 13.142 3.119 4.874 13.141 1.672 4.867
15.570 1.236 3.903 15.442 -0.069 4.125 16.602 -0.938 2.458 14.989
0.245 4.631 14.196 -0.843 6.030 14.276 -1.326 6.149 15.037 -2.097
6.687 14.456 -0.476 6.547 13.080 -2.233 7.317 13.544 -2.849 7.015
11.967 -1.436 6.400 10.848 -1.990 6.988 10.479 -2.831 6.088 9.906
-0.931 5.458 10.138 0.236 4.870 11.036 0.359 5.348 9.061 0.976
6.264 8.204 0.440 6.847 6.908 0.745 6.651 6.170 1.881 7.223
5.353 2.043 5.919 6.503 2.492 7.724 6.258 -0.051 8.055 6.911
-1.133 8.733 6.350 -1.758 7.611 8.067 -1.612 6.702 8.669 -0.786
5.497 12.581 -3.234 4.638 13.239 -3.368 5.025 11.257 -2.596 4.310
11.477 -1.804 4.656 10.519 -3.308 6.275 12.280 -4.409 5.500 11.979
-5.810 6.364 12.522 -6.897 4.057 12.355 -5.719 5.563 10.356 -5.787
6.837 9.674 -5.782 7.436 10.068 -6.603 7.467 9.951 -4.937 6.766
8.167 -5.821 7.761 7.861 -6.145 6.504 7.554 -4.615 5.631 6.478
-4.777 6.220 5.635 -5.138 5.021 6.057 -3.527 4.106 6.700 -2.808
3.608 7.602 -3.130 3.891 6.267 -1.574 4.715 5.106 -1.568 5.004
4.142 -0.516 4.515 4.139 0.665 4.582 3.329 1.265 3.918 4.915
0.912 5.856 3.187 -0.783 6.516 3.177 -1.969 7.128 2.312 -2.178
6.395 4.053 -2.954 5.469 5.061 -2.730 5.695 7.669 -6.892 5.247
8.530 -7.387 4.698 7.045 -5.861 4.027 7.780 -5.415 4.036 6.300
-6.302 6.368 6.732 -7.720 5.550 5.973 -8.849 6.582 5.781 -9.942
4.331 6.706 -9.122 5.266 4.482 -8.287 6.422 3.642 -7.976 6.834
3.307 -8.927 7.227 4.172 -7.468 6.223 2.406 -7.160 7.140 1.824
-7.064 5.807 2.717 -5.804 4.922 1.688 -5.382 5.386 0.743 -5.100
3.934 2.191 -4.410 3.273 3.374 -4.641 3.520 4.012 -5.477 2.365
3.792 -3.708 1.766 4.689 -3.761 2.078 2.984 -2.554 1.209 3.341
-1.563 1.125 2.680 -0.804 0.926 4.310 -1.588 2.654 1.745 -2.437
3.575 1.372 -3.359 4.078 0.289 -3.206 5.147 1.421 -7.789 4.562
1.830 -8.612 4.233 1.287 -6.629 3.351 1.886 -6.854 3.759 0.307
-6.695 5.620 0.160 -8.139 4.776 -1.001 -8.970 5.711 -1.794 -9.809
3.665 -0.325 -9.682 4.158 -1.876 -7.788 4.978 -2.583 -6.908 5.395
-3.398 -7.500 5.655 -1.855 -6.462 4.075 -3.253 -5.851 4.682 -3.799
-5.129 3.376 -2.299 -5.105 2.271 -2.958 -4.528 2.665 -3.615 -3.754
1.200 -2.083 -3.942 0.428 -1.130 -4.568 0.584 -0.845 -5.597 -0.549
-0.668 -3.770 -0.294 -1.290 -2.560 -0.904 -1.194 -1.308 -1.738 -0.395
-0.930 -0.582 -2.276 -0.463 -1.073 -2.348 0.417 0.405 -3.189 -0.691
0.553 -4.133 0.233 0.000 -3.916 1.050 1.278 -4.836 0.195 1.149
-3.204 -1.793 0.678 -2.259 -2.673 3.007 -4.239 -6.431 2.978 -4.109
-7.513 1.745 -3.835 -5.695 1.231 -3.128 -6.347 1.071 -4.637 -5.392
3.558 -5.559 -6.197 2.949 -6.949 -6.940 4.048 -7.972 -7.050 2.134
-6.605 -8.171 1.874 -7.337 -5.751 2.303 -8.072 -4.613 2.716 -9.045
-4.877 3.071 -7.405 -4.221 1.185 -8.106 -3.657 1.455 -8.611 -2.729
0.719 -6.794 -3.149 -0.595 -7.059 -2.616 -0.459 -7.269 -1.555 -1.569
-5.911 -2.717 -1.953 -5.277 -3.908 -1.474 -5.617 -4.814 -2.706 -4.166
-3.893 -3.096 -3.488 -5.112 -2.388 -3.796 -5.882 -4.142 -3.596 -5.397
-2.890 -2.418 -5.073 -3.279 -3.710 -2.610 -3.853 -2.655 -2.364 -3.060
-4.496 -1.499 -3.537 -4.394 -0.615 -2.289 -5.636 -1.524 -2.217 -6.303
-0.486 0.063 -8.894 -4.264 -0.116 -8.555 -5.284 -1.077 -8.377 -3.324
-1.978 -8.223 -3.917 -1.278 -9.037 -2.480 0.371 -10.257 -4.228 -0.593
-11.354 -4.820 0.079 -12.700 -4.841 -1.183 -10.713 -6.023 -1.937 -11.506
-3.790 -1.740 -12.236 -2.559 -1.580 -13.279 -2.833 -0.907 -11.860 -1.965
-3.016 -12.099 -1.663 -2.792 -12.519 -0.683 -3.421 -10.733 -1.550 -4.818
-10.809 -1.504 -5.148 -11.301 -0.590 -5.449 -9.460 -1.810 -5.080 -8.844
-2.978 -4.581 -9.502 -3.674 -5.454 -7.568 -3.266 -5.240 -6.915 -4.711
-4.323 -6.330 -4.783 -5.088 -7.838 -5.270 -6.075 -6.322 -5.085 -6.304
-6.782 -2.261 -6.582 -5.626 -2.350 -6.767 -7.527 -1.201 -7.474 -7.216
-0.551 -6.426 -8.865 -0.985 -6.880 -9.330 0.074 -4.217 -12.944 -2.197
-3.963 -13.567 -3.055 -5.225 -11.836 -2.501 -4.949 -11.507 -3.503 -6.286
-12.085 -2.492 -4.600 -13.782 -1.169 -5.537 -15.022 -1.437 -5.507 -15.813
-0.218 -5.103 -15.712 -2.659 -6.967 -14.462 -1.691 -7.844 -14.186 -0.632
-8.186 -15.177 -0.332 -7.342 -13.780 0.246 -9.126 -13.394 -0.910 -9.757
-13.485 -0.026 -8.845 -12.009 -1.212 -9.871 -11.639 -2.035 -10.772 -11.576
-1.425 -9.649 -10.282 -2.684 -9.063 -9.949 -3.918 -8.593 -10.729 -4.499
-9.018 -8.720 -4.183 -9.731 -8.152 -3.117 -10.120 -6.761 -2.902 -9.749
-5.760 -3.523 -10.911 -6.577 -1.833 -11.199 -5.655 -1.538 -11.276 -7.588
-0.977 -11.880 -7.265 0.103 -12.351 -6.377 0.205 -12.146 -7.941 0.804
-10.912 -8.839 -1.076 -10.167 -9.090 -2.210 -9.968 -14.023 -2.005 -9.416
-14.775 -2.569 -9.988 -12.821 -2.973 -9.081 -12.874 -3.576 -10.880 -12.918
-3.591 -11.202 -14.386 -1.375 -12.550 -14.731 -2.205 -13.386 -15.643 -1.400
-12.206 -15.181 -3.526 -13.148 -13.322 -2.198 -13.675 -12.608 -1.042 -14.302
-13.218 -0.391 -12.803 -12.230 -0.507 -14.692 -11.474 -1.384 -14.956 -11.060
-0.411 -13.986 -10.432 -2.137 -14.905 -9.848 -2.973 -15.665 -9.401 -2.332
-14.236 -8.925 -3.804 -13.363 -9.167 -4.803 -13.144 -10.160 -5.168 -12.971
-8.133 -5.504 -13.520 -7.034 -4.800 -13.412 -5.586 -5.009 -12.679 -4.917
-5.793 -14.098 -4.854 -4.045 -13.957 -3.865 -4.187 -14.871 -5.426 -3.070
-15.703 -4.635 -2.406 -15.707 -3.638 -2.572 -16.551 -5.072 -2.073 -15.053
-6.762 -2.949 -14.273 -7.541 -3.779 -15.945 -11.803 -2.266 -15.952 -12.857
-2.545 -15.653 -11.101 -3.576 -14.896 -11.561 -4.211 -16.561 -10.769 -4.080
-17.233 -11.403 -1.758 -17.885 -11.343 -2.460
#!/usr/bin/tcsh
# Example script to run a MD trajectories whose coordinate
# format is GROMOS. All irrelevent variables are commented out.
# Read the manuals for further explanation.
setenv status On # setenv Program progress - On
setenv dna_or_protein dna #Is it DNA or Protein?
# setenv nresidue #Number of residues- relevent for protein.
# setenv chain #Irrelevent for DNA. For Protein, the chain
# #name
# setenv start_residue #In case of proteins, the residue to start the
# #analysis from.
# setenv end_residue #In case of proteins, the residue to end the
# #analysis at.
setenv nucleotide 10 #Number of nucleotide units/strand.
setenv sequence "CCAACGTTGG"
#setenv sequence "CGAACG" #One letter nucleotide sequence for the left
#hand strand only
setenv insert_into_bdna no #Read the manual for insertion into
setenv insert_into_adna no #straight helical axes.
setenv frequency 1 #Choose every frequency'th structure for
# #analysis.
# setenv distance_tolerance #Used in determining connectivities.
# #If you get errors about missing
# #connectivities, increase this slightly,
# # setenv upto .4. If problems exist even after
# #that, check your coordinates.
setenv include_first no #Show the values for the first structure in
# #the graphics.
setenv draw_structures plate #Create strn.pdb files with just the plates.
setenv drawing_type dynamics #Use dynamics data structure.
setenv user_notes "Vickie Structure Results" #Some notes.
setenv coord_format pdb #The file is in PDB format.
setenv additional_file_format Dynamics #Additional file is in COR format.
setenv do_plotting yes # Don't do the graphics now.
# #I wish to edit the winlim.dat from the
# #dials_and_windows before producing graphics.
setenv layout 0 #Use the layout0.dat file.
#setenv device #set the device to Post, for example.
setenv start_time 0. #First structure corresponds to 20 ps
setenv time_interval 2. #Every subsequent structure is after 20ps.
# setenv canonical_dials #All these default to YES.
# setenv canonical_helics #Which means, All canonical "Dials" and
# setenv canonical_junctions #"Windows" will be drawn.
# setenv minmax_dials #All these default to NO.
# setenv minmax_helics
# setenv minmax_junctions #Don't show average or min and max values.
# setenv average_dials
# setenv average_helics
# setenv average_junctions #Upto here.
# setenv labels #I don't want any labels.
# setenv patterns #No patterns also.
# setenv canonical_patterns #No to this one too.
# setenv draw_dials #The next three default to yes, meaning
# setenv draw_helics #draw all the dials and windows.
# setenv draw_junctions #Set them to "No" if you want.
setenv for007 ../../r/massaged.pdb #I wish to analyze the Dickerson crystal
setenv for008 ../md5.21st #Dynamics results.
dials_and_windows #I have the defined PATH to include the
Attachment Converted: "c:\eudora\attach\massaged.pdb"
Received on Wed Jun 07 2000 - 08:54:34 PDT