应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 请教一个fortran小程序编译出错的问题，谢谢
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 2580  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

li_clifff

银虫 (正式写手)

[求助] 请教一个fortran小程序编译出错的问题，谢谢

如下程序编译不成功，不知道问题出在哪里，mesh.dat文件格式是第一行是numnp, numel的数值，第二行到最后是mesh的各个四边形的具体坐标。错误信息是：Error: A specification expression object must be a dummy argument, a COMMON block object, or an object accessible through host or use association   [NUMNP]  integer*4 numnp, numel   
还有G:\Program Files\Microsoft Visual Studio\MyProjects\mesh\pmesh.for(22) : Error: An automatic object is invalid in a main program.   [X] dimension x(2,numnp),ia(numnp)，请多多指教，谢谢                        


  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                           
      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)
         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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-08-24 11:01:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-08-24 16:26:52
ben_ladeng: 专家考核存档 2011-09-28 11:39:50
你主程中的 numnp 和 numel 必须赋值，否则必须使用动态数组，那么你的 x, ia, ien, mat, in 等数组就需要用 allocatable 之类的关键词重新声明，然后在调用子程序 sizemesh 之后，用 allocated 来分配内存。
赞一下(2人) 回复此楼
2楼2011-08-24 11:13:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
li_clifff(金币+20): 确实是这样，动态数组，学习了，谢谢 2011-08-24 11:36:10
余泽成(金币+3): 辛苦了！ 2011-08-24 16:27:14
dubo(金币+5, 程序强帖+1): 多谢应猪 2011-08-27 16:42:39
ben_ladeng: 专家考核存档 2011-09-28 11:40:00
上面错了一点儿，不是用 allocated，而是用 allocate 来分配内存。下面是修改后的程序：
CODE:
      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
赞一下(3人) 回复此楼
3楼2011-08-24 11:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
★
dubo(金币+1): 欢迎讨论 2011-08-27 16:41:30
但是我真看到有的f77程序主程序不用allocatable，直接调用子程序生成的numnp, numel就可以，就是按我那样开始写的，不知道自己按那样的方式怎么行不通了
赞一下(1人) 回复此楼
4楼2011-08-24 12:53:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★ ★
余泽成(金币+1): 谢谢参与应助！ 2011-08-24 16:27:27
dubo(金币+1): 欢迎讨论 2011-08-27 16:41:43
ben_ladeng: 专家考核存档 2011-09-28 11:40:16
引用回帖:
4楼: Originally posted by li_clifff at 2011-08-24 12:53:08:
但是我真看到有的f77程序主程序不用allocatable，直接调用子程序生成的numnp, numel就可以，就是按我那样开始写的，不知道自己按那样的方式怎么行不通了

在主程序中不行，在子程序中有可能可行，但不是你这种写法……

也就是说，如果声明变量时，numnp, numel 有值就行，否则不行……
赞一下(3人) 回复此楼
5楼2011-08-24 14:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

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

如果在引用call meshsize之前对numnp赋值，或者打开一个文件，读取参数，对numnp赋值，那么是不是就没有问题了呢
赞一下 回复此楼
6楼2011-09-28 11:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
问题是，如果 numnp 没有值，主程序中相关的变量就不能成功声明，所以还是不行……
赞一下 回复此楼
7楼2011-09-28 11:23:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
引用回帖:
7楼: Originally posted by snoopyzhao at 2011-09-28 11:23:24:
问题是，如果 numnp 没有值，主程序中相关的变量就不能成功声明，所以还是不行……

恩，明白，那只有include 一个参数文件，如 Include 'Parameters.inc',
把numnp这样的参数 赋值在里面应该可以吧？
赞一下 回复此楼
8楼2011-09-28 11:40:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_clifff

银虫 (正式写手)
我发现早期f77程序，对数组没有动态赋值，都是先定义一个大的数组比如
x（nnpmax）,nnpmax在parameter.inc文件里，而实际上子程序里的大小nnp是变化的，但只要不超过nnpmax就行，这样的方式也不错。
赞一下 回复此楼
9楼2011-09-28 11:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
引用回帖:
8楼: Originally posted by li_clifff at 2011-09-28 11:40:48:
恩，明白，那只有include 一个参数文件，如 Include 'Parameters.inc',
把numnp这样的参数 赋值在里面应该可以吧？

是这样的，呵呵……
赞一下 回复此楼
10楼2011-09-28 12:34:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li_clifff 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分调剂 +7 Faune 2026-03-30 7/350 2026-03-30 13:26 by 帕尔马拉特
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +13 丹青奶盖 2026-03-26 15/750 2026-03-30 12:35 by fangnagu
[考研] 310求调剂 +9 争取九点睡 2026-03-30 9/450 2026-03-30 11:44 by Wang200018
[考研] 085600材料与化工调剂 +6 kikiki7 2026-03-30 6/300 2026-03-30 10:57 by 星空星月
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +11 辰沐5211314 2026-03-26 11/550 2026-03-30 09:00 by 无际的草原
[考研] 318一志愿吉林大学生物与医药 求调剂 +5 笃行致远. 2026-03-28 5/250 2026-03-30 06:56 by ilovexiaobin
[考博] 26申博自荐 +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 343求调剂 +5 爱羁绊 2026-03-28 5/250 2026-03-28 20:53 by 唐沐儿
[考研] 320分，材料与化工专业，求调剂 +9 一定上岸aaa 2026-03-27 13/650 2026-03-28 15:00 by 神马都不懂
[考研] 292求调剂 +14 鹅鹅鹅额额额额� 2026-03-25 15/750 2026-03-28 08:45 by WYUMater
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 085600，材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 0703化学338求调剂！ +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 各位老师您好：本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
信息提示
关闭
请填处理意见
关闭