<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><div><font face="arial" size="2">Dear All I am working on ZnO(hcp structure). I want to calculate thermal properties&nbsp;</font></div><div><font face="arial" size="2">using QHA.I am using fedora core 9 (gfortran version 4.3.0).&nbsp;</font></div><div><font face="arial" size="2">for thermal properties following errors occures.</font></div><div><font face="arial" size="2"><br></font></div><div><font face="arial" size="2">Would any body tell me how I can cover this problem</font></div><div><font face="arial" size="2">&nbsp;</font></div><div><font face="arial" size="2">[zafar@localhost ZnO]$ ./Run_Me&nbsp;</font></div><div><font face="arial" size="2">****** input tetrahedra for BZ-integration ******</font></div><div><font face="arial" size="2">&nbsp; &nbsp;6 &nbsp; 3 &nbsp;12</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.000000
 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.577350 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp;-0.333330 &nbsp;0.577350 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.000000 &nbsp;0.803030</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.577350 &nbsp;0.803030</font></div><div><font face="arial" size="2">&nbsp;-0.333330 &nbsp;0.577350 &nbsp;0.803030</font></div><div><font face="arial" size="2">****** input tetrahedra for BZ-integration ******</font></div><div><font face="arial" size="2">********************** generate_tetra *********************</font></div><div><font face="arial" size="2">&nbsp; NT0= &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 3 &nbsp;NTETMX= &nbsp; &nbsp; &nbsp; &nbsp;1728</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp;-0.3333 &nbsp;-0.3333</font></div><div><font face="arial"
 size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.8030</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.0000 &nbsp;-0.3333</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.8030 &nbsp; 0.8030 &nbsp; 0.8030</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000
 &nbsp; 0.0000 &nbsp;-0.3333 &nbsp; 0.0000</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.8030 &nbsp; 0.8030 &nbsp; 0.0000</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; total volume of BZ is = &nbsp; 0.0773 &nbsp;omg48= &nbsp;25.8830</font></div><div><font face="arial" size="2">****************** end of generate_tetra ******************</font></div><div><font face="arial" size="2">&nbsp;Recalculating omega(q) from C(R)</font></div><div><font face="arial" size="2">STOP Partial_DOS finished</font></div><div><font face="arial" size="2">&nbsp;natoms== &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 4</font></div><div><font face="arial"
 size="2">&nbsp;irec==== &nbsp; &nbsp; &nbsp; &nbsp; 168</font></div><div><font face="arial" size="2">&nbsp; 0.75000000000000000 &nbsp; &nbsp;&nbsp;</font></div><div><font face="arial" size="2">&nbsp;Zn &nbsp;Zn &nbsp;</font></div><div><font face="arial" size="2">It seems you have imaginary frequences.\ &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</font></div><div><font face="arial" size="2">&nbsp;Hopefully you know what you are doing</font></div><div><font face="arial" size="2">&nbsp;nstep==== &nbsp; &nbsp; &nbsp; &nbsp; 734</font></div><div><font face="arial" size="2">&nbsp; 112.0761 &nbsp;112.0767 &nbsp;136.8345</font></div><div><font face="arial" size="2">&nbsp; 120.2794 &nbsp;129.5796 &nbsp;132.5383</font></div><div><font face="arial" size="2">&nbsp; 123.3896 &nbsp;128.8859 &nbsp;131.0670</font></div><div><font face="arial" size="2">&nbsp; 126.1556 &nbsp;126.1563 &nbsp;131.4016</font></div><div><font face="arial"
 size="2">&nbsp; 123.0497 &nbsp;123.8439 &nbsp;148.2783</font></div><div><font face="arial" size="2">&nbsp; 125.4183 &nbsp;128.9576 &nbsp;144.2191</font></div><div><font face="arial" size="2">&nbsp; 126.6441 &nbsp;130.6791 &nbsp;142.6102</font></div><div><font face="arial" size="2">&nbsp; 126.1474 &nbsp;131.6514 &nbsp;141.8672</font></div><div><font face="arial" size="2">&nbsp; 127.5920 &nbsp;134.0565 &nbsp;139.4667</font></div><div><font face="arial" size="2">&nbsp; 128.0304 &nbsp;136.7692 &nbsp;136.7699</font></div><div><font face="arial" size="2">&nbsp; 116.8975 &nbsp;118.1902 &nbsp;151.5177</font></div><div><font face="arial" size="2">&nbsp;E_min= &nbsp; 0.0000000000000000 &nbsp; &nbsp; &nbsp; &nbsp; E_max= &nbsp; 530.33650000000000 &nbsp; &nbsp;&nbsp;</font></div><div><font face="arial" size="2">&nbsp;nstep==== &nbsp; &nbsp; &nbsp; &nbsp; 734</font></div><div><font face="arial" size="2">****** input tetrahedra for BZ-integration
 ******</font></div><div><font face="arial" size="2">&nbsp; &nbsp;6 &nbsp; 3 &nbsp;12</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.000000 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.577350 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp;-0.333330 &nbsp;0.577350 &nbsp;0.000000</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.000000 &nbsp;0.803030</font></div><div><font face="arial" size="2">&nbsp; 0.000000 &nbsp;0.577350 &nbsp;0.803030</font></div><div><font face="arial" size="2">&nbsp;-0.333330 &nbsp;0.577350 &nbsp;0.803030</font></div><div><font face="arial" size="2">****** input tetrahedra for BZ-integration ******</font></div><div><font face="arial" size="2">********************** generate_tetra *********************</font></div><div><font face="arial" size="2">&nbsp; NT0= &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 3 &nbsp;NTETMX= &nbsp; &nbsp; &nbsp;
 &nbsp;1728</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp;-0.3333 &nbsp;-0.3333</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.8030</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.0000 &nbsp;-0.3333</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.8030 &nbsp; 0.8030 &nbsp; 0.8030</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp;
 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.0000 &nbsp;-0.3333 &nbsp; 0.0000</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.5774 &nbsp; 0.5774 &nbsp; 0.5774</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.0000 &nbsp; 0.8030 &nbsp; 0.8030 &nbsp; 0.0000</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000 &nbsp; 1.0000</font></div><div><font face="arial" size="2">&nbsp; volume of tetrahedron = &nbsp;0.02576</font></div><div><font face="arial" size="2">&nbsp; total volume of BZ is = &nbsp; 0.0773 &nbsp;omg48= &nbsp;25.8830</font></div><div><font face="arial" size="2">****************** end of generate_tetra ******************</font></div><div><font face="arial" size="2">&nbsp;before integration: &nbsp;E_min= &nbsp; 0.0000000000000000 &nbsp; &nbsp;
 &nbsp; &nbsp; E_max= &nbsp; 530.33650000000000 &nbsp; &nbsp;&nbsp;</font></div><div><font face="arial" size="2">&nbsp;Freq,Tot_DOS,DOS== &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp; &nbsp; 0.000000 &nbsp; &nbsp;</font></div><div><font face="arial" size="2">At line 107 of file Integration.f (unit = 32, file = '')</font></div><div><font face="arial" size="2">Fortran runtime error: File already opened in another unit</font></div><div><font face="arial" size="2">&nbsp;ndiv from file === &nbsp; &nbsp; &nbsp; &nbsp; 734</font></div><div><font face="arial" size="2">&nbsp;ndiv=== &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 1</font></div><div><font face="arial" size="2">&nbsp;norm_partial== &nbsp; 0.0000000000000000 &nbsp; &nbsp;&nbsp;</font></div><div><font face="arial" size="2">&nbsp;ndiv from file === &nbsp; &nbsp; &nbsp; &nbsp;
 734</font></div><div><font face="arial" size="2">&nbsp;ndiv=== &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 1</font></div><div><font face="arial" size="2">&nbsp;norm_partial== &nbsp; 0.0000000000000000 &nbsp; &nbsp;&nbsp;</font></div><div><font face="arial" size="2">At line 58 of file Atom_projected_properties.f90 (unit = 9, file = 'projected.DOS')</font></div><div><font face="arial" size="2">Fortran runtime error: End of file</font></div><div><font face="arial" size="2">At line 71 of file Mean_square_displacement.f90 (unit = 1, file = 'projected_DOS.O')</font></div><div><font face="arial" size="2">Fortran runtime error: End of file</font></div><div><font face="arial" size="2">mv: cannot stat `Displacements': No such file or directory</font></div><div><font face="arial" size="2">At line 58 of file Atom_projected_properties.f90 (unit = 9, file = 'projected.DOS')</font></div><div><font face="arial" size="2">Fortran runtime error: End of file</font></div><div><font
 face="arial" size="2">At line 71 of file Mean_square_displacement.f90 (unit = 1, file = 'projected_DOS.O')</font></div><div><font face="arial" size="2">Fortran runtime error: End of file</font></div><div><font face="arial" size="2">mv: cannot stat `Displacements': No such file or directory</font></div><div><font face="arial" size="2">&nbsp;ndiv from file === &nbsp; &nbsp; &nbsp; &nbsp; 734</font></div><div><font face="arial" size="2">&nbsp;ndiv=== &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 0</font></div><div><font face="arial" size="2">&nbsp; &nbsp; 0.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 100.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 200.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp;
 &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 300.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 400.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 500.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 600.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 700.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 800.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp;
 &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp; 900.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1000.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1100.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1200.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1300.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1400.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp;
 &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1500.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1600.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1700.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1800.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;1900.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2000.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp;
 &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2100.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2200.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2300.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2400.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2500.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2600.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp;
 &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2700.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2800.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;2900.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">&nbsp;3000.00 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000 &nbsp; &nbsp;0.00000000</font></div><div><font face="arial" size="2">Phonon DOS and Quasiharmonic calculations have finished.</font></div><div><font face="arial" size="2">Now you can analyse these data using Gnuplot or xmgrace</font></div><div style="font-family: arial; font-size: 10pt; "><br></div><br><font face="arial" size="2">Muhammad
 Zafar</font><br><font face="arial" size="2">
PhD Scholar</font><br><font face="arial" size="2">
Department of Physics</font><br><font face="arial" size="2">
The Islamia University of Bahawalpur,Pakistan</font></td></tr></table>