sim_power_dens function

Simulate power densities