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

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (6041)
>考研 (1553)
>導(dǎo)師招生 (1385)
>文獻(xiàn)求助 (242)
>蟲友互識 (189)
>碩博家園 (164)
>休閑灌水 (162)
>考博 (130)
>博后之家 (80)
>論文投稿 (68)
>公派出國 (65)
>招聘信息布告欄 (55)
>基金申請 (46)
>教師之家 (38)
>有機資源 (25)
>找工作 (25)

返回列表

【懸賞金幣】回答本帖問題，作者量子海洋將贈送您 50 個金幣

當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

量子海洋

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 94.4
帖子: 66
在線: 108.7小時
蟲號: 4158109
注冊: 2015-10-19
專業(yè): 理論和計算化學(xué)

[求助] 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ù)此樓

» 猜你喜歡

一志愿中科大能動297求調(diào)劑，本科川大已經(jīng)有3人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有17人回復(fù)
化工專碩348，一志愿985求調(diào)劑已經(jīng)有10人回復(fù)
一志愿東北大學(xué)化學(xué)314分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿天津大學(xué)085600 319分材料與化工金屬方向求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
環(huán)境工程學(xué)碩288求助調(diào)劑已經(jīng)有5人回復(fù)
0856材料工程，初試313調(diào)劑已經(jīng)有4人回復(fù)
085600 英一數(shù)二272求調(diào)劑已經(jīng)有37人回復(fù)

1樓 2021-07-28 11:36:37

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

get-it

新蟲 (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 1825.8
散金: 26
紅花: 4
帖子: 273
在線: 448.3小時
蟲號: 3212202
注冊: 2014-05-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

抱歉了，我也不太熟。
https://sourceforge.net/p/rings- ... il.com/#msg33203770

我看到rings官網(wǎng)有人問了同樣的問題，不過似乎沒有回答。

贊一下

回復(fù)此樓

5樓2021-07-29 21:10:50

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

查看全部 10 個回答

量子海洋

新蟲 (小有名氣)

發(fā)自小木蟲IOS客戶端

2樓2021-07-28 11:45:54

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

get-it

新蟲 (小有名氣)

應(yīng)助: 23 (小學(xué)生)
金幣: 1825.8
散金: 26
紅花: 4
帖子: 273
在線: 448.3小時
蟲號: 3212202
注冊: 2014-05-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

At line 230 of file input.f90 (unit = 10, file = 'options')
它叫你檢查options 230行

贊一下

回復(fù)此樓

3樓2021-07-28 16:43:58

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

量子海洋

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 94.4
帖子: 66
在線: 108.7小時
蟲號: 4158109
注冊: 2015-10-19
專業(yè): 理論和計算化學(xué)

引用回帖:

3樓: Originally posted by get-it at 2021-07-28 16:43:58
At line 230 of file input.f90 (unit = 10, file = 'options')
它叫你檢查options 230行

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

贊一下

回復(fù)此樓

4樓2021-07-29 20:11:26

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

查看全部 10 個回答

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

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

[考研] 085600 英一數(shù)二272求調(diào)劑 5+6	vida_a 2026-03-01	37/1850	2026-03-03 14:33 by vida_a
[考研] 085600材料與化工調(diào)劑 280分 +9	yyqqhh 2026-03-03	9/450	2026-03-03 13:53 by NUAAZXWS
[考研] 接收調(diào)劑 +12	津萌津萌 2026-03-02	20/1000	2026-03-03 13:29 by 考研考研考研考?/a>
[考研] 080500材料科學(xué)與工程 +4	202114020319 2026-03-03	4/200	2026-03-03 11:55 by ms629
[考研] 289求調(diào)劑 +5	BrightLL 2026-03-02	7/350	2026-03-03 11:22 by BrightLL
[考研] 材料325求調(diào)劑 30+3	mariusuki 2026-03-02	7/350	2026-03-03 11:14 by EBSD
[基金申請] 成果系統(tǒng)訪問量大，請15分鐘后再嘗試。由此給您造成的不便，敬請諒解。 +8	xhuama 2026-03-02	9/450	2026-03-03 08:13 by kudofaye
[考研] 271求調(diào)劑 +4	Ricardo1113 2026-03-02	4/200	2026-03-03 08:00 by 無際的草原
[考研] 化工京區(qū)271求調(diào)劑 +7	11ing 2026-03-02	7/350	2026-03-03 07:30 by 利好利好.
[考研] 298求調(diào)劑 +7	axyz3 2026-02-28	8/400	2026-03-03 06:11 by tgxtgxtgx9
[考研] 288求調(diào)劑 +3	少71.8 2026-03-02	5/250	2026-03-03 06:01 by tgxtgxtgx9
[考研] 290分材料工程085601求調(diào)劑數(shù)二英一 +8	llx0610 2026-03-02	9/450	2026-03-02 22:09 by 無際的草原
[考研] 085600求調(diào)劑 +4	LRZZZZZZ 2026-03-02	5/250	2026-03-02 20:33 by 云旗yunqi
[考研] 285求調(diào)劑 +9	滿頭大汗的學(xué)生 2026-02-28	9/450	2026-03-02 20:29 by hypershenger
[考研] 275求調(diào)劑 +7	明遠(yuǎn)求學(xué) 2026-03-01	7/350	2026-03-02 19:22 by zhukairuo
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5	半截的詩0927 2026-03-02	5/250	2026-03-02 18:37 by 明亮9527
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6	禮堂丁真258 2026-02-28	9/450	2026-03-02 12:04 by 52hz~~
[考研] 材料調(diào)劑 +6	愛擦汗的可樂冰 2026-02-28	7/350	2026-03-02 10:42 by Jy?
[考研] 化工299分求調(diào)劑一志愿985落榜 +5	嘻嘻(^ω^) 2026-03-01	5/250	2026-03-01 19:47 by 無際的草原
[考研] 328求調(diào)劑 +3	aaadim 2026-03-01	5/250	2026-03-01 17:29 by njzyff

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

24小時熱門版塊排行榜

量子海洋

[求助] R.I.N.G.S. code運行報錯

» 猜你喜歡

get-it

量子海洋

get-it

量子海洋

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