<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-2022-jp">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"MS ゴシック";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:游ゴシック;
        panose-1:2 11 4 0 0 0 0 0 0 0;}
@font-face
        {font-family:"MS Pゴシック";
        panose-1:2 11 6 0 7 2 5 8 2 4;}
@font-face
        {font-family:"\@游ゴシック";
        panose-1:2 11 4 0 0 0 0 0 0 0;}
@font-face
        {font-family:"\@MS Pゴシック";}
@font-face
        {font-family:"\@MS ゴシック";
        panose-1:2 11 6 9 7 2 5 8 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0mm;
        text-align:justify;
        text-justify:inter-ideograph;
        font-size:10.5pt;
        font-family:游ゴシック;}
span.18
        {mso-style-type:personal-compose;
        font-family:游ゴシック;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;
        font-family:游ゴシック;}
/* Page Definitions */
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:99.25pt 30.0mm 30.0mm 30.0mm;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026">
<v:textbox inset="5.85pt,.7pt,5.85pt,.7pt" />
</o:shapedefaults></xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="JA" link="#0563C1" vlink="#954F72" style="word-wrap:break-word;text-justify-trim:punctuation">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Dear all,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">I have been trying to run a scf calculation with PWSCF v.6.6,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">using the use_wannier = .ture. and U_projection_type =
</span><span style="font-size:11.0pt">‘<span lang="EN-US">ortho-atomic</span>’<span lang="EN-US"><o:p></o:p></span></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">in order to get satwfc files for pw4gww.x.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x crashes with an error complaining that an array is not being allocated:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">___________________________________________________________________________________<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">forrtl: severe (151): allocatable array is already allocated<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Image              PC                Routine            Line        Source            
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               000000000366B39B  Unknown               Unknown  Unknown<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               0000000000B6CC7E  orthouwfc_.a               78  orthoatwfc.f90<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               000000000143D50D  wfcinit_.a                 52  wfcinit.f90<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               00000000008F7375  init_run_.a               119  init_run.f90<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               0000000000F5FFC1  run_pwscf_                128  run_pwscf.f90<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               000000000041629B  MAIN__                    106  pwscf.f90<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               00000000004160E2  Unknown               Unknown  Unknown<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">libc-2.17.so       00002AE187D863D5  __libc_start_main     Unknown  Unknown<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">pw.x               0000000000415FE9  Unknown               Unknown  Unknown<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">___________________________________________________________________________________<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The input file is as follows:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">______________________________________________________<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&control<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  prefix        = './',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  outdir        = './',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  verbosity     = 'high',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  calculation   = 'scf',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  restart_mode  = 'from_scratch',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  wf_collect    = .True.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  pseudo_dir    = '/home/local/pseudo_dojo/nc-sr-04'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  use_wannier    = .true.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&system<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  ibrav         = 5,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  celldm(1)     = 10.32799684<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  celldm(4)     = 0.55996291<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  nat           = 10,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  nspin         = 2,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  tot_magnetization = 0.0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  starting_magnetization(1) = 0.6,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  starting_magnetization(2) = -0.6,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  starting_magnetization(3) = 0.0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  lda_plus_u    = .true.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  Hubbard_U(1)  = 5.3,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  Hubbard_U(2)  = 5.3,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  Hubbard_U(3)  = 0.0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  U_projection_type  = 'ortho-atomic',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  ntyp          = 3,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  input_dft          = 'pbe',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  ecutwfc       = 25.0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  ecutrho       = 255.0,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  nosym         = .true.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  noinv         = .true.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  tot_charge    = 0.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&electrons<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  conv_thr            = 1d-6,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  mixing_beta         = 0.7,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  mixing_ndim         = 8,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  mixing_mode         = 'plain',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  electron_maxstep    = 100,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  diagonalization     = 'david',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  diago_david_ndim    = 10,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  startingwfc            = 'random',<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&wannier_ac<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  plot_wannier = .false.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  use_energy_int = .false.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  print_wannier_coeff = .true.,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  nwan = 50<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_SPECIES<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  FE1 55.8452 Fe.upf<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  FE2 55.8452 Fe.upf<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">  O 15.9994 O.upf<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_POSITIONS (crystal)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">FE1      0.144585600   0.144585693   0.144584311<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.749973638   0.053142820   0.446718330<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">FE2      0.855580129   0.855582133   0.855581453<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.249973175   0.946718596   0.553142099<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.446718682   0.749973096   0.053142783<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.553141069   0.249973464   0.946719545<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.053143257   0.446718378   0.749972513<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O        0.946718496   0.553142032   0.249973540<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">FE2      0.355581790   0.355580001   0.355580136<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">FE1      0.644584162   0.644583784   0.644585288<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">K_POINTS automatic<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">4 4 4 1 1 1<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">______________________________________________________<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Best Regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">--<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Yuki Ishikawa<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Business Development Division<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">TODA KOGYO CORP.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
</body>
</html>