Fortran - Programacion de multiplicacion de matrices complejas en subrutina

Buenas Tardes comunidad,

Antes que nada, gracias por cualquier comentario o sugerencia.

Mi duda es la siguiente, estoy trabajando con fortran 90, y estoy haciendo el producto de matrices complejas en una subrutina. La cuestion es que todas las matrices son complejas, asi que para verificar que el producto es correcto, determino la parte imaginaria. La cuestion es que al hacer la tercera multiplicacion no la realiza.
el producto es el siguiente,

! ---------------------------------------------------------------------------
matD = matmul(matFG1,matFSRL);
PartD1 = matmul(matFSLR,matD);
matD = matmul(matFG4,PartD1)
! ---------------------------------------------------------------------------

donde matFG1, matFSRL matFSLR, matFG4 son matrices complejas que son leidas desde otra subrutima.

La cuestion es que para verificar que es correcto, las escribo en un archivo con el mismo formato, pero al realizar la tercera multiplicacion y respectivamente escribir el archivo varios elementos de la matriz son (asterisco) ********.

Mi pregunta un poco abusiva porque sucede esto y como lo puedo corregir?

Hola Mauricio.
Tu problema es que no usas el formato adecuado en los datos de salida. Para que quede más claro pongo un ejemplo:
la salida de este programa será:
5.5
****
justamente tu problema. La función WRITE(*,'(F4.1)') dice que escriba en konsola un flotante de 4 dígitos y un decimal. Siempre debes considerar el punto y el signo negativo, por ejemplo en el caso anterior si escribes WRITE(*,'(F5.1)') y la salida será:
5.5
-15.5

Saludos