Alternative small scale meteorology input to a chemical transport ...

More documents

Recommendations

Info

C 3: winter dec, jan, feb C 5: spring mar, apr, may ! season index 4 not used in PBLPKG C 1: summer jun, jul, aug C 2: autumn sep, oct, nov C C** local indices INTEGER ISESN ! season # INTEGER IP1, JP1, L ! loop indices INTEGER C, R, K INTEGER ILU, IZ, IW C C ------------------------------------------------------------------ C For "MAXK" layer meteorology C new method for the estimation of surface fluxes C a la Byun (Nov 1990, Jan 1997) C-- NDEP: number of mm5 layers averaged into the lowest ctm layer C-- to compute deposition velocities. C C** local arrays C 2-D/3-D REAL DTMPV (NCOLS_X,NROWS_X) ! virt pot. temp. difference REAL FCORIS (NCOLS_X,NROWS_X) ! Coriolis factor REAL TSURF (NCOLS_X,NROWS_X) ! surface temperature used REAL QSURF (NCOLS_X,NROWS_X) ! mixing ratio at surface REAL QRATIO (NCOLS_X,NROWS_X) ! mixing ratio at layer 1 c REAL RADIAT (NCOLS_X,NROWS_X) ! radiation REAL WSPEED (NCOLS_X,NROWS_X,MAXK) ! wind speed REAL UUSTC (NCOLS_X,NROWS_X) ! SPEED * USTAR = constant REAL VPTEMP (NCOLS_X,NROWS_X) ! virt. pot. temperature C REAL PPSAVE (NCOLS_X,NROWS_X) ! array for saving PBL power REAL RRSAVE (NCOLS_X,NROWS_X) ! array for saving PBL power REAL Z0B1 (NCOLS_X,NROWS_X) ! temp. array for roughness REAL Z0 (NCOLS_X,NROWS_X,LURADM) ! landuse dependent roughness REAL BRATIO (NCOLS_X,NROWS_X) ! Estimate of Bowen Ratio C 1-D C INTEGER ILANDU (MAXI,MAXJ) ! landuse from calmet REAL WSTR (MAXI,MAXJ) ! wstar from calmet COMMON ILANDU ! directly from calmet COMMON WSTR ! directly from calmet REAL RIBK ( MAXK ) REAL VPTEMPK( MAXK ) REAL ZL(MAXK), ZR(MAXK+1), UL(MAXK), DTHDZ(MAXK) REAL ZTOP(MAXK), TL(MAXK) REAL ZF(0:MAXK), QL(MAXK), PL(MAXK) C C** local scalar variables REAL UNS, VNS ! wind component 100
C REAL VSAT, VPRESS ! saturation and current vepor pressures REAL DENSW ! density of water vapor REAL DENSD ! density of dry air REAL PSURF ! surface pressure REAL VPTSURF ! surface virtual potential temperature REAL RHUMID ! relative humidity REAL THREF ! reference potential temperature REAL QVREF ! reference vapor mixing ratio REAL ZF1 ! height of top of layer 1 at a point REAL ZMEAN1 ! height of first layer middle at a cell REAL ZRUF ! surface roughness length at a cell REAL HFXMAX ! maximum heat flux value (used for QA) REAL UST ! friction velocity at a cell REAL RIB ! bulk Richardson number REAL THSTAR ! temperature scale at a cell REAL QVSTAR ! vapor mixing ratio scale at a cell REAL ARGCH ! an argument scratch variable REAL ZTEMP ! height scratch variable REAL DTEMP ! temperature difference at a point REAL PTEMP ! atmospheric pressure at a point REAL RADIAT ! radiation REAL THETA1, THETA2 ! scratch potential temperatures REAL U1, U2 ! scratch wind speeds REAL ALATDEG ! latitude in degree REAL ANG1 ! an angle scratch variable REAL ALBEDO ! rough estimation of earth's albedo REAL ARG ! tmporary variable for a function argument REAL ZENITH ! zenith angle in degree REAL COSZEN ! cos(zenith) REAL DRIB REAL CPAIR ! REAL FINT ! scratch variable for interpolation REAL ZSL ! surface layer height REAL HPBL ! PBL height REAL HTMAX ! max height of each layer REAL TSFC, PSFC ! scratch for temp and press at surface REAL UHT, THT, QHT ! scratch for wind, temp, m. ratio at ZTEMP INTEGER IFLAG ! scratch flag REAL ZU1, ZU2, Z1 C integration limits for aerodynamic resistance computation C REAL SMALL, AIRDENS DATA & SMALL / 1.E-04 /, & AIRDENS / 1.2250 / ! kg/m**3 101
Page 1 and 2:
Alternative Small Scale Meteorology
Page 3 and 4:
TABLE OF CONTENTS LIST OF TABLES...
Page 5 and 6:
LIST OF TABLES 1. Transformation ch
Page 7 and 8:
11. Inverse Monin-Obukhov length on
Page 9 and 10:
CHAPTER 1. INTRODUCTION In recent y
Page 11 and 12:
techniques (Otte, 1999). It must be
Page 13 and 14:
2.1. Historical remarks CHAPTER 2.
Page 15 and 16:
treat independently. Proper modelin
Page 17 and 18:
2.4. Use of a meteorological model
Page 19 and 20:
spectrum of spatial and temporal sc
Page 21:
CHAPTER 3. APPROACH 3.1. Alternativ
Page 41:
educes eventual inconsistencies in
Page 71 and 72: 11. Scire, J. S., F. R. Robe, M. E.
Page 73 and 74: C** Set up includes for grid & univ
Page 75 and 76: SDATE = YYYY*1000 + DDD ! YYYYDDD S
Page 77 and 78: IF ( .NOT. COLFLAG ) GOTO 1010 IF (
Page 79 and 80: JJ = J0 + JW ! coarse y-grid index
Page 81 and 82: C ENDDO ENDDO X3JACOBF(C,R,0) = 1.0
Page 83 and 84: c JJ = J0 + JW - 1 JJ = J0 + JW COL
Page 85 and 86: ENDIF CALL BILIN2D (NDX,NCOLS_X,NRO
Page 87 and 88: C SUBROUTINE READCM C**************
Page 89 and 90: C NEW RECORD -- #3 - ADDITIONAL RUN
Page 91 and 92: ! 0=NO TURBULENT KINETIC ENERGY ! F
Page 93 and 94: C C --- READ THE 2-D METEOROLOGICAL
Page 95 and 96: C----------------------------------
Page 97 and 98: C NDATHR - INTEGER - DATE AND TIME
Page 99 and 100: C ---------------------------------
Page 101 and 102: X_CROSS(I,J) = X0 + FLOAT(I)*RESOLN
Page 103 and 104: C Compute actual specific humidity
Page 105 and 106: C**********************************
Page 107: C C********************************
Page 111 and 112: RAUS( A1, B1, CKUST ) = PRO * LOG(
Page 113 and 114: C** compute tstar and bulk Richards
Page 115 and 116: ENDDO DO C = 1, NCOLS_X DO R = 1, N
Page 117 and 118: C 640 CONTINUE C C** compute cell a
Page 119 and 120: APPENDIX B: INPUT FILE TO CALMET 4
Page 121 and 122: Starting date: Year (IBYR) -- No de
Page 123 and 124: Type of unformatted output file: (I
Page 125 and 126: !END! (IPR2) (0=no, 1=yes) Default:
Page 127 and 128: (0 = NO, 1 = YES) Layer-dependent b
Page 129 and 130: Multiplicative scaling factor for e
Page 131 and 132: !END! ! XG1 = 0. ! X Grid line 2 de
Page 133 and 134: (will use mixing ht MNMDAV,HAFANG s
Page 135: EDUCATION VITA Krassimira Ilieva La
show all

Alternative small scale meteorology input to a chemical transport ...

Create successful ePaper yourself

Delete template?

Save as template?