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 / IWALLS.FOR < prev    next >
Text File  |  1992-05-08  |  7KB  |  209 lines
  1.       SUBROUTINE IWALLS(NRM,IZW,KON,NZONC,TWIA,AWIA,IDBG) 
  2. C 
  3. C  TEMPERATURES OF INTERIOR WALLS AND LVG SPC-TO-GARAGE WALLS 
  4. C  TWIA AND AWIA ARE USED TO GET THE AREA WEIGHTED PARTITION WALL SURF
  5. C                TEMP FOR USE BY THE THERMOSTAT MODEL 
  6. C 
  7. C - - - PARAMETER CONSTANTS 
  8. C 
  9. CMDK FTR
  10. CMDK NKONST
  11. CMDK NWL
  12. CMDK NWN
  13. CMDK NZN
  14. CMDK NZW
  15. CMDK SIGMA
  16. C 
  17. C - - - COMMON BLOCKS 
  18. CMDK BLKQS
  19. CMDK CNSTRK
  20. CMDK ENCBK1
  21. CMDK ENCBK2
  22. CMDK ENCBLK
  23. CMDK IZPART
  24. CMDK MRTBLK
  25. CMDK OWETHR
  26. CMDK SURFAR
  27. CMDK TCNBLK
  28. CMDK TCONCB
  29. CMDK TEMPB
  30. CMDK TIMEB
  31.       LOGICAL NZY 
  32. C 
  33.       IDBGSV=IDBG 
  34.       K=KON 
  35.       NLM=NLMP(K) 
  36.       NZON=NZONC-8
  37.       NZY=.FALSE. 
  38.       IF(IDBG.NE.0)WRITE(60,506)TIME,NZON,NRM,K,KNSTAT(NRM)
  39.       IF(NZON.EQ.NRM.AND.K.EQ.KNSTAT(NRM))NZY=.TRUE.
  40.       AWL=AWALLA(IZW) 
  41.       IF(NZY)AWIA=AWIA+AWL
  42.       T1=TMP(1,IZW) 
  43.       T2=TMP(2,IZW) 
  44.       T3=TMP(3,IZW) 
  45.       IF(NZON.NE.NRM.AND.NZONC.NE.6)THEN
  46. C  PARTITION WALL BETWEEN ZONES 
  47.           IZWO=IZSET(IZW) 
  48.           NRMO=NZON 
  49. C - - -NOW CHECK TO SEE IF WALL TEMPS HAVE ALREADY BEEN CALCULATED IN 
  50. C        OTHER ZONE (IF CALC ORDER IS KEPT THE SAME AS INPUT ORDER, 
  51. C        THEN IZW FOR BACKSIDE SURFACE WILL BE GREATER THAN IZW FOR 
  52. C        SURFACE IN THIS ZONE 
  53.           IF(IZWO.LT.IZW)THEN 
  54. C - - -TEMPS FOR ENCLOSURE HAVE ALREADY BEEN CALC IN ZONE NZON
  55. C          SO RESET TEMPS AND RETURN
  56.               IF(NLM.EQ.1)THEN
  57.                   TMPN(1,IZW)=TMPN(3,IZWO)
  58.                   TMPN(2,IZW)=TMPN(2,IZWO)
  59.                   TMPN(3,IZW)=TMPN(1,IZWO)
  60.                   QRMRT(IZW)=QRMRT(IZWO)
  61.                   IF(NZY)TWIA=TMPN(3,IZW)*AWL+TWIA
  62.                   IF(IDBG.NE.0)WRITE(60,504)NRM,IZW,TMPN(1,IZW), 
  63.      +              TMPN(2,IZW),TMPN(3,IZW) 
  64.                   GO TO 200 
  65.                   ENDIF 
  66.               IF(NLM.EQ.3)THEN
  67.                   TMPN(1,IZW)=TMPN(4,IZWO)
  68.                   TMPN(2,IZW)=TMPN(3,IZWO)
  69.                   TMPN(3,IZW)=TMPN(2,IZWO)
  70.                   TMPN(4,IZW)=TMPN(1,IZWO)
  71.                   QRMRT(IZW)=QRMRT(IZWO)
  72.                   IF(NZY)TWIA=TMPN(4,IZW)*AWL+TWIA
  73.                   IF(IDBG.NE.0)WRITE(60,505)NRM,IZW,TMPN(1,IZW), 
  74.      +              TMPN(2,IZW),TMPN(3,IZW),TMPN(4,IZW) 
  75.                   GO TO 200 
  76.                   ENDIF 
  77.               ENDIF 
  78.       ELSE
  79. C              IMMERSED PARTITION WALLS AND GARAG WALL!      
  80.           IZWO=IZW
  81.           NRMO=NRM
  82.           ENDIF 
  83.       TC1=TCON(1,K) 
  84.       TC2=TCON(2,K) 
  85.       TC3=TCON(3,K) 
  86.       IF(NZONC.NE.6)THEN
  87.          QPRADO=HRI(IZWO)*(TMR(NRMO)-T1)+QSI(IZWO)
  88.       ELSE
  89. C  RADIATION AND CONVECTION FROM GARAGE-SIDE OF LVG SPC-TO-GARAGE WALL
  90.          TMN=0.5*(TGARAG+T1)
  91.          QPRADO=4.*WALLFO(IZW)*SIGMA*(TMN+FTR)**3 *(TGARAG-T1)
  92.          HIO=HC(NRM,IZW,5,.FALSE.,T1,TGARAG)
  93.          QCGAR=HIO*AWALLA(IZW)*(TGARAG-T1)
  94.          QLOSWL(NRM)=QLOSWL(NRM)-QCGAR-AWALLA(IZW)*QPRADO
  95.       ENDIF
  96.       QCOVAI=QCSURF(IZW)/AWALLA(IZW)
  97.       IF(NZONC.NE.6)THEN
  98.          QCOVAO=QCSURF(IZWO)/AWALLA(IZW)
  99.       ELSE
  100.          QCOVAO=QCGAR/AWALLA(IZW) 
  101.          ENDIF
  102.       IF(NLM.EQ.1)THEN
  103. C  CHECK FOR LIGHT WALL ELEMENT 
  104.           IF(NKSTAB.NE.0)THEN 
  105.               DO 10 IK=1,NKSTAB 
  106.    10         IF(K.EQ.KONSTB(IK))GO TO 20 
  107.               GO TO 40
  108.    20         NRAT=NRATTM(IK) 
  109.               XRAT=1./NRAT
  110.               TS1=T1
  111.               TS2=T2
  112.               TS3=T3
  113.               TCS1=XRAT*TC1 
  114.               TCS2=XRAT*TC2 
  115.               TCS3=XRAT*TC3 
  116.               QPRDO=0.
  117.               QPRDI=0.
  118.               DO 30 NT=1,NRAT 
  119.               IF(NT.NE.1)THEN
  120.                 QPRADO=HRI(IZWO)*(TMR(NRMO)-TS1)+QSI(IZWO) 
  121.               ENDIF
  122.               TS1N=TCS3*(TS2-TS1)+TCS2*(QCOVAO+QPRADO)+TS1
  123.               TS2N=0.5*TCS1*(TS1-2.*TS2+TS3)+TS2
  124.               QPRADI=HRI(IZW)*(TMR(NRM)-TS3)+QSI(IZW) 
  125.               TS3N=TCS3*(TS2-TS3)+TCS2*(QCOVAI+QPRADI)+TS3
  126.               TS1=TS1N
  127.               TS2=TS2N
  128.               TS3=TS3N
  129.               QPRDO=QPRDO+QPRADO
  130.               QPRDI=QPRDI+QPRADI
  131.    30         CONTINUE
  132.               QPRADO=QPRDO/NRAT 
  133.               QPRADI=QPRDI/NRAT 
  134.               TMPN(1,IZW)=TS1 
  135.               TMPN(2,IZW)=TS2 
  136.               TMPN(3,IZW)=TS3 
  137.               IF(NZY)TWIA=TMPN(3,IZW)*AWL+TWIA
  138.               GO TO 50
  139.               ENDIF 
  140.    40     TMPN(1,IZW)=TC1*(T2-T1)+TC2*(QCOVAO+QPRADO)+T1
  141.           TMPN(2,IZW)=0.5*TC1*(T1-2.*T2+T3)+T2
  142.           QPRADI=HRI(IZW)*(TMR(NRM)-T3)+QSI(IZW)
  143.           TMPN(3,IZW)=TC3*(T2-T3)+TC2*(QCOVAI+QPRADI)+T3
  144.           IF(NZY)TWIA=TMPN(3,IZW)*AWL+TWIA
  145.           ENDIF 
  146.       IF(NLM.EQ.3)THEN
  147.           TC4=TCON(4,K) 
  148.           TC5=TCON(5,K) 
  149.           TC6=TCON(6,K) 
  150.           TC7=TCON(7,K) 
  151.           TC8=TCON(8,K) 
  152.           T4=TMP(4,IZW) 
  153.           TMPN(1,IZW)=TC1*(T2-T1)+TC2*(QCOVAO+QPRADO)+T1
  154.           TMPN(2,IZW)=TC3*(T1-T2)+TC4*(T3-T2)+T2
  155.           TMPN(3,IZW)=TC5*(T2-T3)+TC6*(T4-T3)+T3
  156.           QPRADI=HRI(IZW)*(TMR(NRM)-T4)+QSI(IZW)
  157.           TMPN(4,IZW)=TC7*(T3-T4)+TC8*(QCOVAI+QPRADI)+T4
  158.           IF(NZY)TWIA=TMPN(4,IZW)*AWL+TWIA
  159.           ENDIF 
  160.    50 CONTINUE
  161.       QRMRT(IZW)=QPRADI*AWALLA(IZW) 
  162.       IF(NZONC.NE.6)THEN
  163.          QRMRT(IZWO)=QPRADO*AWALLA(IZW) 
  164.       ELSE
  165.          QRGAR=QPRADO*AWALLA(IZW) 
  166.          ENDIF
  167.       IF(IDBG.NE.0)THEN 
  168.           WRITE(60,500)NRM,IZW,KON,NZON
  169.           IF(NLM.EQ.1)THEN
  170.               WRITE(60,501)T1,TMPN(1,IZW),T2,TMPN(2,IZW),T3,TMPN(3,IZW), 
  171.      +        TC1,TC2,TC3 
  172.               ENDIF 
  173.           IF(NLM.EQ.3)THEN
  174.               WRITE(60,502)T1,TMPN(1,IZW),T2,TMPN(2,IZW),T3,TMPN(3,IZW), 
  175.      +        T4,TMPN(4,IZW),TC1,TC2,TC3,TC4,TC5,TC6,TC7,TC8
  176.               ENDIF 
  177.           IF(NZONC.NE.6)THEN
  178.              WRITE(60,503)QCSURF(IZW),QCSURF(IZWO),QRMRT(IZW), 
  179.      +                   QRMRT(IZWO),TMR(NRM),TMR(NRMO),QSI(IZW), 
  180.      +                   QSI(IZWO)
  181.           ELSE
  182.              WRITE(60,503)QCSURF(IZW),QCGAR,QRMRT(IZWO),QRGAR, 
  183.      +                   TMR(NRM),TGARAG,QSI(IZW),0.
  184.              ENDIF
  185.           ENDIF
  186.   200 CONTINUE
  187.       IDBG=IDBGSV 
  188.       RETURN
  189.   500 FORMAT(10X,'IWALLS:NRM,IZW,KON,NZONC= ',4I5)
  190.   501 FORMAT(10X,'T1,T1N= ',2G13.5/ 
  191.      +10X,'T2,T2N= ',2G13.5/
  192.      +10X,'T3,T3N= ',2G13.5/
  193.      +10X,'TC1= ',G13.5,'  TC2= ',G13.5,'  TC3= ',G13.5)
  194.   502 FORMAT(10X,'T1,T1N= ',2G13.5/ 
  195.      +10X,'T2,T2N= ',2G13.5/
  196.      +10X,'T3,T3N= ',2G13.5/
  197.      +10X,'T4,T4N= ',2G13.5/
  198.      +10X,'TC1= ',G13.5,'  TC2= ',G13.5,'  TC3= ',G13.5,' TC4= ',G13.5/ 
  199.      +10X,'TC5= ',G13.5,'  TC6= ',G13.5,'  TC7= ',G13.5,' TC8= ',G13.5) 
  200.   503 FORMAT(10X,'QCSURF(AT INNER AND OUTER SURF)= ',2G13.5/
  201.      +10X,'QRMRT(AT INNER AND OUTER SURF)= ',2G13.5/
  202.      +10X,'TMR(IN ZONE AND CONNECTING ZONE)= ',2G13.5/
  203.      +10X,'QSI(AT INNER AND OUTER SURF)= ',2G13.5)
  204.   504 FORMAT(10X,'NRM,IZW= ',2I5,' TMPN(1-3)= ',3G13.5) 
  205.   505 FORMAT(10X,'NRM,IZW= ',2I5,' TMPN(1-4)= ',4G13.5) 
  206.   506 FORMAT(10X,'IWALLS:TIME= ',F12.5/13X,'LOOK FOR TSTAT;NZON', 
  207.      +',NRM,K,KNSTAT(NRM)=', 4I5) 
  208.       END 
  209.