<div dir="ltr">Hello! I am facing a problem when i use starting_charge command in 5.1 version. <div>The problem is:</div><div>  Â  Â Error in routine Â read_namelists (19):<br>  Â  Â  reading namelist system<br></div><div><br></div><div>I am no sure starting_charge is implemented in this old version... I search for an old pw.x input information, but i cannot found it.</div><div><br></div><div>Someone knows if is implemented and how can i use it?</div><div><br></div><div>My input is:</div><div>&CONTROL<br>  Â  Â  Â  Â  Â  Â  Â  Â calculation = 'scf',<br>  Â  Â  Â  Â  Â  Â  Â  restart_mode = 'from_scratch',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  prefix = 'cu2o',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  outdir = './temp/',<br></div><div>  Â  Â  Â  Â  Â  Â  Â  Â  pseudo_dir = '.',<br>  Â  Â  Â  Â  Â  Â  Â etot_conv_thr = 1.0D-6, Â  Â  Â  Â  Â  Â  Â  Â  Â !1.0D-4 default<br>  Â  Â  Â  Â  Â  Â  Â forc_conv_thr = 1.0D-4, Â  Â  Â  Â  Â  Â  Â  Â  Â !1.0D-3 default<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â tstress = .true.,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â tprnfor = .true.,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nstep = 400,<br>/<br>&SYSTEM<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ibrav = 0,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â celldm(1) = 7.975234267 ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nat = 47 ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ntyp = 3,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ecutwfc = 50 ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ecutrho = 300 Â ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â input_dft = 'pbesol',<br>  Â  Â  Â  Â  Â  Â  Â  Â occupations = 'smearing',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  smearing = 'mp',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â degauss = 0.001,<br>  Â  Â  Â  Â  Â  Â  Â  Â  tot_charge = 2.0<br>  Â  Â  Â  Â  Â  Â assume_isolated = 'esm'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  esm_bc = 'bc1'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â esm_w = 0.0<br>  Â  Â  Â  Â  Â  Â  Â  Â  esm_efield = 0.0<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  esm_nfit = 4<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nosym = .true.<br>  Â  Â  Â  Â  starting_charge(1) = 0.5<br>/<br> &ELECTRONS<br>  Â  Â  Â  Â  Â  Â diagonalization = 'david'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  conv_thr = 1.0e-8<br>  Â  Â  Â  Â  Â  Â  Â  Â mixing_beta = 0.3<br>/<br> &IONS<br>  Â  Â  Â  Â  Â  Â  Â  ion_dynamics = 'bfgs'<br> /<br>&CELL<br>  Â  Â  Â  Â  Â  Â  Â cell_dynamics = 'bfgs'<br>  Â  Â  Â  Â  Â  Â  Â  Â cell_dofree = 'all',<br>/<br><br>ATOMIC_SPECIES<br>  Cu1 63.546 Â Cu.pbe-n-van_ak.UPF<br>  Cu 63.546 Â  Cu.pbe-n-van_ak.UPF<br>  O Â 15.9994 Â O.pbesol-rrkjus.UPF<br></div><div><br></div><div>plus kpoints and atomic positions.</div><div><br></div><div>Regards</div></div>