thanks for the reply. Sorry the input file contains only two colums like below:
2995.2000 -11813.2717 3064.7830 -11813.1994 3010.7103 -11813.2796 3017.9634 -11813.2844 3029.5311 -11813.2865 . . . We have uncommented the parts that converts the anergy in eV and the lattice parameters in Ang^3. And line 48 in the program is: "read(string(2:length),*,err=3) vol(i),ene(i)" What should I change in it to be able to read my input? Regards,