| 查看: 645 | 回复: 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 |
» 猜你喜欢
材料考研求调剂总分280
已经有3人回复
求调剂
已经有4人回复
材料调剂
已经有10人回复
求调剂到材料
已经有3人回复
化工求调剂!
已经有9人回复
301求调剂
已经有13人回复
308求调剂
已经有13人回复
材料工程310专硕调剂
已经有15人回复
304求调剂
已经有9人回复
081700,311,求调剂
已经有16人回复














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