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 / OWALLS.FOR < prev    next >
Text File  |  1992-05-08  |  6KB  |  203 lines
  1.       SUBROUTINE OWALLS(NRM, IZW, ID, KON, WALLF, FAWLOD,
  2.      +                  HO, TODBR4, IDBG) 
  3. C 
  4. C  TEMPERATURES IN OUTSIDE ENCLOSURE WALLS
  5. C 
  6. C - - - PARAMETER CONSTANTS 
  7. C 
  8. CMDK NZN
  9. CMDK NIWL
  10. CMDK NKONST
  11. CMDK NWL
  12. CMDK NWN
  13. CMDK NZW
  14. CMDK SIGMA      
  15. C 
  16. C - - - COMMON BLOCKS 
  17. CMDK BLKBSF
  18. CMDK BLKMMU
  19. CMDK BLKQS
  20. CMDK CNSTRK
  21. CMDK ENCBLK
  22. CMDK ENCBK1
  23. CMDK ENCBK2
  24. CMDK IWLS
  25. CMDK IZWQ
  26. CMDK MRTBLK
  27. CMDK NIJB
  28. CMDK OWETHR
  29. CMDK SOILB
  30. CMDK SOLARB
  31. CMDK SURFAR
  32. CMDK TCNBLK
  33. CMDK TEMPB
  34. CMDK TIMEB
  35. CMDK TSTATC
  36. CMDK UABAS
  37. C 
  38.       DATA FTR/459.7/,NQTM/0/
  39.       
  40.       IF(IGRND.GT.0.AND.NQTM.EQ.0.AND.ID.EQ.9)THEN
  41. C         INITIALIZE BELOW GRND BSMT WALLS, IF IGRND > 0     
  42.         NQTM=1
  43. C         ASSUME HEAT FLOW/UNIT AREA IS SAME FOR EACH WALL SEGMENT 
  44. C         SO USE IZW BELOW
  45.         QSURF=QCSURF(IZW)/AWALLA(IZW)        
  46.         DO 2 I=1,4
  47.         IZ=IZWE(I)
  48.         KONF=KONSTA(IZ)
  49.         NLM=NLMP(KONF)
  50.         TWI=TBASM - 5.
  51. C              GET INITIAL VALUE OF INSIDE SURF WALL TEMP        
  52.         CALL TWISET(4,QSURF,TBASM,TWI)
  53.         TMID=0.5*(TSOIL(I)+TWI)
  54.         TMP(1,IZ)=TSOIL(I)
  55.         TMP(2,IZ)=TMID
  56.         TMP(3,IZ)=TWI
  57.         IF(NLM.NE.1)THEN
  58.           TMP(4,IZ)=TWI
  59.           TMP(3,IZ)=TMID+0.5*(TWI-TMID)
  60.         ENDIF
  61.     2   CONTINUE
  62.       ENDIF
  63.       K=KON 
  64.       NLM=NLMP(K) 
  65. C 
  66. C  OUTSIDE WALL TEMPS 
  67.       T1=TMP(1,IZW) 
  68.       T2=TMP(2,IZW) 
  69.       T3=TMP(3,IZW) 
  70.       IF(NLM.EQ.3)T4=TMP(4,IZW) 
  71.       IF(ID.EQ.9)THEN 
  72.          QLSWLN=0.
  73.       ELSE
  74.          T14=(T1+FTR)**4
  75.          FASIWE=FAWLOD*SIGMA*WALLF
  76.          QLSWLN=FASIWE*(T14-TODBR4)*AWALLA(IZW) 
  77.          ENDIF
  78.       IF(ID.LE.4)THEN 
  79. C  NOTE:  NIJA IS SET IN SOLRLD 
  80. C         AND SEE NOTE ON IWALS IN STRRUN OR SOLRLD 
  81.           IF(IWALS.EQ.1)THEN
  82.              ISURF=NIJA(IZW)
  83.           ELSE
  84.              ISURF=NIJA(ID) 
  85.              ENDIF
  86.           TWO=(SOLARL(ISURF)-QLSWLN/AWALLA(IZW))/HO+TODDB 
  87.       ELSE IF(ID.EQ.10)THEN 
  88. C  NO SOLAR RADIATION SHINING ON ABOVE GROUND BASEMENT WALL!      
  89.           TWO=TODDB - QLSWLN/AWALLA(IZW)/HO
  90.           IZWK=IZW
  91.           ENDIF 
  92.       TC1=TCON(1,K) 
  93.       TC2=TCON(2,K) 
  94.       TC3=TCON(3,K) 
  95.       IF(ID.EQ.9)THEN 
  96. C                                BASEMENT WALL BELOW GROUND 
  97.         IHIT=0
  98. C   IZWE IS SET IN CONSLD ABOVE LABEL 134 
  99.         DO 20 I=1,4 
  100.           IF(IZWE(I).EQ.IZW)THEN
  101.             IHIT=1
  102.             IF(IGRND.EQ.0.OR.ICODE.EQ.1.OR.ICODE.EQ.4)THEN
  103.               TMPN(1,IZW)=TSOIL(I)
  104. C  NOTE:  HEAT FLOW TO SOIL WHEN YOU USE TSOIL OPTION W/ IGRND=0
  105. C         WILL NOT BE INCLUDED IN QLOSWL FOR ENERGY BALANCE CALC.
  106. C         THIS HEAT FLOW IS CALC IN ENRBALS FROM AVAIL DATA.
  107.             ELSE
  108.               IJ=I
  109.               IF(I.GT.2)IJ=I+5
  110. C                    CONDWL IS CALC IN SETNOD.FOR ON GROUND.FOR FILE
  111.               QWALBG=CONDWL(I)*(TG(IJ)-T1)
  112.               TMPN(1,IZW)=TC1*(T2-T1)+TC2*QWALBG+T1
  113. C  QLOSWL(I) WAS ZEROED IN STRLP1 
  114.               QLOSWL(NRM)=QLOSWL(NRM)-AWALLA(IZW)*QWALBG
  115.             ENDIF
  116.           ENDIF 
  117.    20   CONTINUE
  118.         IF(IHIT.EQ.0)STOP ' OWALLS: IZW NOT IN IZWE()!'
  119.       ENDIF
  120.       IF(ID.LE.4.OR.ID.EQ.10)THEN 
  121. C            OUTSIDE WALLS OR    BASEMENT WALL ABOVE GROUND 
  122.         QWALOC=HO*(TWO-T1)
  123. C  FOR ID LE 4, QWALOC INCLUDES EFFECT OF SOLAR AND RADIATION LOSS        
  124. C  FOR ID EQ 10, QWALOC INCLUDES RADIATION LOSS
  125.         TMPN(1,IZW)=TC1*(T2-T1)+TC2*QWALOC+T1 
  126.       ENDIF 
  127.       IF(NLM.EQ.1)THEN
  128.           TMPN(2,IZW)=0.5*TC1*(T1-2.*T2+T3)+T2
  129.           IF(ICODE.EQ.1.OR.ICODE.EQ.4)THEN
  130.             QSURF=QCSURF(IZW)/AWALLA(IZW)
  131.             TWI=T3
  132.             CALL TWISET(4,QSURF,TBASM,TWI)
  133.             TMPN(3,IZW)=TWI
  134.           ELSE
  135.             QPRAD=HRI(IZW)*(TMR(NRM)-T3)
  136.             QSURF=QCSURF(IZW)/AWALLA(IZW)+QPRAD
  137. C  QCSURF IS MEAN CONVECTIVE HEAT FLOW CALC IN TAZONE OR TBZONE .
  138. C  IT IS CALC IN TAZONE, TBZONE AND AVERAGED IN ENDLP1.
  139.             TMPN(3,IZW)=TC3*(T2-T3)+TC2*QSURF+T3
  140.           ENDIF
  141.       ELSE
  142.           TC4=TCON(4,K) 
  143.           TC5=TCON(5,K) 
  144.           TC6=TCON(6,K) 
  145.           TC7=TCON(7,K) 
  146.           TC8=TCON(8,K) 
  147.           TMPN(2,IZW)=TC3*(T1-T2)+TC4*(T3-T2)+T2
  148.           TMPN(3,IZW)=TC5*(T2-T3)+TC6*(T4-T3)+T3
  149.           IF(ICODE.EQ.1.OR.ICODE.EQ.4)THEN
  150. C                                 QCSURF IS SET IN BSMTPRE          
  151.             QSURF=QCSURF(IZW)/AWALLA(IZW)
  152.             TWI=T4
  153.             CALL TWISET(4,QSURF,TBASM,TWI)
  154.             TMPN(4,IZW)=TWI
  155.           ELSE
  156.             QPRAD=HRI(IZW)*(TMR(NRM)-T4)
  157. C                                 QCSURF IS SET IN TBZONE            
  158.             QSURF=QCSURF(IZW)/AWALLA(IZW) + QPRAD
  159.             TMPN(4,IZW)=TC7*(T3-T4)+TC8*QSURF+T4 
  160.           ENDIF
  161.       ENDIF 
  162.       QRMRT(IZW)=QPRAD*AWALLA(IZW)
  163.       IF(IDBG.NE.0)THEN 
  164.           WRITE(60,500)NRM, IZW,ID,KON, WALLF, FAWLOD,
  165.      +                  HO, TODBR4,         QLSWLN
  166.           IF(NLM.EQ.3)THEN
  167.               WRITE(60,501)T1,TMPN(1,IZW),T2,TMPN(2,IZW),T3, 
  168.      +        TMPN(3,IZW),T4,TMPN(4,IZW),AWALLA(IZW), 
  169.      +        ISURF,IT,SOLARL(ISURF),DSOLAR(IT,ISURF),
  170.      +        TODDB,TWO,TA(NRM),HO,
  171.      +        TC1,TC2,TC3,TC4,TC5,TC6,TC7,TC8 
  172.           ELSE IF(NLM.EQ.1)THEN 
  173.               WRITE(60,503)T1,TMPN(1,IZW),T2,TMPN(2,IZW),T3,TMPN(3,IZW), 
  174.      +                    AWALLA(IZW),TODDB,TA(NRM),HO,TC1,TC2,TC3 
  175.               ENDIF 
  176.           WRITE(60,502)QCSURF(IZW),QRMRT(IZW),TMR(NRM) 
  177.           ENDIF 
  178. C     IDBG=0
  179.       RETURN
  180.   500 FORMAT(15X,'OWALLS:NRM,IZW,ID,KON= ',4I5/ 
  181.      +15X,'WALLF,FAWLOD,HO,TODBR4= ',4G13.5/
  182.      +15X,'    QLOSWL(IZW)= ', G13.5) 
  183.   501 FORMAT(15X,'T1,T1N= ',2G13.5/ 
  184.      +15X,'T2,T2N= ',2G13.5/
  185.      +15X,'T3,T3N= ',2G13.5/
  186.      +15X,'T4,T4N= ',2G13.5/
  187.      +15X,'AWALLA(IZW),ISURF,IT,SOLARL,DSOLAR= ',G13.5,2I5,2G13.5/
  188.      +15X,'TODDB,TWO,TA(NRM),H= ',4G13.5/ 
  189.      +15X,'TC1= ',G13.5,' TC2= ',G13.5/ 
  190.      +15X,'TC3= ',G13.5,' TC4= ',G13.5/ 
  191.      +15X,'TC5= ',G13.5,' TC6= ',G13.5/ 
  192.      +15X,'TC7= ',G13.5,' TC8= ',G13.5) 
  193.   502 FORMAT(15X,'QCSURF,QRMRT,TMR= ',3G13.5) 
  194.   503 FORMAT(12X,'T1,T1N= ',2G13.5/ 
  195.      +12X,'T2,T2N= ',2G13.5/
  196.      +12X,'T3,T3N= ',2G13.5/
  197.      +12X,'AWALLA(IZW),TODDB,TA(NRM),H= ',4G13.5/ 
  198.      +12X,'TC1= ',G13.5,' TC2= ',G13.5/ 
  199.      +12X,'TC3= ',G13.5)
  200.   504 FORMAT(5X,'OWALLS:CHECK--I,IZW,IZWE,CONDWL(I),T1,TMPN(1,IZW)= '/
  201.      +5X,3I3,3G13.5)
  202.       END 
  203.