1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
| DIMENSION XU(20,20,20),YV(20,20,20),ZW(20,20,20)
pi=3.1416
ni=10; nj=5; nk=5
OPEN( 6, FILE='aa.DAT')
xl=2.*pi; yl=0.5; zl=2.; zpwr=1.
do 100 k=1,nk+1
do 200 j=1,nj+1
do 300 i=1,ni+1
yr=yl*(j-1)/nj
xx=xl*(i-1)/ni
xu(i,j,k)=yr*cos(xx)
yv(i,j,k)=xx*sin(xx)
zw(i,j,k)=zl*((k-1)/nk)**zpwr
300 continue
200 continue
100 continue
do 110 i=1,ni+1
110 write(*,*)xu(i,j,k)
do 120 j=1,nj+1
120 write(*,*) yv(i,j,k)
do 130 k=1,nk+1
130 write(*,*)zw(i,j,k)
WRITE(40,405) xu(i,j,k),yv(i,j,k),zk(i,j,k)
405 FORMAT(3(3X,1E15.7))
stop
end |
Partager