 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1112 | 回復(fù): 4 | ||
[求助]
DLA模型，F(xiàn)ortran 代碼，出現(xiàn)問題
|
|
一個(gè)DLA模型的Fortran代碼，能編譯，能運(yùn)行，但是出來的結(jié)構(gòu)不對(duì)。畫出的圖形除了中間有點(diǎn)外，四個(gè)角上也有點(diǎn)存在，不知道哪里有問題。希望各位蟲友看看，不勝感激！ 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結(jié)果 |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
wup123» 猜你喜歡
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有6人回復(fù)
-
318求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過
已經(jīng)有6人回復(fù)
-
一志愿武理材料工程348求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有6人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有7人回復(fù)
-
303求調(diào)劑
已經(jīng)有6人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Fortran 編譯問題
已經(jīng)有9人回復(fù)
- Fortran的格式化輸入輸出問題
已經(jīng)有14人回復(fù)
- Fortran 95／2003科學(xué)計(jì)算與工程 對(duì)應(yīng)的各章程序源代碼
已經(jīng)有16人回復(fù)
- 有關(guān)fortran的一次而問題，希望大家能幫幫忙，謝謝
已經(jīng)有4人回復(fù)
- 【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧
已經(jīng)有105人回復(fù)
- 【求助】intel fortran for linux的問題
已經(jīng)有4人回復(fù)
linnature
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2.2
- 帖子: 4
- 在線: 3.9小時(shí)
- 蟲號(hào): 1959201
- 注冊(cè): 2012-08-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2012-08-29 09:54:11
3樓2012-09-01 20:06:39
1234liang
木蟲 (正式寫手)
- 應(yīng)助: 70 (初中生)
- 金幣: 1661.4
- 紅花: 8
- 帖子: 629
- 在線: 97.1小時(shí)
- 蟲號(hào): 1538562
- 注冊(cè): 2011-12-15
- 性別: GG
- 專業(yè): 金屬功能材料
4樓2012-10-17 16:18:43
5樓2012-10-17 20:39:04
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wup123 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|