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...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Revised April, 1999 VIII.13<strong>DOE</strong>-<strong>2.1</strong>E Documentation Update: Weather ProcessorC FOR THE 1 HOUR TIME BINA = SIN(DEABC(1))*SIN(STALAT)B = COS(DEABC(1))*COS(STALAT)COSZIN = A*(UL-BL) + B*(SIN(PIOV12*UL)-SIN(PIOV12*BL))/PIOV12COSZAV = COSZIN/(UL-BL)C GET DIRECT NORMAL SOLARSOLDRN = .31721*FLOAT(IDIRH)/COSZAVC PUT LIMITS ON THE SUNRISE AND SUNSET BEAM RADIATIONCALL MAXDIR(COSZAV,SOLCON,DIRMAX)SOLDRN = AMIN1(SOLDRN,DIRMAX)300 CONTINUEC FILL THE DATA ARRAYSIDRY(IH) = IROUND(TDRYF)IWET(IH) = MIN0(IDRY(IH),IROUND(TWETF))IDEW(IH) = MIN0(IWET(IH),IROUND(TDEWF))IPRESS(IH) = IROUND(100.*PRESHG)IWNDSP(IH) = IROUND(WSKNOT)IWNDIR(IH) = IROUND(.0444444*FLOAT(IWINDR))IF (IWNDIR(IH) .EQ. 16) IWNDIR(IH) = 0ICLAMT(IH) = IROUND(10.*SKYCVR)ISOL(IH) = IROUND(SOLHOR)IDN(IH) = IROUND(SOLDRN)ICLTY(IH) = 2ICLTY1(IH) = 2IRN(IH) = 0ISN(IH) = 01000 CONTINUERETURNENDFORTRAN fragment: direct and diffuse solar calculationCC Fortran fragment to calculate direct and diffuseC solar radiation from total solar and to get wet-bulbC and dew point temperatures from dry-bulb, relative humidityC and atmospheric pressure.CC Start with: SOLHOR - total horiz. solar in Btu/hr-areaC TDRYF - dry-bulb temperature in degrees FahrenheitC PRESHG - atmospheric pressure in inches of mercuryC RELHUM - relative humidity in percentCC Define: DTOR - degrees to radians = pi/180C PIOV12 - pi/12C IH - hour of the day (1 - 24)C IDOY - day of year (1 - 365)C STALAT - weather station latitude in radiansC STALAT - station longitude in radiansC XLONG - station longitude in degreesC NTZ - time zone (PST=8, EST=5, Greenwich=0)CC Externals: SUNPRMC MAXDIRC PPWVMSC WBFCThese are all available from WTHPRC, the <strong>DOE</strong>-2Cweather processor.CDIMENSION DEABC(5)CDIRN = 0.DIF = 0.C GET SOLAR CONSTANTS. DEABC(1) IS THE SOLARC DECLINATION ANGLE; DEABC(2) IS THE EQUATIONC OF TIME.CALL SUNPRM(IDOY,DEABC)

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

Saved successfully!

Ooh no, something went wrong!