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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

本帖產(chǎn)生 1 個程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

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值。求一小程序。請指點(diǎn)。謝謝。

回復(fù)此樓

» 猜你喜歡

本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑已經(jīng)有3人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有11人回復(fù)
求材料調(diào)劑已經(jīng)有9人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有10人回復(fù)
271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(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è)作家)

程序強(qiáng)帖: 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): 很詳細(xì)，待樓主評定后獎勵程序強(qiáng)帖 2011-06-18 17:42:21
kathy2008(金幣+10): 2011-06-19 13:00:43
微塵、夢想(程序強(qiáng)帖+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)對之策，請高手指點(diǎn)。

贊一下(1人)

回復(fù)此樓

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

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

snoopyzhao

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

程序強(qiáng)帖: 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的回帖

snoopyzhao

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

程序強(qiáng)帖: 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:35

這樣可能更好一些：

CODE:

      program ei

      implicit none

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

      character(len=256) :: line

      character(len=40) :: fm

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

      open(unit=12, file='nc3h7-r2-sto-eiv.out', status='old')

      open(unit=13, file='nc3h7-r2-sto-pai.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, '2PZ') /= 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, '2PZ') /= 0) then

             line = line(21:)

!            write(*,*) trim(line)

             i=i+1

             n = len_trim(line)/10

             m = mod(len_trim(line),10)

!            write (*,*) m, n

             if (m /= 0) then

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

             else

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

             end if

!            write (*,*) fm

!            write(*,*) j

             read(line,fm) ppx(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

贊一下(1人)

回復(fù)此樓

5樓2011-07-03 22:59:20

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

運(yùn)行后報錯信息如下:
*** glibc detected *** ./a.out: double free or corruption (out): 0x00000000079e2980 ***
======= Backtrace: =========
/lib64/libc.so.6[0x3e64871ce2]
/lib64/libc.so.6(cfree+0x8c)[0x3e6487590c]
./a.out[0x4012c5]
./a.out[0x4019de]
/lib64/libc.so.6(__libc_start_main+0xf4)[0x3e6481d974]
./a.out[0x400a99]

贊一下(1人)

回復(fù)此樓

6樓2011-07-04 07:23:29

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

snoopyzhao

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

程序強(qiáng)帖: 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的回帖

snoopyzhao

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

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

★ ★
jjdg(金幣+2): 辛苦了 2011-07-04 12:46:22

另外，原程序中

CODE:

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

應(yīng)該改為

CODE:

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

贊一下(1人)

回復(fù)此樓

8樓2011-07-04 09:50:20

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

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

返回列表

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

[考研] 085600材料與化工調(diào)劑 324分 +9	llllkkkhh 2026-03-18	9/450	2026-03-19 07:40 by BruceLiu320
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	12/600	2026-03-18 15:27 by 咪嗚喵嗚
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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