[Pw_forum] i/o error in davcio in the phonon calculation

[Jian ZHOU]

Dear all,

I have compiled the pwscf v2.1.4 successfully and also have run the
example02 without any error.  However, when I tried to do another
phonon calculation of Gamma point, it stops during running the ph.x.

I firstly run a scf calculation from scratch and then do a gamma point
phonon calculation, reading the previous output files(The two outdirs
are the same in the two calculation). The following is the screen
output message of ph.x.

----------------------------------------------------------------------------
...
     Atomic displacements:
     There are  17 irreducible representations

     Representation     1      3 modes - To be done

     Representation     2      3 modes - To be done

     Representation     3      2 modes - To be done

     Representation     4      3 modes - To be done

     Representation     5      3 modes - To be done

     Representation     6      3 modes - To be done

     Representation     7      3 modes - To be done

     Representation     8      1 modes - To be done

     Representation     9      3 modes - To be done

     Representation    10      2 modes - To be done

     Representation    11      3 modes - To be done

     Representation    12      2 modes - To be done

     Representation    13      1 modes - To be done

     Representation    14      3 modes - To be done

     Representation    15      3 modes - To be done

     Representation    16      1 modes - To be done

     Representation    17      3 modes - To be done
  IOS =           253

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from davcio : error #        20
     i/o error in davcio
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
[0] MPI Abort by user Aborting program !
[0] Aborting program!
p0_5498:  p4_error: : 0

-----------------------------------------------stop here

There is also a CRASH file with the content:

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         2
     from set_irr : error #         0
     npert > max_irr_dim
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         3
     from set_irr : error #         0
     npert > max_irr_dim
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


-------------------------------------------------------------

I have searched the pw forum and some one suggests that one should
increase the max_irr_dim in the phcom.f90.  However, it seems not work
after I increased the value.

http://www.democritos.it/pipermail/pw_forum/2005-March/002167.html


Thank you in advance.

Best wishes,

Jian