read(12, *) ((a(i,j), i=1,5),j=1,5)

Originally posted by david830wu at 2011-07-17 18:45:35:
再进一步问一下，也是我真实遇到的，就是GAMBIT的网格划分文件，数据文件是十六进制应该怎么办呢？
MATRIX.dat
1   3f  a6    21   4
a cc        f3      2  2
3    d0    b2    46   4
……
这样的文件用 ...

program read_hex

implicit none
character(len=80) :: line="3 5  F6F  1d 0 0"
character(len=80) :: buff, hex_str, fm
integer, dimension(3) :: hex
integer :: str_len, i, first, last

buff = adjustl(line)
str_len = len_trim(buff)

!! remove the extra space in the string
do i=1,str_len-1
   if (buff(i:i) == ' ' .and. buff(i+1:i+1) == ' ') buff(i:) = buff(i+1:)
end do

!! extract the hexadecimal constant
first = 1
do i=1,3
   first = index(buff(first:), ' ') + first - 1
   last = index(buff(first+1:), ' ') + first
   hex_str = buff(first+1:last-1)
   write(fm, '(a,i0,a)') '(Z', len_trim(hex_str), ')'
   read(hex_str, fm) hex(i)
   first = last
end do

!! hexadecimal to decimal constrant, and output...
write(*,*) hex

end program read_hex