program prep2d
*
      implicit double precision (a-h)
      implicit double precision (o-z)
      parameter (maxelm=15000,mxnode=8000,mxprof=1200000) 
      double precision x(mxnode),y(mxnode)
      integer node(maxelm,3),nd(mxnode),itype(maxelm),jdiag(mxnode)
      integer renum(mxnode+1),lvl(mxnode),lvls1(mxnode),lvls2(mxnode)
      integer ccstor(mxnode),ndeg(mxnode),ncon(mxnode,20)
    character*1,ch
      common /fe/ nnode,nprof
      common /gra/ nequ,idpth,ideg
*
* read data for element definition and nodal coordinates
*
      open (7,file='base.fpl',form='formatted')
*
      rewind 7
*
    
      read (7,*) nelm,nnode
      read (7,*) (x(i),i=1,nnode)
      read (7,*) (y(i),i=1,nnode)
      read (7,*) ((node(i,j),j=1,3),i=1,nelm)
      read (7,*) (itype(i),i=1,nelm)
    print*,'No of elements = ',nelm
    print*,'No of nodes = ',nnode
    
    end
    
这是代码，很简单，就是有限元里先读文档，然后把文档里的数存一下。但是运行就出错，见截图。高手求解答。。。