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 / STRDAY.FOR < prev    next >
Text File  |  1992-05-11  |  17KB  |  569 lines
  1.       SUBROUTINE STRDAY(ATOUT, C, CDAY, DQINT, DSOLRH, DTCR, DTIM1,
  2.      + DTIM2, DTOB, DWHM, DWINT, DWSPED, ENDT, HDUM, HUMON, ICYCCP,
  3.      + ICYC1, ICYC6, IDAYM, IDUM, INC, IREP, IRH, ISWING, KDAYM,
  4.      + MPRINT, NROUT, NZONG, ODDBM, ODHUM, QBASVN, QHUMT, QSPT, QWAS,
  5.      + RELHMI, RHDIR, RHIO, RODB, TBAS, TDIR, TDUM, TEMPDV, TIMAX,
  6.      + TIMAXN, TIMIN, TIO, TONMAX, TRO, VENCRV, WATER, WSC, IDBG)
  7. C  COMFTT REMOVED FROM FORMAL PAR LIST ABOVE TO SAVE MEMORY     
  8. C
  9. C-----START OF DAY LOOP
  10. C
  11. C-----CONSTANTS
  12. CMDK CPAIR
  13. CMDK BTUKWH
  14. CMDK NHBK
  15. CMDK NHRO
  16. CMDK NHRP
  17. CMDK NZN
  18. CMDK NIWL
  19. CMDK NKONST
  20. CMDK NOUTDY
  21. CMDK NWL
  22. CMDK NWN
  23. CMDK NZN7
  24. CMDK NZW
  25. CMDK RHOAIR      
  26. C
  27. C-----COMMON BLOCKS
  28. CMDK ACHBQ
  29. CMDK ACND
  30. CMDK BAR
  31. CMDK BLKCM2
  32. CMDK BLK10
  33. CMDK BLKFMT
  34. CMDK BLKGS1
  35. CMDK BLKGS2
  36. CMDK BLKQGS
  37. CMDK BLKQS
  38. CMDK CNSTRK
  39. CMDK DUCTS1
  40. CMDK DUCTS2
  41. CMDK DUCTS5
  42. CMDK EBALA
  43. CMDK ENCBLK
  44. CMDK ENCBK1
  45. CMDK ENCBK2
  46. CMDK HANDBK
  47. CMDK HUMIDC
  48. CMDK INDIC1
  49. CMDK IWETHR
  50. CMDK IZZQ
  51. CMDK LDPOUT
  52. CMDK LPSTRB
  53. CMDK MZON1
  54. CMDK OUTPUT
  55. CMDK OWETHR
  56. CMDK QUAYCH
  57. CMDK QUAYLE
  58. CMDK SIMA
  59. CMDK SOILB
  60. CMDK SOLARB
  61. CMDK STRUCA
  62. CMDK SURFAR
  63. CMDK TEMP1
  64. CMDK TEMPB
  65. CMDK TIMEB
  66. CMDK TSTB2
  67. CMDK TSTATC
  68. CMDK WIBLK
  69.       CHARACTER*8 CDAY (7)
  70.       INTEGER J10(10),NDEV (3), NROUT(24)
  71.       REAL ATOUT (12, 24), C(5),  DQINT(24)
  72.       REAL DSOLRH(24), DTOB(24), DWHM(24), DWINT(24), DWSPED(24)
  73. C     REAL COMFTT(11,21)
  74.       CHARACTER*7 NL1
  75.       CHARACTER*3 NL2
  76.       REAL ODDBM(25), ODHUM(25), SLR(10)
  77. C
  78. C-----INITIALIZATIONS
  79.       DATA NDEV   /0, 1, -1/
  80.       DATA NL1/'SOLAR('/, NL2/')  '/
  81.       DATA IDAY/0/
  82. C
  83. C
  84. C-----EXTRACT SOLAR DATA AT BEGINNING OF EACH NEW DAY
  85.       IDAY=IDAY+1
  86.       CALL CURS(18)
  87.       CALL SOLRLD(STM,NDYRI,ICC,KM,IDAY,DST,NSEQW,IDBG)
  88. C
  89. C  DATOUT(ISURF,ETC....)
  90. C    ISURF=1,2  SOLAR ON ROOF ABOVE ATTICS 1 AND 2
  91. C        3,4,...2+NROOMS  SOLAR IN WINDOWS OF ZONES 1,2,...NROOMS
  92. C        2+NROOMS+1..... SOLAR ON OUTSIDE ENCLOSURE ELEMENTS
  93.       IQQ2=NIJ
  94.       IF(IQQ2.LT.1)IQQ2=1
  95.       DO 252 ISURF=1, IQQ2
  96.          DO 251 IHR=1, 24
  97.             SOLAR (IHR, ISURF)=DATOUT (ISURF, IHR)
  98. 251         CONTINUE
  99.          SOLAR (25, ISURF)=SOLAR (1, ISURF)
  100. 252   CONTINUE
  101.       SOLARH (25)=SOLARH (1)
  102.       IF (ISKIP.EQ.0 .AND. IGARAGT.GT.0) CALL GARAG (IGARAGT,NZONG,
  103.      +                                               SOLARH(1))
  104.       KDAYM=IDAYM
  105. C
  106. C     DTIME = DTIM1 (2 MIN) HERE
  107.       DTIME=DTIM1
  108. C     WETHR BLOCK 3
  109.       IF (NSEQW.EQ.0 .OR. NSEQW.EQ.1) THEN
  110. C--------CHECK FOR MEAN DAY CALC OR ADJUST NOMINAL DATA SET
  111. C        USED ONLY WITH CARRIER DATA
  112.          IF (IMEAN.LT.0) THEN
  113.             DO 260 J=1, 25
  114.                ODDB (J)=ODDBM (J)
  115. 260            CONTINUE
  116.          ELSE
  117. C-----------DETERMINE OUTDOOR TEMP PROFILE FOR MEAN, MAX, OR MIN DAY
  118. C           CALCULATE WET BULB TEMPERATURE AND SPECIFIC HUMIDITY
  119. C-----------FOR DATA FROM WEATHER FILE
  120.             DO 280 J=1, 25
  121. C--------------NDEV (1) = 0 DEVIATION FOR MEAN DAY CALCULATION
  122. C--------------NDEV (2) = +1 DEVIATION FOR MAX  DAY CALCULATION
  123. C--------------NDEV (3) = -1 DEVIATION FOR MIN  DAY CALCULATION
  124.                ODDB (J)=ODDBM (J)+NDEV (IDAYM)*TEMPDV
  125.                RODB=ODDB (J)
  126.                ODWB=(((C (5)*ODDB (J)+C (4))*ODDB (J)+C (3))*
  127.      1            ODDB (J)+C (2))*ODDB (J)+C (1)
  128. C                   GET SPECIC HUMIDITY, WSC, USING PSY FUNCTIONS     
  129.                PWSS=FPWS(ODWB)
  130.                WSS=FW22(PWSS)
  131.                WSC=FW35(RODB,ODWB,WSS)
  132.                ODHUM (J)=WSC
  133. 280            CONTINUE
  134.          END IF
  135.       ELSE IF (NSEQW.EQ.3) THEN
  136.          ODDB (1)=ODDB (25)
  137.          WSPEED (1)=WSPEED (25)
  138.          ODHUM (1)=ODHUM (25)
  139.          KDAYM=KDY
  140.          DO 285 IHR=2, 25
  141. C 
  142. C  ODDB  - OUTDOOR DRY-BULB TEMP, F
  143. C  ODWB  - OUTDOOR WET-BULB TEMP, F
  144. C  WSPEED- WIND SPEED, MILES/HR
  145. C              PWSS BELOW IS A DUMMY VARIABLE         
  146.             READ(41, *, IOSTAT=IO41) ODDB (IHR), ODWB, WSPEED (IHR),
  147.      1         PWSS, PWSS, KMO, KDY, KHR
  148.             IF(IO41.NE.0)THEN
  149.               WRITE(60,*)' STRDAY: TROUBLE READING TAPE41!'
  150.               STOP ' STRDAY: TROUBLE READING TAPE41!'
  151.               ENDIF
  152.             IF (IHR.LT.25 .AND. KMO.NE.KM) THEN
  153.                WRITE(60, *) 'NSEQW=3 WEATHER OUT OF SEQUENCE',
  154.      1         ' KMO,KM= ',   KMO, KM
  155.                WRITE(*, *) 'NSEQW=3 WEATHER OUT OF SEQUENCE',
  156.      1         ' KMO,KM= ',   KMO, KM
  157.               
  158.                STOP 'STRDAY: NSEQW=3 WEATHER OUT OF SEQUENCE'
  159.             END IF
  160.             RODB=ODDB (IHR)
  161. C                   GET SPECIFIC HUMIDITY, WSC, USING PSY FUNCTIONS            
  162.             PWSS=FPWS(ODWB)
  163.             WSS=FW22(PWSS)
  164.             WSC=FW35(RODB,ODWB,WSS)
  165.             ODHUM (IHR)=WSC
  166. 285         CONTINUE
  167.       END IF
  168.       IF (NSEQW.EQ.0 .OR. NSEQW.EQ.1 .OR. NSEQW.EQ.3) THEN
  169.          IF (INPRT.EQ.1) THEN
  170.             WRITE(60, 2210) KM, CMONTH (KM),
  171.      1         CDAY (IDAYM), WSHADE
  172.             DO 300 IHR=1, 24
  173.                WRITE(60, 2220) IHR, ODDB (IHR), ODHUM (IHR),
  174.      1            (SOLAR (IHR, J), J=1, 6), QINTL (IHR), WINTL (IHR)
  175. 300            CONTINUE
  176.             WRITE(60, *) 'SOLAR GAIN FOR WINDOWS WILL ',
  177.      1         'BE MULTIPLIED BY THE WINDOW SHADE FACTOR, WSHADE=',
  178.      2         WSHADE
  179.          END IF
  180.          TODAVG=0.
  181.          DO 320 IHR=1, 24
  182.             TODAVG=TODAVG+ODDB (IHR)
  183. 320         CONTINUE
  184.          TODAVG=TODAVG/24.
  185. C--------INTERPOLATE WEATHER DATA TO DESIRED TIME INCREMENTS
  186.          DO 325 IHR=1, 24
  187.             IHRP1=IHR+1
  188.             DTOB (IHR)=(ODDB (IHRP1)-ODDB (IHR))*DTIME
  189.             DWHM (IHR)=(ODHUM (IHRP1)-ODHUM (IHR))*DTIME
  190.             DWSPED (IHR)=(WSPEED (IHRP1)-WSPEED (IHR))*DTIME
  191.             DSOLRH (IHR)=(SOLARH (IHRP1)-SOLARH (IHR))*DTIME
  192. 325         CONTINUE
  193.       END IF
  194.       IF (INPRT.NE.0) THEN
  195.         WRITE(60, *) 'HORIZONTAL SOLAR FOR TODAY:'
  196.         WRITE(60, *) (IHR, SOLARH(IHR), IHR=1, 24)
  197.       ENDIF
  198.       IF (WSHADE.LT.1.0) THEN
  199. C--------APPLY WINDOW SHADE FACTOR
  200.          DO 330 IW=3,NROOMS+2
  201.          DO 330 IHR=1, 25
  202.             SOLAR (IHR,IW)=SOLAR (IHR,IW)*WSHADE
  203. 330         CONTINUE
  204.       END IF
  205.       IF (INPRT.NE.0 .AND. NSEQW.EQ.2) THEN
  206.          WRITE(60, 2211) WSHADE, (IHR, (SOLAR (IHR, J), J=1, 6),
  207.      1      QINTL (IHR), WINTL (IHR), IHR=1, 24)
  208.           IF(NIJ.GT.6)THEN
  209.               ISP=6
  210.   332         IST=ISP+1
  211.               ISP=ISP+10
  212.               IF(ISP.GT.NIJ)ISP=NIJ
  213.            WRITE(60,504)
  214.            I=0
  215.            DO 334 J=IST,ISP
  216.            I=I+1
  217.            J10(I)=J
  218.   334      CONTINUE
  219.            WRITE(60,500)(NL1,J10(IK),NL2,IK=1,I)
  220.            DO 337 IHR=1,24
  221.            I=0
  222.            DO 336 J=IST,ISP
  223.            I=I+1
  224.            SLR(I)=SOLAR(IHR,J)
  225.   336      CONTINUE
  226.            WRITE(60,501)(SLR(IK),IK=1,I)
  227.   337      CONTINUE
  228.            IF(ISP.NE.NIJ)GO TO 332
  229.            ENDIF
  230.       END IF
  231. C
  232.       IF (TSOILF (IDAYM).EQ.999.0) THEN
  233.          TODAG1=AFTG1*TODAVG+BFTG1
  234.       ELSE
  235.          TODAG1=TSOILF (IDAYM)
  236.       END IF
  237.       IF (TSOILW (IDAYM).EQ.999.0) THEN
  238.          TODAG2=AWTG2*TODAVG+BWTG2
  239.       ELSE
  240.          TODAG2=TSOILW (IDAYM)
  241.       END IF
  242.       IF(IGRND.EQ.0)THEN
  243.         DO 338 IQW=1,4
  244.   338   TSOIL(IQW)=TODAG2
  245.         TSOIL(5)=TODAG1
  246.         TSOIL(11)=TODAG1
  247.       ENDIF
  248.       DO 340 IHR=1, 24
  249.         DQINT (IHR)=(QINTL (IHR+1)-QINTL (IHR))*DTIME
  250.         DWINT (IHR)=(WINTL (IHR+1)-WINTL (IHR))*DTIME
  251.         DO 340 K=1, NIJ
  252.           DSOLAR (IHR, K)=(SOLAR (IHR+1, K)-SOLAR (IHR, K))*DTIME
  253. 340   CONTINUE
  254. C
  255. C--------CAL MISC. CONSTANTS
  256.       IF (ISKIP.EQ.0) THEN
  257.         AIRMSB=ROMVOL(NZNBAS)*RHOAIR
  258.         IF (IMODE.EQ.1) THEN
  259.           TAIN=TBAS
  260.           CALL ACOND (IDUM, TDUM, TDUM)
  261.         END IF
  262. C
  263. C-----------START OF RUN CALCULATIONS
  264. C-----------INITIALIZE VALUES
  265.         TODDB=TODAVG+TOFFSET
  266.         TGARAG=TODDB
  267.         IF(IGARAGT.EQ.0.AND.NZONG.GT.0)TGARAG=.66*TODDB+.34*TI(1)
  268.         NRDF1=0
  269. C-----------SET INITIAL VALUES FOR DUCT MODEL
  270.         TBAS=TBASI
  271.         IF (TBAS.GT.200.0) TBAS=TI(1)
  272.         TBASA=TBAS
  273.         QFLBAS=0.0
  274.         QBASVN=0.0
  275.         TBASN=TBAS
  276.         TCRW=(TODDB+TI(1))*.5*ICRW
  277.         TCRWA=TCRW
  278.         QFLCRW=0.0
  279.         QCRWVN=0.0
  280.         QGRND=0.0
  281.         QVENTA=0.0
  282.         TCRWN=TCRW
  283. C-----------SET INITIAL VALUES FOR WALL AND ROOF MODELS
  284.         TATTC(1)=(TI(1)+TODDB)*.5
  285.         TATTA(1)=TATTC(1)
  286.         TOTATC=(TODDB+TATTA(1))*.5
  287. C  SET INITIAL TEMPS OF ROOF LAYERS
  288. C   --COMMENTED INITIALIZATIONS WERE USED IN SINGLE-ZONE SP43 MODEL
  289.         TMP(1,IZT+1)=TODDB
  290.         TMP(2,IZT+1)=TOTATC
  291. C        TMP(2,IZT+1)=TATTC(1)
  292.         TMP(3,IZT+1)=TATTA(1)
  293. C        TMP(3,IZT+1)=TATTC(1)
  294.         TMP(4,IZT+1)=TATTA(1)
  295. C        TMP(4,IZT+1)=TI(1)
  296.         IF(NATC.EQ.2)THEN
  297.           TMP(1,IZT+2)=TODDB
  298.           TMP(2,IZT+2)=TOTATC
  299.           TMP(3,IZT+2)=TATTA(1)
  300.           TMP(4,IZT+2)=TATTA(1)
  301.         ENDIF
  302. C  SET INITIAL TEMPS OF ENCLOSURE ELEMENTS
  303.         WRITE(60, 1000)
  304. C
  305.         DO 359 I=1,2
  306.         QCEIL(I)=0.
  307.         QRFSUM(I)=0.
  308.  359    CONTINUE        
  309.         DO 360 J=1,NROOMS
  310.         NRM=NRMA(J)
  311.         TFM(NRM)=TI(NRM)
  312.         EBAL(NRM)=0.
  313.         EBALT(NRM)=0.
  314.         QLOSWL(NRM)=0.
  315.         QLOSFL(NRM)=0.
  316. C                            SET INITIAL TEMPS OF WALL AT THERMOSTAT
  317.         TWSTAT(NRM)=TI(NRM)
  318.         N2=NWALLA(NRM)
  319. C
  320.         DO 360 K2=1,N2
  321.         IZW=NENC(NRM,K2)
  322.         ID=IDEXP(IZW)
  323.         NZONC=NZNC(IZW)
  324.         NZON=NZONC-8
  325.         KON=KONSTA(IZW)
  326.         NLM=NLMP(KON)
  327.         IF(ID.LE.4)THEN
  328. C  OUTSIDE WALLS
  329.           TODTI=(TODDB+TI(NRM))*0.5
  330.           TMP(1,IZW)=TODDB
  331.           TMP(2,IZW)=TODTI
  332.           TMP(3,IZW)=TI(NRM)
  333.           TMP(4,IZW)=TI(NRM)
  334.         ENDIF
  335.         IF(ID.EQ.5)THEN
  336. C  CEILING
  337.           IF(NZONC.EQ.4)IA=1
  338.           IF(NZONC.EQ.5)IA=2
  339.           IF(NZONC.LT.8)THEN
  340.             TWALO=TATTA(IA)
  341.             TMEAN=(TATTA(IA)+TI(NRM))*0.5
  342.           ELSE
  343.             TWALO=TI(NZON)
  344.             TMEAN=(TI(NZON)+TI(NRM))*0.5
  345.           ENDIF
  346.           TMP(1,IZW)=TWALO
  347.           TMP(2,IZW)=TMEAN
  348.           TMP(3,IZW)=TI(NRM)
  349.           TMP(4,IZW)=TI(NRM)
  350.         ENDIF
  351.         IF(ID.EQ.6)THEN
  352. C  FLOOR
  353.           IF(NZONC.EQ.2.OR.NZONC.EQ.7)THEN
  354. C     LVG-SPC FLOOR W/ BSMT BELOW
  355.             TBTI=(TBAS+TI(NRM))*.5
  356.             TMP(1,IZW)=TBAS
  357.             TMP(2,IZW)=TBTI
  358.             TMP(3,IZW)=TI(NRM)
  359.           ENDIF
  360.           IF(NZONC.EQ.3.OR.NZONC.EQ.8)THEN
  361. C    LVG-SPC FLOOR W/ CRAWL SPC BELOW
  362.             TCTI=(TCRW+TI(NRM))*.5
  363.             TMP(1,IZW)=TCRW
  364.             TMP(2,IZW)=TCTI
  365.             TMP(3,IZW)=TI(NRM)
  366.           ENDIF
  367.         ENDIF
  368.         IF(ID.EQ.7)THEN
  369. C  INTERIOR PARTION WALLS OR LVG SPC-GARAGE WALL
  370.           TOUT=TI(NRM)
  371.           IF(NZONC.EQ.6)THEN
  372.             TOUT=TGARAG
  373.           ELSE IF(NZON.GT.0)THEN
  374.             TOUT=TI(NZON)
  375.           ENDIF
  376.           TMID=0.5*(TOUT+TI(NRM))
  377.           TMP(1,IZW)=TOUT
  378.           TMP(2,IZW)=TMID
  379.           TMP(3,IZW)=TI(NRM)
  380.           IF(NLM.NE.1)TMP(4,IZW)=TI(NRM)
  381.           IF(IZW.EQ.20.OR.IZW.EQ.21.OR.IZW.EQ.25)THEN
  382.             WRITE(60,*)' STRDAY: IZW,NRM,TMP(1,IZW)= ',IZW,NRM,
  383.      +                 TMP(1,IZW)
  384.             WRITE(60,*)'    TMP(2-4,IZW)= ',TMP(2,IZW),TMP(3,IZW),
  385.      +                 TMP(4,IZW)
  386.             WRITE(60,*)'    TGARAG,TI= ',TGARAG,TI(NRM)
  387.           ENDIF
  388.         ENDIF
  389.         IF(ID.EQ.8)THEN
  390. C                    BASEMENT FLOOR
  391.           IF(IGRND.EQ.0)THEN
  392.             TMID=0.5*(TODAG1+TBAS)
  393.             TMP(1,IZW)=TODAG1
  394.             TMP(2,IZW)=TMID
  395.             TMP(3,IZW)=TBAS
  396.             IF(NLM.NE.1)THEN
  397.               TMP(3,IZW)=TMID
  398.               TMP(4,IZW)=TBAS
  399.             ENDIF
  400.           ENDIF
  401.         ENDIF
  402.         IF(ID.EQ.9)THEN
  403. C   BELOW-GROUND BSMT WALL
  404.           IF(IGRND.EQ.0)THEN
  405.             TMID=0.5*(TODAG2+TBAS)
  406.             TMP(1,IZW)=TODAG2
  407.             TMP(2,IZW)=TMID
  408.             TMP(3,IZW)=TMID
  409.             TMP(4,IZW)=TBAS
  410.           ENDIF
  411.         ENDIF
  412.         IF(ID.EQ.10)THEN
  413. C   ABOVE-GROUND BSMT WALL
  414.           TINSID=.6*(TBAS-TODDB)+TODDB
  415.           TMID=0.5*(TODDB+TINSID)
  416.           TMP(1,IZW)=TODDB
  417.           TMP(2,IZW)=TMID
  418.           TMP(3,IZW)=TINSID
  419.           IF(NLM.NE.1)THEN
  420.             TMP(3,IZW)=TMID
  421.             TMP(4,IZW)=TINSID
  422.           ENDIF
  423.         ENDIF
  424.         WRITE(60,1010)NRM,IZW,ID,(TMP(I,IZW),I=1,4)
  425.   360   CONTINUE
  426.         WRITE(60,1014)(TMP(I,IZT+1),I=1,4)
  427.         IF(NATC.EQ.2)WRITE(60,1016)(TMP(I,IZT+2),I=1,4)
  428.         TRO=TODDB
  429. C-----------DETERMINE INITIAL STATES OF ALL CONTROLS
  430.         CALL TSTATD (3, IAC, IAC1, IBURN, IEFAN,IDBG)
  431. C--------CALCULATE DUCT MODEL CONFIGURATION CONSTANTS
  432.         DO 472 I=1, NZN7
  433.         TZONE (I)=TI(1)
  434. 472     CONTINUE
  435.         CALL DUCTIN(IDBG)
  436. C            
  437.         IF(IMODE.EQ.1.OR.(IMODE.EQ.2.AND.CFMBP.GT.0.))THEN
  438. C   SET INITIAL ABSOLUTE HUMIDITY IN ZONES
  439.           TIM=AVGWV(TI,NROOMS,RWSDCT)
  440.           IPSY=-3
  441.           CALL PSY (IPSY, TIM, TDUM, TDUM, RELHMI,WIM, HDUM)
  442.         ELSE
  443.           WIM=.01
  444.         ENDIF
  445.         DO 361 JT=1,NROOMS
  446.         IPSY= -3
  447.         CALL PSY(IPSY,TI(JT),TDUM,TDUM,RELHMI,WI(JT),HDUM)
  448.   361   CONTINUE        
  449.         RHI=RELHMI
  450.         RHS=RHI
  451.         RHSN=RHI
  452.         TIO=TI(1)
  453.         RHIO=RHI
  454.         WBAS=WI(NZNBAS)
  455.         WBASN=WBAS
  456.         WAC=WIM
  457.         TPL=TI(1)
  458. C----------SET HUMIDIFIER STATES  
  459.         IF(IMODE.EQ.2.AND.CFMBP.GT.0.)CALL HUMDF1 (TI(1))
  460.         FINFB=0.
  461.         TDIR=1.
  462.         IF (IBURN.EQ.1) TDIR=-1.
  463.         RHDIR=100.
  464.         IF (IBURN.EQ.1) RHDIR=-100.
  465. C--------GO THRU NDPRE DAYS BEFORE TAKING DATA (ONLY ONCE WITH SEQ WTHR)
  466.         IF (NSEQW.EQ.0) IREP=0
  467.       END IF
  468. C  * * * * *    END OF MISC CONSTANTS CALCULATION  (W/ ISKIP = 0)
  469. C-----INITIALIZE COUNTERS TO ZERO
  470. C      DO 431 J=1, 21
  471. C         DO 430 I=1, 11
  472. C            COMFTT (I, J)=0.
  473. C 430         CONTINUE
  474. C 431      CONTINUE
  475.       DO 453 IHR=1, 24
  476.          NROUT (IHR)=0
  477.          DO 440 J=1, 12
  478.             ATOUT (J, IHR)=0.
  479. 440         CONTINUE
  480.       DO 441 J=1,3
  481.       HROT(J,IHR)=0.
  482.       HPRT(J,IHR)=0.
  483. 441   CONTINUE
  484.          DO 450 J=1, NHRP
  485.             HRPRT (J, IHR)=0.
  486. 450         CONTINUE
  487.          DO 451 J=1, NHRO
  488.             HROUT (J, IHR)=0.
  489. 451         CONTINUE
  490. 453      CONTINUE
  491.       DO 452 I=1, NHBK
  492.          HBK (I)=0.0
  493. 452      CONTINUE
  494.       IHUM=(1.-.0001*(RHS-SETHS))
  495.       IF (CFMBP.LE.0.0) IHUM=0
  496.       WATER=0.0
  497.       QSPT=0.0
  498.       QWAS=0.0
  499.       QHUMT=0.0
  500.       HUMON=0.0
  501.       ICYCCP=0
  502.       ICYC1=0
  503.       ICYC6=0
  504.       CLOSE(11)
  505.       OPEN(11,FILE='TAPE11',STATUS='UNKNOWN')
  506.       ISWING=0
  507.       CLOSE(1)
  508.       OPEN(1,FILE='TAPE1',STATUS='UNKNOWN')
  509.       IRH=0
  510.       INC=0
  511.       LINES=0
  512.       TONMAX=0.0
  513.       TIMAX=0.0
  514.       TIMAXN=0.0
  515.       TIMIN=999.0
  516.       MPRINT=1
  517.       IF (NINC.EQ.0) MPRINT=0
  518.       IF ((MPRINT.NE.0) .OR. (IDEBUG.GT.0)) WRITE(60, 2270)
  519.       IF (ISKIP.EQ.0) THEN
  520.          VENCRV=VENCRW*VOLCRW*CPAIR
  521.          DTCR=0.
  522.          IF (ICRW.EQ.1) DTCR=DTIM2/(CPAIR*VOLCRW)
  523.       END IF
  524. C-----CHANGE TIME INCREMENT TO INNER LOOP VALUE
  525.       DTIME=DTIM2
  526.       IF (ISKIP.EQ.0) THEN
  527. C
  528.          TWRT=TBAS
  529.          TWRTB=TBAS
  530.          TWSP=TBAS
  531.          TWSPB=TBAS
  532. C
  533. C--------CALC TIME CONST VALUES
  534.          CALL TSTATD (4, IAC, IAC1, IBURN, IEFAN, IDBG)
  535.          TCRH=0.0
  536.          IF (CFMBP.GT.0.0) TCRH=DTIME/HSTC
  537.          ENDT=HRDY - .5*DTIME
  538.       END IF
  539.       RETURN
  540. C
  541.   500 FORMAT(1X,10(A7,I2,A3))
  542.   501 FORMAT(1X,(10F12.4))
  543.   504 FORMAT(1X///)
  544. 1000  FORMAT (/' INITIAL STRUCTURE TEMPERATURES...')
  545. 1010  FORMAT (1X,'NRM,IZW,ID= ',3I5,' NODE TEMPS(1-4)= ',4G13.5)
  546.  1014 FORMAT(15X,'ROOF TEMPS(ATTIC NO 1)= ',4G13.5)
  547.  1016 FORMAT(15X,'ROOF TEMPS(ATTIC NO 2)= ',4G13.5)
  548. 1020  FORMAT (' ERROR...DUPLICATE GARAGE MODELS SPECIFIED')
  549. C-----INPUT DATA PRINTING
  550. 2210  FORMAT (///'  WEATHER DATA FOR MONTH', I3, '=', A3, 5X, A8,
  551.      1   ' (WSHADE=', F5.2, ')', //,
  552.      2   30X, 4 ('WALL', 8X), 'ROOF', 8X, 'WINDOW', /,
  553.      3   ' HOUR DRY BULB    HUMIDITY    SOLAR(1)    SOLAR(2)    ',
  554.      4   'SOLAR(3)    SOLAR(4)    SOLAR(5)    SOLAR(6)    SEN. LOAD   ',
  555.      5   'LAT. LOAD', /)
  556. 2211  FORMAT (/' SOLAR DATA AND INTERNAL LOADS. (WSHADE=', F5.2, ')', //
  557.      1   8X, 'SOLAR(1)    SOLAR(2)    SOLAR(3)    ',
  558.      2   'SOLAR(4)    SOLAR(5)    SOLAR(6)    SEN. LOAD   LAT. LOAD', /,
  559.      3   (I4, 1X, 8G12.4))
  560. C 2212 FORMAT(/' SOLAR(',I2,')   SOLAR(',I2,')  SOLAR(',I2,')  SOLAR(',
  561. C    +I2,')   SOLAR(',I2,')  SOLAR(',I2,')   SOLAR(',I2,')  SOLAR(',
  562. C    +I2,')   SOLAR(',I2,')'/(2X,10G12.4))
  563. 2220  FORMAT (I4, 1X, 10G12.4)
  564. C-----ITERATION PRINTING
  565. 2270  FORMAT (/// T6, 'TIME', T15, 'TODDB', T28, 'TI(1)', T37, 'RHI',
  566.      1   T48, 'TPL', T56, 'TSUP', T65, 'TRETB', T76, 'TBAS', T85,
  567.      2   'TATTC', T95, 'TFLUE', 3X, 'IAC IBURN IBLOW IHON EXCESS', //)
  568.       END
  569.