应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 请问使用EPW的虫友，在计算epw3.in时参数错误，请帮忙解答
41/1返回列表
查看: 3021  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小莫2016

新虫 (小有名气)

[求助] 请问使用EPW的虫友，在计算epw3.in时参数错误，请帮忙解答

错误报告说我试图分配已经分配的变量“tau”
错误报告：
At line 231 of file ephwann_shuffle.f90
Fortran runtime error: Attempting to allocate already allocated variable 'tau'

Error termination. Backtrace:
#0  0x7f8de11f9ed7 in ???
#1  0x7f8de11faa75 in ???
#2  0x7f8de11fae47 in ???
#3  0x557f17017a58 in ephwann_shuffle_
        at /home/sunxiaomo/下载/qe6.1/qe-6.1/EPW/src/ephwann_shuffle.f90:231
#4  0x557f17005dae in elphon_shuffle_wrap_
        at /home/sunxiaomo/下载/qe6.1/qe-6.1/EPW/src/elphon_shuffle_wrap.f90:766
#5  0x557f16fdea64 in epw
        at /home/sunxiaomo/下载/qe6.1/qe-6.1/EPW/src/epw.f90:153
#6  0x557f16fde41e in main
        at /home/sunxiaomo/下载/qe6.1/qe-6.1/EPW/src/epw.f90:23

问题在于当设置  epwread     =  .true.时就会报错，设置成  epwread     =  .false.时就没问题。这是什么原因？

输入文件：
--
&inputepw
  prefix      = 'Al'
  amass(1)    = 26.98
  outdir      ='./'

  iverbosity  = 0

  elph        = .true.
  epbwrite    = .false.
  epbread     = .true.

  epwwrite    =  .false.
  epwread     =  .true.

  efermi_read = .true.
  fermi_energy = 7.934874

  nbndsub     =  6
  nbndskip    =  0

  wannierize  = .false.
  num_iter    = 300
  iprint      = 2
  dis_win_max = 15
  dis_froz_max= 9
  proj(1)='random'
  ! proj(1)     = 'f=0,0,0:l=-3'   

  !elinterp    = .true.
  !phinterp    = .true.

  !tshuffle2   = .true.
  !tphases     = .false.

  elecselfen  = .true.
  phonselfen  = .false.
  a2f         = .false.

  specfun      = .true.
  wmin_specfun = -13
  wmax_specfun = 12
  nw_specfun   = 500

  parallel_k  = .true.
  parallel_q  = .false.

  fsthick     = 1.36056981 ! eV
  eptemp      = 300 ! K (same as PRB 76, 165108)
  degaussw    = 0.03 ! eV

  dvscf_dir   = '../phonons/save'
  filukk      = './Al.ukk'
!filqf       = 'meshes/path.dat'
  
  nqf1        = 18
  nqf2        = 18
  nqf3        = 18

  nkf1        = 50
  nkf2        = 50
  nkf3        = 50
  
  nk1         = 6
  nk2         = 6
  nk3         = 6

  nq1         = 6
  nq2         = 6
  nq3         = 6
/
      16 cartesian
   0.0000000   0.0000000   0.0000000  0.0092593
  -0.1666667   0.1666667  -0.1666667  0.0740741
  -0.3333333   0.3333333  -0.3333333  0.0740741
   0.5000000  -0.5000000   0.5000000  0.0370370
   0.0000000   0.3333333   0.0000000  0.0555556
  -0.1666667   0.5000000  -0.1666667  0.2222222
   0.6666667  -0.3333333   0.6666667  0.2222222
   0.5000000  -0.1666667   0.5000000  0.2222222
   0.3333333   0.0000000   0.3333333  0.1111111
   0.0000000   0.6666667   0.0000000  0.0555556
   0.8333333  -0.1666667   0.8333333  0.2222222
   0.6666667   0.0000000   0.6666667  0.1111111
   0.0000000  -1.0000000   0.0000000  0.0277778
   0.6666667  -0.3333333   1.0000000  0.2222222
   0.5000000  -0.1666667   0.8333333  0.2222222
  -0.3333333  -1.0000000   0.0000000  0.1111111
回复此楼

» 猜你喜欢
1楼 2018-01-15 22:58:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小莫2016

新虫 (小有名气)
这个怎么调？
回复此楼
2楼2018-01-15 23:49:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小莫2016

新虫 (小有名气)
计算结果差一个数量级是什么原因？
赞一下 回复此楼
3楼2018-01-16 12:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

year@

铜虫 (初入文坛)
引用回帖:
3楼: Originally posted by 小莫2016 at 2018-01-16 12:38:25
计算结果差一个数量级是什么原因？

请问下您  fsthick具体怎么设定的呢
赞一下 回复此楼
4楼2018-10-17 10:28:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小莫2016 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭