| 5 | 1/1 | 返回列表 |
| 查看: 1149 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
alvin510銀蟲 (小有名氣)
|
[求助]
求txt中數(shù)據(jù)重新排序的fortran代碼
|
|||
| 我用MS計算了一個晶體的能帶結(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代碼,萬分感謝。為方便操作,我把原件上傳了上來。 |
銀蟲 (小有名氣)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +14 |
銀蟲 (小有名氣)
|
我把數(shù)據(jù)讀入數(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時能和txt中的k-pont 1對應(yīng)。但當(dāng)nk=2時,讀入的是k-point 7的數(shù)據(jù),而存進(jìn)數(shù)組的序列為nk=2.后面的也跟著對應(yīng)不了。求幫助我解決這個難題。萬分感謝。 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +14 |
|
這里 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) 其他類似的地方也要改 有時間弄這些,花兩天時間學(xué)學(xué)bash腳本,早就弄好了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +5 | 小江小江江江 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑求調(diào)劑 +6 | 121. 2026-04-02 | 6/300 |
|
|
[考研] 286分調(diào)劑 +20 | Faune 2026-03-30 | 22/1100 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +10 | 月月鳥木 2026-04-01 | 10/500 |
|
|
[考研] 329求調(diào)劑 +8 | miaodesi 2026-04-02 | 10/500 |
|
|
[考研] 321求調(diào)劑 +9 | y-yh 2026-04-01 | 10/500 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +11 | @taotao 2026-03-29 | 11/550 |
|
|
[考研] 354求調(diào)劑 +4 | lxb598 2026-03-31 | 5/250 |
|
|
[考研] 化學(xué)工程專碩324分,一志愿中國礦業(yè)大學(xué)求調(diào)劑 +7 | 耿耿1314 2026-04-01 | 7/350 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +7 | 哇呼哼呼哼 2026-03-26 | 7/350 |
|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +14 | 小黃鴨寶 2026-03-30 | 14/700 |
|
|
[考研] 0710生物學(xué)求調(diào)劑! +6 | 敘述文 2026-03-31 | 6/300 |
|
|
[考研] 一志愿a區(qū)211,085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
[考研] 293分求調(diào)劑,外語為俄語 +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|