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

客戶端APP下載

論文輔導(dǎo)

調(diào)劑小程序

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(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ù)此樓

» 猜你喜歡

一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑已經(jīng)有3人回復(fù)
南昌大學(xué)材料專碩311分求調(diào)劑已經(jīng)有6人回復(fù)
316求調(diào)劑已經(jīng)有6人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)
一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！已經(jīng)有3人回復(fù)
310求調(diào)劑已經(jīng)有3人回復(fù)
機械專碩299求調(diào)劑至材料已經(jīng)有4人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有7人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有12人回復(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é)

★
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 個回答

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: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的回帖

查看全部 8 個回答

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

[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿福大288有機化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 0854，計算機類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy

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

24小時熱門版塊排行榜

kathy2008

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

» 猜你喜歡

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

snoopyzhao

snoopyzhao

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

kathy2008

snoopyzhao

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

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

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