1 ! 2 ! 3 REAL FUNCTION STBLU(U,USTR,Z,OBLEN,Z0,VALK) 4 !======================================================================= 5 ! STBLU Module of the AERMOD Dispersion Model 6 ! 7 ! Purpose: To compute the wind speed at the specified height 8 ! for the stable atmosphere 9 ! 10 ! Input: Wind speed at the reference height (U) 11 ! Reference height (REFHT) 12 ! Height at which data are needed (Z) 13 ! Monin-Obukhov length (OBLEN) as a stability parameter 14 ! Surface roughness length (Z0) 15 ! 16 ! Output: Wind speed at the required level (STBLU) 17 ! 18 ! Assumptions: All the input data exist (i.e., are nonmissing) 19 ! 20 ! 21 ! Programmer: Jim Paumier 30 September 1993 22 ! Pacific Environmental Services 23 ! 24 ! Revision history: 25 ! 12/7/94 - Added the PSI(Z0/L) term to the theoretical wind speed 26 ! 27 ! Reference(s): August 6, 1992 Model Coding Abstract and addenda from 28 ! the AERMIC workgroup 29 ! 30 !----------------------------------------------------------------------- 31 ! 32 !---- Variable declarations 33 ! 34 IMPLICIT NONE 35 36 REAL :: U , USTR , Z , OBLEN , Z0 , PSIM , PSIM1 , VALK , PSIMZ0 37 ! 38 ! 39 !---- Data dictionary 40 ! 41 ! VALK = 0.4 (von Karman constant) 42 ! PSIM Correction for stability at Z 43 ! PSIM1 Correction for stability at REFHT 44 ! 45 !---- Data initializations 46 ! 47 !....................................................................... 48 ! 49 PSIM = -17.0*(1.0-EXP(-0.29*Z/OBLEN)) ! 443004 50 51 !jop Compute PSI(Z0/L) 52 PSIMZ0 = -17.0*(1.0-EXP(-0.29*Z0/OBLEN)) 53 54 STBLU = (USTR/VALK)*(ALOG(Z/Z0)-PSIM+PSIMZ0) 55 ! 56 CONTINUE 57 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