Date Posted:09:59:04 05/26/04 Wed Author: Reza Khaksar Fard Subject: IF YOU HAVE PROBLEM WITH DIFFERENTIAL MESH...........(4)
c.....
c..... Subroutine CheckConverg
c.....
subroutine CheckConverg(Converged)
include 'param.inc'
include 'common.inc'
logical Converged
if(Res.lt.Error) Converged = .true.
return
end
c.....
c..... Subroutine PrintResults
c.....
subroutine PrintResults(printcode)
include 'param.inc'
include 'common.inc'
if (printcode==1) then
c..... Print Airfoil
write(1,*) 'title =" Airfoil Plot"'
write(1,*) 'variables =x,y'
write(1,*) 'Zone I=',Imax+1,' J=1 F=Point'
do I = 0,Imax
write(1,*) X(I,0),Y(I,0)
end do
return
Else if (printcode==2) then
c..... Print Initial
write(2,*) 'title =" Initial Plot"'
write(2,*) 'variables =x,y'
write(2,*) 'Zone I=',Imax+1,' J=',Jmax+1,' F=Point'
do J = 0,Jmax
do I = 0,Imax
write(2,*) X(I,J),Y(I,J)
end do
end do
return
Else if (printcode==3) then
c..... Print Final Solution
c..... Print Initial
write(3,*) 'title =" Initial Plot"'
write(3,*) 'variables =x,y'
write(3,*) 'Zone I=',Imax+1,' J=',Jmax+1,' F=Point'
do J = 0,Jmax
do I = 0,Imax
write(3,*) X(I,J),Y(I,J)
end do
end do
return
end if
end
