  We present a complete set of diagnostic tools aimed at reproducing synthetic
non-thermal (synchrotron and/or Inverse Compton, IC) emissivity, integrated
flux energy, polarization and spectral index simulated maps in comparison to
observations. The time dependent relativistic magnetohydrodynamic (RMHD)
equations are solved with a shock capturing code together with the evolution of
the maximum particles energy. Applications to Pulsar Wind Nebulae (PWNe) are
shown.
