1 2 SUBROUTINE REFSVC(HEIGHT,VALUE) 3 !======================================================================= 4 ! REFSVC Module of the AERMOD Dispersion Model 5 ! 6 ! Purpose: To compute the value of sigma-V, the horizontal 7 ! (lateral) dispersion parameter from the surface 8 ! friction velocity and convective scaling velocity 9 ! alone (i.e., no observations of sigma-V) 10 ! 11 ! Input: Mixing height (ZI) 12 ! Friction velocity (USTAR) 13 ! Convective scaling velocity (WSTAR) 14 ! Height at which sigma-V is needed (HEIGHT) 15 ! 16 ! Output: sigmaV at the required level (VALUE) 17 ! 18 ! Assumptions: Mixing height, friction and convective scaling 19 ! velocities are nonmissing 20 ! 21 ! Programmer: Jim Paumier 30 September 1993 22 ! Pacific Environmental Services 23 ! 24 ! Revision history: 25 ! <none> 26 ! 27 ! Reference(s): 28 ! 29 !----------------------------------------------------------------------- 30 ! 31 !---- Variable declarations 32 ! 33 USE MAIN1 34 IMPLICIT NONE 35 CHARACTER MODNAM*12 36 37 SAVE 38 REAL HEIGHT , VALUE , SV2 , ZDCRS , ATZI , SV2DCR , VAL2 39 ! 40 !---- Data dictionary 41 ! SV2 = variance of the horizontal component of the wind 42 ! ZDCRS = value of 1.2 * ZICONV used in profiling SV2 43 ! ATZI = value of SV at ZICONV 44 ! AT1PT2 = the fraction of and above the mixing height 45 ! through which the parameter is changing 46 ! 47 ! AT1PT2 is defined in a PARAMETER statement in MAIN1.INC 48 ! 49 !---- Data initializations 50 ! 51 MODNAM = 'REFSVC' ! 43326 52 ZDCRS = AT1PT2*ZICONV 53 ! 54 !....................................................................... 55 ! 56 SV2 = 0.35*WSTAR**2 57 ! 58 IF ( HEIGHT.LE.ZICONV ) THEN 59 VALUE = SQRT(SV2) ! 9740 60 61 ELSEIF ( HEIGHT.GT.ZICONV .AND. HEIGHT.LE.ZDCRS ) THEN 62 ! COMPUTE sigmaV at 1.2*ZI 63 SV2DCR = AMIN1(SV2,0.25) ! 748 64 ! INTERPOLATE between value of SV2 at ZI and at 1.2*ZI 65 CALL GINTRP(ZICONV,SV2,ZDCRS,SV2DCR,HEIGHT,VAL2) 66 VALUE = SQRT(VAL2) 67 68 ELSE ! requested height is above 1.2*mixing height 69 ATZI = SQRT(SV2) ! 32838 70 VALUE = AMIN1(ATZI,0.5) 71 72 ENDIF 73 ! 74 CONTINUE ! 43326 75 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