国内唯一诱惑网站男人的天堂久久 ,国产xxx在线视频观看,四肖中特期期特期准王中王特色

版塊導(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í)熱門版塊排行榜

返回列表

上一頁(yè)

秋雨子

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時(shí)
蟲(chóng)號(hào): 1027042
注冊(cè): 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

引用回帖:

18樓: Originally posted by asaka at 2011-12-10 16:33:17:
以上作者只考慮閉殼層計(jì)算的情況，你有開(kāi)殼層的計(jì)算嗎？那時(shí)候你需要總的HOMO,LUMO還是分別alpha以及beta的HOMO和LUMO?
他們的程序都只有一點(diǎn)小問(wèn)題，相信LZ可以搞定。
mchen10 的程序少了一個(gè) |，多了一個(gè)奇怪 ...

非常感謝高手的賜教，我的分子群全是閉殼層的，當(dāng)然也有許多別的分子信息需要提取，那些都相對(duì)簡(jiǎn)單一些，自己可以學(xué)著搞定！非常感謝！

贊一下

回復(fù)此樓

To the time to life, rather than to life in time

21樓2011-12-11 20:10:47

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

秋雨子

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時(shí)
蟲(chóng)號(hào): 1027042
注冊(cè): 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

送鮮花一朵

引用回帖:

19樓: Originally posted by mchen10 at 2011-12-11 04:06:32:
[code]#!/bin/bash
# print homo lumo energy of gaussian output
# usage: print_homo_lumo filename(s)
for file in "$@"
do
echo $file
tac $file | tr "\n" "@" | egre ...

謝謝不吝賜教！我會(huì)努力的學(xué)習(xí)的，才接觸腳本，總是覺(jué)得深?yuàn)W難懂，心理的排斥應(yīng)該克服！總是想繞過(guò)難題去走捷徑是不行的！非常感謝！

贊一下

回復(fù)此樓

To the time to life, rather than to life in time

22樓2011-12-11 20:13:09

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

ZDBWHZ

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 7 (幼兒園)
金幣: 1882
紅花: 7
帖子: 736
在線: 228.7小時(shí)
蟲(chóng)號(hào): 655151
注冊(cè): 2008-11-15
專業(yè): 理論和計(jì)算化學(xué)

`ls -l *.log | awk '{print $9}' | sed 's/.log//g' > file_list`
這一行中的9可能是8，需要自己數(shù)一下ls -l *.log出現(xiàn)多少列。

贊一下

回復(fù)此樓

23樓2011-12-12 13:28:55

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

ZDBWHZ

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 7 (幼兒園)
金幣: 1882
紅花: 7
帖子: 736
在線: 228.7小時(shí)
蟲(chóng)號(hào): 655151
注冊(cè): 2008-11-15
專業(yè): 理論和計(jì)算化學(xué)

再來(lái)個(gè)Fortran
PROGRAM HOLU2
!
! Program to extract HOMO and LUMO data from multiple Gaussian output files
! SRK 2011/06/30
! SRK 2011/10/18 Modified to produce 3 output files, bug fix
!
! Input : A file containing a list of filenames, one per line
! Output:  One line per input filename, 3 columns, holding the filename, the HOMO and the LUMO
!       Input filename, HOMO and LUMO are written to file 'HOMOandLUMO.txt'
!       HOMO is written to file 'HOMO.txt'
!       LUMO is written to file 'LUMO.txt'
!
IMPLICIT NONE
! Variable declarations
INTEGER, PARAMETER :: inunit=22, listunit=23 ! Unit numbers for input files and list file
INTEGER, PARAMETER :: hlunit=24, hunit=25, lunit=26 ! Unit numbers for output files
INTEGER :: i       ! Loop variable
INTEGER :: n       ! Number of characters in current line
INTEGER :: stat    ! Holds the current state of the input file, 0=OK
INTEGER :: numfiles ! Number of filenames detected in the input
INTEGER :: blockflag  ! Flag variable to find the correct HOMO-LUMO data
CHARACTER (LEN=100) :: listname ! File containing the list of filenames to be processed
CHARACTER (LEN=100) :: inpfile  ! Current input filename
CHARACTER (LEN=100) :: line    ! Current line in the input file
CHARACTER (LEN=100) :: prevline ! Previous line in the input file
CHARACTER (LEN=9) :: homostring, lumostring ! Strings containing the HOMO/LUMO eigenvalues

!
! Executable code

WRITE (*,*) 'HOLU - a program to extract HOMO and LUMO data'
WRITE (*,*) 'Enter the name of an input file containing a list of filenames (no blank lines!)'
READ (*,'(A)') listname
OPEN(listunit,FILE=listname)
!
! Find out how many files are in the list
numfiles = 0
stat = 0
DO
READ(listunit,'(A)',iostat=stat) line
IF (stat /= 0) EXIT
numfiles = numfiles + 1
END DO
CLOSE(listunit)
WRITE (*,*) 'Detected ',numfiles,' lines in list file'

!
! Open output files
OPEN(hlunit,FILE='HOMOandLUMO.txt')
OPEN(hunit,FILE='HOMO.txt')
OPEN(lunit,FILE='LUMO.txt')

! Now loop over these files
OPEN(listunit,FILE=listname)
DO i = 1, numfiles
READ(listunit,'(A)') inpfile
OPEN(inunit,file=TRIM(inpfile))
blockflag = 0
!
! Search for the needed information in the input file
prevline = ''
DO
   READ(inunit,FMT='(A)',IOSTAT=stat) line ! Read in the next line of the file using '(A)' format
   IF (stat /= 0 ) EXIT
   IF (index(line,'Alpha virt.') > 0 ) THEN ! found a line listing virtual states
      IF (blockflag == 0) blockflag = 1
      IF (blockflag == 2) THEN                   ! we have found the second block
         lumostring = line(30:40)                ! first number on this line
         n = len_trim(prevline)
         homostring = prevline(n-8:n)          ! last number on previous line
         WRITE(*,*) TRIM(inpfile),' ', 'HOMO', ' ', homostring,' ', 'LUMO', ' ', lumostring
         WRITE(hlunit,'(A),(A),(A),(A)') TRIM(inpfile),' ', 'HOMO', ' ', homostring,' ', 'LUMO', ' ', lumostring
         WRITE(hunit,'(A)') TRIM(inpfile),' ', 'HOMO', ' ', homostring
         WRITE(lunit,'(A)') TRIM(inpfile),' ', 'LUMO', ' ', lumostring
         EXIT
      END IF
   ELSE
      IF (blockflag == 1) blockflag = 2
   END IF
   prevline = line
END DO
CLOSE(inunit)
END DO

! Close all files that are still open
CLOSE(listunit)
CLOSE(hlunit)
CLOSE(hunit)
CLOSE(lunit)

!
END PROGRAM HOLU2
這是source code,需要編譯。
編譯后，`ls -l *.log | awk '{print $9}' > file_list
再運(yùn)行編譯后的程序，輸入文件就是file_list。所有的輸出結(jié)果都顯示在屏幕上，可以拷入txt文檔，進(jìn)行后續(xù)的計(jì)算。

贊一下

回復(fù)此樓

24樓2011-12-12 13:37:45

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

秋雨子

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時(shí)
蟲(chóng)號(hào): 1027042
注冊(cè): 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

高手如云啊！學(xué)習(xí)了學(xué)習(xí)了！可以成為經(jīng)典貼了！

贊一下

回復(fù)此樓

To the time to life, rather than to life in time

25樓2011-12-13 23:43:07

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

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

返回列表

上一頁(yè)

普通表情龍兔虎貓

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

[考研] 301求調(diào)劑 +5	yy要上岸呀 2026-03-17	5/250	2026-03-17 20:20 by peike
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 304求調(diào)劑 +8	小熊joy 2026-03-14	8/400	2026-03-17 17:29 by ruiyingmiao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 求調(diào)劑（材料與化工327） +4	愛(ài)吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655

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

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

秋雨子

秋雨子

ZDBWHZ

ZDBWHZ

秋雨子

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