实元k_end,k_start,uz_int中的uz_int并没有在公共块中设置成real*8，而在子程序中相应的哑元S确设置成了real*8，这可能就是不一致的原因。(不知道你把uz_int设置成什么？)再有kxy_integral(k_end,k_start,uz_int)应该是call kxy_integral(k_end,k_start,uz_int)
