Page 348 - Introduction to Computational Fluid Dynamics
P. 348
P2: IWV
P1: ICD/GKJ
CB908/Date
0521853265appc
327
APPENDIX C. 2D CARTESIAN CODE
RETURN 0 521 85326 5 May 25, 2005 11:59
END
C **************************************
SUBROUTINE RESULT
INCLUDE ’COM2D.FOR’
C **************************************
CHARACTER*20 HEADER
JSTEP=-1
WRITE(6,*)’ NITER=’,NITER
HEADER=’ U-VEL’
CALL PRINTK(U,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ V-VEL’
CALL PRINTK(V,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ PRESS’
CALL PRINTK(P,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ OMEGA’
CALL PRINTK(O,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ AMW’
CALL PRINTK(AMW,1,IN,1,JN,HEADER,JSTEP)
C CALCULATE NORALISED EVAPORATION RATE
SUMWAT=0.0
DX=0.0
DO 1 I=2,IB1-1
DX=DX+DXP(I)
1 SUMWAT=SUMWAT+AMW(I,1)*DXP(I)
VBAR=SUMWAT/DX
SUMBR=0.0
DO 2 I=IB2,INM
2 SUMBR=SUMBR+AMW(I,1)*DXP(I)
C DIFFUSION LIMIT
B=(0-1)/(1-OWT)
DL=VISCOS/SC/YC(JB1)*ALOG(1+B)
WRITE(6,*)’ DIFFUSION LIMIT = ’,DL
WRITE(6,*)’ ACTUAL FLUX = ’,VBAR
RR=VBAR/DL
WRITE(6,*)’ RATIO = ’,RR,’ GRM = ’,GRM,’B=’,B
WRITE(6,*)’ SUMWAT = ’,SUMWAT,’ SUMBR = ’,SUMBR
WRITE(6,*)’ OWT = OBR = ’,OWT
IEND=30
JEND=33
OPEN(24,FILE=’EVAP.DAT’)
WRITE(24,*)’TITLE = EVAPORATION’