Page 337 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
316
APPENDIX C. 2D CARTESIAN CODE
C ************************************** May 25, 2005 11:59
SUBROUTINE RESULT
INCLUDE ’COM2D.FOR’
C **************************************
CHARACTER*20 HEADER
JSTEP=-1
WRITE(6,*)’ NITER=’,NITER
HEADER=’ U-VEL’
CALL PRINTK(U,1,IN,2,JNM,HEADER,JSTEP)
HEADER=’ PRESS’
CALL PRINTK(P,1,IN,2,JNM,HEADER,JSTEP)
OPEN(12,FILE=’PORU.DAT’)
WRITE(12,*)’IN = ’,IN
WRITE(12,*)’ X(I) = ’
WRITE(12,500)(X(I),I=1,IN)
WRITE(12,*)’ U- VELOCITY ’
DO 11 NN=1,NITER-1
11 WRITE(12,500)NN,(UU(I,NN),I=1,IN)
CLOSE(12)
OPEN(13,FILE=’PORP.DAT’)
WRITE(13,*)’IN = ’,IN
WRITE(13,*)’ X(I) = ’
WRITE(13,500)(X(I),I=1,IN)
WRITE(13,*)’ PRESSURE’
DO 12 NN=1,NITER-1
12 WRITE(13,500)NN,(PRES(I,NN),I=1,IN)
CLOSE(13)
OPEN(14,FILE=’PORPP.DAT’)
WRITE(14,*)’IN = ’,IN
WRITE(14,*)’ X(I) = ’
WRITE(14,500)(X(I),I=1,IN)
WRITE(14,*)’ PPM’
DO 122 NN=1,NITER-1
122 WRITE(14,500)NN,(PPP(I,NN),I=1,IN)
CLOSE(14)
OPEN(15,FILE=’PORPS.DAT’)
WRITE(15,*)’IN = ’,IN
WRITE(15,*)’ X(I) = ’
WRITE(15,500)(X(I),I=1,IN)
WRITE(15,*)’ PPS’
DO 222 NN=1,NITER-1
222 WRITE(15,500)NN,(PPS(I,NN),I=1,IN)