24小时热门版块排行榜    

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

至尊木虫 (职业作家)

【答案】应助回帖


余泽成(金币+1): 谢谢参与应助! 2012-03-04 01:33:02
引用回帖:
10楼: Originally posted by 夕阳西下 at 2012-02-29 13:15:58:
麻烦告知,怎么申明?

大概这个样子:
CODE:
real, dimension(:,:), allocatable :: a

找本书看吧,从你贴出来的问题来看,你还是需要一本基础的书,让你知道 Fortran 的语法,以及程序结构啥的……
11楼2012-02-29 13:26:39
已阅   回复此楼   关注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个字符以上)
信息提示
请填处理意见