In reactor calculations, a detailed 3-D power density distribution is required for core optimization studies and safety analyses. The general Monte Carlo based neutral particle transport tool, MCNP [1], has the capability to obtain detailed power density distribution in a reactor through its criticality calculation mode (KCODE mode). However, there is no standard tally type in MCNP that is able...