24小时热门版块排行榜    

查看: 1795  |  回复: 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

至尊木虫 (职业作家)

【答案】应助回帖

jjdg(金币+1): 感谢参与 2012-02-29 02:24:36
在声明 A(n,n) 时,n 必须有值,而不能等后面再赋值,除非你动态分配内存……
9楼2012-02-28 18:39:05
已阅   回复此楼   关注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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:08 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:04 by BCdvBAEpt3zm
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 15:04 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 14:29 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 14:24 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
请填处理意见