<div dir="ltr">Dear Giuseppe,<div>Thanks! But it turned out that my problem had a simpler explanation: missing end-of-line character after the final /.</div><div>Best,</div><div>Michal</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 24 May 2019 at 10:02, Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
dear Michal<br>
In the past I sometimes faced with the same problem. I've found (very  <br>
empirically...) that if you break the calculation in two steps the  <br>
problem can be solved. Something like this<br>
<br>
export FILEA="yourprefix"<br>
<br>
export INPFILE=$FILEA-pol.inp<br>
export OUTFILE=$FILEA-pol.out<br>
echo " $FILEA"<br>
echo " $INPFILE"<br>
echo " $OUTFILE"<br>
cat > $INPFILE << EOF<br>
  &inputpp<br>
     prefix  = '$FILEA'<br>
     outdir = '$TMP_DIR/'<br>
     filplot = '$FILEA-pol.dat'<br>
     plot_num= 6<br>
  /<br>
EOF<br>
$PARA_PREFIX $ESPRESSO/pp.x $PARA_POSTFIX < $INPFILE >> $OUTFILE<br>
<br>
export FILE="rho-pol-$FILEA"<br>
export INPFILE=$FILE.inp<br>
export OUTFILE=$FILE.out<br>
echo " $FILE"<br>
echo " $INPFILE"<br>
echo " $OUTFILE"<br>
cat > $INPFILE << EOF<br>
  &inputpp<br>
  /<br>
  &plot<br>
     nfile = 1<br>
     filepp(1) = '$FILEA-pol.dat'<br>
     weight(1) = 1.0<br>
     iflag = 3<br>
     output_format = 6<br>
     fileout = '$FILE.cube'<br>
     e1(1) =1.0, e1(2)=1.0, e1(3) = 0.0,<br>
     e2(1) =0.0, e2(2)=0.0, e2(3) = 1.0,<br>
     nx=56, ny=40<br>
  /<br>
EOF<br>
$PARA_PREFIX $ESPRESSO/pp.x $PARA_POSTFIX < $INPFILE >> $OUTFILE<br>
<br>
You may give it a try<br>
HTH<br>
Giuseppe<br>
<br>
Quoting Michal Krompiec <<a href="mailto:michal.krompiec@gmail.com" target="_blank">michal.krompiec@gmail.com</a>>:<br>
<br>
> Hello,<br>
> I'm struggling with plotting of spin density with pp.x, from a PAW SCF<br>
> calculation. This is my input to pp.x:<br>
> &INPUTPP<br>
>  prefix='Si_slab_mo_40',<br>
>  plot_num=6,<br>
> /<br>
> &PLOT<br>
>  iflag = 3 ,<br>
>  output_format = 6,<br>
>  fileout = 'spin.cube'<br>
>  nx=64<br>
>  ny=64<br>
>  nz=64<br>
> /<br>
> And this is the error I'm getting:<br>
>      Calling punch_plot, plot_num =   6<br>
>      Writing data to file  tmp.pp<br>
>      Message from routine chdens:<br>
>      namelist plot not found or invalid, exiting<br>
><br>
> I'd be very grateful for any advice.<br>
><br>
> Best regards,<br>
> Michal Krompiec<br>
> Merck KGaA and University of Southampton<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div>