ENDIF ! of first C----------------------------------------------------------------------- C start <strong>to</strong> process <strong>input</strong> data arrays below C----------------------------------------------------------------------- C Initialize CALMET 2D-arrays C DO J = 1,MAXJ ! for initialization, it is OK DO I = 1,MAXI F_PSFC (I,J) = 0.0 F_TSFC (I,J) = 0.0 F_RMM (I,J) = 0.0 F_ELEV (I,J) = 0.0 F_ZI (I,J) = 0.0 F_EL (I,J) = 0.0 F_USTAR (I,J) = 0.0 F_WSTAR (I,J) = 0.0 F_DENSS (I,J) = 0.0 F_Z0 (I,J) = 0.0 F_LANDUSE(I,J) = 0.0 ENDDO ENDDO DO I = 1,MAXI DO J = 1,MAXJ DO K = 1,MAXK F_UU (I,J,K) = 0.0 F_VV (I,J,K) = 0.0 F_TT (I,J,K) = 0.0 F_QQ (I,J,K) = 0.0 F_PP (I,J,K) = 0.0 ENDDO DO K = 0,MAXK F_WW (I,J,K) = 0.0 ENDDO ENDDO ENDDO C-------------------------------------------------------------- C --- READ IN U, V, W WIND COMPONENTS FROM CALMET.DAT DO 10 IZ=1,NZ CALL RDR2D(M1UNIT,U(1,1,IZ),XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) CALL RDR2D(M1UNIT,V(1,1,IZ),XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) IF(LCALGRD)CALL RDR2D(M1UNIT,W(1,1,IZ),XBUF,MAXI,MAXJ,NX,NY, & CLABEL,NDATHR) 10 CONTINUE C C --- READ THE 3-D TEMPERATURE FIELDS FROM CALMET.DAT IF(LCALGRD.AND.IRTYPE.EQ.1)THEN DO 12 IZ=1,NZ CALL RDR2D(M1UNIT,ZTEMP(1,1,IZ),XBUF,MAXI,MAXJ,NX,NY, 1 CLABEL,NDATHR) 12 CONTINUE ENDIF 84
C C --- READ THE 2-D METEOROLOGICAL FIELDS FROM CALMET.DAT C - PGT STABILITY CLASS, C - FRICTION VELOCITY (M/S), C - MIXING HEIGHT (M), C - MONIN-OBUKHOV LENGTH (M), C - CONVECTIVE VELOCITY SCALE (M/S), C - PRECIPITATION RATE (MM/HR) C --- (IF RUN TYPE NOT=0,I.E.,ONLY WINDS COMPUTED & STORED ON DISK) IF(IRTYPE.EQ.1)THEN CALL RDI2D(M1UNIT,IPGT,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) CALL RDR2D(M1UNIT,USTR,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) CALL RDR2D(M1UNIT,ZI,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) CALL RDR2D(M1UNIT,EL,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) CALL RDR2D(M1UNIT,WSTR,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR) IF(NPSTA.GT.0)CALL RDR2D(M1UNIT,RMM,XBUF,MAXI,MAXJ,NX,NY,CLABEL & ,NDATHR) C C --- READ THE 1-D METEOROLOGICAL FIELDS FROM CALMET.DAT C - AIR TEMPERATURE (DEG. K), C - AIR DENSITY (KG/M**3), C - SHORT-WAVE SOLAR RADIATION (W/M**2), C - RELATIVE HUMIDITY (PERCENT), C - PRECIPITATION CODE CALL RDR1D(M1UNIT,TEMPK,NSSTA,CLABEL,NDATHR) CALL RDR1D(M1UNIT,RHO,NSSTA,CLABEL,NDATHR) CALL RDR1D(M1UNIT,QSW,NSSTA,CLABEL,NDATHR) CALL RDI1D(M1UNIT,IRH,NSSTA,CLABEL,NDATHR) IF(NPSTA.GT.0)CALL RDI1D(M1UNIT,IPCODE,NSSTA,CLABEL,NDATHR) ENDIF C Compute 2d-pressure, 2d-density, 3d-specific humidity CALL INTCAL(nssta,xssta,yssta,rho,irh,tempk, & elev,zface,ztemp,RHO_GRD,PSFC,PPRES,QQ) C---------------------------------------------------------------------- C FILL UP THE ARRAYS FOR WHICH WE HAVE INFORMATION. C NOTE: IN ORDER TO HAVE CALMET VARIABLES AGREE WITH THE WAY MCIP C HANDLES I,J,K, WE NEED TO FLIP THEM !!!!!!!!! THIS IS DONE C FOR ZTEMP,U,V,W,Z0,ELEV,ZI,EL,USTR,PSFC,RHO_GRD,PPRES C ASSIGN 2D AND 3D PREDICTED FIELDS DO J = 1, METCOL DO I = 1, METROW DO K = 1, METLAY F_QQ (I,J,K) = QQ(J,I,K) ! from INTCAL subroutine F_PP (I,J,K) = PPRES(J,I,K) ! from INTCAL subroutine F_TT (I,J,K) = ZTEMP(J,I,K) ! from calmet.dat ENDDO 85
- 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: ! 0=NO TURBULENT KINETIC ENERGY ! F
- 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
Change language