Originally posted by snoopyzhao at 2011-02-26 10:13:17:
根据上面的意思，重新整理了一下：

[code]
program ex

implicit none
character(len = 128) :: line, fm
integer :: ios, len_line, i, j, ion_num, ion_tmp
real, dimension(5000) :: table
inte ...

Originally posted by 阿黛拉 at 2011-02-26 10:27:25:
您好，分成两组的怎么写？比如1-10一组，11-38一组。三个原子体系的不常见，不过也有用，我保存下来了。

integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 10
integer, parameter :: ion_pos_3 = 28
integer, parameter :: ion_pos_4 = 38

Originally posted by snoopyzhao at 2011-02-26 10:30:51:
改以下四行中数值：

CODE:

integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 10
integer, parameter :: ion_pos_3 = 28
integer, parameter :: ion_pos_4 = 38

按顺 ...

Originally posted by 阿黛拉 at 2011-02-26 12:56:09:
你好，两种程序我都试了一下，出来的都有点问题。 请看附件。比如1  4 之间间距不对，第一组中1 4间距接着应该输出1  3间距，不知道什么原因？

integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 4
integer, parameter :: ion_pos_3 = 5
integer, parameter :: ion_pos_4 = 6

ion  position               nearest neighbor table
   1  0.388  0.481  0.529-   6 2.91   5 2.94   4 2.92   3 2.94
   2  0.596  0.481  0.529-   6 2.94   5 2.92   3 2.94   4 2.90
   3  0.492  0.377  0.529-   5 2.90   6 2.91   1 2.94   2 2.94
   4  0.492  0.585  0.529-   1 2.94   2 2.90   5 2.93   6 2.93
   5  0.492  0.481  0.425-   2 2.90   1 2.94   3 2.94   6 2.94
   6  0.492  0.481  0.633-   2 2.94   1 2.93   3 2.92   5 2.94

  1    4  2.92
  1    3  2.94
  2    3  2.94
  2    4  2.90
  3    1  2.94
  3    2  2.94
  4    1  2.94
  4    2  2.90

  5    6  2.94
  6    5  2.94

  1    6  2.91
  1    5  2.94
  2    6  2.94
  2    5  2.92
  3    5  2.90
  3    6  2.91
  4    5  2.93
  4    6  2.93
  5    2  2.90
  5    1  2.94
  5    3  2.94
  6    2  2.94
  6    1  2.93
  6    3  2.92

  1    6  2.91
  1    5  2.94
  1    4  2.92
  1    3  2.94
  2    6  2.94
  2    5  2.92
  2    3  2.94
  2    4  2.90
  3    5  2.90
  3    6  2.91
  3    1  2.94
  3    2  2.94
  4    1  2.94
  4    2  2.90
  4    5  2.93
  4    6  2.93
  5    2  2.90
  5    1  2.94
  5    3  2.94
  5    6  2.94
  6    2  2.94
  6    1  2.93
  6    3  2.92
  6    5  2.94
   2.9270835   
   23.420002               8   5.8800001               2   40.949997              14
   2.9275002       2.9400001       2.9249997   

Originally posted by snoopyzhao at 2011-02-26 10:30:51:
改以下四行中数值：

CODE:

integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 10
integer, parameter :: ion_pos_3 = 28
integer, parameter :: ion_pos_4 = 38

按顺 ...