1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
| do i = 2,n-1
j=1
l=(j-1)*n+i
l1=(j-1)*n+i+1
l2=(j-1)*n+i-1
l3=j*n+i
A(l,l)=(2.0D0*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/(2.0D0))
A(l,l1)=(-1.0D0)*r(j)*(dr/dz)
A(l,l2)=(-1.0D0)*r(j)*(dr/dz)
A(l,l3)=(-1.0D0)*(r(j)-(dr/2.0D0))*(dz/dr)
if (z(i) <= L_evap) then
vect(l,1)=(q_evap*(r(j)+(dr/2.0D0))*dz)/(Conds)
else if (z(i)>L_evap).and. (z(i)<== L_test)) then
vect(l,1)=0.0D0
else
A(l,l)=(2.0D0*r(j)*(dr/dz))+(r(j)*(dz/dr))-(dz/(2.0D0))+((Coef*(r(j)+(dr/2.0D0))*dz)/(Conds))
vect(l,1)= ((Coef*T_inf*(r(j)+(dr/2.0D0))*dz)/(Conds))
end if
! end if
end do |
Partager