24小时热门版块排行榜

返回列表

xingxiao

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 101.5
帖子: 110
在线: 62.9小时
虫号: 712761
注册: 2009-03-02
专业: 理论和计算化学

[求助] 已算出整个空间的流，下一步想得到在空间中某一位置的流

我已经实现在整个放电空间的粒子流密度，接下来想要得到在电极极板处了流密度，请问怎样实现？？？
下面是我编译的整个放电空间的粒子流密度：
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是一个数组，请问我该怎么办？？？

回复此楼