 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3503 | 回復(fù): 9 | ||
| 【懸賞金幣】回答本帖問題，作者量子海洋將贈送您 50 個金幣 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
量子海洋新蟲 (小有名氣)
|
[求助]
R.I.N.G.S. code運行報錯
|
|
|
用VASP跑完AIMD后，想用RINGS(https://rings-code.sourceforge.net/)處理下數(shù)據(jù)，按照官網(wǎng)的example寫好輸入input（重命名為my-input）和options文件，也吧XDATCAR放在了data文件夾里，在終端中運行 rings my-input 出現(xiàn)報錯“At line 230 of file input.f90 (unit = 10, file = 'options') Fortran runtime error: Bad logical value while reading item 2”，不知道問題出在哪里，怎么解決？ |
回復(fù)此樓» 猜你喜歡
中國林科院林化所（南京）2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名
已經(jīng)有3人回復(fù)
-
306求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料工程，初試313調(diào)劑
已經(jīng)有7人回復(fù)
-
化工專碩調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 英一數(shù)二272求調(diào)劑
已經(jīng)有44人回復(fù)
-
歡迎采礦、地質(zhì)、巖土、計算機、人工智能等專業(yè)的同學(xué)報考
已經(jīng)有10人回復(fù)
-
298求調(diào)劑一志愿中海洋
已經(jīng)有3人回復(fù)
-
環(huán)境調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)306調(diào)劑
已經(jīng)有4人回復(fù)
2樓2021-07-28 11:45:54
量子海洋
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 94.4
- 帖子: 66
- 在線: 108.7小時
- 蟲號: 4158109
- 注冊: 2015-10-19
- 專業(yè): 理論和計算化學(xué)
|
options文件沒有230行，是input.f90文件的第230行，但是我不知道是什么問題，以下是input.f90文件。 ! Copyright ! The University Montpellier 2, a scientific, cultural and professional ! public establishment whose registered office is at ! Place Eugène Bataillon - 34095 Montpellier CEDEX 05 - France. ! N° SIREN 193 410 883, APE Code 803Z ! represented by its President, Danièle HERIN. ! Authors : Dr. Sébastien Le Roux (sebastien.leroux@ipcms.unistra.fr) ! Pr. Philippe Jund (pjund@univ-montp2.fr) ! This software is a computer program whose purpose is to analyse ! the results of molecular dynamics simulations. ! This code is named R.I.N.G.S “Rigorous Investigation of Networks ! Generated using Simulations” and its main feature is the analysis of ! the connectivity using ring statistics. ! This software is governed by the CeCILL-B license under French law and ! abiding by the rules of distribution of free software. ! You can use, modify and/ or redistribute the software under the terms ! of the [CeCILL|CeCILL-B|CeCILL-C] license as circulated by CEA, CNRS ! and INRIA at the following URL: ! "https://www.cecill.info". ! As a counterpart to the access to the source code and rights to copy, ! modify and redistribute granted by the license, users are provided only ! with a limited warranty and the software's author, the holder of the ! economic rights, and the successive licensors have only limited ! liability. ! In this respect, the user's attention is drawn to the risks associated ! with loading, using, modifying and/or developing or reproducing the ! software by the user in light of its specific status of free software, ! that may mean that it is complicated to manipulate, and that also ! therefore means that it is reserved for developers and experienced ! professionals having in-depth computer knowledge. Users are therefore ! encouraged to load and test the software's suitability as regards their ! requirements in conditions enabling the security of their systems and/or ! data to be ensured and, more generally, to use and operate it in the ! same conditions as regards security. ! The fact that you are presently reading this means that you have had ! knowledge of the CeCILL-B license and that you accept its terms. SUBROUTINE READ_INPUT ! ! Read the user input file and the 'options' file ! USE PARAMETERS USE EXTRA LOGICAL :: CALC_R=.false. open (unit=9, file=NOM_FILE, action='read', status='old') read (9, *) read (9, *) read (9, *) read (9, *) SYSTLABEL if (DEBUGINP) write (6, *) SYSTLABEL read (9, *) NA if (DEBUGINP) write (6, *) NA read (9, *) NSP if (DEBUGINP) write (6, *) NSP !allocate (S_LENGTH(NSP), STAT=ERR) !if (ERR .ne. 0) then ! write (6, 001) t_p, "S_LENGTH" ! STOP !endif allocate (LABEL(NSP), STAT=ERR) if (ERR .ne. 0) then write (6, 001) t_p, "LABEL" STOP endif !allocate (M_SS(NSP), STAT=ERR) !if (ERR .ne. 0) then ! write (6, 001) t_p, "M_SS" !endif !allocate (R_DW(NSP), STAT=ERR) !if (ERR .ne. 0) then ! write (6, 001) t_p, "R_DW" !endif do i=1, NSP LABEL(i)="" enddo ! S_LENGTH(i)=0.0d0 ! M_SS(i)=0.0d0 ! R_DW(i)=0.0d0 read (9, *) LABEL!, S_LENGTH(i), M_SS(i), R_DW(i) if (DEBUGINP) then do i=1, NSP write (6, *) LABEL(i)!, S_LENGTH(i), M_SS(i), R_DW(i) enddo endif read (9, *) NS if (DEBUGINP) write (6, *) NS do i=1,3 do j=1, 3 lvect(i,j)=0.0d0 enddo enddo read (9, *) LATTYPE if (LATTYPE .eq. 1) then read (9, *) lvect(1,1), lvect(1,2), lvect(1,3) if (DEBUGINP) write (6, *) lvect(1,1), lvect(1,2), lvect(1,3) read (9, *) lvect(2,1), lvect(2,2), lvect(2,3) if (DEBUGINP) write (6, *) lvect(2,1), lvect(2,2), lvect(2,3) read (9, *) lvect(3,1), lvect(3,2), lvect(3,3) if (DEBUGINP) write (6, *) lvect(3,1), lvect(3,2), lvect(3,3) else read (9, *) modv(1), modv(2), modv(3) if (DEBUGINP) write (6, *) modv(1), modv(2), modv(3) ALPHA=0.0d0 BETA=0.0d0 GAMA=0.0d0 read (9, *) ALPHA, BETA, GAMA if (DEBUGINP) write (6, *) ALPHA, BETA, GAMA endif DELTA_T=0.0d0 read (9, *) DELTA_T if (DEBUGINP) write (6, *) DELTA_T read (9, *) DATA_FORM if (DEBUGINP) write (6, *) DATA_FORM if (DATA_FORM .eq. "MSF") then read (9, 106) MSF_ROOT if (DEBUGINP) write (6, 106) MSF_ROOT read (9, *) FIRST_INDEX if (DEBUGINP) write (6, *) FIRST_INDEX read (9, *) LAST_INDEX if (DEBUGINP) write (6, *) LAST_INDEX read (9, *) NSBG if (NSBG .eq. 0) NSBG=1 if (DEBUGINP) write (6, *) NSBG else read (9, *) NOM_FILE if (DEBUGINP) write (6, *) NOM_FILE if (DATA_FORM.eq."VAS" .or. DATA_FORM.eq."TRJ") then allocate (NBSPBS(NSP), STAT=ERR) if (ERR .ne. 0) then write (6, 001) t_p, "NBSPBS" STOP endif do i=1, NSP read (9, *) LABEL(i), NBSPBS(i) if (DEBUGINP) write (6, *) LABEL(i), NBSPBS(i) enddo else if (DATA_FORM.eq."PDB" .and. NS.gt.1) then write (6, '("Error: for PDB files the number of steps, NS, must be 1")') STOP endif endif read (9, *) NUMBER_OF_DELTA_R NUMBER_OF_DELTA_R=2*NUMBER_OF_DELTA_R if (DEBUGINP) write (6, *) NUMBER_OF_DELTA_R read (9, *) NUMBER_OF_Q_POINTS if (DEBUGINP) write (6, *) NUMBER_OF_Q_POINTS QMAX=0.0d0 read (9, *) QMAX if (DEBUGINP) write (6, *) QMAX SIGMA_LISS=0.0d0 read (9, *) SIGMA_LISS if (DEBUGINP) write (6, *) SIGMA_LISS read (9, *) NUMBER_OF_DELTA_ANG if (DEBUGINP) write (6, *) NUMBER_OF_DELTA_ANG read (9, *) NDV if (DEBUGINP) write (6, *) NDV read (9, *) TAILLE_REF TAILLD=2*TAILLE_REF if (TAILLD .lt. 3) then write (6, '(" Changing option => Max. ring size must be >= 3")') TAILLD=4 TAILLE_REF=2 endif if (DEBUGINP) write (6, *) TAILLE_REF read (9, *) TAILLC if (TAILLC .lt. 2) then write (6, '(" Changing option => Max. chain length must be >= 2")') TAILLC=2 endif if (DEBUGINP) write (6, *) TAILLC read (9, *) NL=NSP*(NSP+1)/2 allocate (GCUT(NSP,NSP), STAT=ERR) if (ERR .ne. 0) then write (6, 001) t_p, "GCUT" STOP endif do i=1, NL read (9, *) A1, A2, DC if (DEBUGINP) write (6, *) A1, " ", A2, DC j=1 do while (LABEL(j) .ne. A1) j=j+1 enddo l=1 do while (LABEL(l) .ne. A2) l=l+1 enddo GCUT(j,l)=0.0d0 GCUT(l,j)=0.0d0 GCUT(j,l)=DC GCUT(l,j)=DC enddo Gr_cutoff=0.0d0 read (9, *) ATOT, Gr_cutoff if (DEBUGINP) write (6, *) ATOT, Gr_cutoff close (9) ! Now read 'options' file open (unit=10, file="options", action="read", status='old') if(t_p .eq. 0) write (6, '(/,"========== Options ==========",/)') read (10, *) read (10, *) read (10, *) read (10, *) Nom1, PBC if (PBC) then if(t_p .eq. 0) write (6, '(" => Periodic Boundary Conditions applied")') else if(t_p .eq. 0) write (6, '(" => No Periodic Boundary Conditions applied")') endif read (10, *) Nom1, FRAC if (FRAC) then if(t_p .eq. 0) write (6, '(" => Atomic positions in fractional coordinates")') else if(t_p .eq. 0) write (6, '(" => Atomic positions in real coordinates")') endif read (10, *) Nom1, CALC_GR if(t_p .eq. 0) then if (CALC_GR) write (6, '(" => Calculate g(r) radial distribution functions")') endif read (10, *) Nom1, CALC_SQ if (CALC_SQ .and. .not.(CALC_GR)) then if(t_p .eq. 0) write (6, '(" Changing option => Calculate g(r) is needed to compute S(q)")') CALC_GR=.true. endif if(t_p .eq. 0) then if (CALC_SQ) write (6, '(" => Calculate S(q) structure factor from FFT[g(r)]")') endif read (10, *) Nom1, CALC_SK if(t_p .eq. 0) then if (CALC_SK) write (6, '(" => Calculate S(k) structure factor from Debye equation")') endif read (10, *) Nom1, CALC_GK if (CALC_GK .and. .not.(CALC_SK)) then if(t_p .eq. 0) write (6, '(" Changing option => Calculate S(k) Debye is needed to compute g(r) from FFT[S(k) Debye]")') CALC_SK=.true. endif if(t_p .eq. 0) then if (CALC_GK) write (6, '(" => Calculate g(r) from FFT[S(k) Debye]")') endif read (10, *) Nom1, CALC_MSD if(t_p .eq. 0) then if (CALC_MSD) write (6, '(" => Calculate Mean Square Displacement MSD")') endif read (10, *) Nom1, MSDEA if(t_p .eq. 0) then if (MSDEA) write (6, '(" => Calculate single atom MSDs")') endif read (10, *) Nom1, CALC_BONDS if(t_p .eq. 0) then if (CALC_BONDS) write (6, '(" => Calculate bond properties")') endif read (10, *) Nom1, CALC_ANG if(t_p .eq. 0) then if (CALC_ANG) write (6, '(" => Calculate angles distribution")') endif ! Chains statistics: read (10, *) Nom1, CALC_CHAINS read (10, *) read (10, *) Nom1, CTLT read (10, *) Nom1, AAAA read (10, *) Nom1, ACAC read (10, *) Nom1, ISOLATED if(t_p .eq. 0) then if (CALC_CHAINS) write (6, '(" => Calculate chains statistics")') if (CALC_CHAINS .and. CTLT .eq. 0)then write (6, '(3x," => Using all atoms as nodes")') elseif (CALC_CHAINS .and. (CTLT.le.NSP .and. CTLT.gt.0)) then write (6, '(3x," => Using species ",I2," as nodes")') CTLT elseif (CALC_CHAINS) then write (6, '(" Problem in ''options'' file => check chains statistics options")') STOP endif if (CALC_CHAINS .and. AAAA .and. ACAC) then write (6, '(" Problem in ''options'' file => check chains statistics options")') write (6, '(" Problem in ''options'' file => AAAA, ABAB or none, not both !")') STOP endif if (CALC_CHAINS .and. AAAA) write (6, '(3x," => Pure AAAA chains")') if (CALC_CHAINS .and. ACAC) write (6, '(3x," => Pure ABAB chains")') if (CALC_CHAINS .and. ISOLATED) write (6, '(3x," => Only 1-(2)n-1 coordinations chains")') if (CALC_CHAINS .and. .not.(CALC_BONDS)) then if(t_p .eq. 0) then write (6, '(" Changing option => Bond properties are needed to compute chains statistics")') endif CALC_BONDS=.true. endif endif read (10, *) ! Rings statistics: read (10, *) Nom1, CALC_R read (10, *) read (10, *) Nom1, LTLT read (10, *) Nom1, ABAB read (10, *) Nom1, CALC_R0 read (10, *) Nom1, CALC_R1 read (10, *) Nom1, CALC_R2 read (10, *) Nom1, CALC_R3 read (10, *) Nom1, CALC_R4 read (10, *) Nom1, CALC_PRINGS read (10, *) Nom1, CALC_STRINGS read (10, *) Nom1, BARYCROUT read (10, *) Nom1, RING_P1 read (10, *) Nom1, RING_P2 read (10, *) Nom1, RING_P3 read (10, *) Nom1, RING_P4 read (10, *) Nom1, RING_P5 if (CALC_R .and. (CALC_R0.or.CALC_R1.or.CALC_R2.or.CALC_R3.or.CALC_R4.or.CALC_PRINGS.or.CALC_STRINGS)) then CALC_RINGS=.true. if (CALC_R1.or.CALC_R2.or.CALC_R3.or.CALC_R4) DOAMPAT=.true. else CALC_RINGS=.false. ABAB=.false. endif if (RING_P1 .or. RING_P2 .or. RING_P3 .or. RING_P4 .or. RING_P5) then CALC_RPROP=.true. else CALC_RPROP=.false. endif if(t_p .eq. 0) then if (CALC_RINGS) write (6, '(" => Calculate rings statistics")') if (CALC_RINGS .and. LTLT .eq. 0)then write (6, '(3x," => Using all atoms as nodes")') elseif (CALC_RINGS .and. (LTLT.le.NSP .and. LTLT.gt.0)) then write (6, '(3x," => Using species ",I2," as nodes")') LTLT elseif (CALC_RINGS) then write (6, '(" Problem in ''options'' file => check rings statistics options")') STOP endif if (CALC_RINGS .and. ABAB) write (6, '(3x," => Pure ABAB rings")') if (CALC_RINGS .and. CALC_R0) write (6, '(3x," => Find all rings in the simulation box")') if (CALC_RINGS .and. CALC_R1) write (6, '(3x," => King''s rings statistics")') if (CALC_RINGS .and. CALC_R2) write (6, '(3x," => Guttman''s rings statistics")') if (CALC_RINGS .and. CALC_R3) write (6, '(3x," => King''s + homopolar bonds rings statistics")') if (CALC_RINGS .and. CALC_R4) write (6, '(3x," => Guttman''s + homopolar bonds rings statistics")') if (CALC_RINGS .and. CALC_PRINGS) write (6, '(3x," => Primitive rings statistics")') if (CALC_RINGS .and. CALC_STRINGS) write (6, '(3x," => Strong rings statistics")') if (CALC_RINGS .and. BARYCROUT) write (6, '(3x," => Output OpenDX rings barycenters positions")') if (CALC_RINGS .and. CALC_RPROP) write (6, '(3x," => Compute detailed rings properties")') if (CALC_RINGS .and. .not.CALC_BONDS) then write (6, '(" Changing option => Bond properties are needed to compute rings statistics")') CALC_BONDS=.true. endif endif read (10, *) read (10, *) Nom1, CALC_VAC if(t_p .eq. 0) then if (CALC_VAC) write (6, '(" => Calculate voids distribution")') endif if (CALC_VAC .and. .not.(CALC_BONDS)) then if(t_p .eq. 0) then write (6, '(" Changing option => Bond properties are needed to evaluate void distribution")') endif CALC_BONDS=.true. endif read (10, *) read (10, *) read (10, *) read (10, *) Nom1, EVOLOUT if(t_p .eq. 0) then if (EVOLOUT) write (6, '(" => Ouptut calculations results for each MD steps")') endif read (10, *) Nom1, DXOUT if(t_p .eq. 0) then if (DXOUT) then write (6, '(" => Ouptut OpenDX .dx file(s)")') write (6, '(" Main file = ''SYSTEM.dx''")') endif endif read (10, *) read (10, *) Nom1, RADOUT if (DXOUT .and. RADOUT) then if(t_p .eq. 0) write (6, '(3x," => Using atomic radius as color map")') IDOUT=.false. endif read (10, *) Nom1, RNGOUT if (DXOUT .and. RNGOUT) then if (LTLT .ne. 0) then if(t_p .eq. 0) then write (6, '(3x," => Atomic ring factor can''t be used as color map")') write (6, '(3x," if rings search options are not changed")') write (6, '(3x," Uing atomic mass as color map instead")') endif else if(t_p .eq. 0) write (6, '(3x," => Using atomic ring factor as color map")') IDOUT=.false. FACTATPNA=.true. endif endif read (10, *) Nom1, DRNGOUT if (DXOUT .and. DRNGOUT) then if (LTLT .ne. 0) then if(t_p .eq. 0) then write (6, '(3x," => Delta ring factor can''t be used as color map")') write (6, '(3x," if rings search options are not changed")') write (6, '(3x," Uing atomic mass as color map instead")') endif else if(t_p .eq. 0) write (6, '(3x," => Using Delta ring factor as color map")') IDOUT=.false. FACTATPNA=.true. endif endif if (DXOUT .and. IDOUT) then if(t_p .eq. 0) write (6, '(3x," => Using atomic mass as color map")') endif read (10, *) Nom1, VACOUT if (DXOUT .and. VACOUT) then if(t_p .eq. 0) write (6, '(3x," => Voids visualization using OpenDX")') if (.not.(CALC_VAC)) then if(t_p .eq. 0) write (6, '(" Changing option => Evaluate voids distribution is needed to visualize voids")') CALC_VAC=.true. if (.not.(CALC_BONDS)) then if(t_p .eq. 0) then write (6, '(" Changing option => Bond properties are needed to evaluate voids distribution")') endif CALC_BONDS=.true. endif endif endif read (10, *) Nom1, DXTETRA if (DXOUT .and. DXTETRA) then if(t_p .eq. 0) write (6, '(3x," => Tetrahedra visualization using OpenDX")') if (.not.(CALC_BONDS)) then if(t_p .eq. 0) then write (6, '(" Changing option => Bond properties are needed to visualize tetrahedra")') endif CALC_BONDS=.true. endif endif read (10, *) Nom1, PATHDX if (PATHDX) then if(t_p .eq. 0) write (6, '(3x," => Trajectories of the ",I4," particles of highest MSD")') NVALT if (.not.CALC_MSD) then if (t_p .eq. 0) write (6, '(" Changing option => MSD has to be computed to output trajectories")') CALC_MSD=.true. endif if (.not.MSDEA) then if (t_p .eq. 0) write (6, '(" Changing option => MSD for single atoms have to be computed to output trajectories")') MSDEA=.true. endif endif read (10, *) read (10, *) Nom1, NOM_OUT if(t_p .eq. 0) write (6, '(/,"======== End Options ========",//)') close (10) !101 FORMAT (I4) !102 FORMAT (I2) !103 FORMAT (I5) !104 FORMAT (f15.10) !105 FORMAT (A3) 106 FORMAT (A13) !107 FORMAT (I10) !109 FORMAT (A17) !111 FORMAT (A2,2x,f9.6) !113 FORMAT (A15) !114 FORMAT (f7.5) 001 FORMAT (/," Allocation Problem -> No space left ?! :: node ",i3," READ_INPUT allocation of ",A15) END SUBROUTINE READ_INPUT |
4樓2021-07-29 20:11:26
量子海洋
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 94.4
- 帖子: 66
- 在線: 108.7小時
- 蟲號: 4158109
- 注冊: 2015-10-19
- 專業(yè): 理論和計算化學(xué)
6樓2021-07-30 10:10:52
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 量子海洋 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +4 | 1一11 2026-03-03 | 4/200 |
|
|
[基金申請] 請問大家，研究風(fēng)險與應(yīng)對措施那里， 大家都怎么寫呢 ？ +3 | cauasen 2026-03-02 | 3/150 |
|
|
|
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[考研] 085600 材料與化工 298 +7 | 小西笑嘻嘻 2026-03-03 | 7/350 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個 2026-02-28 | 11/550 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 調(diào)劑材料學(xué)碩 +4 | 詞凝Y 2026-03-02 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|