12.07.2015 Views

DOE-2 Reference Manual Version 2.1 - DOE2.com

DOE-2 Reference Manual Version 2.1 - DOE2.com

DOE-2 Reference Manual Version 2.1 - DOE2.com

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>DOE</strong>-<strong>2.1</strong>E Documentation Update: Weather ProcessorC GET SOLAR CONSTANT. THIS FORMULA IS FROM KREIDERC AND RABL PAGE 237C BECKMAN, PAGE 7.SOLCON = 435.2*(1. + 0.033*COS(DTOR*360.*FLOAT(IDOY)/365.))C GET HOUR ANGLE AT UPPER AND LOWER HOUR BIN EDGEC THIS IS LOCAL TIME!UL = FLOAT(IH) - 12. + FLOAT(NTZ) + DEABC(2) - XLONG/15.BL = UL - 1.C SUNRISE AND SUNSET HOUR ANGLESSSHA = ACOS(-TAN(STALAT)*TAN(DEABC(1)))/PIOV12SRHA = -SSHAC RESET BIN BOUNDARIES TO ALLOW FOR SUNRISE AND SETIF ((UL .LE. SRHA) .OR. (BL .GE. SSHA)) GO TO 1100IF (SRHA .GT. BL) BL = SRHAIF (SSHA .LT. UL) UL = SSHAIF ((UL-BL) .LT. 0.02) GO TO 1100IF (SOLHOR .EQ. 0) GO TO 1100A = SIN(DEABC(1))*SIN(STALAT)B = COS(DEABC(1))*COS(STALAT)C INTEGRATE SOLAR Z DIREC. COSINE OVER BINCOSZIN = A*(UL-BL) + B*(SIN(PIOV12*UL)-SIN(PIOV12*BL))/PIOV12C AVERAGE COSINE OF THE SOLAR ZENITH ANGLE FOR THE HOURCOSZAV = COSZIN/(UL-BL)C EXTRATERRESTRIAL SOLAR HORIZONTALSOLEXH = SOLCON*COSZINC K sub T IS THE RATIO OF TERRESTRIAL TO EXTRATERRESTRIAL SOLARRKT = AMIN1(SOLHOR/SOLEXH,0.9)C GET DIFFUSE COMPONENT FROM ERBS, KLEIN, AND DUFFIEC CORRELATIONIF (RKT .LE. 0.22) DIF = SOLHOR*(1.-0.09*RKT)RKT2 = RKT*RKTIF ((RKT .GT. 0.22) .AND. (RKT .LE. 0.8)) DIF = SOLHOR*1 (.9511-.1604*RKT+4.388*RKT2-16.638*RKT*RKT2+12.336*RKT2*RKT2)IF (RKT .GT. 0.8) DIF = 0.165*SOLHORC DIRECT HORIZONTALDIRH = AMAX1(0.,SOLHOR-DIF)C DIRECT NORMALDIRN = DIRH/COSZAVC CHECK FOR MAXIMUM DIRECT NORMALCALL MAXDIR(COSZAV,SOLCON,DIRMAX)DIRN = AMIN1(DIRN,DIRMAX)1100 CONTINUEC CALCULATE WET-BULB AND DEW POINT. THE PROCEDURE ISC BASICALLY THAT GIVEN ON PAGE 6.16, ASHRAE FUNDAMENTALSC 1989, SITUATION NO. 3CC SATURATED VAPOR PRESSUREPS = PPWVMS(TDRYF)C PARTIAL PRESSUREPW = .01*RELHUM*PSC HUMIDITY RATIOHUMRAT = .622*PW/(PRESHG-PW)C SPECIFIC ENTHALPYENTH = .24*TDRYF + (1061.+.444*TDRYF)*HUMRATC WET-BULB TEMPERATURETWETF = WBF(ENTH,PRESHG)C DEW POINT TEMPERATUREY = LOG(PW)IF (PW .LE. .1836) THENTDEWF = 71.98 + 24.873*Y + .8927*Y*YELSETDEWF = 79.047 + 30.579*Y + 1.8893*Y*YEND IFRevised April, 1999 VIII.14

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!