1 2 SUBROUTINE PWIDTH(XARG,V1,VN,WIDTH) 3 !*********************************************************************** 4 ! PWIDTH Module of the AMS/EPA Regulatory Model - AERMOD 5 ! 6 ! PURPOSE: Calculates The Effective Area of The Plume for A 7 ! Certain Downwind Distance 8 ! 9 ! PROGRAMMER: Jeff Wang, Roger Brode 10 ! Adapted From Codes By Richard Strelitz, CSC 11 ! 12 ! DATE: July 7, 1993 13 ! 14 ! MODIFIED BY R. Brode, PES, Inc. to move calculation of dispersion 15 ! coefficients to a new ADIS subroutine - 7/21/94 16 ! 17 ! MODIFIED BY R. Brode, PES, Inc. to correct table of GA values 18 ! and extend GA to 79 values - 7/29/94 19 ! 20 ! INPUTS: Downwind Distance 21 ! Crosswind Distance 22 ! Lateral Dispersion Parameter 23 ! Vertical Dispersion Parameter 24 ! Receptor Height Above Ground 25 ! Source Parameter Arrays 26 ! 27 ! OUTPUTS: The Effective Width 28 ! 29 ! CALLED FROM: ZBRENT 30 ! PSIDE 31 ! PLUMEF 32 ! AREAIN 33 !*********************************************************************** 34 35 ! Variable Declarations 36 USE MAIN1 37 IMPLICIT NONE 38 CHARACTER MODNAM*12 39 40 SAVE 41 INTEGER :: ITEMP 42 REAL :: XARG , WIDTH , VN , V1 , TEMP , GA(79) 43 44 ! Variable Initializations 45 ! GA ARE VALUES OF THE CUMULATIVE NORMAL DISTRIBUTION IN 46 ! INCREMENTS OF 0.1 S. 47 DATA GA/0.0 , .0001 , .0001 , .0002 , .0002 , .0003 , .0005 , & 48 & .0007 , .0010 , .0013 , .0019 , .0026 , .0035 , .0047 , & 49 & .0062 , .0082 , .0107 , .0139 , .0179 , .0227 , .0287 , & 50 & .0359 , .0446 , .0548 , .0668 , .0808 , .0968 , .1151 , & 51 & .1357 , .1587 , .1841 , .2119 , .2420 , .2742 , .3085 , & 52 & .3445 , .3821 , .4207 , .4602 , .5000 , .5398 , .5793 , & 53 & .6179 , .6555 , .6915 , .7258 , .7580 , .7881 , .8159 , & 54 & .8413 , .8643 , .8849 , .9032 , .9192 , .9332 , .9452 , & 55 & .9554 , .9641 , .9713 , .9773 , .9821 , .9861 , .9893 , & 56 & .9918 , .9938 , .9953 , .9965 , .9974 , .9981 , .9987 , & 57 & .9990 , .9993 , .9995 , .9997 , .9998 , .9998 , .9999 , & 58 & .9999 , 1.000/ 59 MODNAM = 'PWIDTH' ! 18497K 60 61 IF ( XARG.EQ.0.0 ) THEN 62 SY = 1.0 ! 0 63 VN = V1 64 WIDTH = VN 65 ! Exit Routine 66 GOTO 999 67 ENDIF 68 69 ! Define plume centroid height (CENTER) for use in 70 ! inhomogeniety calculations --- CALL CENTROID 71 CALL CENTROID(XARG) ! 18497K 72 73 ! If the atmosphere is unstable and the stack 74 ! top is below the mixing height, calculate 75 ! the CBL PDF coefficients --- CALL PDF 76 IF ( UNSTAB .AND. (HS.LT.ZI) ) CALL PDF 77 78 ! Determine Effective Plume Height --- CALL HEFF 79 CALL HEFF(XARG) 80 81 ! Iterative average through plume rise layer 82 CALL IBLVAL(XARG) 83 84 ! Calculate lateral dispersion coefficient, SY --- CALL ADISY 85 CALL ADISY(XARG) 86 87 VN = V1/SY 88 TEMP = 10*VN + 40 89 ITEMP = INT(TEMP) 90 WIDTH = 0.0 91 92 IF ( ITEMP.GT.78 ) THEN 93 WIDTH = 1.0000 ! 117050 94 ELSE 95 IF ( ITEMP.GE.1 ) WIDTH = GA(ITEMP) + (TEMP-FLOAT(ITEMP)) & 96 & *(GA(ITEMP+1)-GA(ITEMP)) 97 ENDIF 98 99 999 CONTINUE ! 18497K 100 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