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’
   343   344   345   346   347   348   349   350   351   352   353