1 2 3 SUBROUTINE GRDPT() 4 !======================================================================= 5 ! GRDPT module of the AERMOD Dispersion Model 6 ! 7 ! Purpose: To construct a profile of gridded values of 8 ! potential temperature 9 ! 10 ! Input: Profile of gridded potential temperature gradients 11 ! Temperature at the reference height 12 ! Profile of grid heights 13 ! 14 ! Output: Potential temperature profile at the grid heights. 15 ! 16 ! Assumptions: There is at least one grid level below the reference 17 ! temperature height (which should be satisfied 18 ! because the lowest grid level is 0.5 meters) 19 ! 20 ! Called by: METEXT 21 ! 22 ! Programmer: Jim Paumier 30 Sept 1993 23 ! Pacific Environmental Services 24 ! 25 ! Revision history: 26 ! 12/10/97 - R. Brode, PES, Inc. 27 ! Corrected the order of array indices used for profiling 28 ! potential temperature above the reference height. 29 ! 12/16/94 - J. Paumier, PES, Inc. 30 ! - CALL LOCATE to get the number of levels below the 31 ! temperature reference height, replacing the original 32 ! method which relied on grid heights being every 10 m 33 ! 34 !----------------------------------------------------------------------- 35 ! 36 !---- Variable declarations 37 38 USE MAIN1 39 IMPLICIT NONE 40 CHARACTER MODNAM*12 41 INTEGER :: L , NBELOW 42 REAL :: PTREF 43 44 SAVE 45 46 !---- Data definitions 47 ! 48 ! 49 !---- Data initializations 50 ! 51 ! 52 !....................................................................... 53 ! 54 55 !---- Determine the grid level below the temperature reference 56 ! height (as defined in the scalar file) ---- CALL LOCATE 57 58 CALL LOCATE(GRIDHT,1,MXGLVL,TREFHT,NBELOW) ! 3044 59 60 !---- Compute the potential temperature at the reference level 61 ! using the reference temperature (TA), the reference 62 ! temperature height (TREFHT), and the average stack base 63 ! elevation of all the emission sources (ZBASE) 64 65 PTREF = TA + GOVRCP*(TREFHT+ZBASE) 66 67 !---- Compute the potential temperature at the grid level below 68 ! the temperature reference height 69 70 GRIDPT(NBELOW) = PTREF - 0.5*(GRIDTG(NBELOW+1)+GRIDTG(NBELOW)) & 71 & *(TREFHT-GRIDHT(NBELOW)) 72 73 74 !---- Compute Potential Temp Values for Grid Levels Below Reference Ht. 75 DO L = NBELOW - 1 , 1 , -1 76 77 GRIDPT(L) = GRIDPT(L+1) - 0.5*(GRIDTG(L+1)+GRIDTG(L)) & 78 & *(GRIDHT(L+1)-GRIDHT(L)) 79 80 ENDDO 81 82 83 !---- Compute Potential Temp Values for Grid Levels Above Reference Ht. 84 DO L = NBELOW + 1 , MXGLVL ! 3044 85 86 GRIDPT(L) = GRIDPT(L-1) + 0.5*(GRIDTG(L)+GRIDTG(L-1)) & 87 & *(GRIDHT(L)-GRIDHT(L-1)) 88 89 ENDDO 90 91 CONTINUE ! 3044 92 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