关于自由格式读写数组的问题，求指导

Harry · 发表于 2015-12-7 21:57:56

[Fortran] 纯文本查看 复制代码

01subroutine Stiffb
02    common/fem1/npoin,nelem,nboun,nprop,nnode,nevab,nsvab,ndofn,ndime,nstre
03    common/fem2/props(5,2),coord(50,2),lnods(75,2),ifpre(100),fixed(100),rload(50,2),eload(75,4),&
04         matno(75),stres(75,1),xdisp(100),tdisp(50,2),treac(50,2),astif(100,100),aslod(100),react(100)
05    dimension estif(4,4)
06!
07!  evaluation of member stiffness matrices
08!  for pin-jointed plane frameworks
09!
10    open(unit=10,file='input.txt')
11    open(unit=12,file='output.txt')
12    open(unit=11,file='between.txt')
13    rewind 11
14    do 20 ielem=1,nelem
15    lprop=matno(ielem)
16    young=props(lprop,1)
17    xarea=props(lprop,2)
18    node1=lnods(ielem,1)
19    node2=lnods(ielem,2)
20    xproj=coord(node2,1)-coord(node1,1)
21    yproj=coord(node2,2)-coord(node1,2)
22    eleng=sqrt(xproj**2+yproj**2)
23    sinth=yproj/eleng
24    costh=xproj/eleng
25    fmult=young*xarea/eleng
26    estif(1,1)=fmult*costh**2
27    estif(1,2)=fmult*sinth*costh
28    estif(2,1)=fmult*sinth*costh
29    estif(2,2)=fmult*sinth**2
30    do 10 inode=1,nnode
31    do 10 jnode=1,nnode
32    kount=(-1)**inode*(-1)**jnode
33    do 10 knode=1,nnode
34    do 10 lnode=1,nnode
35    index=(inode-1)*nnode+knode
36    jndex=(jnode-1)*nnode+lnode
37 10 estif(index,jndex)=kount*estif(knode,lnode)
38    write(11) estif                          !问题在此行
39 20 continue
40    return
41  end