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 / SOLRLD.FOR < prev    next >
Text File  |  1992-05-09  |  14KB  |  445 lines
  1.       SUBROUTINE SOLRLD(STM,NDYRI,ICC,KM,IDAY,DST,NSEQW,IDBG)
  2. C     PROGRAM TO CALCULATE MONTHLY SOLAR LOADING ON RES. STRUCTURES 
  3. C  STM      TIME ZONE (SEE ETIME SUBPROGRAM)
  4. C  NDYRI    DAY OF YEAR 
  5. C            =>0, DAY OF YR AT MID POINT IN RUN(USED FOR VALIDATION)
  6. C            =-1, DAY OF YR WILL BE 21ST OF CURRENT MONTH DEFINED IN 
  7. C             MTHS ARRAY. 
  8. C            =-2, DAY OF YR WILL BE CALC FROM MONTH AND DAY ON TAPE19
  9. C  ICC    =1  CLOUD COVER DATA USED TO COMPUTE SOLAR
  10. C         =0  TOTAL SOLAR RADIATION SUPPLIED IN  CC  ARRAY
  11. C  KM       MONTH OF SOLAR CALCULATION
  12. C  IDAY    DAY OF THIS SIMULATION CASE (1....)
  13. C  DST    =0 FOR STANDARD TIME
  14. C         =1 FOR DAYLIGHT SAVINGS TIME
  15. C 
  16. C  INPUTS IN COMMON BLOCKS ENCBLK,WNDBLK,ROOF ARE DEFINED IN READIN 
  17. C 
  18. C  INPUTS IN COMMON BLOCKS: 
  19. C  /SURFAR/  VALUES CALC IN CONSLD SUBPROGRAM:  
  20. C     AWALLA  SURFACE AREA, SQ FT 
  21. C     AWAWA   WINDOW AREA, SQ FT
  22. C     RAREA   TOTAL ROOF AREA OVER ATTICS 1 AND 2, SQ FT
  23. C 
  24. C  /SHDF/  SHDFCT  SHADING FACTOR(0.-1.). INDICES ARE 
  25. C                  (COMPASS DIR, OUTSIDE WALL,   ZONE NO.,HR OF DAY)
  26. C                              , WINDOWS IN WALL, 
  27. C                              , ROOF,
  28. C           COMPASS DIR: 1-E,2-S,3-S!!
  29. C         SHDFCT IS ONLY IMPLEMENTED FOR ONE ZONE AS YET.  MAKE
  30. C         CHANGES IN 3TH DIMENSION LATER IN 4 PLCS WHEN IMPLEMENTED
  31. C         FOR MORE THAN ONE ZONE.
  32. C 
  33. C  /NCR/   IFXDS   COMPASS DIRECTIONS ALWAYS SHADED (1-4) 
  34. C 
  35. C  /SIMA/  REFLG   GROUND REFLECTANCE FOR EACH MONTH (0.-1.)
  36. C
  37. C  /CITINF/ DNLAT  DEG. N. LATTITUDE
  38. C           DLONG  DEG. LONGITUDE 
  39. C           CLEAR  CLEARNESS NO.  (USED IF ICC=1) 
  40. C 
  41. C -- OUTPUTS - - - - - - - - - - - -
  42. C  SOLARH - 24 VALUES OF TOTAL HORIZONTAL SOLAR RADIATION FOR THE DAY
  43. C  DATOUT (I,IHOUR) 
  44. C    I=1  SOLAR LOADING ON ROOF OVER ATTIC NO. 1
  45. C      2  SOLAR LOADING ON ROOF OVER ATTIC NO. 2
  46. C      3 TO (NROOMS+2) TOTAL SOLAR LOAD TRANSMITTED THROUGH WINDOWS 
  47. C                      INTO ZONES NRMA(J) WHERE J=I-2 
  48. C      (NROOMS+3) TO (NROOMS+3+4)SOLAR LOADING ON OUTER SURFACE OF
  49. C                      OUTSIDE WALLS FACING EACH ID COMPASS DIRECTION
  50. C 
  51. C 
  52. C - - - - CONSTANTS 
  53. C   IF IWALS = 0  NIWL = 2 * NROOMS(FOR SOLAR THRU WINDOWS AND ON WALLS)
  54. C                        + 2(FOR 2 ATTIC ROOFS) ... = 8
  55. C   IF IWALS = 1  NIWL = NROOMS(FOR SOLAR THRU WINDOWS) + 
  56. C                        2(FOR 2 ATTIC ROOFS) +
  57. C                        SUM OF NWALLA(I) FOR ALL ZONES I                  
  58. CMDK NZN
  59. CMDK NIWL
  60. CMDK NKONST
  61. CMDK NRF
  62. CMDK NWL
  63. CMDK NWN
  64. CMDK NWZN
  65. CMDK NZW
  66. CMDK NCR
  67. CMDK CITINF
  68. CMDK CNSTRK
  69. CMDK ENCBK1
  70. CMDK ENCBLK
  71. CMDK IWLS
  72. CMDK IZZQ
  73. CMDK LDPOUT
  74. CMDK NIJB
  75. CMDK ROOF
  76. CMDK SHDF
  77. CMDK SIMA
  78. CMDK SURFAR
  79. CMDK WNDBLK
  80. C 
  81.       REAL ITH,IDHC,IDC,IDFC,ITC
  82.       CHARACTER *4 STM
  83.       DIMENSION IDC(4),IDFC(4),DIRECT(4),DIFFUS(4),SFST(4),SFTT(4)
  84.      +,CLSLR(NZN),WSLR(NZW),TRSL(2),SLWL(4),WALS(4),NDYL(12)
  85. C 
  86. C  NDYL IS DAY OF YR AT END OF EACH MONTH 
  87.       DATA NDYL/31,59,90,120,151,181,212,243,273,304,334,365/ 
  88. C 
  89. C  IWALS  OUTPUT OPTION FOR WALL SOLAR LOADING
  90. C  IWALS IS SET IN NAMELIST CONTROL READ IN STRRUN
  91. C     =0  ONLY ONE VALUE FOR EACH OF FOUR COMPASS DIRECTIONS
  92. C     =1  VALUES FOR EACH CONSOLIDATED OUTSIDE-WALL ENCLOSURE ELEMENT 
  93. C            (=1 OPTION WILL REQUIRE TOO MUCH MEMORY FOR A 640K PC)
  94.       CALL CURS(19)
  95.       WRITE(*,513)
  96.       CALL CURS(19)
  97.       WRITE(*,*)'START CALC SOLAR DATA FOR DAY SEQUENCE NO.',
  98.      +          IDAY,' OCCURING IN MONTH ',KM
  99.       WRITE(60,*)'START CALC SOLAR DATA FOR DAY SEQUENCE NO.',
  100.      +          IDAY,' OCCURING IN MONTH ',KM
  101.       IDBGSV=IDBG 
  102.       IDBG=0
  103. C 
  104.       DO 10 I=1,NIWL
  105.       DO 10 J=1,24
  106.    10 DATOUT(I,J)=0.
  107.       DO 12 I=1,NZW 
  108.    12 NIJA(I)=0 
  109.       IF(IWALS.EQ.0)THEN
  110. C  SUM NET WALL AREA IN EACH COMPASS DIRECTION
  111.           DO 18 I=1,4 
  112.           WALS(I)=1.E-13
  113.    18     CONTINUE
  114.           DO 20 J=1,NROOMS
  115.           NRM=NRMA(J) 
  116.           N2=NWALLA(NRM)
  117.           DO 20 K2=1,N2 
  118.           IZW=NENC(NRM,K2)
  119.           ID=IDEXP(IZW) 
  120.           IF(ID.LE.4)WALS(ID)=WALS(ID)+AWALLA(IZW)
  121.    20     CONTINUE
  122.           ENDIF 
  123. C 
  124.       ITIME=-60 
  125. C  READ 24 HRS OF SOLAR OR CLOUD COVER DATA
  126.       NDYRIS=NDYRI
  127.       IF(NSEQW.EQ.3)THEN
  128. C
  129. C NOTE: TAPE 19 WITH SOLARH DATA IS USED TO MINIMIZE MEMORY REQMTS
  130. C       DURING EXECUTION.  IT IS WRITTEN IN WETHR.FOR OR WWEATH.FOR      
  131.         READ(19,*,END=990)KMN,IDYMN
  132.         WRITE(60,*)' SOLRLD FR TAPE19: KMN,IDYMN= ',KMN,IDYMN
  133.         IF(KMN.NE.KM)THEN
  134.           WRITE(*,*)' SOLRLD: KMN FROM TAPE19,KM,IDYMN= ',
  135.      +              KMN,KM,IDYMN
  136.           WRITE(60,*)' SOLRLD: KMN FROM TAPE19,KM,IDYMN= ',
  137.      +               KMN,KM,IDYMN
  138.           STOP ' SOLRLD: KMN FROM TAPE19 NE KM!'
  139.         ENDIF
  140.         IF(NDYRI.EQ. -2)THEN
  141. C  COMPUTE NDYRI FROM DATA ON TAPE19 WHEN NDYRI=-2!        
  142.           IF(KMN.EQ.1)THEN
  143.             NDYRI=IDYMN
  144.           ELSE
  145.             NDYRI=NDYL(KMN-1)+IDYMN
  146.           ENDIF
  147.         ENDIF
  148.       ENDIF
  149.       READ(19,533,END=991)(SOLARH(IHR),IHR=1,24)
  150. C      
  151. C  SET OUTPUT ARRAY SUBSCRIPT FOR HOUR TO ZERO
  152.       IHOUR=0 
  153. C 
  154.    42 ITIME=ITIME+60
  155.       IHOUR=IHOUR+1 
  156.       TIME=FLOAT(ITIME)/60. 
  157.       ISUN=0
  158.       IF(TIME.GT.23.98)GO TO 9996 
  159. C 
  160. C  GET SOLAR DATA FOR THIS HOUR AND DAY 
  161.       CALL SOLR(TIME,NDYRI,IDAY,KM,STM,DNLAT,DLONG, 
  162.      +          DST,ICC,CLEAR,
  163.      +          ISUN,ITH,IDHC,SAZM,SALT,
  164.      +          IDC,IDFC,DIRECT,DIFFUS,SFST,SFTT) 
  165.       SOLARH(IHOUR)=ITH
  166. C 
  167.       DO 60 I=1,NZN 
  168.   60  CLSLR(I)=0. 
  169.       DO 61 IJ=1,4
  170.       SLWL(IJ)=0. 
  171.    61 CONTINUE
  172.       IWAL=0
  173.       DO 62 K=1,NZW 
  174.       WSLR(K)=0.
  175.    62 CONTINUE
  176. C 
  177. C  LOOP FOR EXPOSED WALLS IN EACH ZONE
  178.   
  179.       DO 295 J=1,NROOMS 
  180.       NRM=NRMA(J) 
  181.       N2=NWALLA(NRM)
  182.       DO 295 K2=1,N2
  183.       IZW=NENC(NRM,K2)
  184.       ID=IDEXP(IZW) 
  185. C  SUN CAN ONLY SHINE ON OUTSIDE WALLS
  186. C   NOTE:  ABOVE GROUND BASEMENT WALLS ARE ASSUMED TO BE SHADED 
  187.       IF(ID.GT.4)GO TO 295
  188.       M=ID
  189.       XLEN=XLENA(IZW) 
  190.       HGT=HGTA(IZW) 
  191.       AWALC=XLEN*HGT
  192.       OHANG=OHANGA(IZW) 
  193.       OHNG=OHANG*12.
  194.       HHANG=HHANGA(IZW) 
  195.       NWND=NWNDA(IZW) 
  196.       EPSLM=EPSLR(IZW)
  197.       KON=KONSTA(IZW) 
  198.       WNDSA=0.
  199.       WNDUA=0.
  200.       IF(ISUN.GT.0)GO TO 251
  201. C **************     3RD DIMENSION ON SHDFCT SHOULD BE NRM LATER
  202. C                    WHEN SHDFCT APPLIES TO > 1 ZONE      
  203.       SHDFC=1.
  204.       IF(ID.GT.1)SHDFC=SHDFCT(ID-1,1,1,IHOUR)
  205.       ITC=IDC(ID)*SHDFC+IDFC(ID) 
  206.       IF(NWND.EQ.0)GO TO 241
  207. C 
  208.       DO 240 K4=1,NWND
  209.       NI=NIWND(IZW,K4)
  210.       WHGT=WHGTA(NI)
  211.       WWDTH=WWDTHA(NI)
  212.       REVEAL=REVELA(NI) 
  213.       OHNGH=OHNGHA(NI)
  214.       SFSD=SFSDA(NI)
  215.       WAREA=AWAWA(NI) 
  216. C  IF ANY WALL IS ALWAYS SHADED, GO TO 210
  217.       DO 206 IJ=1,4 
  218.       IF(M.EQ.IFXDS(IJ))GO TO 210 
  219.   206 CONTINUE
  220.       IF (SALT.EQ.0.0 .AND.SAZM.EQ. 0.0) GO TO 210
  221.       WP=REVEAL *SFST(ID)/12. 
  222.       IF (WP.GT.WWDTH) WP=WWDTH 
  223.       IF(WP.LT.0.0)WP=0.0 
  224.       HP=((REVEAL+OHNG)*SFTT(ID)-OHNGH)/12. 
  225.       IF (HP.GT.WHGT) HP=WHGT 
  226.       IF(SFTT(ID).EQ.0.0)HP=0.0 
  227.       IF(HP.LT.0.0)HP=0.0 
  228.       SA=WWDTH*HP+WHGT*WP-WP*HP 
  229.       GO TO 215 
  230. 210   SA=WAREA
  231.       WNDSA=WNDSA+WAREA 
  232.       GO TO 216 
  233. C 
  234. C---- TOTAL SHADED AREA FROM OVERHANG USED FOR FINDING SHADED WALL AREA 
  235. 215   HPP=(OHNG*SFTT(ID)-OHNGH)/12. 
  236.       IF(HPP.LT.0.0)HPP=0.0 
  237.       WNDSA=WNDSA+WWDTH*HPP 
  238. C 
  239. C---- TOTAL UNSHADED WINDOW AREA USED FOR FINDING UNSHADED WALL AREA
  240.       WNDUA=WNDUA+(WAREA-WWDTH*HPP) 
  241. C 
  242. C     CALC INSTANTANEOUS HEAT TRANSMISSION THRU WINDOWS 
  243. C 
  244.   216 SLA=(WAREA-SA)/WAREA
  245.       DIRCT=SLA*DIRECT(M) 
  246.       DIFFS=DIFFUS(M) 
  247.       IF(SFSD.LE.0.)SFSD=1. 
  248. C 
  249. C  TOTAL SOLAR RADIATION TRANSMITTED THRU WINDOW
  250. C  3RD DIMENSION ON SHDFCT SHOULD BE NRM LATER WHEN SHDFCT > 1 ZONE
  251.       SHDFC=1.
  252.       IF(M.GT.1)SHDFC=SHDFCT(M-1,2,1,IHOUR)
  253.       COOL=SFSD*WAREA*(DIRCT*SHDFC+DIFFS) 
  254.       IF(IDBG.NE.0)WRITE(60,503)M,IHOUR,COOL,SFSD,WAREA,SLA
  255.       CLSLR(NRM)=CLSLR(NRM)+COOL
  256.   240 CONTINUE
  257. C 
  258. C  END LOOP FOR WINDOWS IN THIS WALL
  259. C 
  260. C  SOLAR LOADING ON WALL
  261.   241 CONTINUE
  262. C 
  263.       DO 242 I=1,4
  264. C  IF ANY WALL ALWAYS SHADED, GO TO 245 
  265.       ISWA=1
  266.       IF(M.EQ.IFXDS(I))GO TO 245
  267. C  CHECK FOR ALWAYS SHADED CONSTRUCT
  268.       IFX=IFXL(I,KON) 
  269.       ISWA=2
  270.       IF(M.EQ.IFX)GO TO 245 
  271.   242 CONTINUE
  272. C 
  273. C  GET SHADING FROM OVERHANG (ASSUMES SHDFCT'S WERE CALC W/O OVERHANG 
  274.       IF(OHANG.LE.0.0)GO TO 243 
  275.       ISWA=3
  276.       IF (SALT.EQ.0.0 .AND.SAZM.EQ. 0.0) GO TO 245
  277.       HP=(OHNG*SFTT(ID)-HHANG)/12.
  278.       IF(HP.GT.HGT)HP=HGT 
  279.       IF(SFTT(ID).EQ.0.0)HP=0.
  280.       IF(HP.LT.0.0)HP=0.
  281.       SWA=XLEN*HP 
  282.       IF(SWA.LT.0.0)SWA=0.0 
  283.       ISWA=4
  284.       IF(SWA.GE.AWALC) GO TO 245
  285. C 
  286. C---- NET SHADED WALL AREA=TOTAL SHADED WALL AREA - WINDOW AREA 
  287. C---- SHADED FROM OVERHANG
  288. C 
  289.       SWANET=SWA-WNDSA
  290.       IF(SWANET.LT.0.0)GO TO 243
  291. C 
  292. C---- NET UNSHADED WALL AREA=TOTAL WALL AREA-TOTAL SHADED WALL AREA-
  293. C---- TOTAL WINDOW AREA NOT SHADED BY OVERHANG
  294. C 
  295.       UWANET=AWALC-SWA-WNDUA
  296.       GO TO 250 
  297. C 
  298. C---- NO OVERHANG OR SHADED WALL AREA 
  299. 243   SWANET=0.0
  300.       UWANET=AWALC-WNDUA-WNDSA
  301.       GO TO 250 
  302. C 
  303. C---- ENTIRE WALL IS SHADED 
  304. C 
  305. 245   UWANET=0.0
  306.       SWANET=AWALC-WNDUA-WNDSA
  307.       IF(IDBG.NE.0)WRITE(60,506)ISWA,OHANG,SFTT(ID),XLEN,HGT,SWANET
  308. C 
  309. C CALC INSTANTANEOUS SOLAR RAD INTENSITIES FOR SHADED & 
  310. C     UNSHADED AREAS
  311. 250   HIS=IDFC(ID)
  312.       HIUS=ITC
  313.       GO TO 252 
  314. C 
  315. C---- SOLAR LOAD ON WALLS (BTUH)
  316. C 
  317.   251 SLWALL=0. 
  318.       GO TO 253 
  319.   252 CONTINUE
  320.       SLWALL=(HIUS*UWANET+HIS*SWANET)*EPSLM 
  321.       IF(IDBG.NE.0)THEN 
  322.       WRITE(60,504)KM,IZW,M,IHOUR,SLWALL,HIUS,HIS,UWANET,SWANET,EPSLM
  323.           ENDIF 
  324.   253 WSLR(IZW)=SLWALL/AWALLA(IZW)
  325.       SLWL(ID)=SLWL(ID)+SLWALL
  326.       IF(IWALS.EQ.1)THEN
  327.         IWAL=IWAL+1 
  328.         IF(IWAL.GT.(NIWL-2-NROOMS))THEN 
  329.           WRITE(60,502)IWAL,NIWL,NROOMS
  330.           STOP 'SOLRLD:NIWL MUST BE INCR' 
  331.           ENDIF 
  332.         ENDIF
  333.   295 CONTINUE
  334. C 
  335. C ** ** ** END LOOP FOR WALLS IN ZONE AND LOOP FOR ZONES
  336. C 
  337. C  SOLAR LOADING ON ROOF SECTIONS FOR ROOFS ABOVE TWO ATTICS
  338.       TRSL(1)=0.
  339.       TRSL(2)=0.
  340.       IF(ISUN.GT.0)GO TO 364
  341. C 
  342.       DO 360 I=1,NSECTA 
  343.       M=IDEXPR(I) 
  344.       NAT=NATTC(I)
  345. C  3RD DIMENSION ON SHDFCT BELOW SHOULD BE NAT LATER WHEN SHDFCT > 1 
  346. C      ZONE      
  347.       SHDFC=1.
  348.       IF(M.GT.1)SHDFC=SHDFCT(M-1,3,1,IHOUR)
  349.       SLRS=(IDR(I)*SHDFC+IDIFR(I))*SECTAA(I)*EPRSLR(I)
  350.       IF(IDBG.NE.0)THEN 
  351.       WRITE(60,505)I,IDR(I),SHDFC,IDIFR(I),SECTAA(I),EPRSLR(I),SLRS
  352.           ENDIF 
  353.       IF(NAT.EQ.1)TRSL(1)=TRSL(1)+SLRS
  354.       IF(NAT.EQ.2)TRSL(2)=TRSL(2)+SLRS
  355.   360 CONTINUE
  356. C 
  357.   364 TRSL(1)=TRSL(1)/RAREA(1)
  358.       IF(RAREA(2).NE.0)THEN 
  359.           TRSL(2)=TRSL(2)/RAREA(2)
  360.           ENDIF 
  361. C 
  362.       IF(TRSL(1).GT.0.)THEN 
  363.           IF(IDBG.NE.0)WRITE(60,509)1,IHOUR,IDAY,KM,TRSL(1)
  364.   509 FORMAT(1X,'SOLRLD:NIJ,IHR,IDY,KM,TRSL= ',4I5,G13.5) 
  365.           DATOUT(1,IHOUR)=TRSL(1)
  366.           ENDIF 
  367.       IF(TRSL(2).GT.0.)THEN 
  368.           IF(IDBG.NE.0)WRITE(60,510)2,IHOUR,IDAY,KM,TRSL(2)
  369.   510 FORMAT(1X,'SOLRLD:NIJ,IHR,IDY,KM,TRSL= ',4I5,G13.5) 
  370.           DATOUT(2,IHOUR)=TRSL(2)
  371.           ENDIF 
  372.       NIJ=2 
  373.       DO 365 J=1,NROOMS 
  374.       NRM=NRMA(J) 
  375.       NIJ=NIJ+1 
  376.       IF(IDBG.NE.0)WRITE(60,511)NRM,NIJ,IHOUR,IDAY,KM,CLSLR(NRM) 
  377.   511 FORMAT(1X,'SOLRLD:NRM,NIJ,IHR,IDY,KM,CSLR= ',5I5,G13.5) 
  378.       DATOUT(NIJ,IHOUR)=CLSLR(NRM) 
  379.   365 CONTINUE
  380.       IF(IWALS.EQ.1)THEN
  381.           DO 366 J=1,NROOMS 
  382.           NRM=NRMA(J) 
  383.           N2=NWALLA(NRM)
  384. C 
  385.           DO 366 K2=1,N2
  386.           IZW=NENC(NRM,K2)
  387.           ID=IDEXP(IZW) 
  388. C  SUN CAN ONLY SHINE ON OUTSIDE WALLS
  389.           IF(ID.GT.4)GO TO 366
  390.           NIJ=NIJ+1 
  391.           IF(NIJ.GT.NIWL)THEN 
  392.               WRITE(60,501)NIJ,NIWL
  393.               STOP 'SOLRLD:NIWL AT LDPOUT COMMON BLK MUST BE INCR'
  394.               ENDIF 
  395.           NIJA(IZW)=NIJ 
  396.           IF(IDBG.NE.0)WRITE(60,512)NRM,NIJ,IHOUR,IDAY,KM,WSLR(IZW)
  397.   512     FORMAT(1X,'SOLRLD:NRM,NIJ,IHR,IDY,KM,WSLR= ',5I5,G13.5) 
  398. C  SOLAR ON EACH CONSOLIDATED OUTSIDE-WALL ENCLOSURE ELEMENT
  399.           DATOUT(NIJ,IHOUR)=WSLR(IZW)
  400.   366     CONTINUE
  401.       ELSE
  402. C       FIND AVERAGE SOLAR LOADING ON WALLS IN EACH DIRECTION 
  403.           DO 368 ID=1,4 
  404.           SWSLR=SLWL(ID)/WALS(ID) 
  405.           NIJ=NIJ+1 
  406.           IF(NIJ.GT.NIWL)THEN 
  407.               WRITE(60,501)NIJ,NIWL
  408.           STOP 'SOLRLD:NIWL AT LDPOUT COMMON BLK MUST BE INCR'
  409.               ENDIF 
  410.           NIJA(ID)=NIJ
  411.           DATOUT(NIJ,IHOUR)=SWSLR
  412.   368     CONTINUE
  413.           ENDIF 
  414.       GO TO 42
  415. C 
  416.  9996 CONTINUE
  417.       NDYRI=NDYRIS
  418. C  END OF HOURLY LOOP * * * * * * * * * 
  419.       IDBG=IDBGSV 
  420.       CALL CURS(19)
  421.       WRITE(*,513)
  422.       CALL CURS(19)
  423.       WRITE(*,*)'  END CALC SOLAR DATA FOR DAY SEQUENCE NO.',
  424.      +          IDAY,' OCCURING IN MONTH ',KM
  425.       RETURN
  426.   990 WRITE(60,*)' SOLRLD: PROBLEM READING MONTH AND DAY FROM TAPE19'
  427.       STOP ' SOLRLD: PROBLEM READING MONTH AND DAY FROM TAPE19'
  428.   991 WRITE(60,*)' SOLRLD: PROBLEM READING HORIZ SOLAR DATA FROM TAPE19'
  429.       STOP ' SOLRLD: PROBLEM READING HORIZ SOLAR DATA FROM TAPE19'
  430.   500 FORMAT(1X)
  431.   501 FORMAT(1H0,'SOLRLD:NIWL AT LDPOUT COMMON BLK MUST BE INCR', 
  432.      +/'  NIJ,NIWL= ',2I5)
  433.   502 FORMAT(1H0,'OOPS!!SOLRLD:NIWL MUST BE INCREASED. IWAL,NIWL',
  434.      +'NROOMS= ',3I5) 
  435.   503 FORMAT(1X,'ID= ',I3,'  IHOUR= ',I3,'  COOL,SFSD,WAREA,SLA= ', 
  436.      +4G13.5) 
  437.   504 FORMAT(1X,'SOLRLD:KM= ',I3,' IZW= ',I3,' ID= ',I3,'  IHOUR= ',I3/ 
  438.      +' SLWALL,HIUS,HIS,UWANET,SWANET,EPSLM= ',6G13.5)
  439.   505 FORMAT(1X,'I,IDR,SHDFC,IDIFR,SECTA,EPRSLR,SLRS= ',I4,6G13.5)
  440.   506 FORMAT(1X,'ISWA,SFT,XLEN,HGT,SWANET= ',I3,5G13.5) 
  441.   507 FORMAT(1H0,'NOTE!!!! NMONTH MAX IS 5 IN SOLRLD')
  442.   513 FORMAT(80(1H ))
  443.   533 FORMAT(1X,6G13.5/1X,6G13.5/1X,6G13.5/1X,6G13.5)
  444.       END 
  445.