| 查看: 1044 | 回复: 4 | |||
[求助]
DLA模型，Fortran 代码，出现问题
|
|
一个DLA模型的Fortran代码，能编译，能运行，但是出来的结构不对。画出的图形除了中间有点外，四个角上也有点存在，不知道哪里有问题。希望各位虫友看看，不胜感激！ program dla implicit none parameter M=100 parameter N=100 parameter Particles=2000 integer :: A(M,N) integer :: i, j, counter, p real(4):: s, num1, num2, num3 open (30,file="dla2d.dat"  do i =1, M do j=1, N A(i,j)=0 enddo enddo A(50,50)=1 counter=0 do while (counter .LT. Particles) counter=counter+1 write (*,*) 'Current Particles is ', counter call random_number(num1) i=1+int(num1*(N-1)) call random_number(num2) j=1+int(num2*(N-1)) p=0 55 do while (p.EQ.0) call random_number(num3) s = num3 if(s.LE.0.25) then j=j+1 elseif(s.LE.0.50) then j=j-1 elseif(s.LE. 0.75) then i=i+1 else i=i-1 endif if (j.EQ.N+1) then j=1 goto 55 endif if (j.EQ.0) then j=N goto 55 endif if (i .EQ. M+1) then i=1 goto 55 endif if (i.EQ.0) then i=M goto 55 endif if(A(i,j-1).EQ.1.OR.A(i,j+1).EQ.1.OR.A(i+1,j).EQ.1.OR.A(i-1,j).EQ.1) then A(i,j)=1 write(30,*) i, j p=10 endif enddo enddo end program dla  DLA结果 |
回复此楼» 本帖已获得的红花（最新10朵）
wup123» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
Fortran 编译问题
已经有9人回复- Fortran的格式化输入输出问题
已经有14人回复
- Fortran 95／2003科学计算与工程 对应的各章程序源代码
已经有16人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【求助】intel fortran for linux的问题
已经有4人回复
2楼2012-08-29 09:54:11
3楼2012-09-01 20:06:39
1234liang
木虫 (正式写手)
- 应助: 70 (初中生)
- 金币: 1661.4
- 红花: 8
- 帖子: 629
- 在线: 97.1小时
- 虫号: 1538562
- 注册: 2011-12-15
- 性别: GG
- 专业: 金属功能材料
4楼2012-10-17 16:18:43
5楼2012-10-17 20:39:04