24小时热门版块排行榜

返回列表

moyi22334

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 849
帖子: 25
在线: 5.3小时
虫号: 12930850
注册: 2018-11-23

[求助] 想固定x轴作优化，重新编译了VASP，为什么没起作用？

按照网上的方法，修改了constr_cell_relax.F文件，并重新编译了VASP，并没有起作用，可能的原因是什么？
代码如下：
   SUBROUTINE CONSTR_CELL_RELAX(FCELL)
   USE PREC
   REAL(Q) FCELL(3,3), SAVE(3)
   LOGICAL FILFLG
   INTEGER ICELL(3)
   INQUIRE(FILE='OPTCELL',EXIST=FILFLG)
   IF (FILFLG) THEN
      OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD')
      READ(67,"(3I1)"

(ICELL(I),I=1,3)
      CLOSE(67)
      DO I=1,3
      SAVE(I)=FCELL(I,I)
      ENDDO
      FCELL=0.0d0
      DO I=1,3
      IF (ICELL(I)==1) FCELL(I,I)=SAVE(I)
      ENDDO
   ENDIF
!    just one simple example
!    relaxation in x directions only
!    SAVE=FCELL(1,1)
!    FCELL=0 ! F90 style: set the whole array to zero
!    FCELL(1,1)=SAVE
   RETURN
   END SUBROUTINE
方法来自这里：http://blog.sciencenet.cn/blog-671981-731203.html

回复此楼

» 猜你喜欢

1楼 2019-12-03 16:17:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 moyi22334 的主题更新

返回列表