Se presenta un nuevo código paralelo basado en modelos de hidrodinámica relativista especial para describir flujos interactivos. Se considera y verifica un nuevo método numérico altamente preciso. Se describe en detalle una implementación paralela del método mediante la tecnología Coarray Fortran y su eficiencia. La escalabilidad del código es del 92% en un clúster con Intel Xeon 6248R NKS-1P con 192 imágenes de Coarray Fortran. Se consideran diferentes flujos relativistas interactivos como aplicaciones astrofísicas.