1 !---------------------------------------------------------------------- 2 SUBROUTINE WAKE_INI(LDBHR,KST,RURAL,DSBH,DSBW,DSBL,XADJ,YADJ, & 3 & UBLDG,USTACK) 4 !---------------------------------------------------------------------- 5 ! 6 ! --- PRIME Version: 1.0 Level: 970812 WAKE_INI 7 ! D. Strimaitis, EARTH TECH 8 ! Prepared for EPRI under contract WO3527-01 9 ! 10 ! --- PURPOSE: Refreshes variables in /wakedat/ common 11 ! 12 ! --- INPUTS: 13 ! 14 ! LDBHR - logical - Debug output written when .TRUE. 15 ! KST - integer - PG stability class 16 ! RURAL - logical - Denotes rural dispersion when .TRUE. 17 ! DSBH - real - Effective building height (m) 18 ! DSBW - real - Effective building width (m) across flow 19 ! DSBL - real - Effective building length (m) along flow 20 ! XADJ - real - Distance (m) from source to upwind face of bldg 21 ! along flow 22 ! YADJ - real - Distance (m) from source to center of upwind 23 ! face of bldg across flow 24 ! UBLDG - real - Wind speed (m/s) at top of building 25 ! USTACK - real - Wind speed (m/s) at release height 26 ! 27 ! Parameters: 28 ! MXNTR 29 ! 30 ! --- OUTPUT: 31 ! 32 ! Common block /WAKEDAT/ variables: 33 ! HB,WB,XLB,RSCALE,HR,XLR,XLC,XBADJ,YBADJ, 34 ! NWAK, XWAK(mxntr), SZWAK(mxntr), SYWAK(mxntr), 35 ! DRWAK(mxntr), XZVWAK, XYVWAK, UB, URH, 36 ! LRURL, ISTAB 37 ! 38 ! --- WAKE_INI called by: PCALC (HOST subroutine) 39 ! --- WAKE_INI calls: WAKE_SCALES 40 !---------------------------------------------------------------------- 41 ! --- Include parameters 42 INCLUDE 'params.pri' 43 44 ! --- Include common blocks 45 INCLUDE 'wakedat.pri' 46 47 LOGICAL RURAL , LDBHR 48 DATA ZERO/0.0/ 49 50 ! --- Transfer arguments to /wakedat/ variables 51 ISTAB = KST ! 6088 52 LRURL = RURAL 53 HB = DSBH 54 WB = DSBW 55 XLB = DSBL 56 XBADJ = XADJ 57 YBADJ = YADJ 58 UB = UBLDG 59 URH = USTACK 60 61 ! --- Compute wake dimensions and related parameters 62 CALL WAKE_SCALES(LDBHR) 63 64 ! --- Reset contents of sigma arrays for wake region 65 NWAK = 1 66 XWAK(1) = ZERO 67 SZWAK(1) = ZERO 68 SYWAK(1) = ZERO 69 DRWAK(1) = ZERO 70 71 ! --- Reset virtual distances for sigmas beyond wake 72 XZVWAK = ZERO 73 XYVWAK = ZERO 74 75 CONTINUE 76 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