1 2 SUBROUTINE NTRPWS(HTBELO,VBELOW,HTABOV,VABOVE,REQDHT,VALUE) 3 !======================================================================= 4 ! NTRPWS Module of the AERMOD Dispersion Model 5 ! 6 ! Purpose: To compute the wind speed at an intermediate level by 7 ! interpolating between two observed values. 8 ! 9 ! Input: Profile height below the level at which wind speed 10 ! is required (HTBELO) 11 ! Wind speed at HTBELO 12 ! Profile height above the level at which wind speed 13 ! is required (HTABOV) 14 ! Wind speed at HTABOV 15 ! Height at which data are required (REQDHT) 16 ! 17 ! Output: Wind speed at the required level (VALUE) 18 ! 19 ! Called by: 20 ! 21 ! Assumptions: No missing data 22 ! 23 ! 24 ! Programmer: Jim Paumier 30 September 1993 25 ! Pacific Environmental Services 26 ! 27 ! Revision history: 28 ! JOP March 14, 1995 - modified the check on the ratio for a 29 ! positive slope and observed values 30 ! 31 ! Reference(s): Revision dated 6/2/93 to the Interface Model Coding 32 ! Abstract 33 ! 34 !----------------------------------------------------------------------- 35 ! 36 !---- Variable declarations 37 ! 38 IMPLICIT NONE 39 40 REAL :: REFABV , REFBLW , RATIO , VALUE , VABOVE , HTABOV , & 41 & VBELOW , REQDHT , HTBELO , REFREQ , VALINT , REFINT 42 ! 43 !---- Data dictionary 44 ! 45 ! 46 !---- Data initializations 47 ! 48 ! 49 !....................................................................... 50 !---- The computation requires 3 estimates from the reference/theoretical 51 ! profile: one height above, one height below and from the level at 52 ! which the parameter is needed. The ratio of the differences 53 ! [EST(requested ht) - EST(ht below)] / [EST(ht above) - EST(ht below)] 54 ! is applied to the difference between the observed values to obtain 55 ! interpolated value. 56 ! 57 ! COMPUTE the value of the parameter from the reference profile 58 ! at the height below the requested height (REFBLW) --- CALL REFWS 59 ! 60 CALL REFWS(HTBELO,REFBLW) ! 0 61 62 ! COMPUTE the value of the parameter from the reference profile 63 ! at the height above the requested height (REFABV) --- CALL REFWS 64 ! 65 CALL REFWS(HTABOV,REFABV) 66 ! COMPUTE the value of the parameter from the reference profile 67 ! at the requested height (REFREQ) --- CALL REFWS 68 ! 69 CALL REFWS(REQDHT,REFREQ) 70 71 ! 72 ! Linearly interpolate to REQDHT from observed and reference profiles 73 CALL GINTRP(HTBELO,VBELOW,HTABOV,VABOVE,REQDHT,VALINT) 74 CALL GINTRP(HTBELO,REFBLW,HTABOV,REFABV,REQDHT,REFINT) 75 ! REFREQ is value from REFerence profile at REQuired height 76 ! REFINT is value from REFerence profile linearly INTerpolated to req ht 77 ! VALINT is the observed VALue linearly INTerpolated to required height 78 79 RATIO = REFREQ/REFINT 80 VALUE = RATIO*VALINT 81 82 CONTINUE 83 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