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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】如何写FORTRAN程序实现求平均最近邻距离
北京石油化工学院2026年研究生招生接收调剂公告
354/4返回列表上一页1234
查看: 3048  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +4 Mr. Z 2026-04-05 4/200 2026-04-06 00:51 by fmesaito
[考研] 304求调剂 +7 c297914 2026-04-05 8/400 2026-04-05 22:13 by hemengdong
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 304求调剂 +3 luoye0105 2026-04-05 3/150 2026-04-05 18:16 by 土木硕士招生
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 材料科学与工程调剂 +19 深V宿舍吧 2026-03-30 20/1000 2026-04-04 22:13 by hemengdong
[考研] 085400电子信息319求调剂（接受跨专业调剂） +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +11 努力奋斗112 2026-04-04 11/550 2026-04-04 20:51 by 蓝云思雨
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 288求调剂 一志愿哈工大 材料与化工 +39 洛神哥哥 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 抱歉 +5 田洪有 2026-03-30 5/250 2026-04-03 10:24 by linyelide
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
信息提示
关闭
请填处理意见
关闭