看板 Fortran 關於我們 聯絡資訊
截取某一段~ m = 41 n = 41 mp = m-1 np = n-1 ………………… do i = 1, mp do j = 1, np k = j + (i-1)*np ->推測k值可能為粒子編號? x(1, k) = (i-1)*dx + dx/2. ->x值 x(2, k) = (j-1)*dy + dy/2. ->y值 enddo enddo │ ──┼ x值對y值畫圖這部份可以完整畫出1600 do a= 1, mp-1 do b= a+1 , np k = np*(a-1) + (a-1)*a/2 +b -a + 1600 x(1, k) = (a)*dx + dx/2+0.02 ->x值(位移) x(2, k) = (b)*dy + dy/2. ->y值 enddo enddo │ ───┼ 因為想加上排列成倒三角型的粒子 所以就自已寫了一段直接加上去 程式跑出來,總是出現很多0值的點,都是我沒寫到的另一邊三角型的點 不知道大家知不知道我意思囧" 要如何把0值移除掉呢? 這樣粒子都會疊在原點~"~ 請各位給點建議吧!!謝謝>< 不知道能不能PO在這,不能我會自刪:) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.249.162.25 ※ 編輯: yzsungxgcf 來自: 111.249.162.25 (12/15 01:23)
Semisphere:k不要用公式算改用k=k+1,另確保x第2維夠大即可 12/15 09:01
yzsungxgcf:謝謝你的建議,我會試試 12/16 10:07