| 查看: 1719 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
【求助】50金币求大家帮忙改一段程序
|
|||
|
有一组数据如下所示,第一列代表istep变化范围(1-5000),第二列ia(1-30000),第三列ib(1-800),现在我想用一个三维数组H(istep,ia,ib)来标记这些数据,同一行数字定义的H=1,其他的H均为0.也就是说当istep=1时,ia=781,ib=455时H=1。我的主要目的就是能够将istep,ia,ib用一个函数联系起来,因为在程序后面要用这个函数作为进一步的判断条件。我写了一个程序: program exam implicit none integer ia,ib,istep integer ::h(5000,30000,800) integer i h=0 open (10,file='test.out') do i=1,99999999 read(10,*,end=2,err=2) istep,ia,ib h(istep,ia,ib)=1 enddo 2 continue close(10) do ia=1,30000 do ib=1,800 do istep=1,5000 h(istep)=h(istep,ia,ib) enddo do ntel=1,nsteps if(h(ntel)/=0) then 进行后面的程序 ----------------- endif enddo enddo enddo end 但是现在的问题是,定义的三维数组h太大了,程序不能运行, 请大家帮忙看看应该怎么办?先谢了 数据: 1 758 22 1 820 245 1 2050 399 1 2587 232 1 15733 34 1 18193 65 1 20718 256 2 820 245 2 2050 399 2 2587 232 2 20251 125 3 758 22 3 820 245 3 1817 625 3 2050 399 3 16723 555 3 18193 65 3 20718 256 --- ---- ---- |
» 猜你喜欢
有多少人是今天查系统知道结果的?
已经有17人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有8人回复
导师吐槽:我怎么摊上了这么个极品研究生!
已经有7人回复
青基如何破局
已经有7人回复
基金不中,共勉
已经有10人回复
国自然面上复盘~欢迎讨论
已经有11人回复
基金未中,这种答复是模板吗?
已经有6人回复
看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
怎么看青基中了没有啊
已经有5人回复
梦想
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
【急】请大家帮忙看一下这个offer是否合乎要求,新手散光所有金币求助~~
已经有12人回复
帮忙修改一段我的翻译,好难翻啊555
已经有4人回复
也请高人指点,帮忙看看每位评委的意见是什么?(50金币,回帖1金币)
已经有84人回复
帮忙修改一段翻译!谢谢!
已经有3人回复
求高手帮忙改一下英文摘要
已经有12人回复
【求助】撒点金币请大家帮忙出主意,回国都带什么礼物啊
已经有20人回复
急求,30分钟内翻译一段文字(药学),30分钟内完成6金币再追加50金币
已经有1人回复
再求一段英译中,在线等,急,50金币
已经有5人回复
【求助】10金币请帮忙修改一下MATLAB程序
已经有9人回复
» 抢金币啦!回帖就可以得到:
三无人员面上一把中散金币
+2/1910
面上得中
+1/466
中国基,散金币!!!!!!
+5/149
北京化工大学何静教授团队诚聘青年英才和博士后
+1/94
东华大学材料科学与工程学院丁奇老师招收2026级硕士研究生
+1/88
欢迎报考中科院苏州纳米所的氢能、燃料电池方向的博士生!(长期有效)
+1/55
欢迎保送和统考中科院苏州纳米所的氢能、燃料电池方向的研究生!(长期有效)
+1/54
中国科学院苏州纳米技术与纳米仿生研究所 仿生自适应智能系统实验室招收联培生
+1/39
现在生物医药相关有什么兼职可以做吗?
+1/35
国家蛋白质科学中心(北京),招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/17
国家蛋白质科学中心(北京),招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/16
新审查指南堵死了"一案双申"的老路,实用新型和发明只能二选一了?
+1/13
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/6
项目助理招聘启事
+1/5
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/4
今年没中,请各位大神帮我看一下,评审结果,大概几个ABC。
+1/3
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/2
细胞治疗全流程关键检测技术与研究策略
+1/2
重大突破!酶法合成破解生物法DNA合成原料“卡脖子”难题
+1/1
量子物理团队招收博士生与博士后
+1/1
|
以我的第一行数据 1 758 22 为例,1代表istep,758代表ia,22代表ib。 我想用函数将istep,ia,ib联系起来,只有当取我数据中一行上的数字时,定义的函数h=1,其他的h均为0.以第一行为例, istep=1,ia= 758,ib=22时,h=1.或者用多个函数也行,只要能完成后面循环中判断( if(h(ntel)/=0) then)就行。 do ia=1,30000 do ib=1,800 do istep=1,5000 h(istep)=h(istep,ia,ib) enddo do ntel=1,nsteps if(h(ntel)/=0) then 进行后面的程序 ----------------- endif enddo enddo enddo |
5楼2011-03-29 03:05:47
4楼2011-03-28 23:54:35
6楼2011-03-29 15:16:15
|
nstep可以写成5000, 对我的文件(new_total.out)的说明: 文件的格式为默认写入格式。 第一列代表ia,第二列ib,第三列istep(变化范围1-5000) 文件所在网址:http://good.gd/1078236.htm [ Last edited by astringent on 2011-3-30 at 05:47 ] |
7楼2011-03-30 05:44:34









回复此楼