| 查看: 570 | 回复: 4 | ||
牧雪长空铜虫 (正式写手)
|
[求助]
菜鸟紧急求助！！
|
|
正在编一个求均方差位移的城西，各种错误百出。总是运行的时候提示数组越界：array bounds exceed。也不知道哪里出错了。下面是代码： program main implicit none integer i,j,k real*8,dimension(5000,160) :: xb,yb,zb,r,s real*8,dimension(160)::t real*8,dimension(4999)::r2,ms real*8 dx,dy,dz real*8,parameter :: npart=160,nstep=5000 real*8,parameter :: a=9.9222,b=15.9344,c=11.4808,dt=1 open (8,file='xdatcar') do i=1,5000 do j=1,160 read(8,*)xb(i,j),yb(i,j),zb(i,j) xb(i,j)=xb(i,j)*a yb(i,j)=yb(i,j)*b zb(i,j)=zb(i,j)*c end do end do ! to abtain the real coorditate of these atoms, and store them in 2d arrays. do i=1,160 ms=0 do k=1,4999 r2=0 do j=1,k,4999 dx=xb(i+160*(k+j-1),i)-xb(i+160*(j-1),i) dx=dx-a*aint(dx/a) dy=yb(i+160*(k+j-1),i)-yb(i+160*(j-1),i) dx=dy-b*aint(dy/b) dz=zb(i+160*(k+j-1),i)-zb(i+160*(j-1),i) dz=dz-c*aint(dz/c) r2=dx**2+dy**2+dz**2+r2 end do ms=r2/(5000-k) open (9,file='all msd.txt') write(9,*)ms end do ! end up the k loop end do end |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 紧急求助高斯扫势能曲线（同分子中不同基团之间的距离如何定义）
已经有4人回复
- 紧急求助啊 万分火急 我该如何选择？？？？？？
已经有31人回复
- 紧急求助！！！！怎样用matlab编动态断裂应力强度因子的计算程序
已经有7人回复
- 菜鸟紧急求助——怎么在LPSN上找模式菌株
已经有5人回复
- 紧急求助！！！！有关UNSW的免学费申请截止日期及UNSWID，请各位帮忙！谢谢啦~
已经有8人回复
- 求助 nanotechnology 的copyright 什么阶段上传
已经有5人回复
- 【求助】乳液合成（毕业设计），紧急求助！！！
已经有41人回复
- 关于菜鸟生平第一篇文章的紧急求救！！
已经有14人回复
- 【求助】一名菜鸟，紧急求助
已经有13人回复
- [紧急求助]Matlab中的图像拼接
已经有8人回复
- 【讨论】纳米线 纳米棒 纳米晶须 都是怎么区别的？
已经有16人回复
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
2楼2013-02-26 10:04:35
牧雪长空
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1003.8
- 散金: 59
- 帖子: 373
- 在线: 212.7小时
- 虫号: 1465005
- 注册: 2011-10-28
- 专业: 实验地球化学和计算地球化
3楼2013-02-27 09:40:39
xujince
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2516.7
- 红花: 1
- 帖子: 548
- 在线: 376.9小时
- 虫号: 1926989
- 注册: 2012-08-06
- 性别: GG
- 专业: 固体力学
4楼2013-02-27 11:52:50
weiwei43
至尊木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 18316.1
- 红花: 6
- 帖子: 737
- 在线: 94.7小时
- 虫号: 250282
- 注册: 2006-05-13
- 专业: 机械动力学
5楼2013-02-28 10:45:37