Page 338 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
0 521 85326 5
CB908/Date
0521853265appc
317
APPENDIX C. 2D CARTESIAN CODE
CLOSE(15) May 25, 2005 11:59
600 FORMAT(2X,7E10.3)
500 FORMAT(I4,2X,7E10.3)
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 - PFIX CONDITION
10 CONTINUE
SP(2,2)=SP(2,2)+RHO(1,2)*(R(2)*DYP(2))**2/APU(2,2)
SP(INM,2)=SP(INM,2)+RHO(IN,2)*(R(2)*DYP(2))**2/APU(INM,2)
GO TO 1000
C *** FOR U-VEL
20 CONTINUE
DO 21 I=2,INM
TERM=U(I+1,2)/DXMI(I+1)*DYP(2)+U(I-1,2)/DXMI(I)*DYP(2)
SU(I,2)=SU(I,2)+TERM*VISCOS
SP(I,2)=SP(I,2)+VISCOS*DYP(2)*(1./DXMI(I+1)+1./DXMI(I))
21 SP(I,2)=SP(I,2)+VISCOS*RESIST*EPSI*VOL(I,2)
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 P(1,2)=PIN
P(IN,2)=0
PP(IN,2)=0
PP(1,2)=0
C DUMMY VARABLES FOR UU,PRES,PPP,PPS FOR PRINTING