DOE-2 Reference Manual Version 2.1 - DOE2.com

CRevised April, 1999 VIII.12DOE-2.1E Documentation Update: Weather Processor2 IWNDIR(24),ICLTY(24),IRN(24),ISN(24),ICLTY1(24)DIMENSION IRAW(24,12)EQUIVALENCE (IDRY(1),IRAW(1,1))COMMON /LSTIME/ LSTHRS(24)COMMON /REPORC/ IBEGM,IENDMCCOMMON /UNDEF/ IUNDEF,UNDEFCCDIMENSION DEABC(5)C DAY OF YEARIDOY = BEFORE(IMNTH) + IDAYC GET SUN PARAMETERSCALL SUNPRM(IDOY,DEABC)C SOLAR CONSTANTSOLCON = 435.2*(1.+0.033*COS(DTOR*360.*FLOAT(IDOY)/365.))C LOOP OVER HOURS IN THE DAYDO 1000 IH=1,24C READ IN WEATHER DATAREAD (INWTH,9001) IVR,ID,IM,IHR,SKYCVR,IWINDR,WNDSPD,WSVECT,1 IWTHR,DEP,PRESMB,TDRY,RELH,IDIRH,IDIFFH,ILLUMH,IRH,IATMRH9001 FORMAT(I1,3I2,1X,F4.2,1X,I3,2(1X,F4.1),1X,I2,1X,F4.1,1X,F6.1,1 1X,F5.1,1X,F4.2,1X,I4,1X,I4,1X,I7,1X,I4,1X,I4)C CONVERT DRY-BULB FROM CENTIGRADE TO FAHRENHEITTDRYF = 1.8*TDRY + 32.C CALCULATE WET-BULB AND DEW POINTC SATURATED VAPOR PRESSUREPS = PPWVMS(TDRYF)C PARTIAL PRESSUREPW = RELH*PSC CONVERT PRESSURE FROM MILLIBARS TO INCHES OF HGPRESHG = .02953*PRESMBC HUMIDITY RATIOHUMRAT = 0.622*PW/(PRESHG-PW)C SPECIFIC ENTHALPYENTH = 0.24*TDRYF + (1061.+0.444*TDRYF)*HUMRATTWETF = WBF(ENTH,PRESHG)Y = LOG(PW)IF (PW .LE. 0.1836) THENTDEWF = 71.98 + 24.873*Y + 0.8927*Y*YELSETDEWF = 79.047 + 30.579*Y + 1.8893*Y*YEND IFC CONVERT WINDSPEED FROM M/S TO KNOTSWSKNOT = 1.9438*WNDSPDSOLHOR = 0.SOLDRN = 0.CCCCCCCCSET UPPER AND LOWER HOUR ANGLE BIN EDGES FOR SOLAR ZENITHANGLE CALCULATION. HOUR ANGLE IN UNITS OF HOURS.UL = FLOAT(IH) - 12. + FLOAT(NTZ) + DEABC(2) - STALON/PIOV12BL = UL -1.SUNRISE AND SUNSET HOUR ANGLESSSHA = ACOS(-TAN(STALAT)*TAN(DEABC(1)))/PIOV12SRHA = -SSHASKIP IF SUN DOWNIF ((UL .LE. SRHA) .OR. (BL .GE. SSHA)) GO TO 300RESET BIN EDGES AT SUNRISE OR SUNSETIF (SRHA .GT. BL) BL = SRHAIF (SSHA .LT. UL) UL = SSHAIF ((UL-BL) .LT. .02) GO TO 300TOTAL HORIZONTAL SOLAR; CONVERT FROM W/M**2 TOBTU/(FT**2)(HR)SOLHOR = .31721*FLOAT(IDIRH+IDIFFH)GET THE AVERAGE OF THE COSINE OF THE SOLAR ZENITH ANGLE