| 查看: 761 | 回复: 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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
