XDENSWM(COL,ROW,L) = XWVAPOR(COL,ROW,L)*XDENSAM(COL,ROW,L) XCLDWTR(COL,ROW,L) = 0.0 XRANWTR(COL,ROW,L) = 0.0 ENDDO XDENSAF(COL,ROW,0)=XPRESF(COL,ROW,0)/(RDGAS*XTEMPM(COL,ROW,1)) DO L = 0,METLAY XWWIND (COL,ROW,L) = F_WW(II,JJ,L) ENDDO C ------------------------------------------------------ C FILL-UP THE landuse X-ARRAYS by populating C IF(LUTYPE .EQ. 3) THEN C DO JJL = 1,NDX DO IIL = 1,NDX DO ILU = 1,LURADM XLUSE(COL+IIL-1,ROW+JJL-1,ILU) = F_XLUSE(COL,ROW,ILU) ENDDO ENDDO ENDDO C END IF C ENDDO ENDDO C------------------------------------------------------------------- C FILL-UP THE X-ARRAYS BY EITHER INTERPOLATION OR POPULATING IF (NDX .GT. 1) THEN C FIRST, BI-LINEAR INTERPOLATION OF 2-D FIELDS CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XLATC) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XLONC) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XMAPC) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XMAPD) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XTOPO) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XPRSFC) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XTEMPG) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRAINC) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRAINN) IF (.NOT.LCALCPBL .OR. .NOT.LCALCCLD) THEN CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XEMISS) 76
ENDIF CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XGLW) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XGSW) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XHFX) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XQFX) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XPBL) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XZRUF) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCAPG) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XMOL) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRGRND) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRNET) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XFSOIL) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACH) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACM) CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACL) C SECOND, BI-LINEAR INTERPOLATION OF 3-D FIELDS CALL BILIN3D (NDX,NCOLS_X,NROWS_X,0,LURADM, XUSTAR) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XUU) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XVV) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,0,METLAY, XWWIND) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XPRESM) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XTEMPM) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XWVAPOR) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XCLDWTR) CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XRANWTR) ENDIF C----------------------------------------------------------------------- C C CALCULATE DERIVED QUANTITIES SUCH AS LAYER HEIGHTS, DENSITY, etc. c write(*,*) ' calling DERIVED ' c CALL DERIVED DO R = 1, NROWS_X DO C = 1, NCOLS_X XDENSS(C, R) = XDENSAF(C, R, 0) XDLUSE(C, R) = 1. DO ILU = 2, LURADM IF(XLUSE(C, R,ILU).GT.XLUSE(C, R,ILU-1)) XDLUSE(C, R) = ILU ENDDO ENDDO ENDDO C---------------------------------------------------------------------- 77
- 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: c JJ = J0 + JW - 1 JJ = J0 + JW COL
- 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 and 108: C C********************************
- Page 109 and 110: C REAL VSAT, VPRESS ! saturation an
- 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