| 查看: 1429 | 回复: 5 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
alvin510银虫 (小有名气)
|
[求助]
求txt中数据重新排序的fortran代码
|
|||
| 我用MS计算了一个晶体的能带结构，其结果保存在了unitbands.txt文件中。不幸的是，其中的结果没有按自然顺序排列，即原文件是按k-point =1,k-point=7,k-point=13,k-point=19,k-point=25,kpoint=31,k-point=2,k-point=8,k-point=14.....这种顺序排列的。现在，我想要它按k-point=1,k-point=2,k-point=3,k-point=4......这种顺序排列在txt中。求帮我写一段这样的fortran代码，万分感谢。为方便操作，我把原件上传了上来。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : unit_bands.txt
2013-07-25 17:42:52, 38.96 K
» 猜你喜欢
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有5人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：intel visual fortran的格式问题
已经有14人回复
- fortran，带字符的数据读取与存储！
已经有11人回复
- fortran中计算数据写入文本为 NAN的问题
已经有4人回复
- fortran中如何处理字符和数据混合的文件
已经有6人回复
- Fortran数据读入格式的修改
已经有5人回复
- Fortran如何读取数据避开空行？
已经有7人回复
- 求助fortran编程时怎么使用数据库
已经有4人回复
- 【求助】matlab数据导入Fortran
已经有9人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】FORTRAN写文件时每行开头都多了一个空格？
已经有6人回复
- 【求助】FORTRAN读取文本文件
已经有9人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】初学fortran，
已经有15人回复
- 请教fortran导txt数据进origin问题
已经有1人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】fortran中怎么存储不确定长度的数据【已完结】
已经有6人回复
- 【求助】matlab读取fortran输出的txt文件【已解决】
已经有3人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
3楼2013-07-26 17:26:52
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2013-07-25 20:45:07
alvin510
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1784.3
- 散金: 423
- 红花: 4
- 帖子: 288
- 在线: 225.4小时
- 虫号: 1245360
- 注册: 2011-03-26
- 专业: 人工晶体
|
我把数据读入数组时，用的代码是 do nk = 1, nkpts read(fileid,"(A80)"  dummylocat_str = index(dummy,'K-point') read(dummy(locat_str+7  ,*) ndum, kpoints(1,nk), kpoints(2,nk), kpoints(3,nk), weight(nk)do ns = 1, nspins read(fileid,*) dummy do nb = 1, num_bands read(fileid,"(3XF11.8)" eigenvalues(nb, nk, ns)end do enddo enddo close(fileid) 其中，nk就代表了k-point序列. 当nk=1时能和txt中的k-pont 1对应。但当nk=2时，读入的是k-point 7的数据，而存进数组的序列为nk=2.后面的也跟着对应不了。求帮助我解决这个难题。万分感谢。 |
4楼2013-07-26 17:37:39
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
alvin510: 金币+10, ★★★★★最佳答案 2013-07-27 09:22:48
alvin510: 金币+10, ★★★★★最佳答案 2013-07-27 09:22:48
|
这里 read(dummy(locat_str+7,*) ndum, kpoints(1,nk), kpoints(2,nk), kpoints(3,nk), weight(nk) 应该改成 read(dummy(locat_str+7,*) ndum, kpoints(1,ndum), kpoints(2,ndum), kpoints(3,ndum), weight(ndum) 其他类似的地方也要改 有时间弄这些，花两天时间学学bash脚本，早就弄好了。 |
5楼2013-07-26 21:16:23
