| 查看: 1183 | 回复: 0 | ||
[求助]
structure factor
|
|
哪位能帮忙分析分析我计算出来的结构因子为什么在大q时不趋近于1？我调整过qx、qy、dq的值，可是没什么作用，大伙给点意见，谢谢！ 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 !!!!!!!!水分子的坐标 xx(mm) = X_o yy(mm) = Y_o zz(mm) = Z_o mm=mm+1 !!!!!!!水分子的数目 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 |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
-
能否申诉？
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
GSAS精修完之后怎么得到拟合峰强度或者积分峰面积？
已经有4人回复
- 请教中子散射高手分析
已经有7人回复
- 求助高手帮忙解析cif文件当中的A类C类错误以及精修参数是否合适等
已经有4人回复
- OLEX2解析晶体时check cif
已经有8人回复
- checkcif中出现的A类、B类问题怎么解
已经有4人回复
- 晶体化学（英文版）
已经有208人回复
- SCI生物类期刊 影响因子
已经有14人回复
- 求单斜和斜方S8的cif格式文件
已经有15人回复
- 求助cifcheck！！
已经有4人回复
- 【分享】elk-1.3.2
已经有18人回复
- 【求助】fullprof 精修，指标化-欢迎参与讨论
已经有20人回复
- 【求助】structure factor S(q)的问题
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30