1 2 SUBROUTINE REFSWC(HEIGHT,VALUE) 3 !======================================================================= 4 ! REFSWC 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 convective 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 , SW2 43 ! 44 !---- Data dictionary 45 ! 46 ! SW2 = variance of the vertical velocity 47 ! 48 !---- Data initializations 49 ! 50 MODNAM = 'REFSWC' ! 43326 51 ! 52 !....................................................................... 53 ! 54 IF ( HEIGHT.LE.0.1*ZICONV ) THEN 55 SW2 = 1.6*(HEIGHT/ZICONV)**0.667*WSTAR**2 ! 3836 56 VALUE = SQRT(SW2) 57 58 ELSEIF ( HEIGHT.GT.0.1*ZICONV .AND. HEIGHT.LE.ZICONV ) THEN 59 VALUE = SQRT(0.35*WSTAR**2) ! 5904 60 61 ELSEIF ( HEIGHT.GT.ZICONV ) THEN 62 SW2 = 0.35*WSTAR**2*EXP(-(6.*(HEIGHT-ZICONV)/ZICONV)) ! 33586 63 VALUE = SQRT(SW2) 64 65 ENDIF 66 ! 67 68 CONTINUE ! 43326 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