1

Estoy aprendiendo a programar en FORTRAN muy despacio para poder hacer mi trabajo fin de carrera. Estoy interesado en saber escribir un código que me permita almacenar datos de un fichero (.dat o .txt) a un array multidimensional. He intentado el siguiente código para un fichero que tiene 15 filas y 4 columnas, y quería que lo almacenase en una matriz de 15 x 4 a la que he llamado matrix:

program almacenararray

    implicit none

    integer :: i,j

    integer,dimension(15,4) :: matrix


    open (101,file='/Users/manuelperezescribano/Desktop/PROGRAMAS TFG/localizacionmaximosyminimos/datos.txt',status="old")


    do i=1,15
        do j=1,4
            read(101,*)
            matrix(i,j)=i+j
        end do
    end do

    close (101)


    print*,matrix(i,j)

end program almacenararray

Sin embargo, el resultado no es el esperado: la matriz no es impresa por pantalla, con lo que no sé si se han almacenado los datos del fichero .txt a la matriz o array correctamente. ¿Cómo puedo solucionar el error?

Gracias

1

2 respuestas 2

1

Tienes dos errores.

La manera en la que volcamos lo leído en el archivo a la matriz es:

do i = 1, 15
    do j = 1, 4
       read(101, *) matrix(i, j)
    end do
end do

Pero además la forma en la que imprimimos una matriz es elemento a elemento:

do i = 1, 15
    do j = 1, 4
       print*, matrix(i, j)
    end do
end do

Si quieres darle un formato más bonito podemos imprimir fila a fila.

do i = 1, 15
    print '(4I3)', (matrix(i, j), j = 1, 4)
end do

Con el argumento '(4I3)' estamos indicando que imprimiremos 4 elementos donde cada uno es un entero (I) de 3 caracteres. He escogido 3 caracteres asumiendo que el entero más largo que contenga la matriz será de 2 caracteres. De esta forma habrá un espacio de separación entre las columnas.

0

El codigo que explusio tiene algunos errores:

  1. El estado del archivo esta como status = old, significa que ya existio un archivo en ubicacion.
  2. cuando abres el archivo esta en la ubicacion 101, pero antes de que se escriban los datos en 101 cierras el archivo con close(101).
  3. para que se vea la matriz tiene que tener un formato, recomiendo usar write.
  4. read(101,*) no sabe que lee

codigo.f90

program almacenararray
    implicit none
    integer :: i,j
    integer,dimension(15,4) :: matrix
    open (101,file='archivo.txt',status="unknown")
   do i=1,15
        do j=1,4
            matrix(i,j)=i+j
        end do
    end do
do i=1,15
write(101,20)(matrix(i,j),j=1,4)
end do
20 format(1x,4(1x,i3) )
end program almacenararray

El formato que puse fue 10 format(1x,4(1x,i3) ) , es decir que para write(101,20) tendra 4 valores con un entero de 3 digitos y estos tendran 1 separacion.

La salida del archivo es:

    2   3   4   5
    3   4   5   6
    4   5   6   7
    5   6   7   8
    6   7   8   9
    7   8   9  10
    8   9  10  11
    9  10  11  12
   10  11  12  13
   11  12  13  14
   12  13  14  15
   13  14  15  16
   14  15  16  17
   15  16  17  18
   16  17  18  19

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.