24小时热门版块排行榜    

查看: 3221  |  回复: 3

小莫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
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 17:25 by xxniao123
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 面上合作单位盖章 +3 ssyjh 2026-08-27 3/150 2026-08-27 09:59 by 阿拉贡aragon
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见