| 24小時(shí)熱門(mén)版塊排行榜 |
| 1 | 1/1 | 返回列表 |
| 查看: 1120 | 回復(fù): 0 | ||
解子林新蟲(chóng) (小有名氣)
|
[求助]
structure factor
|
|
哪位能幫忙分析分析我計(jì)算出來(lái)的結(jié)構(gòu)因子為什么在大q時(shí)不趨近于1?我調(diào)整過(guò)qx、qy、dq的值,可是沒(méi)什么作用,大伙給點(diǎn)意見(jiàn),謝謝! do i = 1 , conf_2 mm=1 do j = 1 , 9 read(1,*) enddo do j = 1 , all_n read(1,*) num_o , type_o , X_o , Y_o , Z_o if (type_o==1)then !!!!!!!!水分子的坐標(biāo) xx(mm) = X_o yy(mm) = Y_o zz(mm) = Z_o mm=mm+1 !!!!!!!水分子的數(shù)目 endif end do do nx=1,600 do ny=1,600 do nz=1,600 kr=sqrt((nx*qx)**2+(ny*qy)**2+(nz*qz)**2)!!!!!!!q的模 ! if ((kr>9).and.(kr<62))then bin=1+int(kr/0.5) nhist(bin)=nhist(bin)+1 cossum=0.0 sinsum=0.0 mm1=0 do l = 1 , mm-1 rx=0.1*xx(l) ry=0.1*yy(l) rz=0.1*zz(l) cossum=cossum+cos(nx*qx*rx+ny*qy*ry+nz*qz*rz) sinsum=sinsum+sin(nx*qx*rx+ny*qy*ry+nz*qz*rz) mm1=mm1+1 end do hist(bin) =hist(bin) + (cossum**2+sinsum**2)/real(mm-1) write(*,*) nhist(bin),hist(bin),mm1,mm,i !endif enddo enddo enddo end do !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! do k = 1 ,600 if(hist(k)/=0) then g(k) = real(hist(k))/real(nhist(k)) write(2,*) k*0.5,g(k) else g(k)=0 write(2,*) k*0.5,g(k) endif end do |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +11 | Faune 2026-03-21 | 11/550 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-24 | 3/150 |
|
|
[考研] 收08調(diào)劑生 +6 | komorebi69 2026-03-18 | 6/300 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +3 | gdlf9999 2026-03-24 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開(kāi)頭均可,過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|