c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c cc tfermi routine ccc c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c subroutine tfermi( rhom , tf ,nn ) c implicit real*8 (a-h,o-z) real*8 mue common /nika2/ a , z , ic dimension rhom(nn), tf(nn) c mue = a/z do 100 i = 1, ic rho = rhom( i )*mue rho6 = rho/1.d6 tf(i) = 5.9302d9*( sqrt( 1.d0 + 1.018d0 & *(z/a)**(2.d0/3.d0)*rho6**(2.d0/3.d0) ) -1.d0 ) c write(*,*) 'tf===> ',rhom(i), tf(i) 100 continue return end