53 integer :: jpart,nage,ixave,jyave,kz,kzave,kp
54 integer :: k,k1,k2,ix,ix1,ix2,ixs,jy,jy1,jy2
55 real :: xold,yold,zold,xmean,ymean
61 if ((ioutputforeachrelease.eq.1).and.(mdomainfill.eq.0))
then
67 xmean=(xold+xtra1(jpart))/2.
68 ymean=(yold+ytra1(jpart))/2.
70 ixave=int((xmean*dx+xoutshift)/dxout)
71 jyave=int((ymean*dy+youtshift)/dyout)
73 if (outheight(kz).gt.ztra1(jpart)) goto 16
81 if ((ixave.ge.0).and.(jyave.ge.0).and.(ixave.le.numxgrid-1).and. &
82 (jyave.le.numygrid-1))
then
84 if (outheighthalf(kz).gt.zold) goto 11
86 11 k1=min(numzgrid,kz)
88 if (outheighthalf(kz).gt.ztra1(jpart)) goto 21
90 21 k2=min(numzgrid,kz)
94 flux(5,ixave,jyave,kz,k,kp,nage)= &
95 flux(5,ixave,jyave,kz,k,kp,nage)+ &
99 flux(6,ixave,jyave,kz,k,kp,nage)= &
100 flux(6,ixave,jyave,kz,k,kp,nage)+ &
110 if ((kzave.le.numzgrid).and.(jyave.ge.0).and. &
111 (jyave.le.numygrid-1))
then
115 if (abs(xold-xtra1(jpart)).lt.
real(nx)/2.) then
116 ix1=int((xold*dx+xoutshift)/dxout+0.5)
117 ix2=int((xtra1(jpart)*dx+xoutshift)/dxout+0.5)
120 if ((ix.ge.0).and.(ix.le.numxgrid-1))
then
121 flux(1,ix,jyave,kzave,k,kp,nage)= &
122 flux(1,ix,jyave,kzave,k,kp,nage) &
127 if ((ix.ge.0).and.(ix.le.numxgrid-1))
then
128 flux(2,ix,jyave,kzave,k,kp,nage)= &
129 flux(2,ix,jyave,kzave,k,kp,nage) &
140 ixs=int(((
real(nxmin1)-1.e5)*dx+xoutshift)/dxout)
141 if ((ixs.ge.0).and.(ixs.le.numxgrid-1))
then
142 if (xold.gt.xtra1(jpart))
then
144 flux(1,ixs,jyave,kzave,k,kp,nage)= &
145 flux(1,ixs,jyave,kzave,k,kp,nage) &
150 flux(2,ixs,jyave,kzave,k,kp,nage)= &
151 flux(2,ixs,jyave,kzave,k,kp,nage) &
163 if ((kzave.le.numzgrid).and.(ixave.ge.0).and. &
164 (ixave.le.numxgrid-1))
then
165 jy1=int((yold*dy+youtshift)/dyout+0.5)
166 jy2=int((ytra1(jpart)*dy+youtshift)/dyout+0.5)
170 if ((jy.ge.0).and.(jy.le.numygrid-1))
then
171 flux(3,ixave,jy,kzave,k,kp,nage)= &
172 flux(3,ixave,jy,kzave,k,kp,nage) &
177 if ((jy.ge.0).and.(jy.le.numygrid-1))
then
178 flux(4,ixave,jy,kzave,k,kp,nage)= &
179 flux(4,ixave,jy,kzave,k,kp,nage) &
subroutine calcfluxes(nage, jpart, xold, yold, zold)