Page 368 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
347
APPENDIX C. 2D CARTESIAN CODE
DXW=X(IN)-X(INM) May 25, 2005 11:59
DXE=X(2)-X(1)
DYS=Y(JN)-Y(JNM)
DYN=Y(2)-Y(1)
C WEST
SUMW=FINTW(APU,I,J)
AW(I,J)=FINTW(RHO,I,J)*(R(J)*DYP(J))**2/SUMW*(1-LW)*LB
IF(NTAGW(I,J).EQ.17)THEN
JJ=J
IF(IPERIOD.EQ.1)JJ=JN-J+1
RHOW=(DXW*RHO(2,J)+DXE*RHO(INM,JJ))/(DXE+DXW)
SUMW=(DXW*APU(2,J)+DXE*APU(INM,JJ))/(DXE+DXW)
AW(I,J)=RHOW*(R(J)*DYP(J))**2/SUMW*2*LB
ENDIF
C EAST
SUME=FINTE(APU,I,J)
AE(I,J)=FINTE(RHO,I,J)*(R(J)*DYP(J))**2/SUME*(1-LE)*LB
IF(NTAGE(I,J).EQ.27)THEN
JJ=J
IF(IPERIOD.EQ.1)JJ=JN-J+1
RHOE=(DXW*RHO(2,JJ)+DXE*RHO(INM,J))/(DXE+DXW)
SUME=(DXW*APU(2,JJ)+DXE*APU(INM,J))/(DXE+DXW)
AE(I,J)=RHOE*(R(J)*DYP(J))**2/SUME*2.0*LB
ENDIF
C SOUTH
SUMS=FINTS(APV,I,J)
AS(I,J)=FINTS(RHO,I,J)*(RC(J)*DXP(I))**2/SUMS*(1-LS)*LB
IF(NTAGS(I,J).EQ.37)THEN
II=I
IF(JPERIOD.EQ.1)II=IN-I+1
RHOS=(DYS*RHO(I,2)+DYN*RHO(II,JNM))/(DYN+DYS)
SUMS=(DYS*APV(I,2)+DYN*APV(II,JNM))/(DYN+DYS)
AS(I,J)=RHOS*(RC(J)*DXP(I))**2/SUMS*2.0*LB
ENDIF
C NORTH
SUMN=FINTN(APV,I,J)
AN(I,J)=FINTN(RHO,I,J)*(RC(J+1)*DXP(I))**2/SUMN*(1-LN)*LB
IF(NTAGN(I,J).EQ.47)THEN
II=I
IF(JPERIOD.EQ.1)II=IN-I+1
RHON=(DYS*RHO(II,2)+DYN*RHO(I,JNM))/(DYN+DYS)
SUMN=(DYS*APV(II,2)+DYN*APV(I,JNM))/(DYN+DYS)