10/11/01 RK with Newton added 06/01/31 power absorption carrected in wrcalc.f (thanks to Dr. Idei)