应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 程序运行不了问题在哪
21/1返回列表
查看: 1449  |  回复: 11
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 夕阳西下 的 4 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

夕阳西下

金虫 (著名写手)

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

下面这个程序有错误，怎么改呢？
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
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2012-02-05 22:36:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)
★
夕阳西下(金币+1):谢谢参与
余泽成(金币+2): 谢谢参与应助！ 2012-02-21 00:27:51
我感觉就程序编写来说，似乎有两个问题：
1.  不能用变量I 直接定义数组ne的大小。
    应写成(按fortran90风格)
    integer,parameter:: J
      real:: ne(J)
2.  Hx(0),Hx(I)那两行，你似乎把ne数组的下标丢了。
赞一下(1人) 回复此楼
4楼2012-02-06 16:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夕阳西下 的主题更新
21/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭