1 2 SUBROUTINE REFSWM(HEIGHT,VALUE) 3 !======================================================================= 4 ! REFSWM Module of the AERMOD Dispersion Model 5 ! 6 ! Purpose: To compute the value of sigma-W, the vertical 7 ! dispersion parameter from the surface friction 8 ! velocity and convective scaling velocity alone 9 ! for the mechanical boundary layer. 10 ! (i.e., no observations of sigma-W) 11 ! 12 ! Input: 13 ! 14 ! 15 ! 16 ! 17 ! Output: Value of the parameter at the required level 18 ! 19 ! Assumptions: Mixing height, friction and convective scaling 20 ! velocities are nonmissing 21 ! 22 ! Programmer: Jim Paumier 30 September 1993 23 ! Pacific Environmental Services 24 ! 25 ! Revision history: 26 ! Jim Paumier, PES 5 July 1995 27 ! - modified the coefficient on the u* term for CBL; 28 ! - made the value 1.2 a parameter in MAIN1; 29 ! - modified the form of sigma-w at all height for SBL; 30 ! 31 ! Reference(s): 32 ! 33 !----------------------------------------------------------------------- 34 ! 35 !---- Variable declarations 36 ! 37 USE MAIN1 38 IMPLICIT NONE 39 CHARACTER MODNAM*12 40 41 SAVE 42 REAL HEIGHT , VALUE , SWR , SWBL 43 ! 44 !---- Data dictionary 45 ! 46 ! SWR = variance of the vertical velocity due to residual turbulence 47 ! SWBL = variance of the vertical velocity in the boundary layer 48 ! 49 !---- Data initializations 50 ! 51 MODNAM = 'REFSWM' ! 264828 52 ! 53 !....................................................................... 54 ! 55 ! 56 57 SWR = SWRMAX*MIN(1.0,HEIGHT/ZI) 58 59 IF ( HEIGHT.LT.ZI ) THEN 60 SWBL = 1.3*USTAR*SQRT(1.0-HEIGHT/ZI) ! 77090 61 ELSE 62 SWBL = 0.0 ! 187738 63 ENDIF 64 65 VALUE = SQRT(SWR*SWR+SWBL*SWBL) ! 264828 66 67 68 CONTINUE 69 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