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 / TSTAT.FOR < prev    next >
Text File  |  1992-05-08  |  18KB  |  531 lines
  1.       SUBROUTINE TSTAT(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  2. C
  3. C  THERMOSTAT MODEL
  4. C    IACZP - A/C ON-OFF SIGNAL IN EACH ZONE
  5. C    IBURN
  6. C
  7. C  I - ZONE NO.
  8. C  IHONEY  = 0  CARRIER THERMOSTAT
  9. C           1  HONEYWELL ROUND THERMOSTAT
  10. C           2  T8200 HONEYWELL THERMOSTAT(GAS FURNACE)
  11. C           3  T8300 HONEYWELL THERMOSTAT(HEAT PUMP W/ SETBACK)
  12. C
  13. C   MODE = 1  AUTO HTG AND CLG
  14. C          2  CLG ONLY
  15. C          3  HTG ONLY
  16. C          4  NON AUTOMATIC HTG AND CLG
  17. C
  18. C   TIME  MUST BE TIME OF DAY(HRS) AND IS USED FOR SETBACK CALC
  19. C  REVISION:  ADDED T8300 MODEL ON 8/14/85
  20. C           (INTEGRAL ACTION ON HTG SETBACK RECOVERY NOT YET MODELED)
  21. C  1/11/91 RDF REDUCED MEMORY REQMT BY REVISING SETBACK CODE AND
  22. C              ELIMINATING RUNNING AVG OPTION WHICH WAS NEVER USED
  23. C  3/31/92 RDF REVISED THE MULTIZONE VERSION FROM TRNSYS-LIKE W/
  24. C              CONSTANTS IN CON ARRAY, REAL AND INTEG VAR IN VAR AND 
  25. C              IVAR ARRAY TO SIMPLY DIMENSIONING ALL TERMS IN SINGLE-
  26. C              ZONE TSTAT.  NOW TSTATS IN EACH ZONE CAN BE UNIQUE 
  27. C              AND ANY LATENT ERROR SHOULD BE ELIMINATED.
  28. C             
  29.       PARAMETER (SMALL=1.E-37)
  30. CMDK NWL
  31. CMDK NZN
  32. CMDK NZW
  33. CMDK BIDRUN
  34. CMDK BLKTS1
  35. CMDK BLSTBK
  36. CMDK BLKTST
  37. CMDK COLLEC
  38. CMDK DEBG
  39. CMDK ENCBK1
  40. CMDK TIMEB
  41. CMDK TSTATC
  42.       DIMENSION ANTC1(NZN),ANTH1(NZN),ANTH2(NZN),CNINTG(NZN),
  43.      + DIFFC1(NZN),DIFFH1(NZN),DIFFH2(NZN),IACZP(NZN),
  44.      + IACN1(NZN),IBRNP(NZN),IBUR1N(NZN),IBUR2N(NZN),IBURN1Z(NZN),
  45.      + IBURN2(NZN),IHONEY(NZN),SBKCOF(NZN),SBKCON(NZN),
  46.      + SBKHOF(NZN),SBKHOF2(NZN),SBKHON(NZN),SBKHON2(NZN),SETBC(NZN),
  47.      + SETBH(NZN),SETBH2(NZN),SETC1(NZN),SETH1(NZN),SETH2(NZN),
  48.      + SETHON2(NZN),TAC1(NZN),TAC1N(NZN),TAH1(NZN),TAH2(NZN),
  49.      + TATC1(NZN),TATH1(NZN),TATH2(NZN),TCTAC1(NZN),
  50.      + TCTAH1(NZN),TCTAH2(NZN),TCTC1(NZN),TCTH1(NZN),TCTH2(NZN),
  51.      + TCTWC1(NZN),TCTWH1(NZN),TCTWH2(NZN),TSC1(NZN),TSC1N(NZN),
  52.      + TSCT1(NZN),TSH1(NZN),TSH1N(NZN),TSH2(NZN),TSH2N(NZN),
  53.      + TSHT1(NZN),TSHT2(NZN),TSWC1(NZN),TSWC1N(NZN),TSWH1(NZN),
  54.      + TSWH1N(NZN),TSWH2(NZN),TSWH2N(NZN),TWCT1(NZN),TWHT1(NZN),
  55.      + TWHT2(NZN),WALCPG(NZN)
  56.       DIMENSION ICNT(NZN),AIRCPG(NZN),XINTG(NZN),IAC1(NZN),IACQ(NZN)
  57. C
  58.       LOGICAL NDEBUG
  59. C   ---   ____  ----   PARAMETERS IN NAMELIST
  60. C  IHONEY  - THERMOSTAT TYPE(SEE ABOVE)
  61. C  WALCPG  - WALL COUPLING FRACTION
  62. C  IZNCOLL - ZONE IN WHICH TO COLLECT TSTAT SENSOR TEMP IN CDATA(64)
  63. C  CNINTG   INTEGRATOR CONSTANT USED TO MINIMIZE DROOP IN T8300 TSTAT
  64. C  TWCT1,TWHT1,TWHT2  WALL TIME CONST,HR
  65. C  TATC1,TATH1,TATH2   ANTICIPATOR TIME CONST,HR
  66. C  TAH1,TAH2   INITIAL VALUES OF ANTICIPATOR TEMP,F
  67. C ***   TSXTN  - SENSOR TIME CONSTANT (FOR AIR SENSOR IN T8200 AND 
  68. C                T8300 TSTAT,HR
  69. C ***   ANTXN  - ANTICIPATION
  70. C ***   SETXN  - SET POINT
  71. C ***   DIFFXN - DIFFERENTIAL
  72. C ***   SETBX  - AMOUNT OF SET-BACK(POSITIVE FOR SETBACK,NEG FOR SETUP)
  73. C ***   SBKXON - SET-BACK START TIME (0. TO 24. HR)
  74. C ***   SBKXOF - SET-BACK STOP TIME (0. TO 24. HR)
  75.       NAMELIST/TSTATN/
  76.      + IHONEY,WALCPG,IZNCOLL,CNINTG,
  77.      + TSCT1,TWCT1,TATC1,ANTC1,TAC1,
  78.      + DIFFC1,SETC1,SETBC,SBKCON,SBKCOF,
  79.      + TSHT1,TWHT1,TATH1,ANTH1,TAH1,
  80.      + DIFFH1,SETH1,SETBH,SBKHON,SBKHOF,
  81.      + TSHT2,TWHT2,TATH2,ANTH2,TAH2,
  82.      + DIFFH2,SETH2,SETBH2,SBKHON2,SBKHOF2
  83. C     
  84.       DATA WALCPG,TWCT1,TATC1,TWHT1,TATH1,TWHT2,TATH2/7*NZN*0./
  85.       DATA TAC1N/NZN*0./
  86.       DATA IBRNI/1/,IBLW/0/
  87.       DATA NPRMAX/500/
  88.       DATA TAH1,TAH2/2*NZN*0./
  89.       DATA CNINTG/NZN*0./
  90.       DATA IRESETBK/0/
  91.       NDEBUG=.FALSE.
  92.       NPR=0
  93.       OPEN(15,FILE='TAPE15',STATUS='OLD',IOSTAT=IO15)
  94.       IF(IO15.NE.0)THEN
  95.         WRITE(60,*) ' TSTAT: CANT OPEN TAPE15 WITH HONEYW NAMELIST'
  96.         STOP ' TSTAT: CANT OPEN TAPE15'
  97.       ENDIF
  98.       READ(15,TSTATN)
  99.       WRITE(60,TSTATN)     
  100.       CLOSE(15)
  101.       DO 2 J=1,NZN
  102.       ICNT(J)=0
  103.       IF(CNINTG(J).LT.0.)THEN
  104.         CNINTG(J)=ABS(CNINTG(J))
  105.         ICNT(J)=1
  106.       ENDIF
  107. C  - - -  INPUT CNINTG(J) AS NEG IF WANT TO BYPASS INTEG ON DIF1
  108. C          (USED ON RAMP UP TEST TO CHECK OPERATION OF TSTAT ON CLG)
  109.       IF (SETC1(J) .LT. SETH1(J)+DIFFH1(J)) THEN
  110.         WRITE(60,3020)
  111.         STOP ' ***OVERLAP BETWEEN A/C & H/P THERMOSTATS - RUN STOPPED'
  112.       ENDIF
  113.       WRITE(60,*)' THERMOSTAT INPUT DATA FOR ZONE ',J
  114.       WRITE(60,2170) TSHT1(J),ANTH1(J),SETH1(J),DIFFH1(J),SETBH(J),
  115.      +               SBKHON(J),SBKHOF(J)
  116.       WRITE(60,2180) TSHT2(J),ANTH2(J),SETH2(J),DIFFH2(J),SETBH2(J),
  117.      +               SBKHON2(J),SBKHOF2(J)
  118.       WRITE(60,2190) IHONEY,WALCPG,TWCT1(J),TWHT1(J),TWHT2(J),
  119.      +               TATC1(J),TATH1(J),TATH2(J)
  120.       AIRCPG(J)=1.-WALCPG(J)
  121.       IF(ITHERM.EQ.0)AIRCPG(J)=1.
  122.     2 CONTINUE      
  123.       IMODE=1
  124.       IF(MODE.EQ.3)IMODE=2
  125.       RETURN
  126. C 
  127.       ENTRY TSTAT1(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  128. C *** RESET IMODE FOR SPECIFIED HTG AND CLG MONTHS, IF ANY
  129.       NPR=0
  130.       IF ((KM.GE.MONCS).AND.(KM.LE.MONCE)) IMODE=1
  131.       IF ((KM.GE.MONHS).OR.(KM.LE.MONHE)) IMODE=2
  132.       RETURN
  133. C
  134.       ENTRY TSTAT2(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  135. C *** DETERMINE INITIAL STATES OF ALL CONTROLS
  136. C
  137.       TT=TI
  138.       IF(ITHERM.GT.0)THEN
  139.         IF(WALCPG(I).EQ.0.)THEN
  140.           TT=.5*(TI+TWALL)
  141.         ELSE
  142.           TT=TI*AIRCPG(I)+TWALL*WALCPG(I)
  143.         ENDIF
  144.       ENDIF
  145.       TSC1(I)=TT
  146.       TSH1(I)=TT
  147.       TSH2(I)=TT
  148.       TAC1(I)=0.
  149.       TDIF1=0.
  150.       IF(IHONEY(I).LT.2)GO TO 146
  151.       TIA=TI*AIRCPG(I)
  152.       TTPW=TWALL*WALCPG(I)
  153.       TSC1(I)=TIA
  154.       TSH1(I)=TIA
  155.       TSH2(I)=TIA
  156.       TSWC1(I)=TTPW
  157.       TSWH1(I)=TTPW
  158.       TSWH2(I)=TTPW
  159.       TSC1N(I)=TSC1(I)
  160.       TSH1N(I)=TSH1(I)
  161.       TSH2N(I)=TSH2(I)
  162.       TSWC1N(I)=TSWC1(I)
  163.       TSWH1N(I)=TSWH1(I)
  164.       TSWH2N(I)=TSWH2(I)
  165.   146 IAC1(I)=1-IACZP(I)
  166.       IACNZ(I)=IACZP(I)
  167.       IACN1(I)=1-IACNZ(I)
  168.       IF(MODE.EQ.2) THEN
  169.         IBRNP(I)=0
  170.         IBURNNZ(I)=0
  171.         IBUR1N(I)=0
  172.         IBUR2N(I)=0
  173.         IBURN1Z(I)=0
  174.         IBURN2(I)=0
  175.       ELSE
  176.         IBRNP(I)=IBRNI
  177.         IBURNNZ(I)=IBRNI
  178.         IBUR1N(I)=IBRNI
  179.         IBUR2N(I)=IBRNI
  180.         IBURN1Z(I)=IBRNI
  181.         IBURN2(I)=IBRNI
  182.       ENDIF
  183.       IEFAN=IBLW
  184.       CALL SETBAK(TIME,IMODE,SETH1(I),SETH2(I),SETBH(I),SBKHON(I),
  185.      +            SBKHOF(I),SETBH2(I),SBKHON2(I),SBKHOF2(I),TSSH1,
  186.      +            TSSH2,SETC1(I),SETBC(I),SBKCON(I),SBKCOF(I),TSSC1)
  187.       GO TO (390,390,400,410), MODE
  188.   390 TSUM=TSC1(I)
  189.       IF(IHONEY(I).EQ.1)TSUM=TSC1(I)+TAC1(I)
  190.       IF(IHONEY(I).GE.2)THEN
  191.         TSENS=TSC1N(I)+TSWC1N(I)
  192.         TDIF=TSENS-TSSC1
  193.         TSM=TDIF+TAC1N(I)
  194.         GO TO 392
  195.       ENDIF
  196.       TSM=TSUM-TSSC1
  197.   392 IF(TSM.GE.DIFFC1(I))IACZP(I)=1
  198.       IF(TSM.LE.0.)IACZP(I)=0
  199.       IAC1(I)=1-IACZP(I)
  200.       IF(MODE.GT.1) GO TO 420
  201.  400  CONTINUE
  202. C ***   FURNACE CANNOT BE ON IF A/C IS ON . . .
  203.       TSUM1=TSH1(I)
  204.       IF(IHONEY(I).EQ.1)TSUM1=TSH1(I)+TAH1(I)
  205.       IF(IHONEY(I).GE.2)THEN
  206.         TSENS=TSH1N(I)+TSWH1N(I)
  207.         TDIF=TSENS-TSSH1
  208.         TSM=TDIF+TAH1N(I)
  209.         GO TO 402
  210.       ENDIF
  211.       TSM=TSUM1-TSSH1
  212.   402 IF(TSM.GE.DIFFH1(I))IBURN1Z(I)=0
  213.       IF(TSM.LE.0.)IBURN1Z(I)=IAC1(I)
  214.       TSUM2=TSH2(I)
  215.       IF(IHONEY(I).EQ.1)TSUM2=TSH2(I)+TAH2(I)
  216.       IF(IHONEY(I).GE.2)THEN
  217.         TSENS=TSH2N(I)+TSWH2N(I)
  218.         TDIF=TSENS-TSSH2
  219.         TSM=TDIF+TAH2N(I)
  220.         GO TO 404
  221.       ENDIF
  222.       TSM=TSUM2-TSSH2
  223.   404 IF(TSM.GE.DIFFH2(I))IBURN2(I)=0
  224.       IF(TSM.LE.0.)IBURN2(I)=IAC1(I)
  225.       IACP= (1-IBURN1Z(I))*IBURN2(I)
  226.       IBRNP(I)= IBURN1Z(I)+IACP
  227.       GO TO 420
  228. 410   IF((KM.GE.MONCS) .AND. (KM.LE.MONCE)) GO TO 390
  229.       IF((KM.GE.MONHS) .OR. (KM.LE.MONHE)) GO TO 400
  230.   420 CONTINUE
  231.       RETURN
  232. C
  233.       ENTRY TSTAT3(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  234. C  ***  COMPUTE TIME CONSTANT TERMS
  235.       TCTC1(I) = DTIME/TSCT1(I)
  236.       TCTH1(I) = DTIME/TSHT1(I)
  237.       TCTH2(I)= DTIME/TSHT2(I)
  238.       IF(IHONEY(I)-1)430,422,422
  239.   422 TCTAC1(I)=DTIME/TATC1(I)
  240.       TCTAH1(I)=DTIME/TATH1(I)
  241.       TCTAH2(I)=DTIME/TATH2(I)
  242.       IF(IHONEY(I).EQ.1)GO TO 430
  243.       TCTWC1(I)=DTIME/TWCT1(I)
  244.       TCTWH1(I)=DTIME/TWHT1(I)
  245.       TCTWH2(I)=DTIME/TWHT2(I)
  246.       XINTG(I)=CNINTG(I)*DTIME
  247.   430 CONTINUE
  248.       RETURN
  249. C
  250.       ENTRY TSTAT4(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  251. C
  252. C FOR RUNNING 121486 BELOW!!!!!(SEE ALSO TSTATD)
  253.       IF(IDRUN.EQ.'1214'.AND.IRESETBK.EQ.0)THEN
  254.         IF(TIME.GT.12.00)THEN
  255. C  SET BACK IN ZONES 1 AND 3 AFTER NOON OF 1ST DAY OF 121486 RUN!
  256.           IRESETBK=1        
  257.           SBKHON(1)=0.
  258.           SBKHON(3)=0.
  259.           SETBH(1)=40.
  260.           SETBH(3)=40.
  261.           SBKHOF(1)=24.
  262.           SBKHOF(3)=24.
  263.           SETBH2(1)=40.
  264.           SETBH2(3)=40.
  265.           SBKHOF2(1)=24.
  266.           SBKHOF2(3)=24.  
  267.         ENDIF
  268.       ENDIF
  269. C *** FIND NEW SET POINTS FOR ALL THERMOSTATS
  270.       CALL SETBAK(TIME,IMODE,SETH1(I),SETH2(I),SETBH(I),SBKHON(I),
  271.      +            SBKHOF(I),SETBH2(I),SETHON2(I),SBKHOF2(I),TSSH1,
  272.      +            TSSH2,SETC1(I),SETBC(I),SBKCON(I),SBKCOF(I),TSSC1)
  273. C      
  274.       IAC1(I)=1-IACZP(I)
  275.       TT = TI
  276.       IF(ITHERM.GT.0.)THEN
  277.         IF(WALCPG(I).EQ.0.)THEN
  278.           TT=.5*(TI+TWALL)
  279.         ELSE
  280.           TT=TI*AIRCPG(I)+TWALL*WALCPG(I)
  281.         ENDIF
  282.       ENDIF
  283.       IF(IMODE.GT.1)GOTO 690
  284. C  ****  COOLING
  285.       IF(IHONEY(I).EQ.0)THEN
  286. C             CARRIER THERMOSTAT      
  287.         TSC1N(I) = TCTC1(I)*(TT-TSC1(I)+ANTC1(I)*IAC1(I))+TSC1(I)
  288.         TSUM=TSC1N(I)
  289.         TSM=TSUM-TSSC1
  290.       ELSEIF(IHONEY(I).EQ.1)THEN
  291. C             HONEYWELL ROUND THERMOSTAT      
  292.         TSC1N(I)=TCTC1(I)*(TT-TSC1(I))+TSC1(I)
  293.         TAC1N(I)=TCTAC1(I)*(ANTC1(I)*IACZP(I)-TAC1(I))+TAC1(I)
  294.         TSUM=TSC1N(I)+TAC1N(I)
  295.         TSM=TSUM-TSSC1
  296.       ELSE
  297. C              HONEYWELL T8200 AND T8300 THERMOSTATS
  298.         TSC1N(I)=TCTC1(I)*(TI*AIRCPG(I)-TSC1(I))+TSC1(I)
  299.         TSWC1N(I)=TCTWC1(I)*(TWALL*WALCPG(I)-TSWC1(I))+TSWC1(I)
  300.         TAC1N(I)=TCTAC1(I)*(ANTC1(I)*IACZP(I)-TAC1(I))+TAC1(I)
  301.         TSENS=TSC1N(I)+TSWC1N(I)
  302.         TDIF=TSENS-TSSC1
  303. C  INTEGRATE TEMP DIF TO MINIMIZE DROOP BY INCREASING CYCLE RATE
  304.         IF(IHONEY(I).EQ.2)THEN
  305.           TDIF1=0.
  306.         ELSE
  307.           IF(ICNT(I).EQ.0)THEN
  308.             TDIF1=TDIF*XINTG(I)+TDIF1
  309.             IF(TDIF1.GT.3.) TDIF1=3.
  310.             IF(TDIF1.LT.0.) TDIF1=0.
  311.           ELSE
  312.             TDIF1=0.
  313.           ENDIF
  314.         ENDIF
  315.         TSM=TDIF1+TAC1N(I)+TDIF
  316.       ENDIF
  317.       IF(ICOLL.NE.0.AND.I.EQ.IZNCOLL)CDATA(64) = TSM
  318.       IF(TSM.GE.DIFFC1(I)) THEN
  319.         IACNZ(I)=1
  320.         IBURNNZ(I)=0
  321.       ENDIF
  322.       IF(TSM.LE.0.)IACNZ(I)=0
  323.       IEFAN=IACNZ(I)
  324.       IF(NDEBUG) WRITE(60,510) TDIF,TAC1N(I),TDIF1
  325.       IF(MODE.GT.1)GO TO 710
  326.       IACN1(I)=1-IACNZ(I)
  327. C  ***  HEATING
  328.   690 IF(IHONEY(I).EQ.0)THEN
  329. C              CARRIER THERMOSTAT EQUATIONS 
  330.         TSH1N(I) = TCTH1(I)*(TT-TSH1(I)+ANTH1(I)*IBURN1Z(I))+TSH1(I)
  331.         TSH2N(I)=TCTH2(I)*(TT-TSH2(I)+ANTH2(I)*IBURN2(I))+TSH2(I)
  332.         TSUM1=TSH1N(I)
  333.         TSUM2=TSH2N(I)
  334.         TSM1=TSUM1-TSSH1
  335.       ELSE
  336.         TAH1N(I)=TCTAH1(I)*(ANTH1(I)*IBURN1Z(I)-TAH1(I))+TAH1(I)
  337.         IF(TAH1N(I).LT.0.)TAH1N(I)=0.
  338.         TAH2N(I)=TCTAH2(I)*(ANTH2(I)*IBURN2(I)-TAH2(I))+TAH2(I)
  339.         IF(TAH2N(I).LT.0.)TAH2N(I)=0.
  340.         IF(IHONEY(I).EQ.1)THEN
  341. C              HONEYWELL ROUND THERMOSTAT
  342.           TSH1N(I)=TCTH1(I)*(TT-TSH1(I))+TSH1(I)
  343.           TSH2N(I)=TCTH2(I)*(TT-TSH2(I))+TSH2(I)
  344.           TSUM1=TSH1N(I)+TAH1N(I)
  345.           TSUM2=TSH2N(I)+TAH2N(I)
  346.           TSM1=TSUM1-TSSH1
  347.           TSM2=TSUM2-TSSH2
  348.         ELSE
  349. C             HONEYWELL T8200 AND T8300 THERMOSTATS
  350.           TSH1N(I)=TCTH1(I)*(TI*AIRCPG(I)-TSH1(I))+TSH1(I)
  351.           IF(TSH1N(I).LT.SMALL)TSH1N(I)=0.
  352.           TSWH1N(I)=TCTWH1(I)*(TWALL*WALCPG(I)-TSWH1(I))+TSWH1(I)
  353.           IF(TSWH1N(I).LT.SMALL)TSWH1N(I)=0.
  354.           TSH2N(I)=TCTH2(I)*(TI*AIRCPG(I)-TSH2(I))+TSH2(I)
  355.           IF(TSH2N(I).LT.SMALL)TSH2N(I)=0.
  356.           TSWH2N(I)=TCTWH2(I)*(TWALL*WALCPG(I)-TSWH2(I))+TSWH2(I)
  357.           IF(TSWH2N(I).LT.SMALL)TSWH2N(I)=0.
  358.           TSENS=TSH1N(I)+TSWH1N(I)
  359.           TDIF=TSENS-TSSH1
  360.           IF(IHONEY(I).EQ.2)THEN
  361.             TDIF1=0.
  362.           ELSE
  363.             IF(ICNT(I).EQ.0)THEN
  364.               TDIF1=TDIF*XINTG(I)+TDIF1
  365.               IF(TDIF1.GT.3.) TDIF1=3.
  366.               IF(TDIF1.LT.0.) TDIF1=0.
  367.             ELSE
  368.               TDIF1=0.
  369.             ENDIF
  370.           ENDIF
  371.           TSM1=TDIF1+TAH1N(I) + TDIF
  372.           TSENS=TSH2N(I)+TSWH2N(I)
  373.           TDIF=TSENS-TSSH2
  374.           IF(IHONEY(I).EQ.2)THEN
  375.             TDIF1=0.
  376.           ELSE
  377.             IF(ICNT(I).EQ.0)THEN
  378.               TDIF1=TDIF*XINTG(I)+TDIF1
  379.               IF(TDIF1.GT.3.) TDIF1=3.
  380.               IF(TDIF1.LT.0.) TDIF1=0.
  381.             ELSE
  382.               TDIF1=0.
  383.             ENDIF
  384.           ENDIF
  385.           TSM2=TDIF1+TAH2N(I) + TDIF
  386.         ENDIF
  387.       ENDIF
  388. C - - - SET STATE OF THERMOSTAT SWITCH
  389.       IF(TSM1.GE.DIFFH1(I))IBUR1N(I)=0
  390.       IF(TSM1.LE.0.)IBUR1N(I)=IACN1(I)
  391.       IF(ICOLL.NE.0.AND.IZNCOLL.EQ.I) CDATA(64)=TSM1
  392. C      
  393.       IF(IHONEY(I).LT.2)TSM2=TSUM2-TSSH2
  394.       IF(TSM2.GE.DIFFH2(I))IBUR2N(I)=0
  395.       IF(TSM2.LE.0.)IBUR2N(I)=IACN1(I)
  396. C      
  397. 710   IACP=(1-IBUR1N(I))*IBUR2N(I)
  398.       IBURNNZ(I)=IACP+IBUR1N(I)
  399.       IF(NDEBUG)THEN
  400.         WRITE(60,*)' TIME,IBURNN,TSSH1= ',TIME,IBURNNZ(I),TSSH1
  401.         WRITE(60,*)'      TAH1N,TSM1,DIFFH1= ',TAH1N(I),TSM1,DIFFH1(I)
  402.         WRITE(60,*)'      TSENS,TI,TWALL= ',TSENS,TI,TWALL
  403.         WRITE(60,*)' '
  404.       ENDIF
  405.       RETURN
  406. C
  407.       ENTRY TSTAT5(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  408. C  ***  RESET VALUES
  409.       TSC1(I) = TSC1N(I)
  410.       TSH1(I)=TSH1N(I)
  411.       TSH2(I)=TSH2N(I)
  412.       TAC1(I)=TAC1N(I)
  413.       TSWC1(I)=TSWC1N(I)
  414.       TAH1(I)=TAH1N(I)
  415.       TAH2(I)=TAH2N(I)
  416.       TSWH1(I)=TSWH1N(I)
  417.       TSWH2(I)=TSWH2N(I)
  418.       IBURN1Z(I)=IBUR1N(I)
  419.       IBURN2(I)=IBUR2N(I)
  420.       IBRNP(I)=IBURNNZ(I)
  421.       IACQ(I)=IACZP(I)
  422.       IACZP(I)=IACNZ(I)
  423.       IF(.NOT.NDEBUG)RETURN
  424.       GO TO 720
  425. C      
  426.       ENTRY TSTAT6(IACZP,IBURN1Z,IBRNP,IEFAN,TI,TWALL,I)
  427. C *** WRITE VALUES OF THERMOSTAT PARAMETERS      
  428.   720 NPR=NPR+1
  429.       IF(NPR.GT.NPRMAX)RETURN
  430.       IF(NPR.LE.1.AND.MODE.NE.2)WRITE(60,3001)
  431.       IF(NPR.LE.1.AND.MODE.EQ.2)WRITE(60,3002)
  432.       IF(MODE.NE.2)THEN
  433.       WRITE(60,3000)TIME,TI,TWALL,TSH1(I),TSH2(I),TAH1(I),TAH2(I),
  434.      +               TSWH1(I),TSWH2(I),IBURN1Z(I),IBURN2(I),IBRNP(I)
  435.       ELSE
  436.         WRITE(60,3003)TIME,TI,TWALL,TSC1(I),TAC1(I),TSWC1(I),
  437.      +                 IACQ(I),IACNZ(I)
  438.       ENDIF
  439.       RETURN
  440.   500 FORMAT(1X,'TAC1N= ',G13.5)
  441.   510 FORMAT(1X,'TSTAT - TDIF,TAC1N,TDIF1 ',3E11.3)
  442.  2160 FORMAT(1H0,17X,'THERMOSTAT SETTING',11X,'SET BACK'/14X,
  443.      +       'T/C    ANTC   SET   DIFF   AMT    ON    OFF'/
  444.      +       ' COOLING   ',F9.5,F5.2,F7.1,F6.3,F7.1,2F7.2)
  445.  2170 FORMAT(' HEAT #1   ',F9.5,F5.2,F7.1,F6.3,F7.1,2F7.2)
  446.  2180 FORMAT (' HEAT #2   ',F9.5,F5.2,F7.1,F6.3,F7.1,2F7.2)
  447.  2190 FORMAT(1X,' IHONEY= ',I5,' WALL COUPLING= ',G13.5/
  448.      +1X,'WALL SENSOR TIME CONST(CLG,HTG1,HTG2)= ',3G13.5/
  449.      +1X,'ANTICIPATOR TIME CONST(CLG,HTG1,HTG2)= ',3G13.5)
  450.  3000 FORMAT(1X,F10.3,8(1X,F6.2),3(3X,I4))
  451.  3001 FORMAT(1H1,3X,4HTIME,6X,2HTI,4X,5HTWALL,2X,4HTSH1,3X,4HTSH2,
  452.      +3X,4HTAH1,3X,4HTAH2,3X,5HTSWH1,2X,5HTSWH2,2X,6HIBURN1,1X,
  453.      +6HIBURN2,1X,5HIBURN)
  454.  3002 FORMAT(1H1,3X,4HTIME,6X,2HTI,4X,5HTWALL,2X,4HTSC1,3X,4HTAC1,
  455.      +3X,5HTSWC1,7X,3HIAC,9X,
  456.      +4HIACN)
  457.  3003 FORMAT(1X,F10.3,5(1X,F6.2),2(3X,I4))
  458.  3020 FORMAT(' ***OVERLAP BETWEEN A/C & H/P THERMOSTATS - RUN STOPPED')
  459.       END
  460.       SUBROUTINE SETBAK(TIME,IMODE,SETH1,SETH2,SETBH,SBKHON,SBKHOF,
  461.      +                  SETBH2,SBKHON2,SBKHOF2,
  462.      +                  TSSH1,TSSH2,SETC1,SETBC,SBKCON,SBKCOF,TSSC1)
  463. C  SETBACK/SETUP CONTROL MODEL     
  464. C           NOTE:: THE CODE BELOW FROM MULTIZONE SP43 TSTAT
  465. C                  ELIMINATES THE NEED FOR LARGE SETBHN,SETBCN ARRAYS      
  466.       IF(IMODE.EQ.1)THEN
  467. C   **********  COOLING  ***********      
  468.         IF(SBKCON.EQ.SBKCOF)THEN
  469. C                                   ***NO SETBACK 
  470.           TSSC1=SETC1 
  471.         ELSE
  472. C                                   ***SETBACK SOMETIME DURING THE DAY
  473.           SETBCN=0. 
  474.           IF(SBKCON.GT.SBKCOF)THEN
  475. C                                   ***SETBACK STARTS BEFORE MIDNIGHT 
  476.             IF(TIME.GT.SBKCON)SETBCN=SETBC
  477.             IF(TIME.LE.SBKCOF)SETBCN=SETBC
  478.           ELSE
  479. C                                   ***SETBACK STARTS AFTER MIDNIGHT
  480.             IF(SBKCON.LE.TIME.AND.TIME.LE.SBKCOF)SETBCN=SETBC 
  481.           ENDIF 
  482.           TSSC1=SETC1-SETBCN
  483.         ENDIF 
  484.       ELSE
  485. C        ********* HEATING **********      
  486.         ICODE=0
  487.         IF(SBKHON.NE.SBKHOF.AND.SBKHON2.EQ.SBKHOF2)THEN
  488.           ICODE=1
  489.         ELSEIF(SBKHON.EQ.SBKHOF.AND.SBKHON2.NE.SBKHOF2)THEN
  490.           ICODE=2
  491.         ELSE
  492. C              BOTH SETBACKS AVAILABLE!        
  493.           ICODE=3
  494.         ENDIF
  495.         IF(ICODE.EQ.0)THEN
  496.           TSSH1=SETH1
  497.           TSSH2=SETH2
  498.           GO TO 30
  499.         ENDIF
  500.         IF(ICODE.EQ.2)GO TO 10
  501. C              1ST SETBACK!!!!!!!!!!!        
  502. C                                   ***SETBACK SOMETIME DURING THE DAY
  503.         SETBHN=0. 
  504.         IF(SBKHON.GT.SBKHOF)THEN
  505. C                                   ***SETBACK STARTS BEFORE MIDNIGHT 
  506.           IF(TIME.GT.SBKHON)SETBHN=SETBH
  507.           IF(TIME.LE.SBKHOF)SETBHN=SETBH
  508.         ELSE
  509. C                                   ***SETBACK STARTS AFTER MIDNIGHT
  510.           IF(SBKHON.LE.TIME.AND.TIME.LE.SBKHOF)SETBHN=SETBH 
  511.         ENDIF 
  512.         IF(ICODE.EQ.1)GO TO 20
  513. C              2ND SETBACK!!!!!!!!!!!          
  514.    10   CONTINUE          
  515. C                                   ***SETBACK SOMETIME DURING THE DAY
  516.         IF(SBKHON2.GT.SBKHOF2)THEN
  517. C                                   ***SETBACK STARTS BEFORE MIDNIGHT 
  518.           IF(TIME.GT.SBKHON2)SETBHN=SETBH2
  519.           IF(TIME.LE.SBKHOF2)SETBHN=SETBH2
  520.         ELSE
  521. C                                   ***SETBACK STARTS AFTER MIDNIGHT
  522.           IF(SBKHON2.LE.TIME.AND.TIME.LE.SBKHOF2)SETBHN=SETBH2 
  523.         ENDIF 
  524.    20   CONTINUE          
  525.         TSSH1=SETH1-SETBHN
  526.         TSSH2=SETH2-SETBHN
  527.       ENDIF
  528.    30 CONTINUE      
  529.       RETURN
  530.       END
  531.