PROGRAM pmesh
c
c-----------------------------------------------------------------
c     purpose : to eliminate the string from existed mesh data
c               warning : the mesh has to be made out of quadrangles only
c     in :c
c     out :
c             x     : array containing node coordinates
c             ia    : nodes reference number
c             ien   : list of element nodes
c             mat   : material ref. numbers for all elements
c             in    : element sides ref. number
cc-----------------------------------------------------------------------
c
      implicit real*8 (a-h,o-z)
      real xreal,yreal
      integer*4 numnp, numel                           
      real, allocatable :: x(:,:)
      integer*4, allocatable :: ien(:,:), ia(:), mat(:), in(:,:)
!     dimension x(2,numnp),ia(numnp)
!     dimension ien(4,numel),mat(numel),in(4,numel)
   
      open(57,file='pmesh1.dat')
      open(58,file='pmesh2.dat')
      
         call sizemesh(numnp, numel)
         allocate(x(2,numnp),ia(numnp),ien(4,numel),mat(numel),
     +             in(4,numel))
         open(21,file='mesh.dat')
          read(21,*)
          do 300 i=1,numnp                  
          read(21,2100)idum,xreal,yreal,ia(i)
              x(1,i) = xreal
              x(2,i) = yreal
          write(57,*)idum,xreal,yreal,ia(i)
300    continue
          do 500 n=1,numel
          read(21,2300)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
          write(58,*)ndum,(ien(i,n),i=1,4),(in(i,n),i=1,4)
500    continue
          close(21)

2100  format(i5,2x,e15.5,2x,e12.5,2x,i2)
2300  format(5(i5,2x),4(i2,2x))
          end
c**************************************************************


      subroutine sizemesh(numnp,numel)

        integer numnp, numel

      open(21,file='mesh.dat')
      read(21,*) numnp,numel
      close(21)
      return
      end

4楼: Originally posted by li_clifff at 2011-08-24 12:53:08:
但是我真看到有的f77程序主程序不用allocatable，直接调用子程序生成的numnp, numel就可以，就是按我那样开始写的，不知道自己按那样的方式怎么行不通了

2楼: Originally posted by snoopyzhao at 2011-08-24 11:13:08:
你主程中的 numnp 和 numel 必须赋值，否则必须使用动态数组，那么你的 x, ia, ien, mat, in 等数组就需要用 allocatable 之类的关键词重新声明，然后在调用子程序 sizemesh 之后，用 allocated 来分配内存。

7楼: Originally posted by snoopyzhao at 2011-09-28 11:23:24:
问题是，如果 numnp 没有值，主程序中相关的变量就不能成功声明，所以还是不行……

8楼: Originally posted by li_clifff at 2011-09-28 11:40:48:
恩，明白，那只有include 一个参数文件，如 Include 'Parameters.inc',
把numnp这样的参数 赋值在里面应该可以吧？