| 查看: 1117 | 回复: 4 | ||
wup123版主
|
[求助]
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» 猜你喜欢
286求调剂
已经有7人回复
-
290求调剂
已经有5人回复
-
085600，320分求调剂
已经有6人回复
-
292分，材料与化工，申请调剂
已经有15人回复
-
0710生物学求调剂
已经有4人回复
-
一志愿北京科技，085601总分305求调剂
已经有5人回复
-
材料求调剂
已经有4人回复
-
一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业
已经有9人回复
-
材料专硕322分
已经有6人回复
-
070300化学求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Fortran 编译问题
已经有9人回复
- Fortran的格式化输入输出问题
已经有14人回复
- Fortran 95／2003科学计算与工程 对应的各章程序源代码
已经有16人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【求助】intel fortran for linux的问题
已经有4人回复
2楼2012-08-29 09:54:11
wup123
版主
- 应助: 0 (幼儿园)
- 金币: 313.1
- 散金: 27
- 红花: 2
- 帖子: 93
- 在线: 62.5小时
- 虫号: 728617
- 注册: 2009-03-22
- 性别: MM
- 专业: 理论和计算化学
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
wup123
管理员
- 应助: 0 (幼儿园)
- 金币: 313.1
- 散金: 27
- 红花: 2
- 帖子: 93
- 在线: 62.5小时
- 虫号: 728617
- 注册: 2009-03-22
- 性别: MM
- 专业: 理论和计算化学
5楼2012-10-17 20:39:04
