24小时热门版块排行榜    

查看: 1027  |  回复: 1
【悬赏金币】回答本帖问题,作者修炼天木将赠送您 15 个金币

修炼天木

新虫 (初入文坛)

[求助] 关于Fortran编程遇到的问题

我是一个初学者,有大神的话可以私聊,我将程序发给您,请您帮忙解答。
受篇幅所限,这个只是一部分程序,如果需要更多,烦请私聊我

program caphere

dimension porerho(3000),porethta(3000),porephi(3000),porerad(3000)

real inclrho(8000),inclthta(8000),inclphi(8000),inclrad(8000)

integer nsurf(100),nsepd(100)

dimension ashsum(100),ashsum2(100),bshsum(100),bshsum2(100)

real az(100),surfbin(100,15),sepdbin(100,15)

common /inclusion/ inclrho, inclthta, inclphi, inclrad

common /pores/ porerho, porethta, porephi, porerad

common /dimensions/ radchar, radpore, radincl

common /parameters/ iupper, ilower, issed, delz

common /bins/ az, surfbin, sepdbin, nsurf, neepd

common /sums/ ashsum, ashsum2, bshsum, bshsum2

open(5, file='csin')

open(7, file='csout1')

open(8, file='csout2')

read(5,*)nruns,mruns

read(5,*)iseed

read(5,*)radchar, delz

do 10 ik=1,nruns

read(5,*)radpore, radincl

read(5,*)fvpore, fvincl

npores=fvpore*(radchar/radpore)**3

nincls=fvincl*(radchar/radincl)**3

650 format(20x,'mruns=',i8)

800 format(14x'char radius=',f16.1)

325 format(/10x'number of pores=',i8/5x,'number of inclusions=',i8)

600 format(//5x,45('*')

$ /5x,'**** results from the continuous model ****'

$ /5x,45('*')//)

750 format(20x,'iseed=',i12)

do 5 i=1,100

nsurf(i)=0 nsepd(i)=0 ashsum(i)=0. ashsum2(i)=0. bshsum(i)=0. bshsum2(i)=0. do 6 j=1,15 surfbin(i,j)=0 sepdbin(i,j)=0

6 continue

5 continue

do 15 j=1,mruns

vfrpore=0 vfrincl=0 npores=fvpore*(radchar/radpore)**3 write(6,*)'Starting the run no.',j,'......'

call distribn(j,npores,nincls,vfrpore,vfrincl)

if(npores.ne.0)then call sort(1,npores)

endif call sort(2,nincls) call reaction(iz,npores,nincls)

call results(vfrpore,vfrincl,nincls) 这一行报错

write(6,*)'Finished the run no.',j,'......'

15 continue

call results(iz,mruns,vfrpore,vfrincl,nincls)

10 continue stop end

错误如下、 error #6631: A non-optional actual argument must be present when invoking a procedure with an explicit interface. [NINCLS]

error #6631: A non-optional actual argument must be present when invoking a procedure with an explicit interface. [VFRINCL]

error #6633: The type of the actual argument differs from the type of the dummy argument. [NINCLS]

error #6633: The type of the actual argument differs from the type of the dummy argument. [VFRINCL]

error #6633: The type of the actual argument differs from the type of the dummy argument. [VFRPORE]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

修炼天木

新虫 (初入文坛)

2楼2019-06-02 08:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 修炼天木 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +11 王凯12 2024-06-13 22/1100 2024-06-14 19:45 by yifeng11
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +3 babu2015 2024-06-14 3/150 2024-06-14 19:34 by 星火12
[基金申请] 2024国社科通讯评审 +6 qsd10086 2024-06-13 7/350 2024-06-14 18:02 by thesuna
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +7 味道很好啊 2024-06-13 7/350 2024-06-14 17:40 by zeolitess
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +7 安稳22123 2024-06-13 8/400 2024-06-14 14:43 by bear2007
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 投稿kbs被拒后,系统提供的推荐的其他期刊有用吗? +3 chenrui2015 2024-06-10 3/150 2024-06-14 10:30 by Tr.sjx1997
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[教师之家] 请各位老师谈谈身边非升即走的例子 +10 Lucas1121 2024-06-09 25/1250 2024-06-13 21:52 by Lucas1121
[基金申请] 中者为王-观某校级项目为领导瓜分有感 +24 zqwcr 2024-06-08 36/1800 2024-06-13 14:53 by yuting
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[硕博家园] 交友 +5 xkkl 2024-06-08 6/300 2024-06-12 14:47 by xkkl
[教师之家] 难上的课安排给我,杂事安排给我。评优时我就没份儿了,即便我做得不比绝大多数人差 +15 河西夜郎 2024-06-08 15/750 2024-06-12 12:51 by oooooo?o
[论文投稿] 给大家分享一下扎心的投稿经历 +10 lizhengke06 2024-06-08 13/650 2024-06-11 20:37 by lizhengke06
[硕博家园] 求助 +6 LYWwrz 2024-06-09 9/450 2024-06-11 13:12 by powerhours
[考博] 希望能25博士入学,可提前一年做科研助理 +3 干饭版小太阳 2024-06-09 7/350 2024-06-10 00:05 by 干饭版小太阳
信息提示
请填处理意见