| 查看: 613 | 回复: 0 | |||
[求助]
想固定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 |
» 猜你喜欢
求助:我三月中下旬出站,青基依托单位怎么办?
已经有6人回复
北京211副教授,35岁,想重新出发,去国外做博后,怎么样?
已经有8人回复
Cas 72-43-5需要30g,定制合成,能接单的留言
已经有7人回复
磺酰氟产物,毕不了业了!
已经有5人回复
论文终于录用啦!满足毕业条件了
已经有25人回复
2026年机械制造与材料应用国际会议 (ICMMMA 2026)
已经有3人回复
自荐读博
已经有3人回复
不自信的我
已经有5人回复
投稿Elsevier的杂志(返修),总是在选择OA和subscription界面被踢皮球
已经有8人回复









(ICELL(I),I=1,3)
回复此楼