Power fix it

A long effort has been put towards fixing a bug with the power in the code. Now the code clearly identifies what is heat/power Q(:) and what is heat/power density Dens(:). With this clearly defined the code is refactored to use the correct term at all times.