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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

本帖產(chǎn)生 1 個程序強帖，點擊這里進行查看

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

kathy2008

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2161.5
紅花: 2
帖子: 328
在線: 213.2小時
蟲號: 558502
注冊: 2008-05-15
性別: GG
專業(yè): 材料物理化學(xué)

[求助] 如何從高斯輸出文件快速提出 pai 軌道信息。

如題。從高斯輸出文件提出了eigenvector那一部分出來，即附件1�，F(xiàn)在需要得到 pai 軌道信息。即附件2。附件2 對應(yīng)于附件1的32號，35號，38號，39號，40號，41號軌道（占據(jù)軌道），42號一直到47號（非占據(jù)軌道）的2Px值。求一小程序。請指點。謝謝。

回復(fù)此樓

» 猜你喜歡

279分求調(diào)劑一志愿211 已經(jīng)有16人回復(fù)
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）已經(jīng)有5人回復(fù)
上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境）已經(jīng)有3人回復(fù)
急發(fā)核心期刊論文已經(jīng)有5人回復(fù)
291 求調(diào)劑已經(jīng)有5人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有7人回復(fù)
316求調(diào)劑已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有4人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有7人回復(fù)

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

用高斯計算瞬態(tài)吸收時分子軌道后出現(xiàn)A，B是什么意思已經(jīng)有3人回復(fù)
請教，關(guān)于用gaussview看分子軌道已經(jīng)有8人回復(fù)
【求助】用gaussview 讀取windows版g03計算出來的chk文件時能打開卻不能看MO 軌道已經(jīng)有4人回復(fù)
【求助】如何在以下高斯03輸出文件中得到二階超極化率已經(jīng)有15人回復(fù)
【求助】高斯中怎么計算軌道磁矩已經(jīng)有6人回復(fù)

1樓 2011-06-17 19:54:10

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
jjdg(金幣+1): 感謝參與 2011-07-04 12:46:36

我才注意到，你把我在二樓給出的程序中的

CODE:

read(line,fm) px(i,(j+1):(j+n))

改成了

CODE:

read(line,fm) ppx(i,(j+1):(j+n))

害得我弄了半天才知道為啥出來的結(jié)果總是不對。

贊一下(1人)

回復(fù)此樓

7樓2011-07-04 09:48:55

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

查看全部 8 個回答

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★ ★
ben_ladeng(金幣+2): 很詳細，待樓主評定后獎勵程序強帖 2011-06-18 17:42:21
kathy2008(金幣+10): 2011-06-19 13:00:43
微塵、夢想(程序強帖+1): 2011-06-19 17:04:28

大概這個樣子。只是需要手工輸入軌道號（這樣可能靈活一些），每次輸入一個軌道序號，回車，輸入 0 則結(jié)束整個程序……

CODE:

program ei

real, dimension(:,:), allocatable :: px,ppx

character(len=256) :: line

character(len=40) :: fm

integer :: nrow, ncol, i, j, k, ios

open(unit=12, file='eigenvector.out', status='old')

open(unit=13, file='2px.out', status='new')

do

   read(12,'(a)', iostat=ios) line

   if (ios /= 0) exit

   if (index(line,'EIGENVALUES') /= 0) then

      nrow=0

      ncol=0

      do

         read(12,'(a)', iostat=ios) line

         if (ios /= 0) exit

         if (line(1:4) == '    ') exit

         ncol=ncol+1

         if (index(line, '2PX') /= 0) nrow=nrow+1

      end do

      exit

   end if

end do

!write (*,*) nrow, ncol

rewind (12)

allocate(px(nrow,ncol),ppx(nrow,ncol))

i=0

j=0

do

   read(12,'(a)', iostat=ios) line

   if (ios /= 0) exit

   if (i == nrow) then

       i=0

       j=j+n

   end if

   if (index(line, '2PX') /= 0) then

       line = line(21:)

!      write(*,*) trim(line)

       i=i+1

       n = len_trim(line)/10

       write(fm,'(a,i0,a)') '(', n, 'f10.5)'

!      write(*,*) j

       read(line,fm) px(i,(j+1):(j+n))

   end if

end do

k=0

do

   write(*,*) 'please input a number between 1 and ', nrow, 'end the program by 0.'

   read(*,*) i

   if(i==0) exit

   k=k+1

   ppx(:,k) = px(:,i)

end do

!write(*,*) k/5, mod(k,5)

if (k>=5) then

   do j=1,k/5

      do i=1,nrow

         write(13,'(5f10.5)') ppx(i,(j-1)*5+1:j*5)

      end do

      write(13,*)

   end do

end if

if (mod(k,5) /=0) then

   write(fm,'(a,i0,a)') '(', mod(k,5), 'f10.5)'

   do i=1,nrow

      write(13, fm) ppx(i,(k/5*5+1):k)

   end do

end if

end program ei

贊一下(2人)

回復(fù)此樓

2樓2011-06-18 16:41:08

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

kathy2008

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2161.5
紅花: 2
帖子: 328
在線: 213.2小時
蟲號: 558502
注冊: 2008-05-15
性別: GG
專業(yè): 材料物理化學(xué)

★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-07-31 13:36:05

引用回帖:

Originally posted by snoopyzhao at 2011-06-18 16:41:08:
大概這個樣子。只是需要手工輸入軌道號（這樣可能靈活一些），每次輸入一個軌道序號，回車，輸入 0 則結(jié)束整個程序……

[code]
program ei
real, dimension(:,

, allocatable :: px,ppx
character(len=256 ...

利用該程序提取π軌道信息，報錯。信息如下
At line 48 of file eigen-pai-nc3h7-r2.f
Fortran runtime error: Bad value during floating point read

google也沒有找到應(yīng)對之策，請高手指點。

贊一下(1人)

回復(fù)此樓

3樓2011-07-03 15:06:57

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-07-31 13:36:26

你的 .out 文件是咋生成的，這次的這個文件比上次的文件每一行前面多了一個空格……

所以，你把程序中：

CODE:

line = line(21:)

改成

CODE:

line = line(22:)

就可以了……

贊一下(1人)

回復(fù)此樓

4樓2011-07-03 22:28:17

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

查看全部 8 個回答

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

[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	16/800	2026-03-23 17:45 by 材料學(xué)257求調(diào)劑
[考研] 08工學(xué)調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學(xué)習(xí)學(xué)習(xí)使?/a>
[考研] 336求調(diào)劑 +3	收到VS 2026-03-20	3/150	2026-03-23 13:23 by Iveryant
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6	幸運哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 一志愿武理材料工程348求調(diào)劑 +5	￣^￣゜汗 2026-03-19	7/350	2026-03-22 19:44 by 公瑾逍遙
[考研] 求調(diào)劑 +7	Auroracx 2026-03-22	7/350	2026-03-22 12:38 by 素顏傾城1988
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨

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

24小時熱門版塊排行榜

kathy2008

[求助] 如何從高斯輸出文件快速提出 pai 軌道信息。

» 猜你喜歡

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

snoopyzhao

snoopyzhao

【答案】應(yīng)助回帖

kathy2008

snoopyzhao

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

[求助] 如何從高斯輸出文件快速提出 pai 軌道信息。

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