real K(1:8922,1:8922),a(1:89103,1:3)
open(unit=1,file='filestiffness.txt',status='old')
open(unit=2,file='kk.txt',status='new')
read(1,*)((a(i,j),j=1,3),i=1,89103)
do 10 i=1,8922
do 20 j=1,8922
do 30 n=1,89103
if (a(n,1).eq.i.and.a(n,2).eq.j) then
K(i,j)=a(n,3)
end if
30 continue
20 continue
10 continue
!100 format(1x,8922(e19.12,1x))
write(2,*) ((K(i,j),j=1,8922),i=1,8922)
end
这个程序的功能是将数组a的第一列和第二列作为行号和列号给第三列在k中定位.
filestiffness.txt内容如下:
10 9 2558
6 4 5285
3 3 5582
…………
……………
[ 本帖最后由 yxb345 于 2007-5-22 16:39 编辑 ] |