program tawon
c***********************************************************************
c algebraic grid generation code for elliptic circle mesh
c programmed by prof. m.s. kim on august 24, 1999
c***********************************************************************
parameter (idim=300,jdim=150)
common /mesh/ x(idim,jdim),y(idim,jdim)
common /input/ fnx,fny,radi1,strch,arad,brad
open (unit=1,file='tawon.in')
open (unit=2,file='tawon.echo')
open (unit=3,file='tawon.grd')
open (unit=4,file='tawon.plt')
call datain
call dataecho
call meshgen
call gridout
stop
end
c-----------------------------------------------------------------------
subroutine datain
parameter (idim=300,jdim=150)
common /mesh/ x(idim,jdim),y(idim,jdim)
common /input/ fnx,fny,radi1,strch,arad,brad
read(1,100)
100 format(1x)
read(1,200) fnx,fny,radi1,strch,arad,brad
200 format(6f10.2)
return
end
c-----------------------------------------------------------------------
subroutine dataecho
parameter (idim=300,jdim=150)
common /mesh/ x(idim,jdim),y(idim,jdim)
common /input/ fnx,fny,radi1,strch,arad,brad
write(2,100) fnx,fny,radi1,arad,brad,strch
100 format(5x,'fnx =',f10.2/
1 5x,'fny =',f10.2/
2 5x,'rad1 =',f10.2/
3 5x,'arad =',f10.2/
4 5x,'brad =',f10.2/
5 5x,'stretch=',f10.2)
return
end
c-----------------------------------------------------------------------
subroutine meshgen
parameter (idim=300,jdim=150)
common /mesh/ x(idim,jdim),y(idim,jdim)
common /input/ fnx,fny,radi1,strch,arad,brad
pai = 4.*atan(1.)
theta = (pai/2.)/(fnx-1.0)
imax = fnx
jmax = fny
c.... discretize the grid points on the inner and outer boundary
do 100 i=1,imax
x(i,1) = radi1*cos(pai -theta*(i-1))
y(i,1) = radi1*sin(pai -theta*(i-1))
aa = cos(pai -theta*(i-1)) / arad
aa = aa*aa
bb = sin(pai -theta*(i-1)) / brad
bb = bb*bb
radi2 = 1./sqrt(aa +bb)
x(i,jmax) = radi2*cos(pai -theta*(i-1))
y(i,jmax) = radi2*sin(pai -theta*(i-1))
100 continue
c.... fill-in the inside grid points
do 200 j=2,jmax-1
do 200 i=1,imax-1
xleng = x(i,jmax) - x(i,1)
yleng = y(i,jmax) - y(i,1)
coeff1 = 1.0 -strch**(j-1)
coeff2 = 1.0 -strch**(jmax-1)
x(i,j) = x(i,1) + xleng*coeff1/coeff2
y(i,j) = y(i,1) + yleng*coeff1/coeff2
200 continue
c.... setting the imax boundary grid points (set y=0)
do 300 j=1,jmax
yleng = y(imax,jmax) - y(imax,1)
coeff1 = 1.0 -strch**(j-1)
coeff2 = 1.0 -strch**(jmax-1)
x(imax,j) = x(imax,1)
y(imax,j) = y(imax,1) + yleng*coeff1/coeff2
300 continue
return
end
c-----------------------------------------------------------------------
subroutine gridout
parameter (idim=300,jdim=150)
common /mesh/ x(idim,jdim),y(idim,jdim)
common /input/ fnx,fny,radi1,strch,arad,brad
imax = fnx
jmax = fny
c.... to calculate cl,cd, cm coefficients (airfoil geometry)
scal = 1.0
chord= 2.0*radi1
xm = 0.
ym = 0.
xle = 0.
xitl = 0.
xitu = 0.
write(3,100) fnx,fny,scal,chord,xm,ym,xle,xitl,xitu
write(3,100) ((x(i,j),i=1,imax),j=1,jmax)
write(3,100) ((y(i,j),i=1,imax),j=1,jmax)
100 format(10e18.11)
c.... prepare grid file for tecplot
write(4,140) imax,jmax
140 format(4x,'title = "mesh shape"'/
1 4x,'variables = "x" "y" '/
2 4x,'zone i=',i3,',',1x,'j=',i3,',',1x,'f=point')
write(4,150) ((x(i,j), y(i,j), i=1,imax), j=1,jmax)
150 format(2e13.5)
return
end
에러메시지
--------------------Configuration: tawon - Win32 Debug--------------------
Compiling Fortran...
C:\자료\tawon.for
f90: Severe: No such file or directory
... file is 'C:\AU·a\tawon.for'
Error executing df.exe.
tawon.exe - 1 error(s), 0 warning(s)
이렇게 나오는데..
제가 짠건아니고, 프로그램을 돌려보라고 해서 실행할려고 하니까 안되네요
그래프도 나온다고 하는데....그건 잘 모르겠고요....아시는 분은 점검좀 부탁합니다.
타원격자 모형입니다....
카페 게시글
포트란 Q&A
에러 메시지좀 봐주세요...
북마크
번역하기
공유하기
기능 더보기
다음검색
첫댓글 한글폴더이름 때문입니다.
현성님 빙고~ ^^