 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2829 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
astringent銅蟲 (著名寫手)
|
[交流]
【求助】在編譯時出錯 已有5人參與
|
||
|
when i execute the program, the following error appear: At line 111 of file readmd.f90 (unit = 20, file = 'TRAJ1') Fortran runtime error: I/O past end of record on unformatted file How can solve this, please help me, thanks. |
回復(fù)此樓» 猜你喜歡
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有17人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有7人回復(fù)
-
申博26年
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有13人回復(fù)
-
288求調(diào)劑
已經(jīng)有14人回復(fù)
-
307求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- msi2lmp導(dǎo)入MS模型進入lammps時出錯“bash: msi2lmp.exe: command not found”
已經(jīng)有0人回復(fù)
- VASP5.2編譯出錯求助 懇請指點！不勝感激～
已經(jīng)有2人回復(fù)
- 【求助】緊急求助FORTRAN運行問題
已經(jīng)有2人回復(fù)
- 【求助】Siesta安裝報錯
已經(jīng)有3人回復(fù)
- 【求助】VASP計算出錯--新年第一帖
已經(jīng)有1人回復(fù)
- 【求助】molden編譯出錯，求指點
已經(jīng)有1人回復(fù)
- 【求助】科研需要，請人編兩個關(guān)于一維radial distribution function matlab小程序
已經(jīng)有1人回復(fù)
- 【求助】請教一個賦值的問題【已完結(jié)】
已經(jīng)有3人回復(fù)
- 【求助】編譯Abinit6.4.1出錯
已經(jīng)有5人回復(fù)
- 【求助】linux系統(tǒng)下編制程序后，在鏈接生成可執(zhí)行文件時出錯怎么回事？
已經(jīng)有1人回復(fù)
- 【整理】《提問的智慧》——獻給那些不會提問的人
已經(jīng)有78人回復(fù)
- 【求助】FORTRAN中增加計算區(qū)域后出錯，
已經(jīng)有5人回復(fù)
- 【求助】VASP編譯出錯
已經(jīng)有5人回復(fù)
- 【求助】vasp安裝問題
已經(jīng)有27人回復(fù)
- 【求助】mkl+ifort+openmpi+siesta3.并行安裝出錯
已經(jīng)有8人回復(fù)
- 【neweroica個人文集】最新版vaspmo (v0.1) ――可視化VASP分子軌道
已經(jīng)有103人回復(fù)
astringent
銅蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 2210.7
- 散金: 1233
- 紅花: 4
- 帖子: 1180
- 在線: 354.8小時
- 蟲號: 900839
- 注冊: 2009-11-12
- 專業(yè): 理論和計算化學(xué)
|
astringent 新蟲 (小有名氣) 程序強帖: 0 金幣: 72.5 帖子: 112 在線: 14.9小時 蟲號: 900839 注冊: 2009-11-12 狀態(tài): 在線 下面是各個子程序 subroutine findhbond & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) implicit none include 'constants.h' integer io,ih,iat,nattot,ihbnum double precision xxx(nattot),yyy(nattot),zzz(nattot) double precision cell(3) double precision theta,thetamax,roamax,rhamax double precision dx0,dy0,dz0,r0 double precision dx1,dy1,dz1,r1 double precision dx2,dy2,dz2,r2 thetamax = 30.d0 roamax = 3.5d0 rhamax = 2.45d0 ihbnum = 0 call bondcalculation(io,ih,xxx,yyy,zzz,nattot,cell,dx0,dy0,dz0,r0) call bondcalculation(io,iat,xxx,yyy,zzz,nattot,cell,dx1,dy1,dz1,r1) if (r1.lt.roamax) then theta = dacos(dx0*dx1+dy0*dy1+dz0*dz1) if(theta.lt.thetamax) then call bondcalculation & (ih,iat,xxx,yyy,zzz,nattot,cell,dx2,dy2,dz2,r2) if (r2.lt.rhamax) then ihbnum = 1 endif endif endif return end subroutine test & (io,ih,ia,xxx,yyy,zzz,nattot,cell,atomname,resname,ihbnum) implicit none include 'constants.h' integer io,ih,ia,nattot,ihbnum double precision xxx(nattot),yyy(nattot),zzz(nattot) double precision cell(*) character*4 atomname(nattot),resname(nattot) integer iat ia = 0 ihbnum = 0 do iat=1, nattot if((('resname'.eq.'DC5').or.('resname'.eq.'DC')).and.& ('atomname'.eq.'N1'.or.'atomname'.eq.'N4'.or.'atomname'.eq.'N3'& &.or.'atomname'.eq.'O2')) then call findhbond & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) elseif (('resname'.eq.'DG3'.or.'resname'.eq.'DG').and.& &('atomname'.eq.'N9'.or.'atomname'.eq.'N7'.or.'atomname'.eq.& &'O6'.or.'atomname'.eq.'N1'.or.'atomname'.eq.'N2'& &.or.'atomname'.eq.'N3')) then call findhbond & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) elseif (('resname'.eq.'DA').and.('atomname'.eq.'N9'.or.& 'atomname'.eq.'N7'.or.'atomname'.eq. & 'N6'.or.'atomname'.eq.'N1'& &.or.'atomname'.eq.'N3')) then call findhbond & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) elseif (('resname'.eq.'DT').and.('atomname'.eq.'N1'.or. & 'atomname'.eq.'O4'.or.'atomname' & .eq.'N3'.or.'atomname'.eq.'O2')) then call findhbond & (io,ih,iat,xxx,yyy,zzz,nattot,cell,ihbnum) endif if (ihbnum.eq.1) then ia=iat endif enddo return end subroutine bondcalculation & (i1,i2,xxx,yyy,zzz,nattot,cell,dx,dy,dz,roh) implicit none integer i1,i2,nattot double precision xxx(nattot),yyy(nattot),zzz(nattot) double precision cell(3),dx,dy,dz double precision xo,yo,zo,xh,yh,zh,roh xo = xxx(i1) yo = yyy(i1) zo = zzz(i1) xh = xxx(i2) yh = yyy(i2) zh = zzz(i2) dx = xh - xo dy = yh - yo dz = zh - zo dx = dx - anint(dx/cell(1))*cell(1) dy = dy - anint(dy/cell(2))*cell(2) dz = dz - anint(dz/cell(3))*cell(3) roh = dsqrt(dx**2+dy**2+dz**2) dx = dx / roh dy = dy / roh dz = dz / roh return end subroutine readmd & (mdtype,nattot,nin,atomname,xxx,yyy,zzz,cell) implicit none integer mdtype,nattot,nin double precision :: xxx(nattot),yyy(nattot),zzz(nattot) real*4 :: x_dcd(nattot),y_dcd(nattot),z_dcd(nattot) double precision cell(3) character*4 atomname(nattot) integer iatm,i double precision xtable(6) integer ::npdb=70 if (mdtype.eq.0) then ! original TINKER read(nin,*) do iatm=1, nattot read(nin,'(i6,2x,a2,x,3f12.6)')i,atomname(iatm),xxx(iatm),& yyy(iatm),zzz(iatm) enddo elseif (mdtype.eq.1) then ! modified TINKER read(nin,*) read(nin,*) do iatm=1, nattot read(nin,'(a2,9x,6f12.4)')atomname(iatm),xxx(iatm), & yyy(iatm), zzz(iatm) enddo elseif ((mdtype.eq.2).or.(mdtype.eq.3)) then ! DLPOLY read(nin,*) read(nin,*) read(nin,*) read(nin,*) do iatm=1, nattot read(nin,'(a2,x,6f10.5)')atomname(iatm),xxx(iatm), & yyy(iatm), zzz(iatm) enddo elseif (mdtype.eq.4) then ! Orsay - FX COUDERT read(nin,*) read(nin,*) do iatm=1, nattot read(nin,'(2x,a2,6x,3f14.5)')atomname(iatm),xxx(iatm),& yyy(iatm),zzz(iatm) enddo elseif (mdtype.eq.5) then ! xyz read(nin,*) read(nin,*) do iatm=1, nattot read(nin,'(a2,8x,3f12.4)')atomname(iatm),xxx(iatm),& yyy(iatm),zzz(iatm) enddo elseif (mdtype.eq.6) then ! dcd binary read(nin) (xtable(i),i=1,6) read(nin) (x_dcd(i),i=1,nattot) read(nin) (y_dcd(i),i=1,nattot) read(nin) (z_dcd(i),i=1,nattot) cell(1)=xtable(1) cell(2)=xtable(3) cell(3)=xtable(6) OPEN(npdb,file='system.pdb',status='old') call read_pdb(atomname,npdb,nattot) CLOSE(npdb) ! send from single precision to double do i=1,nattot xxx(i)=x_dcd(i) yyy(i)=y_dcd(i) zzz(i)=z_dcd(i) enddo endif return end subroutine readmdinit(mdtype,nin) implicit none integer mdtype,nin ! local variables -------------- integer ntitle,nset,istart,icntrl(6),jcntrl(9) character*80 title(10) character*4 hdr double precision delta integer i,nfr,dfr,natms if (mdtype.eq.3) then read(nin,*) read(nin,*) elseif(mdtype.eq.6) then read(nin) hdr,nfr,istart,dfr,(icntrl(i),i=1,6),delta,(jcntrl(i) & ,i=1,9) read(nin) ntitle,(title(i),i=1,ntitle) read(nin) natms endif return end subroutine read_pdb(atomname,n,nattot) integer ::n,nattot character(len=4),dimension(nattot),intent(inout) ::atomname character(len=4),dimension(nattot) ::resname character(len=4),dimension(nattot) ::itmp integer,dimension(nattot) ::resindex,label !--- local variables ------------------- integer :: i do i=1,nattot read(n,1000,err=20,end=20) itmp,label(i),atomname(i),resname(i),& resindex(i) enddo 20 CONTINUE 1000 format(3x,a4,6x,i5,a4,x,a4,3x,i5) |
4樓2010-11-16 03:22:17
2樓2010-11-15 22:02:38
astringent
銅蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 2210.7
- 散金: 1233
- 紅花: 4
- 帖子: 1180
- 在線: 354.8小時
- 蟲號: 900839
- 注冊: 2009-11-12
- 專業(yè): 理論和計算化學(xué)
|
這個是源文件 open(nin,file='TRAJ1',status='old',form & ='unformatted') open(nout,file='csbw.out') do istep=1, nstep call readmdinit (mdtype,nin) call readmd & (mdtype,nattot,nin,atomname,xxx,yyy,zzz,cell) ! loop over all waters do i=1, nattot if (('atomname'.eq.'O').and.('resname'.eq.'WAT')) then io = i ih1 = io+1 ih2 = io+2 call test & (io,ih1,ia1,xxx,yyy,zzz,nattot,cell,atomname,resname,ihbnum) if (ihbnum.eq.1) then acc(j)=ia1 call test & (io,ih2,ia2,xxx,yyy,zzz,nattot,cell, & atomname,resname,ihbnum) if (ihbnum.eq.1) then if (resindex(ia1).eq.resindex(ia2)+1 & .or. resindex(ia1).eq. resindex(ia2)-1) then write(*,*) 'not find a cswb' else write(*,*) 'find a cswb' nhb=nhb+1 endif endif endif endif enddo sum=sum+nhb do j=1,nattot write (nout,"('acceptor is',i2)"  acc(j)enddo enddo avgnhb=sum/nstep write (nout,"('avergaed number is',i2)" avgnhbclose(nout) end |
3樓2010-11-16 03:01:44
5樓2010-11-17 02:03:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +8 | 冷笙123 2026-03-17 | 8/400 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|