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