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 / STRCTL.FOR < prev    next >
Text File  |  1992-05-09  |  29KB  |  710 lines
  1.       SUBROUTINE STRCTL(IDBG,NZONG) 
  2. C 
  3. C-----READS STRUCTURE INFORMATION FOR RESIDENTIAL SIMULATION PROGRAM
  4. C  NOTES:  1. IN BASMT ZONE, INPUT VALUES FOR THE CEILING AS IF IT
  5. C             WERE AN INTERIOR PARTITION BETWEEN ZONES.  ENTER BASE-
  6. C             MENT ENCLOSURE DATA AFTER OTHER ZONE DATA 
  7. C-----CONSTANTS 
  8. CMDK NKONST
  9. CMDK NRF
  10. CMDK NWL
  11. CMDK NWN
  12. CMDK NWZN
  13. CMDK NZN
  14. CMDK NZW
  15. C 
  16. C-----COMMON BLOCKS 
  17. CMDK CNSTRK
  18. CMDK BLKFMI
  19. CMDK ENCBK1
  20. CMDK ENCBK2
  21. CMDK ENCBLK
  22. CMDK DUCTS4
  23. CMDK IRDFQ
  24. CMDK IWLS
  25. CMDK IZPART
  26. CMDK IZZQ
  27. CMDK MZON1
  28. CMDK NCR
  29. CMDK QUAYCH
  30. CMDK ROOF
  31. CMDK STRUCA
  32. CMDK SURFAR
  33. CMDK WNDBLK
  34. C 
  35.       CHARACTER*1 IDEXPI
  36.       CHARACTER*1 IDXP (10), IFXDI (4)
  37.       CHARACTER*4 NHEADE(10)
  38. C 
  39. C-----INITIALIZATIONS 
  40.       DATA IDXP   /'N', 'E', 'S', 'W', 'C', 'F', 'I','B','G','H'/ 
  41.       NLEAKS=0
  42. C         SET SHADING ALWAYS FOR CONSTRUCT 11 (SEE BLKDAT)      
  43.       IFXL(1,11)=4
  44. C ZERO NWALLA,NRMA,AND IZSET ARRAYS 
  45.       CALL ZERVI(NZN,NWALLA)
  46.       CALL ZERVI(NZN,NRMA)
  47.       CALL ZERVI(NZW,IZSET) 
  48. C  NOTE:MULTIPLE INTERIOR WALLS THAT ARENT SEPARATING ZONES AND 
  49. C       HAVE SAME CONSTRUCT NO. SHOULD BE ENTERED CONSECUTIVELY 
  50.       WRITE(60,524)
  51.       CALL FLUSH(6)
  52. C 
  53. C  ---READ DATA THAT IS TO BE LOADED INTO STRUCTURE FILE
  54. C   NIZ IS DIMENSION OF IZWSS AND IZWSV SET IN CMDECKS(6*NZN)
  55.       CALL ZERVI(NIZ,IZWSS) 
  56.       NIZWS=0 
  57.       IZW=0 
  58.       NW=0
  59.       IBAS=0
  60.       ICRW=0
  61. C      
  62.       OPEN(10,FILE='TAPE10',STATUS='OLD',IOSTAT=IO10)
  63.       IF(IO10.NE.0)THEN
  64.         WRITE(60,*)' STRCTL: CANT OPEN TAPE10 WITH STRUCTURE DATA'
  65.         STOP ' STRCTL: CANT OPEN TAPE10 WITH STRUCTURE DATA'
  66.         ENDIF
  67. C 
  68. C     IFXDI   ARRAY OF UP TO FOUR COMPASS DIRECTIONS THAT ARE ALWAYS
  69. C             SHADED
  70. C     WROTA   ANGLE THAT STRUCTURE IS ROTATED FROM THE CARDINAL COMPASS 
  71. C             DIRECTIONS (+ CLOCKWISE, - COUNTER-CLOCKWISE) 
  72.       READ(10, 500) IFXDI, WROTA 
  73. C 
  74. C     CONVERT IFXDI TO NUMERIC VALUE IN IFXDS ARRAY 
  75.       NSQ=0 
  76.       DO 81 I=1, 4
  77.          IFXDS (I)=0
  78.          DO 80 J=1, 4 
  79.             IF (IFXDI (I).EQ.IDXP (J)) THEN 
  80.                IFXDS (I)=J
  81.                NSQ=NSQ+1
  82.             END IF
  83. 80          CONTINUE
  84. 81       CONTINUE 
  85. C 
  86.       IF (NSQ.NE.0) WRITE(60, 503) (IFXDS (I), I=1, NSQ) 
  87.       WRITE(60, 504) WROTA 
  88.       CALL FLUSH(6)
  89.       IAREA=0 
  90. C 
  91. C     NROOMS  NUMBER OF ZONES 
  92. C     NSECTA  NUMBER OF ROOF SECTIONS 
  93. C    NZNBAS  ZONE NUMBER OF BASEMENT (IE. NRM OF BASEMENT ZONE) 
  94. C    NATC  NUMBER OF ATTIC ZONES
  95. C 
  96.       READ(10, *) NROOMS, NSECTA, NZNBAS, NATC 
  97.       WRITE(60, 505) NROOMS, NSECTA, NZNBAS, NATC
  98.       CALL FLUSH(6)
  99.       IF (NROOMS.GT.NZN) STOP 'STRCTL: NROOMS TOO LARGE, INCREMENT NZN' 
  100.       IF (NSECTA.GT.NRF) STOP 'STRCTL: NSECT TOO LARGE, INCREMENT NRF'
  101. C-----READ ZONE DATA
  102.       IZT=0 
  103.       DO 130 J=1, NROOMS
  104. C 
  105. C        NRM     ZONE NUMBER (CAN BE ANY ORDER (BEST TO USE 1,2,..5,
  106. C                HOWEVER), BUT KEEP MAXIMUM VALUE AS LOW AS POSSIBLE 
  107. C        NWALLI  NUMBER OF ENCLOSURE ELEMENTS (OUTSIDE WALLS, INTERIOR
  108. C                WALLS, FLOOR, CEILING, ETC.) IN THIS ZONE
  109. C NOTE: NWALLA MAY GET MODIFIED WITH VALUES FROM AN INTERIOR PARTITION
  110. C       BETWEEN ZONES SET IN ANOTHER ZONE .
  111. C        CLHT    CEILING HEIGHT IN THIS ZONE, FT
  112. C        KNSTAT   CONSTRUCT THAT HAS TSTAT MOUNTED ON IN THIS ZONE
  113.          READ(10, *) NRM, NWALLI, CLHT (NRM), KNSTAT(NRM)
  114.          N2I=1
  115. C   NWALLA(NRM) CAN BE > 0 IF MIRROR IMAGE WAS SET IN THIS ZONE 
  116.          IF(NWALLA(NRM).GT.0)N2I=NWALLA(NRM)+1
  117.          NWALLA(NRM)=NWALLA(NRM)+NWALLI 
  118.          WRITE(60,529) 
  119.          WRITE(60, 506) NRM, NWALLI, CLHT (NRM), KNSTAT(NRM) 
  120.          WRITE(60,530) 
  121.          CALL FLUSH(6)
  122.          NRMA (NRM)=NRM 
  123.          N2=NWALLA (NRM)
  124. C 
  125. C NOTE!!!!!!! 
  126. C  INTERIOR PARTITION WALLS BETWEEN ZONES SHOULD ONLY BE ENTERED
  127. C  IN 1ST ZONE ( A MIRROR IMAGE IS SET UP IN 2ND ZONE)
  128. C 
  129.          DO 128 K2=N2I, N2
  130.             IZW=IZW+1 
  131.             IF (IZW.GT.NZW) STOP 'STRCTL: NZW MUST BE INCREMENTED'
  132.             NENC (NRM, K2)=IZW
  133. C 
  134. C    IDEXPI:                   IDEXP
  135. C       N  NORTH OUTSIDE WALL  - 1 
  136. C       E  EAST OUTSIDE WALL   - 2 
  137. C       S  SOUTH OUTSIDE WALL  - 3 
  138. C       W  WEST OUTSIDE WALL   - 4 
  139. C       C  CEILING INSIDE      - 5 
  140. C       F  FLOOR INSIDE        - 6 
  141. C       I  INTERIOR INSIDE WALL- 7 
  142. C       B  BSMT FLOOR          - 8 (SET NZNC=0)
  143. C       G  BELOW GRND BSMT WALL- 9 (SET NZNC=0)
  144. C       H  ABOVE GRND BSMT WALL- 10(SET NZNC=0)
  145. C    NHEADE - UP TO 40 CHAR DESCRIPTION OF THIS WALL
  146.             READ(10, 501) IDEXPI, NHEADE 
  147. C           CONVERT IDEXPI TO NUMERIC VALUE 
  148.             DO 90 I=1, 10 
  149.                IF (IDEXPI.EQ.IDXP (I)) GO TO 92 
  150. 90             CONTINUE 
  151.             WRITE(60, *)' STRCTL: IDEXPI BAD; NRM,K2,IDEXPI= ',
  152.      +            NRM,K2,IDEXPI
  153.             CALL FLUSH(6)
  154.             STOP 'STRCTL: IDEXPI BAD' 
  155. C 
  156. 92          IDX=I 
  157.             WRITE(60, 518) NHEADE
  158.             WRITE(60, 507) IDX, IZW
  159.             CALL FLUSH(6)
  160.             IDEXP (IZW)=IDX 
  161. C 
  162. C           NZNC  CONNECTING SPACE NUMBER:  
  163. C              0 FOR BSMT FLR,BLW GRND BSMT WALL, AND ABV GRND BSMT WALL
  164. C              1 OUTDOORS     4 ATTIC NUMBER 1       7 BASEMENT W/DUCTS 
  165. C              2 BASEMENT     5 ATTIC NUMBER 2       8 CRAWLSPC W/DUCTS 
  166. C              3 CRAWLSPACE   6 ATTACHED GARAGE 
  167. C              9-13  FIVE ZONES IN LIVING SPACE 
  168. C           KONSTA  CONSTRUCT NUMBER (SEE BLOCK DATA SUBPROGRAM)
  169. C                   IF KONSTA < 0, CONSTRUCT IN BLOCK DATA WILL BE 
  170. C                   REPLACED.
  171. C           HGTA    HEIGHT OF ENCLOSURE ELEMENT, FT 
  172. C           XLENA   LENGTH OF ENCLOSURE ELEMENT, FT 
  173. C  NOTE 1. ENTER TWO ELEMENTS FOR BASEMENT FLOOR(ID=8).  THE SMALLER
  174. C          ELEMENT WILL REPRESENT THE FLOOR NEAR THE OUTSIDE WALL AND 
  175. C          SHOULD BE ENTERED SECOND.
  176. C  NOTE 2. ENTER FOUR ELEMENTS FOR THE BASEMENT WALL BELOW GROUND(ID=9) 
  177. C          IN ORDER OF DEPTH.  HEIGHT DZN OF GROUND NODES USED IN SETNOD
  178. C          OF SOIL OR GROUND TEMPERATURE ROUTINES WILL BE THE SAME AS 
  179. C          THE HEIGHT OF THESE ELEMENTS.
  180. C  NOTE 3. FOR FLOOR ABOVE BSMT, ENTER ONLY ONE ELEMENT W/O DUCTS BELOW 
  181. C          AND ONLY ONE ELEMENT W/ DUCTS BELOW.  IF MORE THAN ONE ELEMENT 
  182. C          IS ENTERED FOR THESE, FLOOR LOSSES WONT GET CALC PROPERLY).
  183. C          ENTER THESE ELEMENTS IN ZONE ABOVE BSMT AND DONT REPEAT
  184. C          IN BASEMENT ZONE FOR BASEMENT CEILING.  BSMT CEILING WILL
  185. C          GET GENERATED FROM LVG SPC FLOOR INFO. 
  186. C 
  187.             READ(10, *) NZNC (IZW), KONSTA (IZW), HGTA (IZW),
  188.      1         XLENA (IZW)
  189.             WRITE(60, 508) NZNC (IZW), KONSTA (IZW), HGTA (IZW), 
  190.      1         XLENA (IZW)
  191.             CALL FLUSH(6)
  192. C  OVERRIDE IDEXPI IF <= 4 AND NZONC=6 (GARAGE) SO THAT IWALLS IS USED TO 
  193. C  CALC WALL TEMPS OF WALL BETWEEN LVG SPC AND GARAGE 
  194.          IF(IDX.LE.4.AND.NZNC(IZW).EQ.6)THEN
  195.             IDXGAR=IDX 
  196.             IDX=7 
  197.             IDEXP(IZW)=IDX
  198.             WRITE(60,532)IZW,IDXGAR,IDX
  199.             ENDIF 
  200.          IF(NZNC(IZW).EQ.2.OR.NZNC(IZW).EQ.7)IBAS=1 
  201.          IF(NZNC(IZW).EQ.3.OR.NZNC(IZW).EQ.8)ICRW=1 
  202. C 
  203. C           CHECK TO SEE IF CONSTRUCT IS AVAILABLE; IF NOT, READ
  204. C           CONSTRUCT DATA
  205.          IF (KONSTA (IZW).GT.NCONS.OR.KONSTA(IZW).LT.0) THEN
  206.             IF(KONSTA(IZW).GT.0)THEN
  207.                NCONS=NCONS+1
  208.                KONSTA (IZW)=NCONS 
  209.                N=NCONS
  210.             ELSE
  211.                KONSTA(IZW)=IABS(KONSTA(IZW))
  212.                N=KONSTA(IZW)
  213.                WRITE(60,3060)N 
  214.                ENDIF
  215.                CALL FLUSH(6)
  216. C 
  217. C              NLMP   NUMBER OF LUMPS IN THIS CONSTRUCT (1 OR 3)
  218. C              IFXL   = (1, 2, 3, 4) FOR ALWAYS SHADING (N, E, S, W) DIR
  219. C                     = 0 FOR DEACTIVATE
  220.             WRITE(60,*)' TRYING TO READ NEW CONSTRUCT DATA!!!!'
  221.             CALL FLUSH(6)
  222.             READ(10, *) NLMP (N), (IFXL (I, N), I=1, 4)
  223.             WRITE(60, 509) NLMP (N), (IFXL (I, N), I=1, 4) 
  224.             CALL FLUSH(6)
  225. C 
  226. C              WK   THERMAL CONDUCTIVITY OF LUMP MATERIAL, BTU/HR-FT-F
  227. C              WX   THICKNESS OF LUMP LAYER, FT 
  228. C              WRHO DENSITY OF LUMP MATERIAL, LBM/FT3 
  229. C              WCP  SPECIFIC HEAT OF LUMP MATERIAL, BTU/LBM-F 
  230.             READ(10, *) WK (1, N), WX (1, N), WRHO (1, N), 
  231.      1            WCP (1, N)
  232.             WRITE(60, 510) 1, WK (1, N), WX (1, N),
  233.      1            WRHO (1, N), WCP (1, N) 
  234.             CALL FLUSH(6)
  235.             IF (NLMP (N).EQ.3) THEN 
  236.                READ(10, *) WK (2, N), WX (2, N), WRHO (2, N),
  237.      1               WCP (2, N) 
  238.                WRITE(60, 510) 2, WK (2, N), WX (2, N), 
  239.      1               WRHO (2, N), WCP (2, N)
  240.                CALL FLUSH(6)
  241.                READ(10, *) WK (3, N), WX (3, N), WRHO (3, N),
  242.      1               WCP (3, N) 
  243.                WRITE(60, 510) 3, WK (3, N), WX (3, N), 
  244.      1               WRHO (3, N), WCP (3, N)
  245.                CALL FLUSH(6)
  246.                END IF 
  247.             END IF
  248. C 
  249. C           LEAKAGE PARAMETERS
  250. C           DON'T READ LEAKAGE IF ENCLOSURE ELEMENT IS INTERIOR WALL
  251. C           WITHIN THE SAME ZONE (IR=0) OR CEILING OR FLOOR
  252.             IR=0
  253. C 
  254. C           IF THE ENCLOSURE ELEMENT IS AN INTERIOR WALL BETWEEN ZONES
  255. C           , AN ABOVE-GRND BSMT WALL, OR AN OUTSIDE WALL (IR=1)
  256. C   INTERIOR WALL BETWEEN ZONES 
  257. C                   NZON IS CONNECTING ZONE
  258.               NZON=NZNC(IZW)-8
  259.               IF(IDX.EQ.7.AND.NRM.NE.NZON)IR=1
  260. C   BASEMENT WALL ABOVE GROUND
  261.               IF(IDX.EQ.10)IR=1 
  262. C   OUTSIDE WALL
  263.               IF(IDX.LE.4)IR=1
  264. C              
  265. C              WDLEAK - WIDTH OF LEAKAGE OPENING, IN.
  266. C              HGTLEAK - HEIGHT OF LEAKAGE OPENING, IN.
  267. C              HNPL - HEIGHT FROM CENTER OF OPENING UP TO NEUTRAL
  268. C                     PRESSURE LEVEL IN BUILDING (USED ON OUTSIDE
  269. C                     WALLS WHEN WIND VELOCITY IS ZERO.--SEE P.23.4
  270. C                     OF 1989 ASHRAE HANDBOOK OF FUNDAMENTALS)
  271. C                    
  272.              IF(IR.EQ.1)THEN
  273.                READ(10, *) WDLEAK (IZW), HGTLEAK (IZW), HNPL (IZW)
  274.                WRITE(60, 512) WDLEAK (IZW), HGTLEAK (IZW), HNPL (IZW)
  275.                CALL FLUSH(6)
  276.              END IF 
  277. C 
  278. C           WALLFI   LONGWAVE EMISSIVITY INSIDE SURFACE OF
  279. C                    ENCLOSURE ELEMENT
  280. C           WALLFO   LONGWAVE EMISSIVITY AT OUTER SURFACE OF
  281. C                    ENCLOSURE ELEMENT
  282. C           EPSLR    SHORTWAVE (SOLAR) ABSORPTIVITY AT OUTER
  283. C                    SURFACE OF ENCLOSURE ELEMENT 
  284.             IF(IDX.LE.4)THEN
  285.               READ(10, *) WALLFI (IZW), WALLFO (IZW), EPSLR(IZW) 
  286.               WRITE(60, 531) WALLFI (IZW), WALLFO (IZW),EPSLR(IZW) 
  287.               CALL FLUSH(6)
  288.             ELSE
  289.               READ(10, *) WALLFI(IZW), WALLFO(IZW) 
  290.               WRITE(60, 513) WALLFI (IZW), WALLFO (IZW)
  291.               CALL FLUSH(6)
  292.               ENDIF 
  293.             WRITE(60, 520) IDX 
  294.             CALL FLUSH(6)
  295. C 
  296. C           NWNDA   NUMBER OF WINDOWS FOR THIS ENCLOSURE ELEMENT
  297. C           OHANGA  ROOF OVERHANG ABOVE THESE WINDOWS, FT 
  298. C           HHANGA  HEIGHT TO OVERHANG FROM TOP OF OUTSIDE WALL, IN.
  299.             IF (IDX.LE.4) THEN
  300.                READ(10, *) NWNDA (IZW), OHANGA (IZW), HHANGA (IZW) 
  301.                WRITE(60, 514) NWNDA (IZW), OHANGA (IZW), 
  302.      1            HHANGA (IZW)
  303.                CALL FLUSH(6)
  304.                N4=NWNDA (IZW) 
  305.                IF (N4.NE.0) THEN
  306. C                 READ WINDOW DATA
  307.                   DO 94 K4=1, N4
  308.                      NW=NW+1
  309.                      NIWND (IZW, K4)=NW 
  310. C 
  311. C                    WHGTA   HEIGHT OF TRANSPARENT WINDOW, FT 
  312. C                    WWDTHA  WIDTH OF TRANSPARENT WINDOW REDUCED TO 
  313. C                            ACCOUNT FOR ANY BLOCKAGE, FT 
  314. C                    REVELA  DEPTH OF REVEAL OR SETBACK OF WINDOW FROM
  315. C                            OUTER SURFACE, IN. 
  316. C                    OHNGHA  DISTANCE FROM TOP OF WINDOW AND ROOF 
  317. C                            OVERHANG, IN.
  318. C                    SFSDA   WINDOW SHADING COEFFICIENT (ASHRAE 
  319. C                            HANDBOOK OF FUNDAMENTALS)
  320.                      READ(10, *) WHGTA (NW), WWDTHA (NW),
  321.      1                  REVELA (NW), OHNGHA (NW), SFSDA (NW)
  322.                      WRITE(60, 515) NW, WHGTA (NW), WWDTHA (NW), 
  323.      1                  REVELA (NW), OHNGHA (NW), SFSDA (NW)
  324.                      CALL FLUSH(6)
  325. 94                   CONTINUE 
  326.                END IF 
  327.             END IF
  328.       KONI=KONSTA(IZW)
  329.       NZONC=NZNC(IZW) 
  330.       NZON=NZONC-8
  331.       IBZ=0 
  332. C   IBZ = 1, IF THIS ENCLOSURE ELE IS AN INTERIOR WALL BETWN ZONES      
  333.       IF(IDX.EQ.7.AND.NRM.NE.NZON.AND.NZONC.NE.6)IBZ=1
  334.       IF(IDX.EQ.6)THEN
  335.           IF(NZONC.EQ.2.OR.NZONC.EQ.7)THEN
  336. C   IBZ = 2, IF THIS ENCLOSURE ELE IS FLOOR BETWN LVG ZONE AND BSMT          
  337.               IBZ=2 
  338.               NZON=NZNBAS 
  339.               ENDIF 
  340.           ENDIF 
  341.       IF(IBZ.NE.0)THEN
  342. C - - INTERIOR PARTITION WALL OR FLOOR BETWEEN ZONES
  343. C  SET VALUES FOR OTHER SIDE OF WALL OR FLOOR 
  344.           N2=NWALLA(NZON)+1 
  345.           NWALLA(NZON)=N2 
  346.           IF(NRMA(NZON).EQ.0)NRMA(NZON)=NZON
  347.           IZWX=IZW
  348.           IZW=IZW+1 
  349.           IF(IZW.GT.NZW)STOP 'STRCTL:NZW MUST BE INCREASED' 
  350.           NENC(NZON,N2)=IZW 
  351.           IF(IBZ.EQ.1)THEN
  352.               WRITE(60,3040)IZWX,NRM,IZW,NZON
  353.               CALL FLUSH(6)
  354. C  IDQ IS THE IDEXP OF THE CONNECTING ZONE IN WHICH A MIRROR IMAGE
  355. C  IS SET              
  356.               IDQ=7 
  357.               ENDIF 
  358.           IF(IBZ.EQ.2)THEN
  359.               WRITE(60,3050)IZWX,NRM,IZW,NZON
  360.               CALL FLUSH(6)
  361.               IDQ=5 
  362.               ENDIF 
  363.           IDEXP(IZW)=IDQ
  364.           NZNC(IZW)=NRM+8 
  365. C - - LOOK FOR MIRROR IMAGE OF THIS CONSTRUCT; IF NOT AVAILABLE, THEN 
  366. C         SET UP A NEW CONSTRUCT WHICH IS A MIRROR IMAGE
  367.           CALL KONSET(KONI,KONO)
  368.           KONSTA(IZW)=KONO
  369.           HGTA(IZW)=HGTA(IZWX)
  370.           XLENA(IZW)=XLENA(IZWX)
  371. C          WDLEAK(IZW)=WDLEAK(IZWX)
  372. C          HGTLEAK(IZW)=HGTLEAK(IZWX)
  373. C          HNPL(IZW)=HNPL(IZWX)
  374.           WALLFI(IZW)=WALLFO(IZWX)
  375.           WALLFO(IZW)=WALLFI(IZWX)
  376.           AWALLA(IZW)=AWALLA(IZWX)
  377. C  SET IZW ON BACKSIDE OF PARTITION WALL BETWEEN ZONES
  378.           NIZWS=NIZWS+1 
  379.           IF(NIZWS.GT.NIZ)THEN
  380.             WRITE(60,528)2*NIZ 
  381.            STOP 'INCR DIM OF IZWSS AND IZWSV IN CMDECKS & NIZ IN STRCTL' 
  382.             ENDIF 
  383. C  IZSET IS IZW ON THIS SIDE OF INTERIOR PART WALL BETWN ZONES
  384. C  IZWSS IS IZW ON BACK SIDE OF INTERIOR PART WALL BETWEEN ZONES            
  385.           IZSET(NIZWS)=IZWX 
  386.           IZWSS(NIZWS)=IZW
  387.           NIZWS = NIZWS + 1
  388.           IF(NIZWS.GT.NIZ)THEN
  389.             WRITE(60,528)2*NIZ 
  390.            STOP 'INCR DIM OF IZWSS AND IZWSV IN CMDECKS & NIZ IN STRCTL' 
  391.           ENDIF 
  392.           IZSET(NIZWS)=IZW
  393.           IZWSS(NIZWS)=IZWX 
  394.           ENDIF 
  395. 128         CONTINUE
  396.          WRITE(60, 521)
  397. C 
  398. C        CHECK FOR CONSOLIDATION OF ENCLOSURE ELEMENTS IN ZONE NRM
  399. C        IZT IS IZW AT END OF ENCLOSURE ELEMENTS FOR PREVIOUS ZONE
  400.          CALL CONSLD (NRM, IAREA, 0,IDBG) 
  401. C         
  402. C   SET LEAKAGE CONNECTION DATA
  403.          N2=NWALLA(NRM)
  404.          DO 129 K2=1,N2
  405.          IZW=NENC(NRM,K2)
  406.          NZONC=NZNC(IZW)
  407. C              SET THE ZONE THAT CONNECTS TO AN ATTACHED GARAG         
  408.          IF(NZONC.EQ.6)THEN
  409.            NZONG=NRM
  410.          ENDIF
  411.          IF(WDLEAK(IZW).GT.0.)THEN
  412.            NLEAKS=NLEAKS+1
  413.            IZWLEAK(NLEAKS)=IZW
  414.          ENDIF
  415.   129    CONTINUE         
  416. C        IZT IS IZW AT END OF ENCLOSURE ELEMENTS FOR THIS ZONE
  417.           IZW=IZT 
  418.   130    CONTINUE 
  419. C   NOW RESET IZSET ARRAY WITH ORDER OF IZWSS VALUES
  420.       DO 131 IQ=1,NIZWS 
  421.   131 IZWSV(IQ)=IZSET(IQ) 
  422.       CALL ZERVI(NZW,IZSET) 
  423.       DO 132 IQ=1,NIZWS 
  424.       IZW=IZWSS(IQ) 
  425.   132 IZSET(IZW)=IZWSV(IQ)
  426.       WRITE(60,*)'  ***  NOW PRINT IZW INDICES OF MIRROR-IMAGE ELEMENTS'
  427.       DO 133 IZW=1,NZW
  428.       IF(IZSET(IZW).NE.0)THEN
  429.         WRITE(60,*)' IZW,IZSET(IZW)= ',IZW,IZSET(IZW)
  430.       ENDIF
  431.   133 CONTINUE
  432. C 
  433. C     READ ROOF DATA
  434.       DO 140 K5=1, NSECTA 
  435. C        READ COMPASS DIRECTION THAT ROOF FACES 
  436. C 
  437. C        IDEXPI:  
  438. C           N  NORTH FACING SECTION- 1
  439. C           E  EAST FACING SECTION - 2
  440. C           S  SOUTH FACING SECTION- 3
  441. C           W  WEST FACING SECTION - 4
  442. C  NHEADE - UP TO 40 CHAR DESCRIPTION OF ROOF ELEMENT
  443.          READ(10,501) IDEXPI,NHEADE
  444.          WRITE(60, 519) IDEXPI,NHEADE
  445.          DO 138 I=1, 4
  446.             IF (IDEXPI.EQ.IDXP (I)) GO TO 139 
  447. 138         CONTINUE
  448.          STOP 'STRCTL: IDEXP OF ROOF SECTION IS BAD'
  449.   139      IDEXPR (K5)=I
  450.          WRITE(60, 516) K5, IDEXPR (K5)
  451. C 
  452. C        NATTC   ATTIC NUMBER BELOW ROOF (1 OR 2) 
  453. C        KONSTR  CONSTRUCT NUMBER (SEE BLOCK DATA SUBPROGRAM) 
  454. C        SECTAA  SURFACE AREA OF THIS ROOF SECTION, SQ FT 
  455. C        PAA     PITCH ANGLE BETWEEN ROOF LINE AND HORIZONTAL, DEG
  456. C        ROOFI   LONGWAVE EMISSIVITY/ABSORPTIVITY FOR UNDERSIDE OF ROOF 
  457. C        ROOFO   LONGWAVE EMISSIVITY/ABSORPTIVITY FOR OUTER SURF OF ROOF
  458. C       EPRSLR   ABSORPTIVITY FOR SOLAR RADIATION AT OUTER SURFACE OF 
  459. C                ROOF 
  460.          READ(10, *) NATTC (K5), KONSTR (K5), SECTAA (K5), PAAI, 
  461.      1      ROOFI (K5), ROOFO (K5), EPRSLR(K5)
  462.          WRITE(60, 517) NATTC (K5), KONSTR (K5), SECTAA (K5), PAAI,
  463.      1      ROOFI (K5), ROOFO (K5),EPRSLR(K5) 
  464.          IF(NATTC(K5).EQ.1)THEN 
  465.             KONRF(1)=KONSTR(K5) 
  466.             ROFOM(1)=ROOFO(K5)
  467.             ROFIM(1)=ROOFI(K5)
  468.             ENDIF 
  469.          IF(NATTC(K5).EQ.2)THEN 
  470.             KONRF(2)=KONSTR(K5) 
  471.             ROFOM(2)=ROOFO(K5)
  472.             ROFIM(2)=ROOFI(K5)
  473.             ENDIF 
  474.          PAA (K5)=90.0-PAAI 
  475. 140      CONTINUE 
  476. C 
  477. C     CALCULATE WINDOW AND ROOF AREAS 
  478.       IAREA=1 
  479.       CALL CONSLD (1, IAREA, 0,IDBG)
  480. C  NRMVNT IS NUMBER OF LVG SPACE ZONES(OTHER THAN BASEMENT) 
  481. C    THAT VENT GOES THROUGH.
  482. C 
  483. C  HIR - CONVECTIVE HEAT TR COEF AT INSIDE SURF OF ROOF,BTU/HR-FT2-F
  484. C  HAC - CONVECTIVE HEAT TR COEF AT ATTIC SIDE OF CEILING, BTU/HR-FT2-F
  485.       READ(10, *)HIR,HAC 
  486. C  IF IHCOLD = 0, H'S WILL BE CALC IN SUB: HC FROM STILL AIR VALUES 
  487. C                 AS A FUNCT OF HT FLOW DIR. AND SURF. ORIENTATION
  488. C                 IF HIBLWR>0., IT WILL BE USED AT WALLS WHEN BLOWER
  489. C                 IS ON(TARP USED A VALUE OF 1.09--SEEMS HIGH RDF)
  490. C                 ,BUT INDEP. OF BLWR STATE AND SURF-TO-AIR DT
  491. C  IF IHCOLD = 1, H'S WILL BE CALC IN SUB: HC AS FUNCT OF BLOWER
  492. C                 STATE AND SURF-TO-AIR DELTA T.  HIC,HIG,ETC ARE 
  493. C                 MULTIPLIERS ON BUILT IN FUNCTIONS 
  494. C  IF IHCOLD = 2 , THEN ALL H'S (HIG,HIW,HIC,ETC.) WILL BE
  495. C                  USED AS IS.
  496. C  NOTE:::: UNITS FOR THE H VALUES BELOW ARE BTU/HR-FT2-F
  497. C          NOTE: SET HIC,H1F,HIW,H3B,HIG,H3F,H3FB = 1. OR OTHER
  498. C                VALUE OF MULTIPLIER IF IHCOLD < 2 .  IF IHCOLD
  499. C                =2, THE VALUES HIC,H1F,...WILL BE USED AS IS.
  500. C  HIG   - 
  501. C  HIW   - CONVECTIVE HEAT TR COEF AT INSIDE SURF OF EXTERIOR WALLS
  502. C  HIC   - CONVECTIVE HEAT TR COEF AT INSIDE SURF OF CEILING
  503. C  H1F   - CONVECTIVE HEAT TR COEF AT BASEMENT CEILING
  504. C  H3F   - CONVECTIVE HEAT TR COEF AT LIVING-SPACE SIDE OF LVG SPC FLOOR
  505. C  H3B   - CONVECTIVE HEAT TR COEF AT INSIDE SURF OF BASEMENT WALLS
  506. C  H3FB  - CONVECTIVE HEAT TR COEF AT INSIDE SURF OF BASEMENT FLOOR
  507. C  VOLCRW - VOLUME OF CRAWL SPACE, FT3
  508. C  UWGR  - OVERALL HEAT TR COEF FOR AN UNDERGROUND WALL EXPOSED TO 
  509. C          LIVING SPACE(DO NOT INCL FILM COEF ON SOIL-SIDE OF WALL),
  510. C          BTU/HR-FT2-F
  511. C  USLFL - OVERALL HEAT TR COEF FOR SLAB FLOOR EXPOSED TO LIVING SPACE,
  512. C          BTU/HR-FT2-F
  513. C  UWIND - OVERALL HEAT TR COEF FOR WINDOW GLASS SURFACE, BTU/HR-FT2-F
  514. C  HAD   - CONVECTIVE HEAT TR COEF INSIDE DUCTS,BTU/HR-FT2-F
  515. C  HOD   - CONVECTIVE HEAT TR COEF OUTSIDE DUCTS, BTU/HR-FT2-F
  516. C  ASUP  - TOTAL SURFACE AREA OF SUPPLY DUCT, FT2
  517. C  WSUP  - MASS OF SUPPLY DUCT, LBM
  518. C  ARET  - TOTAL SURFACE AREA OF RETURN DUCT, FT2
  519. C  WRET  - MASS OF RETURN DUCT, LBM
  520. C  FSDIB - FRACTION OF SUPPLY DUCT IN BASEMENT OR OUTSIDE
  521. C  FRDIB - FRACTION OF RETURN DUCT IN BASEMENT OR OUTSIDE
  522. C  NRMVNT- NO. OF LIVING SPACE ZONES ABOVE THE BASEMENT THAT VENT
  523. C          GOES THROUGH AND RECEIVE VENT HEAT LOSS
  524. C  NATVNT- ATTIC NO. THAT VENT GOES THROUGH AND RECEIVES VENT HEAT LOSS
  525.       READ(10, *) IHCOLD, HIG, HIW,
  526.      +    HIC, H1F, H3F,        H3B,  H3FB, 
  527.      +   VOLCRW,               UWGR, USLFL, UWIND,
  528.      +   HAD, HOD, ASUP, WSUP, ARET, WRET, FSDIB, FRDIB,
  529.      +   NRMVNT, NATVNT 
  530. C  ATCVOL - NATC VALUES OF ATTIC VOLUMES, FT3
  531.       READ(10,*)(ATCVOL(IA),IA=1,NATC)
  532. C  ROMVOL - NROOM VALUES OF VOLUME OF AIR IN EACH ZONE, FT3     
  533.       READ(10,*)(ROMVOL(NRM),NRM=1,NROOMS)
  534.       IF(NRMVNT.LE.0.OR.NRMVNT.GT.2)THEN
  535.          STOP 'STRCTL:NRMVNT MUST BE 1 OR 2'
  536.          ENDIF
  537. C  IRMVNT(I) IS LIVING SPACE ZONE NO. THAT VENT GOES THROUGH
  538. C  FRCVNT(I) IS FRACTION OF LVG SPC VENT HEAT LOSS TO ZONE IRMVNT(I)
  539. C 
  540.       READ(10,*)(IRMVNT(IJ),FRCVNT(IJ),IJ=1,NRMVNT) 
  541.       DO 142 J=1,NROOMS 
  542.       NRM=NRMA(J) 
  543. C  NRMBGW - NO. OF BELOW-GROUND WALLS IN THIS ZONE
  544.       READ(10,*) NRMBGW(NRM)
  545. C  AWGR - AREA OF BELOW-GROUND WALLS IN THIS ZONE      
  546.       IF(NRMBGW(NRM).GT.0)READ(10,*)AWGR(NRM) 
  547.       IF(NRMBGW(NRM).EQ.0)THEN
  548.           WRITE(60,522)NRM
  549.       ELSE
  550.           WRITE(60,3020)NRMBGW(I),AWGR(I)
  551.           ENDIF 
  552.   142 CONTINUE
  553.       DO 144 J=1,NROOMS 
  554.       NRM=NRMA(J) 
  555. C  NRMSLB - NO. OF SLAB FLOORS IN THIS ZONE      
  556.       READ(10,*) NRMSLB(NRM)
  557. C  ASLFL - AREA OF SLAB FLOOR IN THIS ZONE, FT2      
  558.       IF(NRMSLB(NRM).GT.0)READ(10,*)ASLFL(NRM)
  559.       IF(NRMSLB(NRM).EQ.0)THEN
  560.           WRITE(60,523)NRM
  561.       ELSE
  562.           WRITE(60,3030)NRMSLB(I),ASLFL(I) 
  563.       ENDIF 
  564.   144 CONTINUE
  565.       WRITE(60,525)
  566.       WRITE(60, 3000) IHCOLD,HIG,
  567.      +   HIW, HIR, HAC, HIC, H1F, H3F,        H3B,
  568.      +           VOLCRW,               UWGR, USLFL, UWIND,
  569.      +         HAD, HOD, ASUP, WSUP, ARET, WRET, FSDIB, FRDIB, NRMVNT,
  570.      +   NATVNT 
  571.       WRITE(60,3010)(IA,ATCVOL(IA),IA=1,NATC) 
  572.       WRITE(60,3012)(NRM,ROMVOL(NRM),NRM=1,NROOMS)
  573. C  GET FURNITURE DATA 
  574. C  RHOFM - DENSITY OF FURNITURE MASS, LBM/FT3
  575. C  CPFM  - SPECIFIC HEAT OF FURNITURE MASS, BTU/LBM-F
  576.       READ(10,*)RHOFM,CPFM
  577.       WRITE(60,533)RHOFM,CPFM
  578.       DO 146 JT=1,NROOMS
  579.       NRM=NRMA(JT)
  580. C  THKFM - THICKNESS OF FURNITURE ELEMENTS ( USED TO CALC SURFACE AREA
  581. C          OF FURNITURE), INCHES
  582. C  XMFM  - MASS OF FURNITURE, LBM      
  583.       READ(10,*)THKFM(NRM),XMFM(NRM)
  584.       WRITE(60,534)NRM,THKFM(NRM),XMFM(NRM)
  585.   146 CONTINUE
  586.       CLOSE(10)
  587. C 
  588. C-----NORMAL EXIT 
  589.       WRITE(60, 2010)
  590.       IF(IDBG.NE.0)WRITE(60,526)(I,IZSET(I),I=1,NZW) 
  591.       WRITE(60,527)NROOMS,NZN,NIZWS,2*NIZ
  592.       CALL FLUSH(6)
  593.       RETURN
  594. C 
  595. 500   FORMAT (4 (A1, 1X), F10.0)
  596. 501   FORMAT (A1, 10A4) 
  597. 503   FORMAT (//1X,'INPUT DATA FOR ENCLOSURE ELEMENTS - - - - ', /, 
  598.      +   ' COMPASS DIRECTIONS THAT ARE ALWAYS SHADED ARE = ', 4I5)
  599. 504   FORMAT (' ANGLE HOUSE IS ROTATED FROM CARDINAL COMPASS POINTS= ', 
  600.      +   G13.5) 
  601. 505   FORMAT (' NUMBER OF ZONES= ', I5, /,
  602.      +   ' NUMBER OF ROOF SECTIONS= ', I5/
  603.      +   ' BASEMENT IS ZONE  ',I5/
  604.      +   ' NUMBER OF ATTIC ZONES= ',I5) 
  605. 506   FORMAT (/1X,'ZONE NUMBER = ', I5,
  606.      +   '  NUMBER OF ENCLOSURE ELEMENTS = ', I5, 
  607.      +   ' CEILING HEIGHT = ', G13.5, ' FT'/
  608.      + 10X, 'CONSTRUCT THAT TSTAT IS MOUNTED ON = ',I5) 
  609. 507   FORMAT (' ENCLOSURE ID= ', I5, ' ENCLOSURE NUMBER= ', I5) 
  610. 508   FORMAT (' THIS ENCLOSURE ELEMENT IS CONNECTED TO SPACE', I5, /, 
  611.      +   ' CONSTRUCT NUMBER  = ', I5, /,
  612.      +   ' HEIGHT AND LENGTH = ', 2G13.5, ' FT')
  613. 509   FORMAT ('  NEW CONSTRUCT NUMBER !!!', /,
  614.      +   10X, 'NUMBER OF LUMPS= ', I5,
  615.      +   '  DIRECTIONS ALWAYS SHADED ARE = ', 4I3)
  616. 510   FORMAT (10X, 'FOR LUMP NUMBER', I3, 
  617.      +   '  K, DX, RHO, CP= ', 4G13.5)
  618. 512   FORMAT (10X, 'WIDTH OF LEAKAGE OPENING = ', G13.5, ' IN', /, 
  619.      +   10X, 'HEIGHT OF LEAKAGE OPENING = ', G13.5,' IN',/,
  620.      +   10X, 'HEIGHT FROM CENTER OF OPENING UP TO NEUTRAL PRESSURE',
  621.      +        ' LEVEL = ',G13.5,' IN')
  622. 513   FORMAT (' EMISSIVITY/ABSORPTIVITY OF INNER AND OUTER SURFACE',
  623.      1   ' OF ENCLOSURE= ', 2G13.5) 
  624. 514   FORMAT (' NUMBER OF WINDOWS= ', I5, 
  625.      1   ' ROOF OVERHANG= ', G13.5, 'FT', 
  626.      2   3X, 'HGT TO OVERHANG FROM TOP OF ENCLOSURE ELEMENT = ', G13.5, 
  627.      3   'IN.') 
  628. 515   FORMAT (' WINDOW INDEX NUMBER = ', I5,
  629.      1   ' HEIGHT AND WIDTH OF WINDOW= ', 2G13.5, ' FT', /, 
  630.      2   10X, 
  631.      3   'REVEAL AND HEIGHT FROM TOP OF WINDOW TO ROOF OVERHANG = ',
  632.      4   2G13.5, ' IN.', /, 
  633.      5   10X, 'WINDOW SHADING COEFFICIENT = ', G13.5) 
  634. 516   FORMAT (' ROOF SECTION INDEX NUMBER = ', I5, /, 
  635.      1   10X, 'ORIENTATION ID = ', I5)
  636. 517   FORMAT (' ATTIC NUMBER THIS SECTION IS CONNECTED TO= ', I5, /,
  637.      1   10X, 'CONSTRUCT NUMBER=', I5, ' AREA= ', G13.5, ' SQ FT', /, 
  638.      2   ' PITCH ANGLE FROM HORIZONTAL= ', G13.5, /,
  639.      3   10X, 'LONGWAVE EMISSIVITY/ABSORPTIVITY OF INNER AND OUTER'/
  640.      4   10X,'SURFACE= ',2G13.5/10X,'SHORTWAVE (SOLAR) ABSORPTIVITY', 
  641.      5   ' AT OUTER SURFACE= ',G13.5) 
  642. 518   FORMAT (/1X,' ****', 10A4)
  643. 519   FORMAT (/1X,' ID FOR ROOF SECTION= ', A4,5X,10A4) 
  644. 520   FORMAT (' IDX= ---- ', I4)
  645. 521   FORMAT ('  - - - - -  END OF ENCLOSURE INPUT DATA FOR THIS ZONE', 
  646.      1   /, 
  647.      2   '     NOW CONSOLIDATE ENCLOSURE ELEMENTS FOR THIS ZONE, ', 
  648.      3   ' IF POSSIBLE', /) 
  649.   522 FORMAT(/' NO BELOW GROUND WALLS IN LVG SPC ZONE ',I5)
  650.   523 FORMAT(/' NO SLAB FLOOR IN LVG SPC ZONE ',I5)
  651.   524 FORMAT(/1X,'*--*--*NOTE:MULTIPLE INTERIOR WALLS THAT ARENT'/
  652.      +12X,'SEPARATING ZONES AND HAVE THE SAME CONSTRUCT NO. SHOULD'/
  653.      +12X,'BE ENTERED CONSECUTIVELY SO THEY WILL GET CONSOLIDATED', 
  654.      +' PROPERLY')
  655.   525 FORMAT('  NOTE:::HIG,HIW,HIR,HAC,HIC,H1F,H3F, AND H3B'/ 
  656.      +' ARE MULTIPLIERS ON BUILT IN NAT. CONV.CORRELATION'/ 
  657.      +'  AND SO SHOULD BE INPUT = OR CLOSE TO 1.0. IF IRF'/ 
  658.      +'  =0, HIR AND HAC ARE USED AS IS (INCL. RAD.COMPONENT)IN'/ 
  659.      +'   THE ATTIC FOR H') 
  660.   526 FORMAT((1X,'STRCTL:I,IZSET(I)= ',2I4))
  661.   527 FORMAT(/1X,5(1H*),' NO. OF ZONES= ',I3, 
  662.      +' VS MAX ZONES ALLOWED AS COMPILED= ',I3/ 
  663.      +1X,5(1H*),'2 X NO. OF MIRROR-IMAGE ENCLOSURES= ',I3,
  664.      +' VS MAX ALLOWED AS COMPILED(2 X NIZ)= ',I3)
  665.   528 FORMAT(/1X,'STRCTL: INCR DIM OF IZWSS AND IZWSV', 
  666.      +' IN CMDECKS AND XX OF NIZ=XX*NZN IN STRCTL'/ 
  667.      +15X,'2 X NIZ IN STRCTL = ',I4)
  668.   529 FORMAT(/1X, 40(2H- ))
  669.   530 FORMAT(1X,40(2H- )) 
  670.   531 FORMAT (' LONGWAVE EMISS/ABSORPT OF INNER AND OUTER SURFACE', 
  671.      1   ' OF ENCLOSURE= ', 2G13.5/ 
  672.      2' SHORTWAVE (SOLAR) ABSORPTIVITY OF OUTER SURFACE= ',G13.5) 
  673.   532 FORMAT(1X,'ENCLOSURE ',I5,': CHANGED ID FROM ',I4,' TO ', 
  674.      +I4,' SO THAT IWALLS '/
  675.      +'  WILL BE USED TO CALC TEMPS OF LVG SPC/GARAGE WALL')
  676.   533 FORMAT(1X,1H0,'FURNITURE DATA'/ 
  677.      +' DENSITY, LB/FT3= ',G13.5/ 
  678.      +' SPEC HEAT, BTU/LB-F= ',G13.5) 
  679.   534 FORMAT((1X,'ZONE NO.= ',I4,' THICKNESS, INCHES = ',G13.5, 
  680.      +' MASS IN ZONE, LB = ',G13.5/)) 
  681.   
  682. 2010  FORMAT (' ********** END OF STRCTL INPUT **********') 
  683. 3000  FORMAT ('0', T2, 'IHCOLD= ',I3, T2, 'HIG   =', G13.5, 
  684.      3   T2, 'HIW   =', G13.5, T24, 'HIR   =', G13.5, 
  685.      4   T46, 'HAC   =', G13.5, T68, 'HIC   =', G13.5,
  686.      5   T90, 'H1F   =', G13.5, T112, 'H3F   =', G13.5, /,
  687.      6                         T24, 'H3B   =', G13.5, 
  688.      8                          T112, 'VOLCRW=', G13.5, /,
  689.      A   T46, 'UWGR  =', G13.5, T68, 'USLFL =', G13.5,
  690.      B   T90, 'UWIND =', G13.5, 
  691.      D   T46, 'HAD   =', G13.5, T68, 'HOD   =', G13.5,
  692.      E   T90, 'ASUP  =', G13.5, T112, 'WSUP  =', G13.5, /,
  693.      F   T2, 'ARET  =', G13.5, T24, 'WRET  =', G13.5, 
  694.      G   T46, 'FSDIB =', G13.5, T68, 'FRDIB =', G13.5/
  695.      +   T2, 'NRMVNT=',I7,6X, T24, 'NATVNT=',I7,6X) 
  696.  3010 FORMAT(1X,'IA,ATCVOL= ',2(I4,G13.5))
  697.  3012 FORMAT(1X,'NRM,ROMVOL= ',5(I4,G13.5)) 
  698.  3020 FORMAT((1X,'ZONE NO. W/ BELOW-GRND WALL= ',I4,' WALL AREA= ', 
  699.      +       G13.5/)) 
  700.  3030 FORMAT((1X,'ZONE NO. W/ SLAB FLOOR= ',I4,' SLAB FLOOR ',
  701.      +       'AREA= ',G13.5/))
  702.  3040 FORMAT(/' MIRROR-IMAGE OF INTERIOR ENCLOSURE ELE NO. ',I4,
  703.      +' IN ZONE ',I4/'  SET AS ELE NO. ',I4, ' IN ZONE ',I4)
  704.  3050 FORMAT(/' MIRROR-IMAGE OF FLOOR ELE NO. ',I4,' IN ZONE ', 
  705.      +I4/'  SET AS INTERIOR ENCLOSURE ELE NO. ',I4,' FOR CEILING ', 
  706.      +'ELE IN ZONE ',I4)
  707.  3060 FORMAT(/'YOU HAVE REPLACED CONSTRUCT NO. ',I4, 
  708.      +' IN BLKDAT') 
  709.       END 
  710.