Page 344 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
0 521 85326 5
CB908/Date
0521853265appc
323
APPENDIX C. 2D CARTESIAN CODE
WRITE(12,*)’VARIABLES = XX YY UU VV TT ’ May 25, 2005 11:59
WRITE(12,*)’ZONE T = ZONE1, I = ’,ITOT,’, J = ’,JN,’ ,F = BLOCK’
DO 111 J=1,JN
111 WRITE(12,*)(X(I),I=IBEG,IEND)
DO 112 J=1,JN
112 WRITE(12,*)(Y(J),I=IBEG,IEND)
DO 113 J=1,JN
113 WRITE(12,*)(U(I,J),I=IBEG,IEND)
DO 114 J=1,JN
114 WRITE(12,*)(V(I,J),I=IBEG,IEND)
DO 115 J=1,JN
115 WRITE(12,*)(T(I,J),I=IBEG,IEND)
CLOSE(12)
RETURN
END
C **************************************
SUBROUTINE ADSORB(NN)
INCLUDE ’COM2D.FOR’
C **************************************
N=NN
GO TO (10,20,30,40,50,60,70,80,90),N
C *** FOR PRESSURE CORRECTION
10 GO TO 1000
C *** FOR U-VEL
20 GO TO 1000
C *** FOR V-VEL
30 GO TO 1000
C *** FOR W-VEL
40 GO TO 1000
C *** FOR K. ENERGY
50 GO TO 1000
C *** FOR DISSIPATION
60 GO TO 1000
C *** FOR TEMPERATURE
70 GO TO 1000
C *** FOR FLUID PROPERTIES
80 GO TO 1000
C *** CALLED FROM BOUNDP
90 CONTINUE
1000 CONTINUE
RETURN
END