Cooling function

The routine W_HD.f prompts for values of the density, , and kinetic temperature, T(K).  These parameters should lie in the ranges

and The cooling function W is expressed per HD molecule and calculated in erg s-1 (10-7 W).