How to output several arrays?

[tocarnas]

I am writing a program which computes the LU decomposition of a matrix, with partial pivoting, and I would like the function to output several (2 or 3) matrices without running the program several times, which is a waste of time since it gets me everything I want in one run. Is there a way of doing this?

 

[mikrom]

Inseveral examples I posted here, there are subroutines print_matrix, e.g. here:

http://www.tek-tips.com/viewthread.cfm?qid=1702652

or search in this forum

http://www.tek-tips.com/search.cfm?...t=&h=mikrom&searchSite=1&searchArea=1&d=1&s=1

 

[salgerman]

program pm
    real :: matrix(2,4)
    character(len=20) :: fmt
    matrix(1,:) = [ 1, 2, 3, 4]
    matrix(2,:) = [ 5, 6, 7, 8]
    write(fmt,"('('i2,'(f8.2))')") size(matrix,2)
    write(*,fmt) transpose(matrix)
end program pm