IF YOU HAVE PROBLEM WITH DIFFERENTIAL MESH...........(4) -- CFD

[ Show ]

Support VoyForums

[ Shrink ]

VoyForums Announcement: Programming and providing support for this service has been a labor of love since 1997. We are one of the few services online who values our users' privacy, and have never sold your information. We have even fought hard to defend your privacy in legal cases; however, we've done it with almost no financial support -- paying out of pocket to continue providing the service. Due to the issues imposed on us by advertisers, we also stopped hosting most ads on the forums many years ago. We hope you appreciate our efforts.

Show your support by donating any amount. (Note: We are still technically a for-profit company, so your contribution is not tax-deductible.) PayPal Acct: Feedback:

Donate to VoyForums (PayPal):

[ Login ] [ Contact Forum Admin ] [ Main index ] [ Post a new message ] [ Search | Check update time ]

[ Next Thread | Previous Thread | Next Message | Previous Message ]

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

[ Next Thread | Previous Thread | Next Message | Previous Message ]