/*theavy time dependent heavy element ionization balance */ void theavy(long int );