c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c cc Calculate the Energy loss rate due to ccc c cc Bremsstruhlng process ccc c cc Using the Fitting Formulae ccc c cc By Hayashi Hiroshi ccc c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c cc main routine ccc c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c parameter ( nn = 10000 ) implicit real*8 (a-h,o-z) integer fityp dimension qb(nn) , tem(nn), rho(nn), tf(nn) , fityp(nn) & ,gamma(nn) common /nika2/ a , z , ic c call init1 call init2 ( rho , tem , nn ) call tfermi( rho , tf , nn ) call sftype( rho , tem , tf , fityp ,gamma, nn ) call qbcalc( rho , tem , tf , fityp , qb ,gamma, nn ) call output( rho , tem , fityp , qb ,gamma, nn ) stop end