版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

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

kathy2008

木蟲 (正式寫手)

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

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

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

回復此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復
281求調(diào)劑（0805）已經(jīng)有16人回復
304求調(diào)劑已經(jīng)有6人回復
材料工程專碩調(diào)劑已經(jīng)有6人回復
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復
求材料調(diào)劑已經(jīng)有8人回復
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復
一志愿華中科技大學，080502，354分求調(diào)劑已經(jīng)有4人回復
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有6人回復

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

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

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

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

snoopyzhao

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

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

【答案】應助回帖

★ ★
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 則結束整個程序……

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人)

回復此樓

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

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

kathy2008

木蟲 (正式寫手)

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

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

引用回帖:

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

[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也沒有找到應對之策，請高手指點。

贊一下(1人)

回復此樓

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

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

snoopyzhao

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

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

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

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

所以，你把程序中：

CODE:

line = line(21:)

改成

CODE:

line = line(22:)

就可以了……

贊一下(1人)

回復此樓

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

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

snoopyzhao

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

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

★
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人)

回復此樓

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

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

運行后報錯信息如下:
*** 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人)

回復此樓

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

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

snoopyzhao

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

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

★
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))

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

贊一下(1人)

回復此樓

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

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

snoopyzhao

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

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

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

另外，原程序中

CODE:

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

應該改為

CODE:

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

贊一下(1人)

回復此樓

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

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

相關版塊跳轉我要訂閱樓主 kathy2008 的主題更新

返回列表

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

[考研] 一志愿天大材料與化工（085600）總分338 +4	蔡大美女 2026-03-13	4/200	2026-03-18 22:47 by walc
[考研] 一志愿華中科技大學，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +5	小材化本科 2026-03-18	5/250	2026-03-18 21:53 by 無懈可擊111
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 0703化學求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 268求調(diào)劑 +8	一定有學上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

24小時熱門版塊排行榜

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

» 猜你喜歡

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

【答案】應助回帖

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