版塊導(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)劑公告

返回列表

damao4361556

鐵蟲 (正式寫手)

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

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

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

自帶的 plotband.x 等感覺不好用。

回復(fù)此樓

» 猜你喜歡

材料科學(xué)與工程求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有5人回復(fù)
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）已經(jīng)有6人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有21人回復(fù)
26申博自薦已經(jīng)有5人回復(fù)
296求調(diào)劑已經(jīng)有5人回復(fù)
總分293求調(diào)劑已經(jīng)有10人回復(fù)
材料專碩調(diào)劑已經(jīng)有3人回復(fù)
調(diào)劑310 已經(jīng)有12人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(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畫能帶圖已經(jīng)有14人回復(fù)
電子態(tài)密度怎么畫圖�。� 已經(jīng)有4人回復(fù)
【求助】ms能帶圖已經(jīng)有11人回復(fù)
【求助】band.dat用origin畫能帶圖已經(jīng)有5人回復(fù)
【求助】pwscf的聲子態(tài)密度圖和色散曲線圖能用origin畫出來嗎？已經(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的回帖

司馬星空

銀蟲 (初入文壇)

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

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

dos算完的結(jié)果是可以直接導(dǎo)入origin的，band算完的程序要用一個(gè)小程序轉(zhuǎn)換下。下面這個(gè)程序是我從小木蟲找到的。我改了下，有兩個(gè)輸入文件一個(gè)band.dat,還有一個(gè)ev寫費(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的回帖

damao4361556

鐵蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 38.8
紅花: 3
帖子: 479
在線: 435.8小時(shí)
蟲號(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è)程序是我從小木蟲找到的。我改了下，有兩個(gè)輸入文件一個(gè)band.dat,還有一個(gè)ev寫費(fèi)米能級(jí)。
program prog
real, allocatable ...

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

請(qǐng)問為什么我這個(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的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 damao4361556 的主題更新

返回列表

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

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

[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 0856求調(diào)劑 +13	zhn03 2026-03-25	14/700	2026-03-29 08:13 by fmesaito
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 289求調(diào)劑 +13	新時(shí)代材料 2026-03-27	13/650	2026-03-29 01:16 by 544594351
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9	天天好運(yùn)來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +3	蘿卜燉青菜 2026-03-28	3/150	2026-03-28 14:11 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 291求調(diào)劑 +15	hhhhxn.. 2026-03-23	21/1050	2026-03-28 11:26 by self2008
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 070300化學(xué)求調(diào)劑 +4	起個(gè)名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 283求調(diào)劑（080500） +4	A child 2026-03-27	4/200	2026-03-27 15:34 by XPU李慶
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 材料考研求調(diào)劑 +3	Dendel 2026-03-23	6/300	2026-03-26 17:51 by fmesaito
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

damao4361556

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

» 猜你喜歡

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

司馬星空

damao4361556

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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