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 / GARAG.FOR < prev    next >
Text File  |  1992-04-15  |  4KB  |  105 lines
  1.       SUBROUTINE GARAG(IGARAGT,NZON,SOLARH) 
  2. C 
  3. C  COMPUTE AIR TEMP IN AN ATTACHED GARAGE 
  4. C   APPROACH: INNER SURF OF OUTSIDE WALLS AND CEILING DRIVEN TOWARD 
  5. C             AVG OF SOLAIR AND GARAGE AIR TEMPS WITH A TIME CONST EQN
  6. C   IGARAGT:  SEE STRRUN
  7. C   NZON - ZONE CONNECTING TO ATTACH GARAGE (1 - 5)
  8. C - - CONSTANTS 
  9. CMDK NZN
  10. CMDK OWETHR
  11. CMDK TEMP1
  12. CMDK TIMEB
  13.       DIMENSION B(3),DTG(3),BB(3) 
  14. C      
  15. C INPUTS IN NAMELIST GAR:    
  16. C  ALFSOL  AREA WEIGHTED MEAN SOLAR ABSORPTIVITY OF OUTSIDE SURF OF ROOF
  17. C          AND OUTSIDE WALLS OF GARAGE(USED IN SOLAIR TEMP CALC)
  18. C  AWALI   AREA OF WALL BETWEEN GARAGE AND LIVING SPACE, FT SQ
  19. C  AWALO   AREA OF OUTSIDE GARAGE WALLS, FT SQ
  20. C  HCEIL   HEAT TRANSFER COEF AT INSIDE SURF OF CEILING(CONV AND RAD),BTU/HR-FTF
  21. C  HWALO   HEAT TRANSFER COEF AT INSIDE SURF OF OUTSIDE WALLS(CONV AND RAD),
  22. C          BTU/HR-FT SQ-F 
  23. C  TCWALO  TIME CONSTANT FOR CHANGE IN TEMP OF INSIDE SURF OF OUTSIDE WALL,HR 
  24. C  TCCEIL  TIME CONSTANT FOR CHANGE IN TEMP OF INSIDE SURF OF CEILING,HR
  25. C  UWALI   OVERALL HEAT TRANSFER COEF FOR HEAT FLOW THRU WALL BETWEEN GARAGE
  26. C          AND LIVING SPACE, BTU/HR-FT SQ-F 
  27. C  VOLGAR  VOLUME OF GARAGE, FT SQ
  28.       NAMELIST/GAR/ALFSOL,AWALI,AWALO,HCEIL,HWALO,
  29.      +             TCWALO,TCCEIL,UWALI,VOLGAR
  30. C  VALUES IN DATA STMTS SET FOR HOUSE A
  31. C     (SINGLE STORY HOUSE WITH SINGLE-CAR ATTACHED GARAGE ON SOUTH SIDE)
  32.       DATA FTR/459.7/ 
  33.       DATA HWALO,AWALO/.22,400./
  34.       DATA HCEIL,ACEIL/.37,308./
  35.       DATA UWALI,AWALI/.31,176./
  36.       DATA VOLGAR/2464./
  37.       DATA TCWALO,TCCEIL/2.,3./ 
  38.       DATA ALFSOL/.65/
  39. C  ACHG1  REF. AIR CHANGE DUE TO REF. TEMP DIFFERENCE TGAR1-TOD1
  40. C  ACHG2  REF. AIR CHANGE DUE TO REF. WIND SPEED OF WS2 MPH 
  41.       DATA ACHG1,TOD1,TGAR1/0.5,20.,50./
  42.       DATA ACHG2,WS2/0.75,20./
  43. C  FITTED COEF BELOW
  44.       DATA B/.06976,.6911,28.04/
  45.       DATA NTIM/0/
  46. C 
  47.       IF(IGARAGT.EQ.2)THEN
  48.        OPEN(17,FILE='TAPE17',STATUS='OLD',IOSTAT=IO17)
  49.        IF(IO17.NE.0)THEN
  50.          WRITE(60,*) ' GARAG: CANT OPEN TAPE17 WITH NAMELIST GAR'
  51.          STOP ' GARAG: CANT OPEN TAPE17 WITH NAMELIST GAR'
  52.          END IF
  53.        READ(17,GAR,END=99)
  54.        CLOSE(17)
  55.        WRITE(60,GAR)
  56.        ENDIF
  57.       CAPGAR=0.075*0.24*VOLGAR
  58.       BB(1)=DTIME*HCEIL*ACEIL/CAPGAR
  59.       BB(2)=DTIME*HWALO*AWALO/CAPGAR
  60.       BB(3)=DTIME*UWALI*AWALI/CAPGAR
  61.       B4=ACHG1*ACHG1/(1./(TOD1+FTR)-1./(TGAR1+FTR)) 
  62.       B5=(ACHG2*ACHG2-ACHG1*ACHG1)*(TOD1+FTR)/(WS2*WS2) 
  63.       WRITE(60,500) (BB(I),I=1,3),B4,B5 
  64.       TCTWLO=DTIME*B(1)/TCWALO
  65.       TCTCLG=DTIME*B(1)/TCCEIL
  66.       TGARAG=.66*TI(NZON)+.34*TODDB
  67.       RETURN
  68. C 
  69.       ENTRY GARAGA(IGARAGT,NZON,SOLARH)
  70.       IF(NTIM.EQ.0)THEN 
  71.           TCEIL=TGARAG
  72.           TWALO=TCEIL 
  73.           NTIM=1
  74.           ENDIF 
  75.       IF(WSPED.GT.0.4)U=SQRT(WSPED)/.6584 
  76.       IF(WSPED.LE.0.4)U=1.
  77.       TSLAIR=(TODDB+ALFSOL*SOLARH/U+TGARAG)*0.5 
  78.       TCEIL=TCTCLG*(TSLAIR-TCEIL)+TCEIL 
  79.       TWALO=TCTWLO*(TSLAIR-TWALO)+TWALO 
  80.       DTG(1)=BB(1)*(TCEIL-TGARAG)+BB(2)*(TWALO-TGARAG)
  81.       DTG(2)=B(2)*BB(3)*(TI(NZON)-TGARAG)
  82.       TOD11=1./(TODDB+FTR)
  83.       SQTRM=B4*(TOD11-1./(TGARAG+FTR))+B5*B(3)*WSPED*WSPED*TOD11
  84.       IF(SQTRM.LT.0.)ACHG=0.
  85.       IF(SQTRM.LT.0.)GO TO 48 
  86.       ACHG=SQRT(SQTRM)
  87.    48 CONTINUE
  88.       DTG(3)=DTIME*ACHG*(TODDB-TGARAG)
  89. C 
  90.       SUM=0.
  91.       DO 50 K=1,3 
  92.    50 SUM=SUM+DTG(K)
  93. C      TGARSV=TGARAG 
  94.       TGARAG=TGARAG+SUM 
  95. C      NQQ=NQQ+1 
  96. C      IF(NQQ.GT.360.AND.NQQ.LT.381)WRITE(60,501) WSPED,SOLARH,TODDB,
  97. C     +TGARSV,(DTG(IK),IK=1,3),TGARAG,ACHG 
  98.       RETURN
  99.    99 STOP ' GARAG: EOF ON TAPE17'      
  100.   500 FORMAT(1X,'GARAG: BB S= ',3G13.5,4HB4= ,G13.5,4HB5= ,G13.5) 
  101.   501 FORMAT(1X,'WSPED,SOLARH,TODDB= ',3G13.5/
  102.      +1X,'TGARAG LAST= ',G13.5,8HDTG(I)= ,3G13.5/ 
  103.      +1X,'TGARAG=      ',G13.5,6HACHG= ,G13.5)
  104.       END 
  105.