home *** CD-ROM | disk | FTP | other *** search
/ The World of Computer Software / World_Of_Computer_Software-02-387-Vol-3of3.iso / h / house_ii.zip / FOR / PLENUM.FOR < prev    next >
Text File  |  1992-05-11  |  4KB  |  113 lines
  1.       SUBROUTINE PLENUM 
  2. C  PLENUM MODEL FOR OUTLET FROM FURNACE OR AIR HANDLER TO ENTRY TO DUCTS
  3. C    EXTRACTED FROM SUBPROGRAM GASF SO THAT GASF ISNT CALLED
  4. C    IN COOLING MODE
  5. CMDK BLK83
  6. CMDK BLK85
  7. CMDK BLK87
  8. CMDK BLK89
  9. CMDK BLKCM2
  10. CMDK BLKGS1
  11. CMDK BLKGS2
  12. CMDK BLKQGS
  13. CMDK BLKSK
  14. CMDK COLCT
  15. CMDK COLLEC
  16. CMDK DTRQ
  17. CMDK FANBLK
  18. CMDK INDIC1
  19. CMDK PRT1
  20. CMDK TIMEB
  21.       DIMENSION BCON(6),HIPLE(2),TFILM(3)
  22.       NAMELIST /INPGSF/  BCON, CPMACL, DTRJC, DTRPL, EFFHXV,
  23.      + EMISJC, EMISPL, FANTOF, FANTON, HAACL,
  24.      + HIPL2, HOJ, IDV, IELC, IVFA,
  25.      + LINWFL, NACL, NFLT, NHXITR, NINTGF,
  26.      + NTBB, NTHX, NTJJ, PCLOSE, PCTJAC,    
  27.      + POPEN, RCLOSE, RNSLS1, ROPEN, TCMOF1, 
  28.      + TCMOF2, TCYCOF, TCYCON, TDESH, TDESL,
  29.      + TFILM, TFLT, TFLUE, THLOFF,THLON,
  30.      + TIME1,TIME3,TIMOVD, TJSV, TMACSV, 
  31.      + TMFLT, TMPLSV, TMSV, TMV, TSHLIM,
  32.      + TVENTI, TWFMAX, UAHXOF, UAHXON, UAHXV,
  33.      + VOUTSD, WFL1, WFLMIN, WGTHX, ZAJ,
  34.      + ZKBND2, ZRELFI, ZRELFO
  35. C 
  36. C  DTRPL IS TEMP ELEVATION OF PLEN SURF ABOVE TAIN
  37. C  EMISPL IS EMISSIVITY OF PLEN SURFACE 
  38. C  THE ENTIRE NAMELIST INPGSF FROM GASF IS INCLUDED TO ENABLE 
  39. C  READ OF THIS NAMELIST FOR COOLING
  40.       DATA DTMTOL/.01/
  41.       DATA DTMPL/0./
  42.       DATA NTBB/100/
  43. C  IRGSF SHOULD BE 1 FOR COOLING(SET IN LOOP)
  44.       IF(IRGSF.EQ.1)THEN
  45. C      
  46.        OPEN(12,FILE='TAPE12',STATUS='OLD',IOSTAT=IO12)
  47.        IF(IO12.NE.0)THEN
  48.          WRITE(60,*)' GASF: CANT OPEN TAPE12 WITH NAMELIST INPGSF'
  49.          STOP ' GASF: CANT OPEN TAPE12 WITH NAMELIST INPGSF '
  50.          END IF
  51.        READ(12,INPGSF,END=999)
  52.        CLOSE(12)
  53.        WRITE(60,INPGSF)
  54.        ENDIF
  55.       C1=U45*BARPSI/(NCELLS*14.7) 
  56.       WAIRF(1)=CFM(1)*C1
  57.       WAIRF(2)=CFM(2)*C1
  58. C - - - - - - -  SINGLE-SPEED BLOWER SET BELOW
  59.       IBL=1 
  60.       INC=1 
  61.       CONS12= DTIME/(RHOMPL*CPMPL*APLEN)*NCELLS 
  62.       HOPL= HOPL1+EMISPL*.171E-8*4.*(TAIN+DTRPL*.5+459.7)**3
  63.       HIPL1A=HIPL1*APLEN/(.24*WAIRF(1)*NCELLS)
  64.       IF(HIPL1A.GT.173.)HIPL1A=173. 
  65.       HIPLE(1)=1./EXP(HIPL1A) 
  66.       HIPL2A=HIPL2*APLEN/(.24*WAIRF(2)*NCELLS)
  67.       HIPLE(2)=HIPL2A*.5
  68.       HOPLE= HOPL*APLEN/NCELLS
  69.       RETURN
  70. C  *  *  *  *  *  PLENUM EQNS 
  71. C   Q=H*A*LOG MEAN TEMP DIF (BLWR ON) 
  72. C   Q=H*A*AVG DELTA T (BLWR OFF)
  73. C 
  74.       ENTRY PLENM1
  75.       WAIRF(1)=CFM(1)*C1
  76.       WAIRF(2)=CFM(2)*C1
  77.       IEFAN1=1-IEFAN
  78.       NTB=0 
  79.       HIPL=HIPLE(1)*IEFAN+HIPLE(2)*IEFAN1 
  80.   317 TMPL=TMPLSV+0.5*DTMPL 
  81.       IF(IEFAN.EQ.1)TPL=TMPL+HIPL*(TA2-TMPL)
  82.       IF(IEFAN1.EQ.1)THEN 
  83.            TPL=(TA2+HIPL*(2.*TMPL-TA2))/(1.+HIPL) 
  84.            IF(TA2.GT.TMPL.AND.TPL.LT.TMPL)TPL=TA2-HIPL*(TA2-TMPL) 
  85.            IF(TA2.LT.TMPL.AND.TPL.GT.TMPL)TPL=TA2-HIPL*(TA2-TMPL) 
  86.       END IF
  87.       QIPL=CPWAIR*(TA2-TPL) 
  88.       QPLEN=HOPLE*(TMPL-TAIN) 
  89.       QMPL=QIPL-QPLEN 
  90.       DTMPLQ=QMPL*CONS12
  91.       NTB=NTB+1 
  92.       IF(NTBB.LT.50)THEN
  93.         WRITE(60,516)TIME*60.,TMPL,DTMPLQ,DTMPL,TPL 
  94.         NTBB=NTBB+1
  95.         ENDIF
  96.   516 FORMAT(1X,' TIME,TMPL,DTMPLQ,DTMPL,TPL= ',F8.3,4(2X,F6.1))
  97.       IF(ABS(DTMPLQ-DTMPL).LE.DTMTOL)GO TO 318
  98.       IF(NTB.GT.15)GO TO 318
  99.       DTMPL=DTMPLQ
  100.       GO TO 317 
  101.   318 TMPLSV=TMPLSV+DTMPL 
  102. C  - - - -  DTRA IS AIR TEMP RISE INCL BLOWER RISE
  103.       DTRA=TPL-TRETP+DTBLR*IEFAN
  104. C        COLLECT THE PLENUM TEMPERATURE 
  105.       IF(ICOLL.NE.0) CDATA(12)=TPL
  106.       IF(INC.GT.IDEBUG)RETURN 
  107.       WRITE(60,500) DTRA,TA2,TPL,TMPL,QIPL,QPLEN,QMPL 
  108.       INC = INC +1
  109.       RETURN
  110.   999 STOP ' PLENUM: EOF ON TAPE12 WITH NAMELIST INPGSF'
  111.   500 FORMAT(1X,'PLENUM:DTRA,TA2,TPL,TMPL,QIPL,QPLEN,QMPL= ',7G13.5)
  112.       END 
  113.