1 2 !----------------------------------------------------------------------- 3 SUBROUTINE WAKE_TURB(KSTINP,LRURL,XINP,L,RINP,TIZ,TIY) 4 !----------------------------------------------------------------------- 5 ! 6 ! --- PRIME Version: 1.0 Level: 970812 WAKE_TURB 7 ! L. Schulman, D. Strimaitis, EARTH TECH 8 ! Prepared for EPRI under contract WO3527-01 9 ! 10 ! --- PURPOSE: Calculates turbulence intensity as a function of 11 ! location within the wake from modified Weil (1996) 12 ! analytical expressions 13 ! 14 ! --- MODIFIED: For use with the AERMOD model. Added MAIN1 module to 15 ! access global AERMOD data. Modified ambiy and ambiz to 16 ! use AERMOD turbulence intensities. 17 ! R.W. Brode, PES, Inc. - 07/05/01 18 ! 19 ! --- INPUTS: 20 ! kst - integer - PG stability class (1-6) 21 ! lrurl - logical - Rural flag (T=Rural, F=Urban) 22 ! x - real - distance (m) from upwind bldg wall 23 ! L - real - dist (m) of downwind bldg wall from 24 ! upwind bldg wall 25 ! R - real - wake scaling length (m) 26 ! 27 ! Common block /DFSN/ variables: 28 ! wiz0,wiy0,wfz,wfy, 29 ! dua_ua,xdecay,xdecayi, 30 ! rurliz,rurliy,urbniz,urbniy 31 ! 32 ! --- OUTPUT: 33 ! 34 ! tiz - real - turbulence intensity sigw/u 35 ! tiy - real - turbulence intensity sigv/u 36 ! 37 ! --- WAKE_TURB called by: WAKE_DFSN 38 ! --- WAKE_TURB calls : none 39 !---------------------------------------------------------------------- 40 ! 41 42 USE MAIN1 43 44 INCLUDE 'dfsn.pri' 45 46 REAL L , SWFACT 47 LOGICAL LRURL 48 49 ! --- Misc. constants 50 DATA ONE/1.0/ , ZERO/0.0/ 51 52 ! --- Specify ambient turbulence intensities from AERMOD effective parameters 53 AMBIZ = SWEFF/UEFF !3843816 54 AMBIY = SVEFF/UEFF 55 56 ! --- Compute asymptotic turbulence intensity in far wake 57 FARIZ = AMIN1(WIZ0,AMBIZ) 58 FARIY = AMIN1(WIY0,AMBIY) 59 60 ! --- Compute turbulence intensity at position downwind of bldg 61 XML = AMAX1(ZERO,XINP-L) 62 XFAC = ONE/(((XML+RINP)/RINP)**XDECAY-DUA_UA) 63 TIZ = FARIZ*(ONE+((WFZ*WIZ0/FARIZ-ONE)+DUA_UA)*XFAC) 64 TIY = FARIY*(ONE+((WFY*WIY0/FARIY-ONE)+DUA_UA)*XFAC) 65 66 ! --- Interpolate turbulence intensity if over roof of bldg 67 IF ( XINP.LT.L ) THEN 68 XFRAC = XINP/L ! 147816 69 TIZ = AMBIZ + XFRAC*(TIZ-AMBIZ) 70 TIY = AMBIY + XFRAC*(TIY-AMBIY) 71 ENDIF 72 73 CONTINUE !3843816 74 END
HyperKWIC - Version 1.00DD executed at 20:00 on 1 Mar 2018 | Personal or Academic or Evaluation User | Free for Non-Commercial, Non-Government Use