| 查看: 3015 | 回复: 1 | ||
andy101199银虫 (正式写手)
|
[求助]
QE计算声电耦合最后lambda.x出现问题
|
|
QE计算声电耦合在进行最后一步ambda.x出现问题，之前步骤都正常结束。 lambda.in如下（ Program PWSCF v.5.4.0）： 75 0.02 1 ! emax (something more than highest phonon mode in THz), degauss, smearing method 8 ! Number of q-points 0.000000000 0.000000000 0.000000000 1.00 0.000000000 0.000000000 -0.136976104 1.00 0.000000000 0.288675135 0.000000000 6.00 0.000000000 0.288675135 -0.136976104 6.00 0.000000000 -0.577350269 0.000000000 3.00 0.000000000 -0.577350269 -0.136976104 3.00 0.250000000 0.433012702 0.000000000 6.00 0.250000000 0.433012702 -0.136976104 6.00 elph_dir/elph.inp_lambda.1 elph_dir/elph.inp_lambda.2 elph_dir/elph.inp_lambda.3 elph_dir/elph.inp_lambda.4 elph_dir/elph.inp_lambda.5 elph_dir/elph.inp_lambda.6 elph_dir/elph.inp_lambda.7 elph_dir/elph.inp_lambda.8 0.10 错误提示 At line 60 of file lambda.f90 (unit = 5, file = 'stdin') Fortran runtime error: End of file 查看了lambda.f90中第60行附近代码： 58 real*8 mustar, omegalog(20), Tc, x 60 read(5,*) emax, degaussq, ngaussq 61 deltae=emax/(nex-1) 62 read(5,*) nks 63 if (nks.gt.npk) call errore('lambda','too many q-points',npk) 64 sum=0.d0 65 do ik=1,nks 66 read(5,*) q(1,ik), q(2,ik), q(3,ik), wk(ik) 67 sum = sum + wk(ik) 68 end do 69 do ik=1,nks 70 wk(ik)=wk(ik)/sum 71 end do 且能输出lambda.out文件，但是感觉明显是错误的 lambda omega_log T_c -0.15572 2477.084 73200.645 -0.14600 2490.699 88108.999 -0.14491 2488.983 89922.337 -0.14131 2486.484 96426.060 -0.13638 2486.396 ********* -0.13187 2491.111 ********* -0.12821 2495.033 ********* -0.12502 2501.643 ********* -0.12261 2504.104 ********* -0.12064 2512.136 ********* -0.11938 2516.256 ********* -0.11828 2520.488 ********* -0.11777 2523.697 ********* -0.11731 2528.947 ********* -0.11708 2531.623 ********* -0.11712 2534.020 ********* -0.11707 2536.620 ********* -0.11729 2538.828 ********* -0.11730 2541.920 ********* -0.11756 2543.034 ********* -0.11773 2545.929 ********* -0.11785 2545.691 ********* -0.11799 2546.964 ********* -0.11818 2550.281 ********* -0.11844 2550.324 ********* -0.11855 2551.999 ********* -0.11878 2551.735 ********* -0.11905 2553.195 ********* -0.11936 2554.038 ********* -0.11976 2556.991 ********* 请教各位大神，这个问题如何解决？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有3人回复
-
基金未中，这种答复是模板吗？
已经有4人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有5人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
-
面上合作单位盖章
已经有3人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有5人回复
-
基金不中，共勉
已经有7人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
2楼2018-12-19 22:32:34
