 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1147 | 回復(fù): 5 | ||||
alvin510銀蟲 (小有名氣)
|
[求助]
求txt中數(shù)據(jù)重新排序的fortran代碼
|
| 我用MS計(jì)算了一個(gè)晶體的能帶結(jié)構(gòu)，其結(jié)果保存在了unitbands.txt文件中。不幸的是，其中的結(jié)果沒有按自然順序排列，即原文件是按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.....這種順序排列的�，F(xiàn)在，我想要它按k-point=1,k-point=2,k-point=3,k-point=4......這種順序排列在txt中。求幫我寫一段這樣的fortran代碼，萬(wàn)分感謝。為方便操作，我把原件上傳了上來。 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : unit_bands.txt
2013-07-25 17:42:52, 38.96 K
» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有7人回復(fù)
-
抱歉
已經(jīng)有4人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助：intel visual fortran的格式問題
已經(jīng)有14人回復(fù)
- fortran，帶字符的數(shù)據(jù)讀取與存儲(chǔ)！
已經(jīng)有11人回復(fù)
- fortran中計(jì)算數(shù)據(jù)寫入文本為 NAN的問題
已經(jīng)有4人回復(fù)
- fortran中如何處理字符和數(shù)據(jù)混合的文件
已經(jīng)有6人回復(fù)
- Fortran數(shù)據(jù)讀入格式的修改
已經(jīng)有5人回復(fù)
- Fortran如何讀取數(shù)據(jù)避開空行？
已經(jīng)有7人回復(fù)
- 求助fortran編程時(shí)怎么使用數(shù)據(jù)庫(kù)
已經(jīng)有4人回復(fù)
- 【求助】matlab數(shù)據(jù)導(dǎo)入Fortran
已經(jīng)有9人回復(fù)
- 【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離
已經(jīng)有34人回復(fù)
- 【求助】FORTRAN寫文件時(shí)每行開頭都多了一個(gè)空格？
已經(jīng)有6人回復(fù)
- 【求助】FORTRAN讀取文本文件
已經(jīng)有9人回復(fù)
- 【求助】求高手關(guān)于Fortran數(shù)據(jù)讀寫的問題�！疽淹杲Y(jié)】
已經(jīng)有12人回復(fù)
- 【求助】初學(xué)fortran，
已經(jīng)有15人回復(fù)
- 請(qǐng)教fortran導(dǎo)txt數(shù)據(jù)進(jìn)origin問題
已經(jīng)有1人回復(fù)
- 【求助】求指點(diǎn)，在Fortran里面，怎樣實(shí)現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組?
已經(jīng)有5人回復(fù)
- 【求助】請(qǐng)問一下fortran的可視化編程
已經(jīng)有7人回復(fù)
- 【求助】fortran中怎么存儲(chǔ)不確定長(zhǎng)度的數(shù)據(jù)【已完結(jié)】
已經(jīng)有6人回復(fù)
- 【求助】matlab讀取fortran輸出的txt文件【已解決】
已經(jīng)有3人回復(fù)
- 【求助】有關(guān)fortran進(jìn)行數(shù)值模擬
已經(jīng)有10人回復(fù)
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時(shí)
- 蟲號(hào): 1023452
- 注冊(cè): 2010-05-19
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
2樓2013-07-25 20:45:07
alvin510
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1784.3
- 散金: 423
- 紅花: 4
- 帖子: 288
- 在線: 225.4小時(shí)
- 蟲號(hào): 1245360
- 注冊(cè): 2011-03-26
- 專業(yè): 人工晶體
3樓2013-07-26 17:26:52
alvin510
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1784.3
- 散金: 423
- 紅花: 4
- 帖子: 288
- 在線: 225.4小時(shí)
- 蟲號(hào): 1245360
- 注冊(cè): 2011-03-26
- 專業(yè): 人工晶體
|
我把數(shù)據(jù)讀入數(shù)組時(shí)，用的代碼是 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序列. 當(dāng)nk=1時(shí)能和txt中的k-pont 1對(duì)應(yīng)。但當(dāng)nk=2時(shí)，讀入的是k-point 7的數(shù)據(jù)，而存進(jìn)數(shù)組的序列為nk=2.后面的也跟著對(duì)應(yīng)不了。求幫助我解決這個(gè)難題。萬(wàn)分感謝。 |
4樓2013-07-26 17:37:39
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時(shí)
- 蟲號(hào): 1023452
- 注冊(cè): 2010-05-19
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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) 應(yīng)該改成 read(dummy(locat_str+7,*) ndum, kpoints(1,ndum), kpoints(2,ndum), kpoints(3,ndum), weight(ndum) 其他類似的地方也要改 有時(shí)間弄這些，花兩天時(shí)間學(xué)學(xué)bash腳本，早就弄好了。 |
5樓2013-07-26 21:16:23
rsb_long
至尊木蟲 (著名寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 19110.9
- 帖子: 1268
- 在線: 405.3小時(shí)
- 蟲號(hào): 80232
- 注冊(cè): 2005-07-14
- 專業(yè): 計(jì)算機(jī)軟件
6樓2013-07-28 08:54:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 alvin510 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 0703一志愿9，初試成績(jī)：338，四六級(jí)已過，有科研經(jīng)歷，求調(diào)劑！ +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +9 | XX小鄧 2026-03-29 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 10/500 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|