| 查看: 1323 | 回复: 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» 猜你喜欢
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Fortran 编译问题
已经有9人回复
- Fortran的格式化输入输出问题
已经有14人回复
- Fortran 95／2003科学计算与工程 对应的各章程序源代码
已经有16人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有106人回复
- 【求助】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
