- #1
chammakh24
- 2
- 0
please i need help in this code
program wide_boiling
implicit none
integer :: i, j, m
real :: Trosee, Tbulle, P, RR(100), vf(100),
real :: aT1(3), aT6(3), ap1(3), z(3), ap2(3)
real ::, x(3), T, k(3,10)
m=0
T=518.67
m=m+1
do i = 1,3
k(i,m) = (aT1(i)/(T(m)**2))+(aP1(i)*log(P))+aT6(i)+(aP2(i)/P**2)
k(i,m) = exp(k(i,m))
write(6,*) ¢k(i,m) = ¢,¢ ¢,k(i,m),¢ ¢,¢ i = ¢,¢ ¢,i,' ¢,¢m = ¢,¢ ¢,m
end do
100 continue
j=0
vf =0.25
j=j+1
RR(j) = 0
do i = 1,3
RR(j) = RR(j)+(k(i,j)-1)*z(i)/(1+(k(i,j)-1)*vf(j))
end do
write(6,*) ?RR(j)=?,RR(j),? ?,?j=?,j
if(abs(RR(j)).LE.5E-3) then
do i = 1,3
x(i) = z(i)/(1+(k(i,j)-1)*vf(j))
end do
write(6,*) ?x(i)= ?,x(i),? ?,?i = ?,i,? ?,? j = ?,j
else
dRR(j) = 0
do i = 1,3
dRR(j)=dRR(j)+(z(i)*((k(i,m)-1)**2))/((z(i)*((k(i,m)-1)**2))/(1+(k(i,m)*vf(j)))**2)
end do
write(6,*) ?dRR(j)=?,dRR(j),? ?,?j=?,j
vf(j)=vf(j)-(RR(j)/dRR(j))
write(6,*) ?vf(j)=?,vf(j),? ?,?j=?,j
go to 100
end if
end program wide_boiling
program wide_boiling
implicit none
integer :: i, j, m
real :: Trosee, Tbulle, P, RR(100), vf(100),
real :: aT1(3), aT6(3), ap1(3), z(3), ap2(3)
real ::, x(3), T, k(3,10)
m=0
T=518.67
m=m+1
do i = 1,3
k(i,m) = (aT1(i)/(T(m)**2))+(aP1(i)*log(P))+aT6(i)+(aP2(i)/P**2)
k(i,m) = exp(k(i,m))
write(6,*) ¢k(i,m) = ¢,¢ ¢,k(i,m),¢ ¢,¢ i = ¢,¢ ¢,i,' ¢,¢m = ¢,¢ ¢,m
end do
100 continue
j=0
vf =0.25
j=j+1
RR(j) = 0
do i = 1,3
RR(j) = RR(j)+(k(i,j)-1)*z(i)/(1+(k(i,j)-1)*vf(j))
end do
write(6,*) ?RR(j)=?,RR(j),? ?,?j=?,j
if(abs(RR(j)).LE.5E-3) then
do i = 1,3
x(i) = z(i)/(1+(k(i,j)-1)*vf(j))
end do
write(6,*) ?x(i)= ?,x(i),? ?,?i = ?,i,? ?,? j = ?,j
else
dRR(j) = 0
do i = 1,3
dRR(j)=dRR(j)+(z(i)*((k(i,m)-1)**2))/((z(i)*((k(i,m)-1)**2))/(1+(k(i,m)*vf(j)))**2)
end do
write(6,*) ?dRR(j)=?,dRR(j),? ?,?j=?,j
vf(j)=vf(j)-(RR(j)/dRR(j))
write(6,*) ?vf(j)=?,vf(j),? ?,?j=?,j
go to 100
end if
end program wide_boiling