| 查看: 1265 | 回复: 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» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
