24小时热门版块排行榜    

CyRhmU.jpeg
查看: 604  |  回复: 0

moyi22334

新虫 (初入文坛)

[求助] 想固定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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 moyi22334 的主题更新
信息提示
请填处理意见