应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】如何写FORTRAN程序实现求平均最近邻距离
354/4返回列表上一页1234
查看: 3789  |  回复: 34
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

阿黛拉

银虫 (小有名气)
引用回帖:
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 ...

您好，分成两组的怎么写？比如1-10一组，11-38一组。三个原子体系的不常见，不过也有用，我保存下来了。
赞一下 回复此楼
进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。
31楼2011-02-26 10:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by 阿黛拉 at 2011-02-26 10:27:25:
您好，分成两组的怎么写？比如1-10一组，11-38一组。三个原子体系的不常见，不过也有用，我保存下来了。

改以下四行中数值：
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

按顺序分别改成 1、10、11、38。如果只是两原子体系的话，其实不用这么麻烦，把我前面给的那个程序稍改一下就可以了……，但这种情况可能更加通用一些，呵呵……
赞一下(1人) 回复此楼
32楼2011-02-26 10:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿黛拉

银虫 (小有名气)
引用回帖:
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

按顺 ...

你好，两种程序我都试了一下，出来的都有点问题。 请看附件。比如1  4 之间间距不对，第一组中1 4间距接着应该输出1  3间距，不知道什么原因？
赞一下 回复此楼
进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。
33楼2011-02-26 12:56:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by 阿黛拉 at 2011-02-26 12:56:09:
你好，两种程序我都试了一下，出来的都有点问题。 请看附件。比如1  4 之间间距不对，第一组中1 4间距接着应该输出1  3间距，不知道什么原因？

你咋改的程序？我这里很正常啊……
CODE:
integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 4
integer, parameter :: ion_pos_3 = 5
integer, parameter :: ion_pos_4 = 6

输入的数据文件
CODE:
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

输出的相关文件：

fort.1
CODE:
  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

for.2
CODE:
  5    6  2.94
  6    5  2.94

fort.3
CODE:
  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

终端的输出
CODE:
  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   

赞一下(1人) 回复此楼
34楼2011-02-26 13:05:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿黛拉

银虫 (小有名气)
引用回帖:
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

按顺 ...

integer, parameter :: ion_pos_1 = 1
integer, parameter :: ion_pos_2 = 4
integer, parameter :: ion_pos_3 = 5
integer, parameter :: ion_pos_4 = 6
我是按这样改的，我又运行了一下，没有问题了。呵呵 谢谢你了。刚才用你上一个程序运行时出的错，可能是源文件的问题，我又换了一下OUTCAR没有出现间距错误，只是输出不是我想要的。 哈哈 我先用你最后编的吧 。谢谢你不厌其烦的讲解和帮忙。
赞一下 回复此楼
进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。
35楼2011-02-26 13:37:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿黛拉 的主题更新
354/4返回列表上一页1234
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 17:01 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 16:58 by BCdvBAEpt3zm
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 16:45 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 16:32 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:08 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 15:01 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 14:55 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭