作者yzsungxgcf (=。=)
看板Fortran
标题[问题] 改写的问题(新手)
时间Thu Dec 15 01:19:57 2011
截取某一段~
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)
1F:→ Semisphere:k不要用公式算改用k=k+1,另确保x第2维够大即可 12/15 09:01
2F:→ yzsungxgcf:谢谢你的建议,我会试试 12/16 10:07