版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

damao4361556

鐵蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 38.8
紅花: 3
帖子: 479
在線: 435.8小時(shí)
蟲(chóng)號(hào): 917440
注冊(cè): 2009-12-01
專業(yè): 理論和計(jì)算化學(xué)

[交流] PWSCF 算完能帶和 dos 后的數(shù)據(jù) 如何轉(zhuǎn)換成origin 能打開(kāi)的格式？已有1人參與

PWSCF 算完能帶和 dos 后的數(shù)據(jù) 如何轉(zhuǎn)換成origin 能打開(kāi)的格式？
求分享有關(guān)的小代碼

自帶的 plotband.x 等感覺(jué) 不好用。

回復(fù)此樓

» 猜你喜歡

0856求調(diào)劑已經(jīng)有13人回復(fù)
0703化學(xué) 已經(jīng)有11人回復(fù)
329求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有8人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有12人回復(fù)
309求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有4人回復(fù)
330分求調(diào)劑已經(jīng)有3人回復(fù)
299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

如何選擇能帶計(jì)算軟件和入門能帶計(jì)算已經(jīng)有7人回復(fù)
學(xué)習(xí)PWSCF(二)：態(tài)密度的計(jì)算已經(jīng)有7人回復(fù)
pwscf-5.0得到的聲子譜很亂是怎么回事？已經(jīng)有8人回復(fù)
QE-PWscf計(jì)算聲子譜流程。已經(jīng)有10人回復(fù)
Problem - nscf - c_bands: eigenvalues not converged 已經(jīng)有20人回復(fù)
[軟件交流]您在使用哪些免費(fèi)視圖、分析軟件呢？已經(jīng)有23人回復(fù)
pwscf4.3.2畫(huà)能帶圖已經(jīng)有14人回復(fù)
電子態(tài)密度怎么畫(huà)圖�。� 已經(jīng)有4人回復(fù)
【求助】ms能帶圖已經(jīng)有11人回復(fù)
【求助】band.dat用origin畫(huà)能帶圖已經(jīng)有5人回復(fù)
【求助】pwscf的聲子態(tài)密度圖和色散曲線圖能用origin畫(huà)出來(lái)嗎？已經(jīng)有5人回復(fù)
【重點(diǎn)討論】計(jì)算能帶時(shí)如何確定高對(duì)稱點(diǎn)和k-points path? 已經(jīng)有40人回復(fù)

1樓 2015-05-01 06:00:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

damao4361556

鐵蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 38.8
紅花: 3
帖子: 479
在線: 435.8小時(shí)
蟲(chóng)號(hào): 917440
注冊(cè): 2009-12-01
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by 司馬星空 at 2015-05-01 20:58:08
dos算完的結(jié)果是可以直接導(dǎo)入origin的，band算完的程序要用一個(gè)小程序轉(zhuǎn)換下。下面這個(gè)程序是我從小木蟲(chóng)找到的。我改了下，有兩個(gè)輸入文件一個(gè)band.dat,還有一個(gè)ev寫(xiě)費(fèi)米能級(jí)。
program prog
real, allocatable ...

你好非常感謝您的回復(fù)

請(qǐng)問(wèn) 為什么我這個(gè) 編譯總出錯(cuò)呢？是直接編譯 gfortran ***.f 嗎？

[wl_903@hpc-login3 band]$ gfortran Getband.f
Getband.f:1.2:

program prog
  1
Error: Non-numeric character in statement label at (1)
Getband.f:1.2:

program prog
  1
Error: Unclassifiable statement at (1)
Getband.f:2.31:

      real, allocatable :: e(:,
                           1
Error: Missing ')' in statement at or before (1)
Getband.f:3.42:

      real, allocatable :: k(:, ,dk( ,x(
                                       1
Error: Missing ')' in statement at or before (1)
Getband.f:14.18:

      allocate(k(nk,3))
               1
Error: Syntax error in ALLOCATE statement at (1)
Getband.f:15.18:

      allocate(e(nk,nbands))
               1
Error: Syntax error in ALLOCATE statement at (1)
Getband.f:16.19:

      allocate(dk(0:nk))
               1
Error: Syntax error in ALLOCATE statement at (1)
Getband.f:17.18:

      allocate(x(1:nk))
               1
Error: Syntax error in ALLOCATE statement at (1)
Getband.f:20.16:

   read(10,*) (k(i,j),j=1,3)
            1
Error: Expected variable in READ statement at (1)
Getband.f:21.16:

   read(10,*) (e(i,n),n=1,nbands)
            1
Error: Expected variable in READ statement at (1)
Getband.f:26.11:

         dk(0)=0
         1
Error: Unclassifiable statement at (1)
Getband.f:28.72:

      k0=k(n,
                                                                     1
Error: Syntax error in argument list at (1)
Getband.f:30.72:

   a=k(n,-k0
                                                                     1
Error: Syntax error in argument list at (1)
Getband.f:31.12:

   dk(n)=dk(n-1)+sqrt(dot_product(a,a))
         1
Error: Statement function at (1) is recursive
Getband.f:32.72:

   x(n)=dk(n)
                                                                     1
Error: Unexpected STATEMENT FUNCTION statement at (1)
Getband.f:33.72:

   k0=k(n,
                                                                     1
Error: Syntax error in argument list at (1)
Getband.f:38.6:

8888 format (1x, f9.4, 2x ,1000(f12.6,2x))
   1
Error: Bad continuation line at (1)
Getband.f:38.6:

8888 format (1x, f9.4, 2x ,1000(f12.6,2x))
   1
Error: Unclassifiable statement at (1)
Getband.f:39.22:

   end program prog
                  1
Error: Syntax error in END PROGRAM statement at (1)
Error: Unexpected end of file in 'Getband.f'

需要當(dāng)前目錄存在兩個(gè)文件是嗎？一個(gè) 是ev 里面的數(shù)據(jù)是費(fèi)米能級(jí)  另外一個(gè)是 band.dat 里面的內(nèi)容是 QE 算能帶后得到的？

贊一下

回復(fù)此樓

3樓2015-05-04 06:13:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

司馬星空

銀蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 319
散金: 120
紅花: 3
帖子: 36
在線: 323.5小時(shí)
蟲(chóng)號(hào): 2497523
注冊(cè): 2013-06-05
專業(yè): 心理學(xué)史

★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖

dos算完的結(jié)果是可以直接導(dǎo)入origin的，band算完的程序要用一個(gè)小程序轉(zhuǎn)換下。下面這個(gè)程序是我從小木蟲(chóng)找到的。我改了下，有兩個(gè)輸入文件一個(gè)band.dat,還有一個(gè)ev寫(xiě)費(fèi)米能級(jí)。
program prog
real, allocatable :: e(:,

real, allocatable :: k(:,

,dk(

,x(

real, dimension(3) ::k0,a
real ev
character(len=32):: xx, yy, zz,evv
integer :: nbands, i, j, n, nk
open(10, file='band.dat', status='old')
open(22, file='ev', status='old')
read(10,*) xx, yy, nbands, zz, nk
read(22,*) evv
read(evv,*) ev
write(*,*)"nk=",nk
allocate(k(nk,3))
allocate(e(nk,nbands))
allocate(dk(0:nk))
allocate(x(1:nk))
open(11,file='phoband.dat')
   do i=1,nk
   read(10,*) (k(i,j),j=1,3)
   read(10,*) (e(i,n),n=1,nbands)

      enddo

            do n=1,nk
         dk(0)=0
         if (n.eq.1) then
      k0=k(n,

endif
a=k(n,

-k0
   dk(n)=dk(n-1)+sqrt(dot_product(a,a))
   x(n)=dk(n)
   k0=k(n,

   enddo
   do i=1,nk
      write(11, 8888) x(i), (e(i,j)-ev,j=1,nbands)
   enddo
8888 format (1x, f9.4, 2x ,1000(f12.6,2x))
   end program prog

贊一下

回復(fù)此樓

2樓2015-05-01 20:58:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +12	稻妻小編 2026-03-26	12/600	2026-03-29 01:04 by 我是小康
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 材料與化工272求調(diào)劑 +9	阿斯蒂芬2004 2026-03-28	9/450	2026-03-28 15:21 by 1018329917
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 調(diào)劑 +3	好好讀書(shū)。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +9	vv迷 2026-03-22	9/450	2026-03-27 15:59 by 不吃魚(yú)的貓
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 343求調(diào)劑 +4	贈(zèng)我一本書(shū) 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4	我愛(ài)學(xué)電池 2026-03-23	4/200	2026-03-25 00:59 by 1027_324
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt