Page 349 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
328
APPENDIX C. 2D CARTESIAN CODE
WRITE(24,*)’VARIABLES = XX YY UU VV OO ’ May 25, 2005 11:59
WRITE(24,*)’ZONE T = ZONE1, I = ’,IEND
1 ,’,J= ’,JEND,’ ,F = BLOCK’
DO 11 J=1,JEND
11 WRITE(24,*)(X(I),I=1,IEND)
DO 12 J=1,JEND
12 WRITE(24,*)(Y(J),I=1,IEND)
DO 13 J=1,JEND
13 WRITE(24,*)(U(I,J),I=1,IEND)
DO 14 J=1,JEND
14 WRITE(24,*)(V(I,J),I=1,IEND)
DO 15 J=1,JEND
15 WRITE(24,*)(O(I,J),I=1,IEND)
CLOSE(24)
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 DO 31 J=2,JNM
DO 31 I=2,INM
31 SU(I,J)=SU(I,J)+GRM*O(I,J)*VOL(I,J)*(1-NTAG(I,J))
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 FORM BOUNDP - FOR PRESSURE