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 / HBKF.FOR < prev    next >
Text File  |  1992-05-09  |  7KB  |  196 lines
  1.       SUBROUTINE HBKFACT(T) 
  2. C ROUTINE TO CALCULATE HANDBOOK FACTORS AT THE END OF EACH DAY OF 
  3. C SIMULATION. AND TO ACCUMULATE THEM FOR THE TOTAL RUN
  4. C  REF:  SEE REF NO. 3 ON FILE \HII\SCR\S2.DOC AND 
  5. C        ASHRAE 1988 EQUIPMENT HANDBOOK, P. 24.8 TO 24.17
  6. C        FOR MORE INFORMATION ABOUT HANDBOOK FACTORS AND DATA 
  7. C        ON TAPES 25 AND TAPES33-37
  8. C TAPE25 IS THE REFERANCE NET LOAD FILE. WRITE IF HBKRNL=0               
  9. C                                        READ  IF HBKRNL.NE.0 
  10. CMDK NHBK
  11. CMDK HANDBK
  12.       CHARACTER*80 NAME(30) 
  13.       DOUBLE PRECISION T
  14.       LOGICAL FIRST,I25F,I32F,I37F
  15.       REAL HBKOUT(30) 
  16.       REAL MGF,NL,ILF,LMF,LRC(3),CRF(3),CMSI(3) 
  17.       DATA LRC/1.0, 3.0, 5.0/, FIRST/.TRUE./
  18.       DATA I25F,I32F,I37F/.TRUE.,.TRUE.,.TRUE./
  19. 100   IF(FIRST)GOTO 800 
  20.       HBKOUT(01)=HBK(1) 
  21.       FO1=HBK(1)
  22.       HBKOUT(02)=HBK(2) 
  23. C      FO2=HBK(2)
  24.       HBKOUT(03)=HBK(3) 
  25. C      FO3=HBK(3)
  26.       HBKOUT(04)=HBK(4) 
  27.       PEI=HBK(4)
  28.       HBKOUT(05)=HBK(5) 
  29.       AEI=HBK(5)
  30.       TEI=PEI+AEI 
  31.       HBKOUT(06)=TEI
  32.       DENOM=TEI 
  33.       IF(TEI.EQ.0.0)DENOM=1.0E10
  34.       FE=FO1/DENOM
  35.       HBKOUT(07)=FE 
  36.       HBKOUT(08)=FO1
  37.       DI=FO1
  38.       HBKOUT(09)=HBK(6) 
  39.       DO=HBK(6) 
  40.       DENOM=DI
  41.       IF(DI.EQ.0.0)DENOM=1.0E10 
  42.       DE=DO/DENOM 
  43.       HBKOUT(10)=  DE 
  44.       HDE=FE*DE 
  45.       HBKOUT(11)=  HDE
  46.       HBKOUT(12)=HBK(7) 
  47.       THD=HBK(7)
  48.       DENOM=DO
  49.       IF(DO.EQ.0.0)DENOM=1.0E10 
  50.       MGF=THD/DENOM 
  51.       HBKOUT(13)=  MGF
  52.       SE=HDE*MGF
  53.       HBKOUT(14)=  SE 
  54.       RNL=HBKRNL
  55.       HBKOUT(15)=  RNL
  56.       HBKOUT(16) = THD
  57.       NL=THD
  58.       IF(RNL.EQ.0)THEN
  59. C            WRITE DATA FOR RNL ON TAPE25      
  60.         RNL=NL 
  61.         IF(I25F)THEN
  62.           I25F=.FALSE.
  63.           OPEN(25,FILE='TAPE25',STATUS='NEW',IOSTAT=IO25)
  64. C    FOLLOWING THREE LINES ADDED BY RDF   3/11/91 !!!!!          
  65.           IF(IO25.NE.0)THEN
  66.             OPEN(25,FILE='TAPE25',STATUS='OLD')
  67.           ENDIF
  68.         ENDIF
  69.         WRITE(25,*)RNL 
  70.       ELSE
  71. C          READ DATA FOR RNL ON TAPE25      
  72.         IF(I25F)THEN
  73.           I25F=.FALSE.
  74.           OPEN(25,FILE='TAPE25',STATUS='OLD',IOSTAT=IO25)
  75. C   FOLLOWING IF SEGMENT ADDED BY RDF  3/5/91  !!!!!!!          
  76.           IF(IO25.NE.0)THEN
  77.             OPEN(25,FILE='TAPE25',STATUS='NEW')
  78.             RNL=NL
  79.             WRITE(25,*)RNL
  80.             CLOSE(25)
  81.             OPEN(25,FILE='TAPE25',STATUS='OLD')
  82.           ENDIF
  83. C RNL - REFERENCE NET LOAD  (FROM A PREVIOUS RUN OF BASELINE CASE)          
  84.         READ(25,*,END=910)RNL
  85.         ENDIF 
  86.       ENDIF
  87.       SIL=NL-RNL
  88.       HBKOUT(17)=  SIL
  89.       DENOM=NL 
  90.       IF(NL.EQ.0.0)DENOM=1.0E10
  91.       ILF=SIL/DENOM 
  92.       HBKOUT(18)=  ILF
  93.       LMF=1.0-ILF 
  94.       HBKOUT(19)=  LMF
  95.       SI=SE*LMF 
  96.       HBKOUT(20)=  SI 
  97.       DENOM=PEI
  98.       IF(PEI.EQ.0.0)DENOM=1.0E10 
  99.       AER=AEI/DENOM 
  100.       HBKOUT(21)=  AER
  101.       DO 210 I=1,3
  102.       II=18+3*I 
  103.       HBKOUT(II+1)=  LRC(I) 
  104.       CRF(I)=(1.0+AER)/(1.0+LRC(I)*AER) 
  105.       HBKOUT(II+2)=  CRF(I) 
  106.       CMSI(I)=SI*CRF(I) 
  107. 210   HBKOUT(II+3)=  CMSI(I)
  108. C 
  109. C           PRINT OUT RESULTS 
  110. C 
  111.       IF(I32F)THEN
  112.         I32F=.FALSE.
  113.         OPEN(32,FILE='TAPE32',STATUS='NEW')
  114.       ENDIF
  115.       WRITE(32,1000)T,
  116.      +     (I,NAME(I)(1:10),HBKOUT(I),NAME(I)(11:80),I=1,30)
  117. 1000  FORMAT('1HANDBOOK FACTORS ',8X,'(TIME=',F13.7,'):'//
  118.      +(' FACTOR NUMBER ',I2,') ',A10,F20.5,A60) ) 
  119.       IF(I37F)THEN
  120.         I37F=.FALSE.
  121.         OPEN(37,FILE='TAPE37',STATUS='NEW')
  122.       ENDIF
  123.       WRITE(37,1010)(HBK(I),I=1,7)
  124. 1010  FORMAT(7F18.5)
  125. C 
  126.       RETURN
  127. 910   STOP 'SUB.HBK: EOF ON T25'
  128. C 
  129. C-------------------------------------------------------------------- 
  130. C  FIRST TIME THROUGH, SET UP THE OUTPUT FILE 
  131. C 
  132. 800   CONTINUE
  133.       NAME(1)=
  134.      +'FO1..(BTU) FURNACE OUTPUT = ENTHALPY CHANGE ACROSS FURNACE   ' 
  135.       NAME(2)=
  136.      +'FO2..(BTU) FURNACE OUTPUT = INPUT - FURNACE LOSSES           ' 
  137.       NAME(3)=
  138.      +'FO3..(BTU) FURNACE OUTPUT = ENERGY DELIVERED + DUCT LOSS     ' 
  139.       NAME(4)=
  140.      +'PEI..(BTU) PRIMARY ENERGY INPUT                              ' 
  141.       NAME(5)=
  142.      +'AEI..(BTU) AUXILIARY ENERGY INPUT                            ' 
  143.       NAME(6)=
  144.      +'TEI..(BTU) TOTAL ENERGY INPUT                                ' 
  145.       NAME(7)=
  146.      +'FE.....(%) FURNACE EFFICIENCY                                ' 
  147.       NAME(8)=
  148.      +'DI...(BTU) DUCT INPUT = FO1 (ABOVE)                          ' 
  149.       NAME(9)=
  150.      +'DO...(BTU) DUCT OUTPUT                                       ' 
  151.       NAME(10)= 
  152.      +'DE.....(%) DUCT EFFICIENCY                                   ' 
  153.       NAME(11)= 
  154.      +'HDE....(%) HEAT DELIVERY EFFICIENCY                          ' 
  155.       NAME(12)= 
  156.      +'THD..(BTU) TOTAL HEAT DELIVERY (TO CONDITIONED SPACE)        ' 
  157.       NAME(13)= 
  158.      +'MGF....(-) MISCELLANEOUS GAIN FACTOR                         ' 
  159.       NAME(14)= 
  160.      +'SE.....(%) SYSTEM EFFICIENCY                                 ' 
  161.       NAME(15)= 
  162.      +'RNL..(BTU) REFERENCE NET LOAD                                ' 
  163.       NAME(16)= 
  164.      +'NL...(BTU) NET LOAD                                          ' 
  165.       NAME(17)= 
  166.      +'SIL..(BTU) SYSTEM INDUCED LOAD                               ' 
  167.       NAME(18)= 
  168.      +'ILF....(-) INDUCED LOAD FACTOR                               ' 
  169.       NAME(19)= 
  170.      +'LMF....(-) LOAD MODIFICATION FACTOR                          ' 
  171.       NAME(20)= 
  172.      +'SI.....(%) SYSTEM INDEX                                      ' 
  173.       NAME(21)= 
  174.      +'AER....(-) AUXILIARY ENERGY RATIO                            ' 
  175.       NAME(22)=
  176.      +'LRC1.($/$) LOCAL RATIO OF ELECTRICITY TO GAS (BTU BASIS) =1.0' 
  177.       NAME(23)= 
  178.      +'CRF1...(-) COST RATIO FACTOR                                 ' 
  179.       NAME(24)= 
  180.      +'CMSI1..(%) COST MODIFIED SYSTEM INDEX                        ' 
  181.       NAME(25)= 
  182.      +'LRC2.($/$) LOCAL RATIO OF ELECTRICITY TO GAS (BTU BASIS) =3.0' 
  183.       NAME(26)= 
  184.      +'CRF2...(-) COST RATIO FACTOR                                 ' 
  185.       NAME(27)= 
  186.      +'CMSI2..(%) COST MODIFIED SYSTEM INDEX                        ' 
  187.       NAME(28)= 
  188.      +'LRC3.($/$) LOCAL RATIO OF ELECTRICITY TO GAS (BTU BASIS) =5.0' 
  189.       NAME(29)= 
  190.      +'CRF3...(-) COST RATIO FACTOR                                 ' 
  191.       NAME(30)= 
  192.      +'CMSI3..(%) COST MODIFIED SYSTEM INDEX                        ' 
  193.       FIRST=.FALSE. 
  194.       GOTO 100
  195.       END 
  196.