24小时热门版块排行榜

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者夕阳西下的 4 个金币，回帖就立即获得 1 个金币，每人有 1 次机会

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

夕阳西下

金虫 (著名写手)

应助: 138 (高中生)
贵宾: 0.041
金币: 9247.1
帖子: 1268
在线: 1164.4小时
虫号: 220750

[交流] 程序运行不了问题在哪

下面这个程序有错误，怎么改呢？
program Ne
integer :: I,J, N, Ny
real::  ne(J)
real::  Hx(I)
N=1000
Ny=1000
Dx=0.1
Dt=0.1
n0=5E+2
m=2.1
p=2
Te=5
q=0.6E-19
E=2
k=0.514E-9
v=n0*k
Hx(0)=sqrt(8*Te/pm)*ne*(1-0.25)
Hx(I)=-(q/m*v)*E(ne/(1-exp(q*E*Dx/Te)+ne/(1-exp(-q*E*Dx/Te))

open (unit=10,FILE='1.txt')
DO I=2,N-1

  DO J=2,Ny-1
HX(I)=(H1(I+1)-H1(I-1))/(2*DX)
ne(J)=(ne(J+1)-ne(J-1))/(2*Dt)
ne(J)=Hx(I)+R
end do
end do

write (10,"('N=?,Ny=?')" ne(I)
stop
end

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有10人回复
2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-02-05 22:36:24

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

夕阳西下

金虫 (著名写手)

应助: 138 (高中生)
贵宾: 0.041
金币: 9247.1
帖子: 1268
在线: 1164.4小时
虫号: 220750

引用回帖:

2楼: Originally posted by snoopyzhao at 2012-02-05 23:25:40:
你这个是 fortran 的程序吗？
里面有很多变量根本就没有定义，比如 H1
ne 是数组还是普通的变量？
I、J作为 HX 和 NE 的维度的时候没有赋值……

你这个程序用任何一个编译器都编译不过去，根据编译器的提示 ...

主要是ne的问题，ne是个函数，就是根据编译器不知道ne该如何修改，改了很多次。

赞一下

回复此楼

3楼2012-02-06 10:01:51

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

夕阳西下

金虫 (著名写手)

应助: 138 (高中生)
贵宾: 0.041
金币: 9247.1
帖子: 1268
在线: 1164.4小时
虫号: 220750

引用回帖:

4楼: Originally posted by baobiao007 at 2012-02-06 16:53:00:
我感觉就程序编写来说，似乎有两个问题：
1.  不能用变量I 直接定义数组ne的大小。
应写成(按fortran90风格)
integer,parameter:: J
   real:: ne(J)
2.  Hx(0),Hx(I)那两行，你似乎把ne数组的下 ...

这里I是变量，ne是函数

赞一下

回复此楼

5楼2012-02-06 17:24:16

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