| 查看: 699 | 回复: 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 |
» 猜你喜欢
轴手性化合物的确定
已经有5人回复
MSER送审了还被拒稿
已经有6人回复
内心匮乏
已经有7人回复
明天E口面上会评
已经有8人回复
不要再数国自然申请书的 filecode 的分隔符个数了
已经有17人回复
祈祷青基必中
已经有16人回复
生命口会评
已经有12人回复
27届辽宁大学应届毕业生申博
已经有3人回复
关于如何从代码看上不上会
已经有23人回复
chemdraw
已经有6人回复











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