24小时热门版块排行榜    

查看: 1648  |  回复: 13
【悬赏金币】回答本帖问题,作者夕阳西下将赠送您 20 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

夕阳西下

金虫 (著名写手)

[求助] 程序无法运行怎么回事

program L
implicit none
integer,parameter:: nx=100,nt=100,p=2,te=5,E=2
integer:: i,j
real v

real:: ne(x,t)
real:: n0=5E+2,m=2.1,q=0.6E-19,k=0.514E-9,R=50.0,me=1.2E-3
v=n0*k

x=0.0
t=0.0
do i=1,nx
  x=x+i*0.02
  do j=1,nt
  t=t+j*0.01
  ne(0,0)=0.6E-19
ne(x+i*0.02,t+(j+1)*0.01)=-25*te/(me*v)*ne(x+(i+1)*0.02,t+j*0.01)+(1+50*te/(me*v))*ne(x+i*0.02,t+j*0.01)-25*te/(me*v)*ne(x+(i-1)*0.02,t+j*0.01)
  
open (unit=10,file="1.txt"
write (10,*) ne(x,t)

  end do
end do

stop
end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
12楼: Originally posted by 夕阳西下 at 2012-03-01 17:14:07:
你是不是搞错了,我们有这样写啊啊!后面的程序是书上的,我没有改动。

你给我找找看你的程序中哪里有 allocatable 这样的关键字……
13楼2012-03-01 17:54:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

lurencyj

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
余泽成(金币+2): 谢谢参与应助!欢迎常来程序语言版! 2012-03-04 01:31:56
fortran的东东我都忘了,哈哈,只是楼主写程序的时候,风格很差,难读。

可以一步一步自己测试的,不会用debug工具的话,可以用print 屏幕输出; 如果,编译通不过的,那就摸着石头过河,注释掉不确定的,从程序头部开始慢慢往下测。
很女子很弓虽大
2楼2012-02-20 16:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
余泽成(金币+2): 谢谢参与应助! 2012-03-04 01:32:09
你这个程序中有太多的错误,先根据编译器的提示进行修改……

举个例子:
real:: ne(x,t)
中的 x,t 都不是常数,那么,你的 ne 的定义就是错误的
ne(0,0)=0.6E-19
默认情况下数组的下标是从 1,1 开始的……
ne(x+i*0.02,t+(j+1)*0.01)=
下标不能是实数,只能是整数……
等等……
3楼2012-02-23 14:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)

引用回帖:
3脗楼: Originally posted by snoopyzhao at 2012-02-23 14:09:20:
脛茫脮芒赂枚鲁脤脨貌脰脨脫脨脤芦露脿碌脛麓铆脦贸拢卢脧脠赂霉戮脻卤脿脪毛脝梅碌脛脤谩脢戮陆酶脨脨脨脼赂脛隆颅隆颅

戮脵赂枚脌媒脳脫拢潞
real:: ne(x,t)
脰脨碌脛 x,t 露录虏禄脢脟鲁拢脢媒拢卢脛脟脙麓拢卢脛茫碌脛 ne 碌脛露篓脪氓戮脥脢脟麓铆脦贸碌脛
ne(0,0)=0.6E-19
脛卢脠脧脟茅驴枚脧脗脢媒脳茅碌脛脧脗卤锚脢脟麓脫 1,1 驴陋脢录碌脛隆颅隆颅 ...

ne(x,t)脠莽鹿没卤铆脢戮碌脛脢脟潞炉脢媒鹿脴脧碌脛脴拢驴脪虏戮脥脢脟脣碌ne鹿脴脫脷x,t碌脛潞炉脢媒拢卢赂脙脭玫脙麓露篓脪氓ne?
4楼2012-02-23 23:33:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我面上完蛋了 +10 且听虎啸 2026-08-20 11/550 2026-08-25 07:12 by echo8914667
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见