| 查看: 637 | 回复: 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 |
» 猜你喜欢
327求调剂
已经有8人回复
26考研一志愿中国石油大学(华东)305分求调剂
已经有3人回复
伙伴们,祝我生日快乐吧
已经有23人回复
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有13人回复
0856专硕279求调剂
已经有5人回复
本人考085602 化学工程 专硕
已经有7人回复
22408总分284求调剂
已经有3人回复
070305求调剂
已经有4人回复
材料专硕326求调剂
已经有4人回复
294求调剂
已经有4人回复













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