| 查看: 370 | 回复: 1 | |||
[求助]
已算出整个空间的流，下一步想得到在空间中某一位置的流
|
|
我已经实现在整个放电空间的粒子流密度，接下来想要得到在电极极板处了流密度，请问怎样实现？？？ 下面是我编译的整个放电空间的粒子流密度： implicit none integer,intent(in)::Ng real,intent(in)::Sbc,gone,q integer,dimension(1:Ng),intent(in)::ink real,dimension(1:Ng,1:3000,1:7),intent(in)::AN real,dimension(0:Ng),intent(out)::Jf integer::i,j,m real::Vq,qV,qv1,qv2 real,parameter::br=0.064,pi=3.14159 do i=0,Ng 整个空间包含的网格数 Jf(i)=0.0 整个空间的粒子流 end do do i=1,Ng do j=1,ink(i) 在某一个网格里德粒子数目 if(q<0.0) then 如果是负电荷 qV=abs(AN(i,j,7))*AN(i,j,4) 其中AN(i,j,7)是粒子的电荷量， AN(i,j,4)是离子在轴向的速度 else if(q>0.0) then 如果是正电荷 qV=-abs(AN(i,j,7))*AN(i,j,4) end if call particleweighting(AN(i,j,1),Sbc,qV,m,qv1,qv2) 粒子加权到网格点上 if(m==Ng) then 如果在最后一个网格，即右极板处 Jf(m)=Jf(m)+qv1 else Jf(m)=Jf(m)+qv1 Jf(m+1)=Jf(m+1)+qv2 end if end if end do end do Vq=abs(q)*pi*br**2*Sbc 每个网格的体积 Jf(0)=Jf(0)/(0.5*Vq);Jf(Ng)=Jf(Ng)/(0.5*Vq) do i=1,Ng-1 Jf(i)=Jf(i)/Vq end do 现在我想实现粒子在左极板处的流密度，它应该是一个数，而不是数组，但是上述算出来的Jf是一个数组，请问我该怎么办？？？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用什么方法，能输出某一个空间群的对应对称操作矩阵？
已经有7人回复
- 亲和取代反应-空间位阻
已经有13人回复
- 想得到粒径分布函数是用电镜统计法还是粒度分析仪做
已经有9人回复
- 突然停电或者死机之后，想要继续上次的计算的解决方案！
已经有7人回复
2楼2012-03-07 21:41:17
10