主题:[讨论]请问如何写个FORTRAN语言实现平均最近邻距离的计算?
我想请教fortran编程高手如何编写个程序求出最近邻的平均值?谢谢
数据如下(如这个例子1号离子的最近邻有25 21 20 31 35 28 10 12 14,只要能实现诸如2.79 2.80 这样的数字加和求平均即可)
ion position nearest neighbor table
1 0.332 0.445 0.489- 25 2.79 21 2.80 20 2.80 31 2.80 35 2.80 28 2.80 10 2.89 12 2.89 14 2.90
2 0.501 0.393 0.386- 20 2.79 29 2.79 16 2.80 38 2.80 30 2.80 28 2.80 10 2.89 14 2.90 9 2.90
3 0.610 0.411 0.558- 37 2.79 34 2.80 17 2.80 16 2.80 30 2.80 24 2.80 14 2.89 13 2.90 9 2.90
4 0.645 0.608 0.515- 17 2.80 36 2.80 22 2.80 34 2.80 32 2.80 19 2.80 9 2.90 13 2.80 11 2.90
5 0.366 0.641 0.446- 31 2.80 26 2.80 23 2.80 18 2.80 33 2.80 25 2.80 11 2.88 10 2.89 12 2.89
6 0.441 0.463 0.661- 35 2.79 37 2.80 21 2.80 27 2.80 15 2.80 24 2.80 12 2.89 14 2.89 13 2.90
7 0.535 0.590 0.344- 38 2.79 26 2.80 32 2.80 18 2.80 29 2.80 19 2.80 9 2.89 10 2.89 11 2.89
8 0.476 0.660 0.618- 23 2.79 22 2.80 15 2.80 36 2.80 33 2.80 27 2.80 12 2.90 11 2.90 13 2.90
9 0.568 0.502 0.454- 11 2.72 13 2.72 10 2.72 14 2.72 19 2.83 34 2.83 16 2.83 38 2.84 7 2.89 4 2.90 2 2.90 3 2.90
10 0.437 0.518 0.421- 9 2.72 11 2.72 14 2.72 12 2.73 26 2.83 28 2.83 29 2.83 25 2.84 1 2.89 2 2.89 5 2.89 7 2.89
11 0.504 0.619 0.482- 9 2.72 12 2.72 13 2.72 10 2.72 22 2.83 32 2.83 18 2.83 33 2.83 5 2.88 7 2.89 8 2.90 4 2.90
12 0.409 0.551 0.551- 13 2.72 11 2.72 14 2.72 10 2.73 27 2.83 21 2.83 23 2.83 31 2.83 1 2.89 5 2.89 6 2.89 8 2.90
13 0.540 0.535 0.583- 14 2.72 12 2.72 9 2.72 11 2.72 15 2.83 17 2.83 37 2.83 36 2.83 4 2.90 3 2.90 8 2.90 6 2.90
14 0.472 0.434 0.522- 13 2.72 10 2.72 9 2.72 12 2.72 24 2.83 30 2.83 20 2.83 35 2.83 3 2.89 6 2.89 2 2.90 1 2.90
15 0.512 0.569 0.717- 27 2.79 36 2.79 37 2.79 8 2.80 6 2.80 13 2.83
16 0.635 0.381 0.424- 34 2.79 38 2.80 30 2.80 2 2.80 3 2.80 9 2.83
17 0.676 0.519 0.618- 34 2.79 36 2.80 37 2.80 4 2.80 3 2.80 13 2.83
18 0.468 0.708 0.378- 33 2.79 32 2.80 26 2.80 5 2.80 7 2.80 11 2.83
19 0.668 0.571 0.382- 32 2.79 34 2.79 38 2.79 7 2.80 4 2.80 9 2.83
20 0.402 0.328 0.460- 35 2.79 28 2.79 2 2.79 30 2.80 1 2.80 14 2.83
21 0.308 0.481 0.622- 27 2.79 31 2.79 35 2.80 1 2.80 6 2.80 12 2.83
22 0.574 0.725 0.544- 36 2.79 32 2.79 33 2.80 8 2.80 4 2.80 11 2.83
23 0.341 0.672 0.581- 31 2.79 8 2.79 33 2.79 27 2.80 5 2.80 12 2.83
24 0.508 0.344 0.626- 30 2.79 35 2.80 37 2.80 3 2.80 6 2.80 14 2.83
25 0.300 0.534 0.386- 26 2.79 1 2.79 28 2.79 31 2.80 5 2.80 10 2.84
26 0.399 0.604 0.316- 25 2.79 29 2.79 7 2.80 18 2.80 5 2.80 10 2.83
27 0.378 0.585 0.684- 15 2.79 21 2.79 23 2.80 6 2.80 8 2.80 12 2.83
28 0.365 0.414 0.357- 29 2.79 20 2.79 25 2.79 2 2.80 1 2.80 10 2.83
29 0.464 0.484 0.287- 28 2.79 26 2.79 38 2.79 2 2.79 7 2.80 10 2.83
30 0.537 0.311 0.494- 24 2.79 16 2.80 20 2.80 2 2.80 3 2.80 14 2.83
31 0.272 0.568 0.518- 23 2.79 21 2.79 5 2.80 25 2.80 1 2.80 12 2.83
32 0.603 0.692 0.411- 19 2.79 22 2.79 18 2.80 7 2.80 4 2.80 11 2.83
33 0.440 0.742 0.510- 23 2.79 18 2.79 22 2.80 5 2.80 8 2.80 11 2.83
34 0.705 0.485 0.486- 16 2.79 19 2.79 17 2.79 3 2.80 4 2.80 9 2.83
35 0.373 0.361 0.593- 20 2.79 6 2.79 21 2.80 24 2.80 1 2.80 14 2.83
36 0.611 0.639 0.647- 22 2.79 15 2.79 17 2.80 4 2.80 8 2.80 13 2.83
37 0.577 0.449 0.688- 3 2.79 15 2.79 17 2.80 24 2.80 6 2.80 13 2.83
38 0.599 0.468 0.320- 7 2.79 29 2.79 19 2.79 16 2.80 2 2.80 9 2.84
数据如下(如这个例子1号离子的最近邻有25 21 20 31 35 28 10 12 14,只要能实现诸如2.79 2.80 这样的数字加和求平均即可)
ion position nearest neighbor table
1 0.332 0.445 0.489- 25 2.79 21 2.80 20 2.80 31 2.80 35 2.80 28 2.80 10 2.89 12 2.89 14 2.90
2 0.501 0.393 0.386- 20 2.79 29 2.79 16 2.80 38 2.80 30 2.80 28 2.80 10 2.89 14 2.90 9 2.90
3 0.610 0.411 0.558- 37 2.79 34 2.80 17 2.80 16 2.80 30 2.80 24 2.80 14 2.89 13 2.90 9 2.90
4 0.645 0.608 0.515- 17 2.80 36 2.80 22 2.80 34 2.80 32 2.80 19 2.80 9 2.90 13 2.80 11 2.90
5 0.366 0.641 0.446- 31 2.80 26 2.80 23 2.80 18 2.80 33 2.80 25 2.80 11 2.88 10 2.89 12 2.89
6 0.441 0.463 0.661- 35 2.79 37 2.80 21 2.80 27 2.80 15 2.80 24 2.80 12 2.89 14 2.89 13 2.90
7 0.535 0.590 0.344- 38 2.79 26 2.80 32 2.80 18 2.80 29 2.80 19 2.80 9 2.89 10 2.89 11 2.89
8 0.476 0.660 0.618- 23 2.79 22 2.80 15 2.80 36 2.80 33 2.80 27 2.80 12 2.90 11 2.90 13 2.90
9 0.568 0.502 0.454- 11 2.72 13 2.72 10 2.72 14 2.72 19 2.83 34 2.83 16 2.83 38 2.84 7 2.89 4 2.90 2 2.90 3 2.90
10 0.437 0.518 0.421- 9 2.72 11 2.72 14 2.72 12 2.73 26 2.83 28 2.83 29 2.83 25 2.84 1 2.89 2 2.89 5 2.89 7 2.89
11 0.504 0.619 0.482- 9 2.72 12 2.72 13 2.72 10 2.72 22 2.83 32 2.83 18 2.83 33 2.83 5 2.88 7 2.89 8 2.90 4 2.90
12 0.409 0.551 0.551- 13 2.72 11 2.72 14 2.72 10 2.73 27 2.83 21 2.83 23 2.83 31 2.83 1 2.89 5 2.89 6 2.89 8 2.90
13 0.540 0.535 0.583- 14 2.72 12 2.72 9 2.72 11 2.72 15 2.83 17 2.83 37 2.83 36 2.83 4 2.90 3 2.90 8 2.90 6 2.90
14 0.472 0.434 0.522- 13 2.72 10 2.72 9 2.72 12 2.72 24 2.83 30 2.83 20 2.83 35 2.83 3 2.89 6 2.89 2 2.90 1 2.90
15 0.512 0.569 0.717- 27 2.79 36 2.79 37 2.79 8 2.80 6 2.80 13 2.83
16 0.635 0.381 0.424- 34 2.79 38 2.80 30 2.80 2 2.80 3 2.80 9 2.83
17 0.676 0.519 0.618- 34 2.79 36 2.80 37 2.80 4 2.80 3 2.80 13 2.83
18 0.468 0.708 0.378- 33 2.79 32 2.80 26 2.80 5 2.80 7 2.80 11 2.83
19 0.668 0.571 0.382- 32 2.79 34 2.79 38 2.79 7 2.80 4 2.80 9 2.83
20 0.402 0.328 0.460- 35 2.79 28 2.79 2 2.79 30 2.80 1 2.80 14 2.83
21 0.308 0.481 0.622- 27 2.79 31 2.79 35 2.80 1 2.80 6 2.80 12 2.83
22 0.574 0.725 0.544- 36 2.79 32 2.79 33 2.80 8 2.80 4 2.80 11 2.83
23 0.341 0.672 0.581- 31 2.79 8 2.79 33 2.79 27 2.80 5 2.80 12 2.83
24 0.508 0.344 0.626- 30 2.79 35 2.80 37 2.80 3 2.80 6 2.80 14 2.83
25 0.300 0.534 0.386- 26 2.79 1 2.79 28 2.79 31 2.80 5 2.80 10 2.84
26 0.399 0.604 0.316- 25 2.79 29 2.79 7 2.80 18 2.80 5 2.80 10 2.83
27 0.378 0.585 0.684- 15 2.79 21 2.79 23 2.80 6 2.80 8 2.80 12 2.83
28 0.365 0.414 0.357- 29 2.79 20 2.79 25 2.79 2 2.80 1 2.80 10 2.83
29 0.464 0.484 0.287- 28 2.79 26 2.79 38 2.79 2 2.79 7 2.80 10 2.83
30 0.537 0.311 0.494- 24 2.79 16 2.80 20 2.80 2 2.80 3 2.80 14 2.83
31 0.272 0.568 0.518- 23 2.79 21 2.79 5 2.80 25 2.80 1 2.80 12 2.83
32 0.603 0.692 0.411- 19 2.79 22 2.79 18 2.80 7 2.80 4 2.80 11 2.83
33 0.440 0.742 0.510- 23 2.79 18 2.79 22 2.80 5 2.80 8 2.80 11 2.83
34 0.705 0.485 0.486- 16 2.79 19 2.79 17 2.79 3 2.80 4 2.80 9 2.83
35 0.373 0.361 0.593- 20 2.79 6 2.79 21 2.80 24 2.80 1 2.80 14 2.83
36 0.611 0.639 0.647- 22 2.79 15 2.79 17 2.80 4 2.80 8 2.80 13 2.83
37 0.577 0.449 0.688- 3 2.79 15 2.79 17 2.80 24 2.80 6 2.80 13 2.83
38 0.599 0.468 0.320- 7 2.79 29 2.79 19 2.79 16 2.80 2 2.80 9 2.84