| 查看: 3003 | 回复: 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 ********* 请教各位大神，这个问题如何解决？ |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有62人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有29人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
2楼2018-12-19 22:32:34
