1 ! SPAG options set to convert source form only 2 ! Also redirected main output to screen to facilitate validation 3 ! (at line 4094) 4 MODULE MAIN1 5 6 !*********************************************************************** 7 ! MAIN1 8 ! AERMOD Model Data - Parameter, Variable and Array Declarations 9 ! Global Data for All Modules 10 ! 11 !*********************************************************************** 12 13 IMPLICIT NONE 14 EXTERNAL ERF 15 16 17 ! ---------------------------------------------------------------------- 18 ! --- ISC-PRIME Version 1.0 19 ! 20 ! ISCST3 modified to accept additional building data required by 21 ! the PRIME building downwash modules 22 ! 23 ! Changes are denoted in comment fields 24 ! 25 ! Prepared by Earth Tech, Inc 26 ! Prepared for EPRI under contract WO3527-01 27 ! ---------------------------------------------------------------------- 28 29 !*********************************************************************** 30 ! User Specified Model Parameters for Array Dimensions 31 !*********************************************************************** 32 33 ! --- Most array limits for data storage are now allocated at runtime. 34 INTEGER , PARAMETER :: NSEC = 36 , NWSCAT = 6 , NKST = 6 , & 35 & NHR = 24 , NPAIR = 100 , NWET = 2 , & 36 & NHIANN = 10 , NMXPM = 10 , MXPLVL = 50 , & 37 & MXGLVL = 87 38 39 !** NSRC = Max Number of Sources 40 !** NREC = Max Number of Receptors 41 !** NGRP = Max Number of Source Groups 42 !** NAVE = Max Number of Short Term Averaging Periods 43 !** NVAL = Max Number of High Values by Receptor (RECTABLE Keyword) 44 !** NTYP = Max Number of Output Types per Run (CONC, DEPOS, DDEP and WDEP) 45 !** NMAX = Max Number of Overall Maximum Values (MAXTABLE Keyword) 46 !** NSEC = Number of Sectors for Building Dimensions 47 !** NQF = Number of Variable Emission Rate Factors Per Source 48 !** NPDMAX = Max Number of Particle Diameter Categories Per Source 49 !** NWSCAT = Number of Wind Speed Categories 50 !** NKST = Number of Stability Categories 51 !** IXM = Max Number of X-coord (Distance) Values Per Receptor Network 52 !** IYM = Max Number of Y-coord (Direction) Values Per Receptor Network 53 !** NNET = Max Number of Cartesian and/or Polar Receptor Networks 54 !** NHR = Number of Hours in Met Data Loop 55 !** NPAIR = Number of Pairs of TXCONC and IDCONC for TOXXFILE Output 56 !** NWET = Number of phases of precipitation (2 = liquid,frozen) 57 !** NHIANN = Number of high period/annual averages to include in 58 !** the summary page of the output file (formerly controlled 59 !** by NVAL) 60 !** NMXPM = Number of high average H4H 24-hour averages to include 61 !** in the summary table for post-1997 PM10 processing 62 !** MXPLVL = Maximum number of levels in the observed profile file 63 !** MXGLVL = Maximum number of levels in the gridded profiles (0 - 4000 m) 64 !** NARC = Maximum number of Receptor Groupings ('ARCs') for EVALCART 65 66 67 !*********************************************************************** 68 ! Model Constants Specified as Parameters 69 !*********************************************************************** 70 71 ! & AT1PT2 = 1.2, FRACZI = 0.5, FRACSW = 0.05, 72 REAL , PARAMETER :: PI = 3.141593 , TWOPI = 6.283185 , & 73 & RTOFPI = 1.772454 , SRT2PI = 2.506628 , & 74 & G = 9.80616 , VONKAR = 0.4 , & 75 & GOVRCP = 0.00977 , RTOF2 = 1.414214 , & 76 & RTPIBY2 = 1.2533141 , DTORAD = 0.017453293 , & 77 & RTODEG = 57.29578 , DCTODK = 273.16 , & 78 & BETA1 = 0.6 , BETA2 = 0.4 , AT1PT2 = 1.2 , & 79 & STABIY = 0.04 , STABIZ = 0.016 , & 80 & UMINGR = 0.01 , GSIGV = 0.073864 , & 81 & EFOLDH = 0.44 , SVUMIN = 0.05 , SVMIN = 0.2 , & 82 & SWMIN = 0.02 , XVAL = 0.0 , TVAL = 0.01 , & 83 & FFVAL = 0.0050 , SPTGMN = 0.002 , & 84 & ALPHAF = 15.0 , ALPHA1 = 1.0 , BSUBC = 0.5 , & 85 & SZEPS = 0.01 , SZCOEF = 2.15 , HTFACT = 1.0 , & 86 & ALPHAR = 1.4 , LAMDAY = 2.3 , ASUBE = 0.1 , & 87 & REFPOP = 2000000. , DELTRUR = 12.0 , & 88 & CSUBD = 0.15 , RGAS = 8.3145 89 90 !** PI = PI 91 !** TWOPI = 2.*PI 92 !** RTOFPI = SQRT(PI) 93 !** SRT2PI = SQRT(2.*PI) 94 !** RTPIBY2= SQRT(PI/2.) 95 !** BETA1 = Coeff. in the calculation of 'direct' plume rise 96 !** BETA2 = Coeff. in the calculation of buoyancy-induced dispersion 97 !** G = Acceleration Due to Gravity (m/s**2) 98 !** DTORAD = Degrees to Radians Conversion Factor 99 !** RTODEG = Radians to Degrees Conversion Factor 100 !** AT1PT2 = The fraction of the mixed layer and above the mixed 101 ! layer through which a variable changes its value 102 !** FRACZI = The fraction of the value at ZI a parameter has at 103 !** 1.2*ZI, where ZI is the boundary layer height. 104 !** FRACSW = The fractional decrease of sigma_W at ZI 105 !** STABIY = Stability F turbulence intensity, horizontal component 106 !** STABIZ = Stability F turbulence intensity, vertical component 107 !** UMINGR = Minimum value for a gridded wind speed 108 !** GSIGV = Constant used in converting sigma_A to sigma_V 109 !** SVUMIN = Minimum value applied to Sigma_V / U when calculating 110 !** Sigma_Y. 111 !** SVMIN = Minimum value applied to measured Sigma_V values 112 !** SWMIN = Minimum value applied to measured Sigma_W values 113 !** DCTODK = Degrees Celsius to kelvin conversion factor 114 !** VONKAR = von Karman constant 115 !** GOVRCP = Gravity divided by specific heat at constant pressure 116 !** RTOF2 = Square root of 2 117 !** EFOLDH = Constant in computation of dTHETA/dZ in stable atmosphere 118 !** XVAL = Gradient in the miXing layer (unstable atmosphere) 119 !** TVAL = Gradient in the Transition layer (between ZI and 1.2*ZI) 120 !** FFVAL = Gradient in Free flow (above 1.2*ZI) 121 !** SPTGMN = Minimum vert. potential temp. gradient for stable atmosphere 122 !** ALPHAF = Constant in computation of Xm 123 !** ALPHA1 = Constant used in computing sigma_Z for surface layer releases 124 !** BSUBC = Constant used in computing sigma_Z for surface layer releases 125 !** SZEPS = Convergence criterion for convergence of sigma-z in OPTD3 126 !** SZCOEF = Coefficient of sigma-z to define the plume half-width in OPTD3 127 !** HTFACT = Height factor (no longer used) 128 !** ALPHAR = Parameter used in January 18, 1995 indirect source MCA 129 !** LAMDAY = Parameter used in January 18, 1995 indirect source MCA 130 !** ASUBE = Parameter used in January 18, 1995 indirect source MCA 131 !** REFPOP = Reference population for urban option (4/1/96 MCA) 132 !** DELTRUR= Surface cooling in the rural area (4/1/96 MCA) 133 !** RGAS = ideal gas law constant = 8.3145 Pa-m^3/mol-K 134 !** 135 136 137 !*********************************************************************** 138 ! Programmer Specified Model Parameters 139 !*********************************************************************** 140 141 INTEGER , PARAMETER :: IFMAX = 40 , IKN = 85 , ISTRG = 132 , & 142 & IERRN = 170 , ILEN_FLD = 80 143 !*# 144 145 !** IFMAX = Max Number of Fields Per Runstream Record 146 !** IPN = Number of Pathway IDs (Includes '**') 147 !** IKN = Number of Keywords 148 !** ISTRG = Length of Runstream Image Record 149 !** IERRN = Number of Error/Warning/Informational Messages 150 !** ILEN_FLD = Length of Runstream Input Fields. Also used to 151 !** specify length of input filenames and formats. 152 153 154 !*********************************************************************** 155 ! Common Block for Input/Output File Units (Initialized in BLOCK DATA) 156 !*********************************************************************** 157 158 INTEGER :: INUNIT , IOUNIT , MFUNIT , MPUNIT , IERUNT , IERWRT , & 159 & IDPUNT , IDPUN2 , IRSUNT , IEVUNT , ITEVUT , IHREMI , & 160 & IZUNIT , INCUNT , DBGUNT , DBMUNT , ICVUNT , IO3UNT 161 162 !** INUNIT = Input Runstream File Unit (Initialized to 7) 163 !** IOUNIT = Main Printed Output File Unit (Initialized to 8) 164 !** MFUNIT = Input Surface Met Data File Unit (Initialized to 19) 165 !** MPUNIT = Input Profile Met Data File Unit (Initialized to 21) 166 !** IERUNT = Temporary Error/Message File Unit (Initialized to 10) 167 !** IERWRT = Permanent Detailed Error/Message File Unit (Init. to 11) 168 !** IDPUNT = Main SAVEFILE Unit for Re-start Option (Init. to 12) 169 !** IDPUN2 = Secondary SAVEFILE Unit for Re-start Option (Init. to 14) 170 !** IRSUNT = INITFILE Unit for Re-start Option (Initialized to 15) 171 !** IEVUNT = Event File Unit for Use With ISC2EV Model (Init. to 17) 172 !** ITEVUT = Temporary Event File Used to Store High Value Events for 173 !** Summary Tables and for EVENTFIL Option (Initialized to 18) 174 !** IZUNIT = Gridded Terrain File Used (Initialized to 13) 175 !** IHREMI = Hourly Emission Parameters File Unit (Init. to 16) 176 !** INCUNT = INCLUDED File Unit (Initialized to 20) 177 !** DBGUNT = Debug Output File for Calculations (Init. to 24) 178 !** DBMUNT = Debug Output File for Meteorology Profiles (Init. to 25) 179 !** ICVUNT = Cavity Concentration Output File (Initialized to 26) 180 !** IO3UNT = Hourly Ozone Data File for PVMRM and OLM Options (Init. to 29) 181 !** The following file units are declared below: 182 !** ISUNIT = Surface Meteorology File for SCIM'd data (Initialized to 22) 183 !** IPUNIT = Profile Meteorology File for SCIM'd data (Initialized to 23) 184 185 !*********************************************************************** 186 ! This is The Global Variable Definition Block for Runstream Data 187 !*********************************************************************** 188 189 LOGICAL BLINE , INFLD , MARK , ECHO 190 191 CHARACTER PATH*2 , PPATH*2 , KEYWRD*8 , PKEYWD*8 , KEYWD*8 , & 192 & KTYPE*5 , RUNST*1 193 194 CHARACTER(LEN=ILEN_FLD) :: FIELD , INPFIL , OUTFIL , INCFIL 195 CHARACTER(LEN=ISTRG) :: RUNST1 196 197 INTEGER :: LOCB(IFMAX) , LOCE(IFMAX) , IFC , IDC1 , IPNUM , IPPNUM 198 DIMENSION FIELD(IFMAX) , KEYWD(IKN) , RUNST(ISTRG) 199 200 201 !*********************************************************************** 202 ! This is The Global Variable Definition Block for Error Handling 203 !*********************************************************************** 204 205 LOGICAL FATAL , ISTART , IFINIS , RECERR , ERRLST , EOF 206 207 CHARACTER ERRMSG*50 , ERRCOD*3 , VERSN*5 208 CHARACTER(LEN=ILEN_FLD) :: MSGFIL 209 210 DIMENSION ERRMSG(IERRN) , ERRCOD(IERRN) 211 INTEGER :: ILINE , IERROR , IFTL , IWRN , INFO , ICLM , IMSG , & 212 & IHEZ , NFATAL , NWARN , IPAGE 213 REAL :: EXPLIM 214 INTEGER :: ICSTAT(30) , ISSTAT(30) , IRSTAT(30) , IMSTAT(30) , & 215 & IOSTAT(30) , ITSTAT(30) , IESTAT(30) 216 INTEGER :: INCSET , IXYSET , IEVSET , IHLSET=0 , IFGSET 217 218 219 !*********************************************************************** 220 ! This is The Global Variable Definition Block for COntrol Pathway 221 !*********************************************************************** 222 223 LOGICAL :: DFAULT , CONC , DEPOS , DDEP , WDEP , RURAL , URBAN , & 224 & GRDRIS , NOSTD , NOBID , CLMPRO , MSGPRO , PERIOD , & 225 & ANNUAL , MONTH , FLAT , ELEV , FLGPOL , RUN , EVENTS , & 226 & RSTSAV , RSTINP , DAYTAB , MXFILE , PPFILE , PLFILE , & 227 & ANPOST , ANPLOT , STATOK , FSTREC , MULTYR , TXFILE , & 228 & RKFILE , DDPLETE , WDPLETE , FSTCMP , EVONLY , SOCONT , & 229 & DETAIL , NEWMET , SEASONHR , ARDPLETE , DEBUG , METEOR , & 230 & NOCHKD , NOWARN , SCIM , WETSCIM , SCIMHR , & 231 & WETHR , FIRSTWET , TOXICS , SCREEN , URBSTAB , & 232 & PRM_FSTREC , O3FILE 233 !JRA 30 Sept 2005 - following variables were used without being defined 234 ! this could lead to erroneous results 235 ! spotted by Salford FTN95 version 4.8.0 236 LOGICAL :: PM10AVE=.FALSE. , ROMBERG=.FALSE. , PVMRM=.FALSE. , & 237 & OLM=.FALSE. 238 !*# 239 240 CHARACTER ELTYPE*6 , TITLE1*68 , TITLE2*68 , EVPARM*6 , CHRAVE*5 ,& 241 & CHIDEP*4 , SOELEV*6 , REELEV*6 , TGELEV*6 , OUTTYP*5 242 243 CHARACTER(LEN=ILEN_FLD) :: SAVFIL , SAVFL2 , INIFIL , EVFILE , & 244 & DBGFIL , DBMFIL , URBNAM , OZONFL , & 245 & O3FILUNITS , O3VALUNITS , O3FORM 246 247 REAL :: O3CONC , O3BACK , NO2EQUIL 248 LOGICAL :: O3MISS 249 250 INTEGER :: NHIVAL , NMXVAL , NDUMP 251 INTEGER :: NSRC , NREC , NGRP , NQF, & 252 & NPDMAX , NNET , IXM , IYM , NUMEVE , IEVENT , & 253 & NARC , NOLM 254 !JRA 30 Sept 2005 - following variables were used without being defined 255 ! this could lead to erroneous results 256 ! spotted by Salford FTN95 version 4.8.0 257 INTEGER :: NAVE=0 , NVAL=0 , NTYP=0 , NMAX=0 , NEVE=0 258 259 INTEGER , ALLOCATABLE :: KAVE(:) 260 261 LOGICAL , ALLOCATABLE :: EVAL(:) 262 263 ALLOCATABLE ::CHRAVE(:) , CHIDEP(:,:) , OUTTYP(:) 264 CHARACTER(LEN=6) :: MODOPS(18) 265 266 267 268 !*********************************************************************** 269 ! This is The Global Variable Definition Block for SOurce Pathway 270 !*********************************************************************** 271 272 CHARACTER SRCID*8 , SRCTYP*8 , SOPCRD*1 , SOGAS*1 , URBSRC*1 , & 273 & GRPID*8 , EMILBL*40 , OUTLBL*40 , POLLUT*8 , PSOID*8 , & 274 & QFLAG*8 , PERLBL*40 , OLMID*8 275 276 CHARACTER(LEN=ILEN_FLD) :: HRFILE 277 !*# 278 CHARACTER PREVGRPID*8 279 280 LOGICAL LDPART , LWPART , LDGAS , LWGAS 281 282 REAL , ALLOCATABLE :: AXS(:) , AYS(:) , AZS(:) , AQS(:) , AHS(:) ,& 283 & ATS(:) , AVS(:) , ADS(:) , ASYINI(:) , & 284 & ASZINI(:) , ANO2_RATIO(:) 285 REAL , ALLOCATABLE :: ADSBH(:,:) , ADSBW(:,:) , ADSBL(:,:) , & 286 & ADSXADJ(:,:) , ADSYADJ(:,:) 287 288 INTEGER , ALLOCATABLE :: INPD(:) , NDXSTK(:) 289 REAL , ALLOCATABLE :: QFACT(:,:) 290 REAL :: EMICON , HAFLIF , DECOEF , VAIRMS , ZRDEP , VDPHOR 291 REAL , ALLOCATABLE :: EMIFAC(:) , APDIAM(:,:) , APHI(:,:) , & 292 & APDENS(:,:) , AVGRAV(:,:) , ATSTOP(:,:) 293 294 !*---- ISCSTM Modification: allow for hourly emissions - jah 11/3/94 295 REAL :: HRQS , HRTS , HRVS 296 REAL , ALLOCATABLE :: EV_HRQS(:,:) , EV_HRTS(:,:) , EV_HRVS(:,:) 297 INTEGER :: KURHRQ 298 !*---- 299 !*# 300 301 INTEGER , ALLOCATABLE :: IGROUP(:,:) , IGRP_OLM(:,:) 302 ALLOCATABLE ::SRCID(:) , SRCTYP(:) , SOPCRD(:) , SOGAS(:) , & 303 & URBSRC(:) , GRPID(:) , QFLAG(:) , EMILBL(:) , & 304 & OUTLBL(:) , PERLBL(:) , OLMID(:) 305 LOGICAL , ALLOCATABLE :: L_OLMGRP(:) 306 307 308 !** NDXSTK Index of the gridded height immediately below stack top 309 !** KURHRQ Date/Time Stamp for Hourly Emissions Data 310 !** HRFILE Hourly Emissions Data File Name 311 !*# 312 313 !*********************************************************************** 314 ! This is The Global Variable Definition Block for the New Area 315 ! Source Algorithm - 7/7/93 316 ! 317 !* Added XVERT,YVERT - Jayant Hardikar, PES, 7/20/94 318 !*********************************************************************** 319 320 INTEGER , PARAMETER :: NVMAX = 24 , NVMAX2 = NVMAX*2 321 322 !** NVMAX = Maximum Number of Vertices in a Single Area Source + 4 323 !** NVMAX2= NVMAX * 2 324 !** Note that the '+ 4' added to number of vertices for NVMAX refers 325 !** to the maximum number of sides on an AREAPOLY source that a 326 !** straight line can intersect. For complex shapes, '+ 4' may not 327 !** by large enough. R. Brode, PES, Inc. 328 329 LOGICAL LSEG 330 331 INTEGER :: IVERT , NVERT , NSEGS 332 REAL :: UVERT(NVMAX) , VVERT(NVMAX) , VNVERT(NVMAX) , WVERT(NVMAX)& 333 & , UASEGS(NVMAX) , UBSEGS(NVMAX) , XVERT(NVMAX) , & 334 & YVERT(NVMAX) 335 REAL :: SPA(NVMAX,2) 336 REAL , ALLOCATABLE :: AXINIT(:) , AYINIT(:) , AANGLE(:) , & 337 & AXVERT(:,:) , AYVERT(:,:) , RADIUS(:) , & 338 & AXCNTR(:) , AYCNTR(:) 339 INTEGER , ALLOCATABLE :: NVERTS(:) 340 341 342 343 !*********************************************************************** 344 ! This is The Global Variable Definition Block for the New OPENPIT 345 ! Source Algorithm - 7/19/94 346 !*********************************************************************** 347 348 REAL , PARAMETER :: ALPHA = 0.029 349 REAL , ALLOCATABLE :: AALPHA(:) , APDEFF(:) , AVOLUM(:) 350 REAL , ALLOCATABLE :: EFRAC(:) , QPART(:) 351 REAL :: PALPHA , THETA , PDEFF , PDREL , PITFRA , QEFF 352 REAL :: PITLEN , PITWID , PITL , EMIHGT , XEFF , YEFF 353 354 355 !* ALPHA = Proportionality Constant for determining Escape Fraction 356 !* AALPHA = Array of Angles of Long Pit Dimension w.r.t. 357 !* North for Each Source 358 !* PALPHA = Angle of Long Pit Dimension w.r.t. North for 359 !* the Current Source 360 !* THETA = Wind Direction angle w.r.t Long Axis 361 !* of the Pit 362 !* APDEFF = Array of Effective Pit Depths for Each Source 363 !* PDEFF = Effective Pit Depths for Current Source 364 !* PDREL = Relative Pit Depth 365 !* AVOLUM = Array of Volume of the OPENPIT Sources 366 !* EFRAC = Array of Escape Fractions 367 !* QPART = Array of Adjusted Emission Rates 368 !* PITFRA = Fractional Size of the Effective Pit Area 369 !* PITLEN = Length of the Pit 370 !* PITWID = Width of the Pit 371 !* PITL = Along-Wind Length of the Pit 372 !* EMIHGT = Height of Emissions Above Base of Pit 373 !* XEFF = X-dimension of Effective Pit 374 !* YEFF = Y-dimension of Effective Pit 375 376 377 378 !---------------------------------------------------------------------- 379 ! --- COMMON BLOCK /DRYGAS/ -- Dry deposition parameters CALPUFF 380 !---------------------------------------------------------------------- 381 ! 382 REAL , ALLOCATABLE :: PDIFF(:) , PDIFFW(:) , RMOLWT(:) , ALPHAS(:)& 383 & , REACT(:) , HENRY(:) , RCLI(:) , & 384 & FINEMASS(:) , SCF(:) 385 LOGICAL , ALLOCATABLE :: L_METHOD2(:) 386 INTEGER :: ISEAS_GD(12) , ILAND_GD(36) , NCLOUD 387 REAL :: RM , RCUT , QSW , XLAI , VDEPG , USERVD , ZSUBP , & 388 & DELTA_Z , FO , FSEAS2 , FSEAS5 , FRACSAT , LIQCONT , & 389 & DENOM , XNU 390 391 REAL :: WOLD , WNEW , F2 , ESTA 392 393 CHARACTER*40 REFSPE 394 395 LOGICAL LUSERVD 396 ! 397 ! REFSPE - Reference Species (Default is SO2) 398 ! 399 ! --- COMMON BLOCK /DRYGAS/ Variables: 400 ! PDIFF(NSRC) - real - Molecular diffusivity (m**2/s) 401 ! of each pollutant. 402 ! SEE NOTE #1 403 ! PDIFFW(NSRC) - real - Molecular diffusivity in water (m**2/s) 404 ! of each pollutant. 405 ! SEE NOTE #1 406 ! RMOLWT(NSRC) - real - Molecular weight of pollutant (g/mol) 407 ! of each pollutant. 408 ! ALPHAS(NSRC) - real - Solubility enhancement factor due 409 ! to the aqueous phase reactivity of 410 ! the pollutant. 411 ! REACT(NSRC) - real - Reactivity factor for each 412 ! pollutant. 413 ! RM - real - Mesophyll resistance (s/m) 414 ! SEE NOTE #2 415 ! HENRY(NSRC) - real - Henry's law constant (ratio of 416 ! gas to aqueous phase concentration 417 ! of the pollutant). 418 ! RCUT - real - Cuticle resistance (s/m). 419 ! 420 ! NOTE #1: Input units of this variable are cm**2/s. Conversion to m**2/s 421 ! is made internally in the SETUP phase. 422 ! 423 ! NOTE #2: Input units of s/cm are converted to s/m in the SETUP phase. 424 425 426 427 !*********************************************************************** 428 ! This is The Global Variable Definition Block for REceptor Pathway 429 !*********************************************************************** 430 431 LOGICAL ISTA , IEND , NEWID 432 433 CHARACTER :: NETID*8 , NETIDT*8 , PNETID*8='' , NTID*8 , NTTYP*8 , & 434 & RECTYP*2 , PXSOID*8 , PESOID*8 , ARCID*8 435 436 REAL , ALLOCATABLE :: AXR(:) , AYR(:) , AZELEV(:) , AZFLAG(:) , & 437 & AZHILL(:) 438 INTEGER , ALLOCATABLE :: IREF(:) , NDXARC(:) 439 ALLOCATABLE ::NETID(:) , RECTYP(:) , NTID(:) , NTTYP(:) , ARCID(:) 440 INTEGER :: ICOUNT , JCOUNT , IZE , IZH , IZF , IRZE , IRZH , & 441 & IRZF , IRXR , IRYR , IRHZ , IBND , IBELEV , INNET 442 REAL :: XINT , YINT 443 REAL , ALLOCATABLE :: XCOORD(:,:) , YCOORD(:,:) , XORIG(:) , & 444 & YORIG(:) 445 INTEGER , ALLOCATABLE :: NETSTA(:) , NETEND(:) , NUMXPT(:) , & 446 & NUMYPT(:) 447 448 !** AZHILL Hill Height Associated with the Receptor 449 !** HCRIT Critical dividing streamline associated with the receptor 450 451 452 !*********************************************************************** 453 ! This is The Global Variable Definition Block for MEteorology Pathway 454 !*********************************************************************** 455 456 CHARACTER SFNAME*40 , UANAME*40 , ONNAME*40 , ALAT*10 , ALON*10 457 458 CHARACTER(LEN=ILEN_FLD) :: METINP , SCIM_SFCFIL , SCIM_PROFIL , & 459 & PROINP 460 CHARACTER(LEN=105) :: METFRM , PROFRM 461 LOGICAL SCIMOUT 462 463 ! RWB/MJ - allow for SCIM option - May, 1998. 464 INTEGER :: ISDATE , IEDATE , ISYR , ISMN , ISDY , ISHR , IEYR , & 465 & IEMN , IEDY , IEHR , IPROC(366) , ISYEAR , IUYEAR , & 466 & IOYEAR , IDSURF , IDUAIR , IDSITE , ISJDAY , IEJDAY , & 467 & NDAYS , INCRST , ISTRT_CENT , ISTRT_WIND , NREGSTART , & 468 & NREGINT , IFIRSTHR , ISUNIT , IPUNIT , NSKIPTOT , & 469 & NSKIPWET , NSKIPDRY , NSWETCLM , NSDRYCLM , NSWETMSG , & 470 & NSDRYMSG , NWETHR , NWETINT , NWETSTART 471 472 REAL :: UCAT(5) , ZREF , ROTANG , UMIN , VIRTPNT_URB(NKST) , & 473 & VIRTPNT_RUR(NKST) , VP_FACT 474 475 476 !*********************************************************************** 477 ! This is The Global Variable Definition Block for METEXT 478 !*********************************************************************** 479 480 LOGICAL CLMHR , MSGHR , UNSTAB , NEUTRL , STABLE , RUNERR , & 481 & PFLERR , NEWDAY , ENDMON , METHDR , HOURLY 482 483 INTEGER :: KSTMSG 484 INTEGER :: IHOUR , IYEAR , IMONTH , IDAY , KURDAT , ISEAS , & 485 & KHOUR , KYEAR , KMONTH , KDAY , KURPFL , NTOTHRS , & 486 & IPHOUR , IPDATE , IPCODE , KST , IYR , IDAY_OF_WEEK , & 487 & IDAY_OF_WEEK7 , NPLVLS , NTGLVL , IFLAG(MXPLVL) 488 !JRA 30 Sept 2005 - following variables were used without being defined 489 ! this could lead to erroneous results 490 ! spotted by Salford FTN95 version 4.8.0 491 INTEGER :: JDAY=0 492 INTEGER :: FULLDATE 493 REAL :: SFCHF , USTAR , WSTAR , VPTGZI , ZICONV , ZIMECH , & 494 & OBULEN , SFCZ0 , BOWEN , ALBEDO , UREF , WDREF , UREFHT , & 495 & TA , TREFHT , ZI , AFV , BVF , BVPRIM , XLAT , TSIGN , & 496 & ZIRUR , ZIURB , URBWSTR , URBPOP , PRATE , PREC1 , PREC2 ,& 497 & UREF10 , URBZ0 , URBUSTR , URBOBULEN , RURUSTR , & 498 & RUROBULEN , RH , SFCP 499 500 INTEGER :: IKST(NHR) , IAPCODE(NHR) , NACLOUD(NHR) 501 REAL :: APRATE(NHR) , AQSW(NHR) , ARH(NHR) , ASFCP(NHR) 502 REAL :: ASFCHF(NHR) , AUREF(NHR) , AUREFHT(NHR) , ATA(NHR) , & 503 & ATREFHT(NHR) , AWDREF(NHR) , AUSTAR(NHR) , AWSTAR(NHR) , & 504 & AZICONV(NHR) , AZIMECH(NHR) , AOBULEN(NHR) , AVPTGZI(NHR) & 505 & , ASFCZ0(NHR) , ABOWEN(NHR) , AALBEDO(NHR) , AWNEW(NHR) , & 506 & AWOLD(NHR) , AESTA(NHR) , AF2(NHR) , APREC1(NHR) , & 507 & APREC2(NHR) 508 509 INTEGER :: IENDHOUR , IENDDY , IENDMN , NUMYRS , NREMAIN , NDX4ZI 510 511 REAL :: PFLHT(MXPLVL) , PFLWD(MXPLVL) , PFLWS(MXPLVL) , & 512 & PFLTA(MXPLVL) , PFLSA(MXPLVL) , PFLSW(MXPLVL) , & 513 & PFLSV(MXPLVL) , PFLTG(MXPLVL) , PFLTGZ(MXPLVL) 514 REAL :: APFLHT(NHR,MXPLVL) , APFLWD(NHR,MXPLVL) , & 515 & APFLWS(NHR,MXPLVL) , APFLTA(NHR,MXPLVL) , & 516 & APFLSA(NHR,MXPLVL) , APFLSW(NHR,MXPLVL) , & 517 & APFLSV(NHR,MXPLVL) , APFLTG(NHR,MXPLVL) , & 518 & APFLTGZ(NHR,MXPLVL) 519 INTEGER :: AIFLAG(NHR,MXPLVL) , ANPLVLS(NHR) , ANTGLVL(NHR) 520 !--- Add density profile for PRIME 521 !--- Add tubulence dissipation rate (epsilon) profile for PVMRM 522 REAL :: GRIDHT(MXGLVL) , GRIDWD(MXGLVL) , GRIDWS(MXGLVL) , & 523 & GRIDSW(MXGLVL) , GRIDSV(MXGLVL) , GRIDTG(MXGLVL) , & 524 & GRIDPT(MXGLVL) , GRIDRHO(MXGLVL) , GRIDEPS(MXGLVL) 525 REAL :: GRDSWU(MXGLVL) , GRDSVU(MXGLVL) , GRDTGU(MXGLVL) , & 526 & GRDSWR(MXGLVL) , GRDSVR(MXGLVL) , GRDTGR(MXGLVL) , & 527 & GRDPTU(MXGLVL) , GRDPTR(MXGLVL) 528 REAL :: HNPREV , HDPREV , USPREV , HN , HTRANS , SFCLVL , TG4PFL ,& 529 & TG4XTR , THSTAR , SVAVG , SWAVG , UAVG , SVATZI , SWATZI ,& 530 & UATZI , PTATZI , UATHE , SVATHE , SWATHE , UAVH3 , & 531 & SVAVH3 , SWAVH3 , SWRMAX 532 533 !** BVF = Brunt-Vaisala frequency 534 !** BVPRIM = 0.7*BVF 535 !** SFCHF = Surface heat flux (W/sq m) 536 !** USTAR = Surface friction velocity (m/s) 537 !** WSTAR = Convective scaling velocity (m/s) 538 !** VPTGZI = Vertical potential temperature gradient from ZI to ZI+500 539 !** (degrees/m) 540 !** ZICONV = Hourly convective mixing height estimated by AERMET (m) 541 !** ZIMECH = Hourly mechanical mixing height estimated by AERMET (m) 542 !** OBULEN = Monin-Obukhov length (m) 543 !** SFCZ0 = Surface roughness length (m) 544 !** BOWEN = Bowen ratio = sensible heat flux/latent heat flux 545 !** ALBEDO = Albedo at the earth's surface (nondimensional) 546 !** UREF = Reference height wind speed (m/s) 547 !** WDREF = Reference height wind direction (degrees from north) 548 !** UREFHT = Reference height for winds (m) (first nonmissing level 549 !** of wind speed AND direction above 7.0*SFCZ0) 550 !** TA = Ambient temperature at a reference height (kelvin) 551 !** TREFHT = Reference height for temperature (m) (first nonmissing 552 !** level of temperature) 553 !** ZI = The mixing height used by AERMOD after any manipulation 554 !** and massaging (m) 555 !** NPLVLS = Number of levels in the observed hourly profile data 556 !** NTGLVL = Number of levels of observed potential temperature gradient 557 !** IFLAG = Top of profile flag: 1 = top level, 0 = level below top 558 !** PFLHT = Profile height above local ground level (m) 559 !** PFLWD = Profile wind direction (degrees from north) 560 !** PFLWS = Profile wind speed (m/s) 561 !** PFLTA = Profile ambient temperature (kelvins) 562 !** PFLSA = Profile sigma_A (degrees) 563 !** PFLSW = Profile sigma_W (m/s) 564 !** PFLSV = Profile sigma_V (m/s), computed from sigma_A and wind speed 565 !** PFLTG = Profile of Vertical Potential Temperature Gradient (kelvin/m) 566 !** PFLTGZ = Profile of VPTG heights (midpoint of interval) (m) 567 !** GRIDHT = Gridded height (m) 568 !** GRIDWD = Gridded wind direction (degrees from north) 569 !** GRIDWS = Gridded wind speed (m/s) 570 !** GRIDSW = Gridded sigma_W (m/s) 571 !** GRIDSV = Gridded sigma_V (m/s) 572 !** GRIDTG = Gridded vertical potential temperature gradient (deg/m) 573 !** GRIDPT = Gridded potential temperature profile 574 !** GRIDRHO= Gridded density profile 575 !** GRIDEPS= Gridded tubulence dissipation rate (epsilon) profile for PVMRM 576 !** HNPREV = Previous hour's computed smoothed PBL height (m) 577 !** HDPREV = Previous hour's PBL height from AERMET (m) 578 !** USPREV = Previous hour's friction velocity (m/s) 579 !** HN = Current hour's smoothed PBL height (m) 580 !** HTRANS = PBL height at transition hour (heat flux from - to +) (m) 581 !** SFCLVL = Level at which the theoretical temperature gradient 582 !** profile is initialized 583 !** TG4PFL = Potential temperature gradient at 2.0 meters 584 !** TG4XTR = Potential temperature gradient at 100.0 meters 585 !** XLAT = Station latitude, decimal degrees 586 !** TSIGN = Sign used for turning of wind: 1.0 for northern hemis. 587 !** -1.0 for southern hemis. 588 !** NDX4ZI = Index of gridded height immediately below ZI 589 !** SVAVG = Average sigma_V from the surface to ZI (m/s) 590 !** SWAVG = Average sigma_W from the surface to ZI (m/s) 591 !** UAVG = Average wind speed from the surface to ZI (m/s) 592 !** SVATZI = sigma_V at ZI (m/s) 593 !** SWATZI = sigma_W at ZI (m/s) 594 !** UATZI = Wind speed at ZI (m) 595 !** PTATZI = Potential temperature at ZI (kelvin) 596 !** SVATHE = Average sigma_V from the surface to HS for HS > ZI (m/s) 597 !** SWATHE = Average sigma_W from the surface to HS for HS > ZI (m/s) 598 !** UATHE = Average wind speed from the surface to HS for HS > ZI (m/s) 599 !** SVAVH3 = Average sigma_V from the surface to HE3 for penetrated plume 600 !** SWAVH3 = Average sigma_W from the surface to HE3 for penetrated plume 601 !** UAVH3 = Average wind speed from the surface to HE3 for penetrated plume 602 !** SWRMAX = Residual vertical turbulence, average of measured sigma-w 603 !** above ZI or 0.02 * UATZI 604 605 606 !*********************************************************************** 607 ! This is The Global Variable Definition Block for Terrain Grid Pathway 608 !*********************************************************************** 609 610 CHARACTER(LEN=ILEN_FLD) :: TERINP 611 LOGICAL LTGRID 612 613 REAL :: TGX0 , TGY0 , GRDXLL , GRDXUR , GRDYLL , GRDYUR , XYINT 614 615 616 !*********************************************************************** 617 ! This is The Global Variable Definition Block for Calculation 618 !*********************************************************************** 619 620 LOGICAL CALCS , WAKE , WAKESS , BUOYNT , TALL , SQUAT , SSQUAT 621 LOGICAL SURFAC 622 623 DOUBLE PRECISION PHID1 , PHID2 , PHIN1 , PHIN2 624 625 INTEGER :: IREC , ISRC , IGRP , IAVE , ITYP , ISET , NUMREC , & 626 & NUMSRC , NUMGRP , NUMAVE , NUMARC , NUMTYP , NUMYR , & 627 & ICYEAR , NUMURB , NPD , IFVSEC , IUCAT , IOLM , & 628 !JRA 30 Sept 2005 - following variables were used without being defined 629 ! this could lead to erroneous results 630 ! spotted by Salford FTN95 version 4.8.0 631 & NUMOLM=0 632 ! --- PRIME Modification ------------------------------------------- 633 ! ------------------------------------------------------------------ 634 REAL :: XS , YS , ZS , QS , HS , DS , VS , TS , SYINIT , SZINIT , & 635 & XINIT , YINIT , ANGLE , XCNTR , YCNTR , DSBH , DSBW , & 636 & DSBL , XADJ , YADJ , B_SUBS , B_SUBL , RSCALE , D , VD , & 637 & E , WDRAD , WDSIN , WDCOS , ZBASE 638 639 !DEP REAL, ALLOCATABLE :: V(:) 640 !DEPC Accomodate Wet SCIM'ing, M. Jindal, PES, 6/10/98 641 !DEP & ,VDRY(:) 642 REAL , ALLOCATABLE :: PDIAM(:) , PHI(:) , PDENS(:) , VGRAV(:) , & 643 & TSTOP(:) , SCHMIDT(:) , VDEP(:) , WQCOR(:) ,& 644 & DQCOR(:) , PSCVRT(:) , WASHOUT(:) 645 REAL :: WQCORG , GSCVRT , DQCORG , WASHOUTG , VSETL 646 REAL :: XR , YR , X , Y , ZELEV , ZFLAG , ZR , ZEFF , DISTR , & 647 & ZHILL , HCRIT , ZRT , XDIST 648 REAL :: HE , HSP , HEFLAT , HTER , HEMWAK , HEDHH , ZB , ZM , & 649 & HED1 , HED2 , HEN1 , HEN2 , HE3 , HPEN , HED1M , HED2M , & 650 & HEN1M , HEN2M , HE3M , HSBL , QSUBN , QSUB3 , XY , XZ , & 651 & SBID , FM , FB , DTDZ , DHF , DHFAER , DHP , DHP1 , DHP2 ,& 652 & DHP3 , DELT , DHPB , DHPM , XF , XMAX , XFM , XFB , XRAD ,& 653 & WPB , DHCRIT , HTEFF , CENTER , Z4GAMMA , XTR4GAMMA 654 REAL :: HESETL , HE3SETL , HV 655 REAL :: US , SVS , SWS , TGS , TYS , PTS , UP , WDIR , DA , ZLY , & 656 & ZLB , RINIT , CB , CM , QTK , PPF , PSUBS , FHC , SY , & 657 & SYB , SYN , SY3 , SZ , SZUPR , SYAMB , SZAMB , SZAS , & 658 & SZAD1 , SZAD2 , SZAN1 , SZAN2 , SYAN , SZA3 , SZB , SZBD ,& 659 & SZBN , SZ3 , SZD1 , SZD2 , SZN1 , SZN2 , SZEFF , SZSURF , & 660 & SYA3 , SYB3 , SZB3 , VSY3 , VSIGY , VSIGZ , VSYN , VSZD1 ,& 661 & VSZD2 , VSZN1 , VSZN2 , VSZ3 , SZD1M , SZD2M , SZN1M , & 662 & SZN2M , SZ3M , U3 , SV3 , SW3 , TGP 663 DOUBLE PRECISION :: FSUBY , FSUBYD , FSUBYN , FSUBY3 664 REAL :: FSUBZ , FSUBZD , FSUBZN , FSUBZ3 , PHEE , FOPT , CWRAP , & 665 & CLIFT , XMDBG , CWRAPC , CLIFTC , FSUBYC , FSBY3C 666 REAL :: UEFF , SVEFF , SWEFF , TGEFF , UEFFD , SVEFFD , SWEFFD , & 667 & UEFFN , SVEFFN , SWEFFN , UEFF3 , SVEFF3 , SWEFF3 , & 668 & TGEFF3 , EPSEFF , EPSEFFD , EPSEFF3 , XMIXED , XFINAL , & 669 & ZMIDMX 670 REAL :: SKEW , R , ALPHPD , BETAPD , ASUB1 , ASUB2 , BSUB1 , & 671 & BSUB2 , LAMDA1 , LAMDA2 672 REAL :: CHIW , CHIDW , CHINW , CHI3W , CHIL , CHIDL , CHINL , & 673 & CHI3L 674 REAL :: GAMFACT 675 676 !** AZSAVG = Average stack base elevation (m) 677 !** ZBASE = Base elevation used for potential temperature profile (m MSL) 678 !** US = Wind speed at stack height (m/s) 679 !** UP = Stack top wind speed for plume rise computations 680 !** WDIR = Stack top wind direction used for plume transport 681 !** SVS = sigma_V at stack height (m/s) 682 !** SWS = sigma_W at stack height (m/s) 683 !** TGS = Potential temperature gradient at stack top 684 !** PTS = Stack top potential temperature for plume rise 685 !** xxEFF = "effective" value for parameter xx 686 !** FHC = Function of Plume material above HCRIT 687 !** PHEE = "PHI" Term : Fraction of Plume Below Hcrit 688 !** FSUBY = Fy Term (Horizontal Gaussian term) 689 !** FSUBYN = Fy Term (Horizontal Gaussian term) for the 690 !** Indirect Source 691 !** FSUBY3 = Fy Term (Horizontal Gaussian term) for the 692 !** Penetrated Source 693 !** HEDx = Effective Source Heights for Direct Plume, 694 !** x corresponding to each of the 2 distributions 695 !** HENx = Effective Source Heights for Indirect Plume, 696 !** x corresponding to each of the 2 distributions 697 !** HE3 = Effective Source Height for Penetrated Plume 698 !** HEDxM = Effective Source Heights for Direct Plume at Xm, 699 !** x corresponding to each of the 2 distributions 700 !** HENxM = Effective Source Heights for Indirect Plume at Xm, 701 !** x corresponding to each of the 2 distributions 702 !** HE3M = Effective Source Height for Penetrated Plume at Xm 703 704 !** QSUBN = Source Term for Indirect Source 705 !** QSUB3 = Source Term for Penetrated Source 706 !** SKEW = Skewness of the Vertical Velocity 707 !** R = Lagrangian Correlation Coefficient 708 !** ALPHPD = ALPHA Coefficient for the CBL PDF 709 !** BETAPD = BETA Coefficient for the CBL PDF 710 !** ASUB1 = 711 !** ASUB2 = 712 !** BSUB1 = 713 !** BSUB2 = 714 !** LAMDA1 = Relative Frequencies of Updrafts 715 !** LAMDA2 = Relative Frequencies of Downdrafts 716 717 718 !*********************************************************************** 719 ! This is The Global Variable Definition Block for EVent Pathway 720 !*********************************************************************** 721 722 CHARACTER EVNAME*8 , EVGRP*8 723 INTEGER , ALLOCATABLE :: EVAPER(:) , EVDATE(:) , EVJDAY(:) , & 724 & IDXEV(:) 725 726 ALLOCATABLE ::EVNAME(:) , EVGRP(:) 727 728 729 730 731 !*********************************************************************** 732 ! This is The Global Variable Definition Block for OUtput Pathway 733 !*********************************************************************** 734 735 LOGICAL OUTPART 736 737 LOGICAL , ALLOCATABLE :: ANPART(:) , ALLPARTS(:) , ALLPARTG(:) 738 739 CHARACTER(LEN=ILEN_FLD) :: THRFIL , PSTFIL , PLTFIL , ANNPST , & 740 & ANNPLT , THRFRM , PSTFRM , PLTFRM , & 741 & TOXFIL , SEAHRS , RNKFIL , RNKFRM , & 742 & EVLFIL , ANNPART 743 744 INTEGER , ALLOCATABLE :: NHIAVE(:,:) , MAXAVE(:) , IMXVAL(:) , & 745 & IDYTAB(:) , MAXFLE(:,:) , IPSTFL(:,:) , & 746 & IPLTFL(:,:,:) , IANPST(:) , IANPLT(:) , & 747 & INHI(:) , ITOXFL(:) , ISEAHR(:) , & 748 & IRNKFL(:) , IRKVAL(:) , IANPART(:) 749 REAL , ALLOCATABLE :: THRESH(:,:) , TOXTHR(:) 750 INTEGER , ALLOCATABLE :: IMXUNT(:,:) , IPSUNT(:,:) , IPSFRM(:,:) ,& 751 & IPLUNT(:,:,:) , IAPUNT(:) , IANFRM(:) , & 752 & IPPUNT(:) , ITXUNT(:) , ISHUNT(:) , & 753 & IRKUNT(:) , IELUNT(:) , IUPART(:) 754 755 ALLOCATABLE ::THRFIL(:,:) , PSTFIL(:,:) , PLTFIL(:,:,:) , & 756 & ANNPST(:) , ANNPLT(:) , TOXFIL(:) , SEAHRS(:) , & 757 & RNKFIL(:) , EVLFIL(:) , ANNPART(:) 758 759 INTEGER , ALLOCATABLE :: IDCONC(:,:) 760 761 INTEGER :: ITAB , NXTOX , NYTOX , NHOURS , IPAIR 762 763 REAL , ALLOCATABLE :: TXCONC(:,:) 764 765 766 767 !*********************************************************************** 768 ! This is The Global Variable Definition Block for Working Space 769 !*********************************************************************** 770 771 CHARACTER WORKID*8 , DUMMY*8 772 773 INTEGER :: IMIT , INUM , IDUM , INDAVE , INDGRP , INDVAL , ISC , & 774 & IOERRN , NCPP , NRPP , NGPP , NPPX , NPPY 775 REAL :: FNUM , RNUM 776 777 ALLOCATABLE ::WORKID(:) 778 INTEGER , ALLOCATABLE :: IWRK2(:,:) 779 780 ! Declare Temporary Work Arrays for ZELEV and ZFLAG Receptor Data 781 REAL , ALLOCATABLE :: ZETMP1(:) , ZETMP2(:) 782 REAL , ALLOCATABLE :: ZFTMP1(:) , ZFTMP2(:) 783 REAL , ALLOCATABLE :: ZHTMP1(:) , ZHTMP2(:) 784 785 786 SAVE 787 788 !*********************************************************************** 789 ! Formerly MAIN3.INC 790 ! ISCST2 Model Data - Array Names, Array Limits, Named Common Blocks 791 ! Necessary for Model Results 792 ! MODIFIED - 4/17/95 Output CONC/DEPOS in same model run 793 !*********************************************************************** 794 795 796 !*********************************************************************** 797 ! This is The Global Variable Definition Block For The Maximum 798 ! Value, Highest Value, Average Value, Annual Average Value and 799 ! Model Result Arrays. Also Included are Calm/Missing Flag Arrays. 800 !*********************************************************************** 801 802 803 CHARACTER HCLMSG , MCLMSG , HMCLM 804 805 ! Accomodate Wet SCIM'ing, M. Jindal, PES, 6/10/98 806 REAL , ALLOCATABLE :: HRVAL(:) , AVEVAL(:,:,:,:) , HRVALD(:) , & 807 & HRVALJD(:,:) , AERVAL(:) , PRMVAL(:) , & 808 & AERVALD(:) , PRMVALD(:) 809 REAL , ALLOCATABLE :: HIVALU(:,:,:,:,:) , HMAX(:,:,:,:) 810 INTEGER , ALLOCATABLE :: HMLOC(:,:,:,:) , HMDATE(:,:,:,:) , & 811 & NHIDAT(:,:,:,:,:) 812 813 ! Accomodate Wet SCIM'ing, M. Jindal, PES, 6/10/98 814 REAL , ALLOCATABLE :: ANNVAL(:,:,:) , AMXVAL(:,:,:) , & 815 & SHVALS(:,:,:,:,:) , ANNVALD(:,:,:) , & 816 & ANNVALW(:,:,:) , ANNVALJD(:,:,:,:) , & 817 & ANNVALJW(:,:,:,:) 818 INTEGER , ALLOCATABLE :: IMXLOC(:,:,:) 819 ! Accomodate Wet SCIM'ing, M. Jindal, PES, 6/10/98 820 INTEGER :: IANHRS , IANCLM , IANMSG , NSEAHR(4,24) , NSEACM(4,24) & 821 & , IANWET , IWETCLM , IWETMSG 822 REAL , ALLOCATABLE :: RMXVAL(:,:,:,:) 823 INTEGER , ALLOCATABLE :: MXDATE(:,:,:,:) , MXLOCA(:,:,:,:) 824 INTEGER , ALLOCATABLE :: NUMHRS(:) , NUMCLM(:) , NUMMSG(:) 825 ALLOCATABLE ::HCLMSG(:,:,:,:,:) , MCLMSG(:,:,:,:) , HMCLM(:,:,:,:) 826 827 REAL , ALLOCATABLE :: SUMANN(:,:,:) 828 REAL , ALLOCATABLE :: SUMH4H(:,:) , MXPMVAL(:,:) 829 INTEGER , ALLOCATABLE :: MXPMLOC(:,:) 830 831 REAL , ALLOCATABLE :: CHI(:,:,:) , HECNTR(:,:) , HECNTR3(:,:) , & 832 & PPFACT(:) , UEFFS(:,:) , UEFF3S(:,:) , & 833 & FOPTS(:,:) 834 835 REAL , ALLOCATABLE :: ARCMAX(:) , QMAX(:) , DXMAX(:) , UMAX(:) , & 836 & SVMAX(:) , SWMAX(:) , SYMAX(:) , SY3MX(:) , & 837 & U3MAX(:) , HEMAX(:) , ARCCL(:) , SZMAX(:) , & 838 & CHIDMW(:) , CHINMW(:) , CHI3MW(:) , & 839 & CHIDML(:) , CHINML(:) , CHI3ML(:) , & 840 & HSBLMX(:) 841 842 !*********************************************************************** 843 ! This is The Global Variable Definition Block For The 844 ! EVENT Model Result Arrays 845 !*********************************************************************** 846 847 REAL , ALLOCATABLE :: EV_AVEVAL(:) , HRVALS(:,:) , GRPVAL(:) 848 849 REAL :: GRPAVE 850 INTEGER :: EV_NUMHRS , EV_NUMCLM , EV_NUMMSG , ISTAHR , IENDHR 851 852 853 !*********************************************************************** 854 ! 855 ! BLOCK DATA SUBPROGRAM OF THE ISC - Version 2 MODEL 856 ! 857 ! PURPOSE: Initialize Data in COMMON Blocks 858 ! 859 ! MODIFIED: To Include Terrain Grid pathway - 12/15/93 860 ! 861 ! MODIFIED: To Include WET DEPOSITION Arrays - 11/8/93 862 ! 863 ! MODIFIED: To Include New Area Source Arrays - 7/7/93 864 ! 865 ! MODIFIED: For revised DRY DEPOSITION code - 2/15/93 866 ! 867 ! MODIFIED: To Include TOXXFILE Option - 9/29/92 868 ! 869 !*********************************************************************** 870 871 ! Variable Declarations 872 873 874 !*********************************************************************** 875 ! Initialize Model Version Number, VERSN (Year, Julian Day), as a 876 ! Character Variable 877 !*********************************************************************** 878 879 DATA VERSN/'04300'/ 880 881 882 !*********************************************************************** 883 ! Input/Output File Units and Input/Output File Names 884 !*********************************************************************** 885 886 !LF95 Change INUNIT from 5 to 7, and IOUNIT from 6 to 8, for the 887 !LF95 Lahey LF95 compiler (version 5.0). 888 DATA INUNIT/7/ , IOUNIT/8/ , MFUNIT/19/ , MPUNIT/21/ , & 889 & IERUNT/10/ , IERWRT/11/ , IDPUNT/12/ , IZUNIT/13/ , & 890 & IDPUN2/14/ , IRSUNT/15/ , IHREMI/16/ , IEVUNT/17/ , & 891 & ITEVUT/18/ , INCUNT/20/ , ISUNIT/22/ , IPUNIT/23/ , & 892 & DBGUNT/24/ , DBMUNT/25/ , ICVUNT/26/ , IO3UNT/29/ 893 !*# 894 DATA INPFIL/' '/ , OUTFIL/' '/ 895 896 897 !*********************************************************************** 898 ! Initialize Keyword Array 899 !*********************************************************************** 900 901 INTEGER , PRIVATE :: I 902 903 904 ! --- PRIME ------------------------------ 905 ! ---------------------------------------- 906 907 DATA (KEYWD(I),I=1,IKN)/'STARTING' , 'FINISHED' , 'TITLEONE' , & 908 & 'TITLETWO' , 'MODELOPT' , 'AVERTIME' , 'POLLUTID' , & 909 & 'HALFLIFE' , 'DCAYCOEF' , 'DEBUGOPT' , 'ELEVUNIT' , & 910 & 'FLAGPOLE' , 'RUNORNOT' , 'EVENTFIL' , 'SAVEFILE' , & 911 & 'INITFILE' , 'MULTYEAR' , 'ERRORFIL' , 'GASDEPDF' , & 912 & 'GDSEASON' , 'GASDEPVD' , 'GDLANUSE' , 'EVENTFIL' , & 913 & 'URBANOPT' , 'METHOD_2' , 'LOCATION' , 'SRCPARAM' , & 914 & 'BUILDHGT' , 'BUILDWID' , 'BUILDLEN' , 'XBADJ ' , & 915 & 'YBADJ ' , 'EMISFACT' , 'EMISUNIT' , 'PARTDIAM' , & 916 & 'MASSFRAX' , 'PARTDENS' , ' ' , ' ' , & 917 & ' ' , 'CONCUNIT' , 'DEPOUNIT' , 'HOUREMIS' , & 918 & 'GASDEPOS' , 'URBANSRC' , 'EVENTPER' , 'EVENTLOC' , & 919 & 'SRCGROUP' , 'GRIDCART' , 'GRIDPOLR' , 'DISCCART' , & 920 & 'DISCPOLR' , 'SURFFILE' , 'PROFFILE' , 'PROFBASE' , & 921 & ' ' , 'SURFDATA' , 'UAIRDATA' , 'SITEDATA' , & 922 & 'STARTEND' , 'DAYRANGE' , 'WDROTATE' , 'DTHETADZ' , & 923 & 'WINDCATS' , 'RECTABLE' , 'MAXTABLE' , 'DAYTABLE' , & 924 & 'MAXIFILE' , 'POSTFILE' , 'PLOTFILE' , 'TOXXFILE' , & 925 & 'EVENTOUT' , 'INCLUDED' , 'SCIMBYHR' , 'SEASONHR' , & 926 & 'AREAVERT' , 'PARTSIZE' , 'RANKFILE' , 'EVALCART' , & 927 & 'EVALFILE' , 'NO2EQUIL' , 'OZONEVAL' , 'OZONEFIL' , & 928 & 'NO2RATIO' , 'OLMGROUP'/ 929 930 931 !*********************************************************************** 932 ! Initialize Miscellaneous Variables 933 !*********************************************************************** 934 935 !JRA removed space after "-" 936 DATA IPROC/366*1/ , EXPLIM/ -50.0/ 937 DATA UCAT/1.54 , 3.09 , 5.14 , 8.23 , 10.8/ 938 DATA MODOPS/18*' '/ 939 940 941 !*********************************************************************** 942 ! Initialize distance factors used in determining when to switch 943 ! to point source approximation for area sources under the TOXICS 944 ! option. 945 !*********************************************************************** 946 947 ! STAB. CLASS A B C D E F 948 ! *** *** *** *** *** *** 949 DATA VIRTPNT_URB/3.5 , 3.5 , 5.5 , 10.5 , 15.5 , 15.5/ , & 950 & VIRTPNT_RUR/3.5 , 5.5 , 7.5 , 12.5 , 15.5 , 25.5/ 951 952 953 !*********************************************************************** 954 ! Initialize Setup Status Arrays 955 !*********************************************************************** 956 957 DATA ICSTAT/30*0/ , ISSTAT/30*0/ , IRSTAT/30*0/ , IMSTAT/30*0/ , & 958 & IOSTAT/30*0/ , ITSTAT/30*0/ 959 960 961 !*********************************************************************** 962 ! Initialize Gridded Profile Height Array 963 !*********************************************************************** 964 965 DATA GRIDHT/0.0 , 0.5 , 1.0 , 2.0 , 4.0 , 8.0 , 14.0 , 20.0 , & 966 & 30.0 , 40.0 , 50.0 , 60.0 , 70.0 , 80.0 , 90.0 , 100.0 , & 967 & 120.0 , 140.0 , 160.0 , 180.0 , 200.0 , 250.0 , 300.0 , & 968 & 350.0 , 400.0 , 450.0 , 500.0 , 550.0 , 600.0 , 650.0 , & 969 & 700.0 , 750.0 , 800.0 , 850.0 , 900.0 , 950.0 , 1000.0 , & 970 & 1050.0 , 1100.0 , 1150.0 , 1200.0 , 1250.0 , 1300.0 , & 971 & 1350.0 , 1400.0 , 1450.0 , 1500.0 , 1550.0 , 1600.0 , & 972 & 1650.0 , 1700.0 , 1750.0 , 1800.0 , 1850.0 , 1900.0 , & 973 & 1950.0 , 2000.0 , 2100.0 , 2200.0 , 2300.0 , 2400.0 , & 974 & 2500.0 , 2600.0 , 2700.0 , 2800.0 , 2900.0 , 3000.0 , & 975 & 3100.0 , 3200.0 , 3300.0 , 3400.0 , 3500.0 , 3600.0 , & 976 & 3700.0 , 3800.0 , 3900.0 , 4000.0 , 4100.0 , 4200.0 , & 977 & 4300.0 , 4400.0 , 4500.0 , 4600.0 , 4700.0 , 4800.0 , & 978 & 4900.0 , 5000.0/ 979 980 981 !*********************************************************************** 982 ! Initialize Error Code and Message Arrays 983 !*********************************************************************** 984 985 DATA ERRCOD(1)/'100'/ , ERRMSG(1) & 986 & /'Invalid Pathway Specified. The Troubled Pathway is'/ 987 DATA ERRCOD(2)/'105'/ , ERRMSG(2) & 988 & /'Invalid Keyword Specified. The Troubled Keyword is'/ 989 DATA ERRCOD(3)/'110'/ , ERRMSG(3) & 990 & /'Keyword is Not Valid for This Pathway. Keyword is'/ 991 DATA ERRCOD(4)/'115'/ , ERRMSG(4) & 992 & /'STARTING or FINISHED Out of Sequence: Pathway = '/ 993 DATA ERRCOD(5)/'120'/ , ERRMSG(5) & 994 & /'Pathway is Out of Sequence: Pathway = '/ 995 DATA ERRCOD(6)/'125'/ , ERRMSG(6) & 996 & /'Missing FINISHED-Runstream File Incomplete: ISTAT='/ 997 DATA ERRCOD(7)/'130'/ , ERRMSG(7) & 998 & /'Missing Mandatory Keyword. The Missing Keyword is'/ 999 DATA ERRCOD(8)/'135'/ , ERRMSG(8) & 1000 & /'Duplicate Nonrepeatable Keyword Specified:Keyword='/ 1001 DATA ERRCOD(9)/'140'/ , ERRMSG(9) & 1002 & /'Invalid Order of Keyword. The Troubled Keyword is'/ 1003 DATA ERRCOD(10)/'141'/ , ERRMSG(10) & 1004 & /'Conflicting Options: PVMRM and OLM both specified'/ 1005 DATA ERRCOD(11)/'142'/ , ERRMSG(11) & 1006 & /'Following Keyword Invalid Without PVMRM or OLM: '/ 1007 DATA ERRCOD(12)/'143'/ , ERRMSG(12) & 1008 & /'Following Keyword Invalid Without PVMRM Option: '/ 1009 DATA ERRCOD(13)/'144'/ , ERRMSG(13) & 1010 & /'Following Keyword Invalid Without OLM Option: '/ 1011 DATA ERRCOD(14)/'145'/ , ERRMSG(14) & 1012 & /'Conflicting Options: MULTYEAR and Re-Start Option '/ 1013 DATA ERRCOD(15)/'150'/ , ERRMSG(15) & 1014 & /'Conflicting Options: MULTYEAR for Wrong Pollutant '/ 1015 DATA ERRCOD(16)/'152'/ , ERRMSG(16) & 1016 & /'ELEVUNIT card must be first for this Pathway: '/ 1017 DATA ERRCOD(17)/'154'/ , ERRMSG(17) & 1018 & /'Conflicting options: SCIM cannot be used with '/ 1019 DATA ERRCOD(18)/'155'/ , ERRMSG(18) & 1020 & /'Conflicting Decay Keyword. Inputs Ignored for '/ 1021 DATA ERRCOD(19)/'156'/ , ERRMSG(19) & 1022 & /'Option ignored - not valid with SCIM. Option = '/ 1023 DATA ERRCOD(20)/'157'/ , ERRMSG(20) & 1024 & /'Wet SCIM Option Not Operational Yet. Input Ignored'/ 1025 DATA ERRCOD(21)/'158'/ , ERRMSG(21) & 1026 & /'EMISUNIT Keyword Used With More Than 1 Output Type'/ 1027 DATA ERRCOD(22)/'159'/ , ERRMSG(22) & 1028 & /'EMISUNIT Keyword Used With the Following Keyword: '/ 1029 DATA ERRCOD(23)/'160'/ , ERRMSG(23) & 1030 & /'Duplicate ORIG Secondary Keyword for GRIDPOLR: '/ 1031 DATA ERRCOD(24)/'170'/ , ERRMSG(24) & 1032 & /'Invalid Secondary Keyword for Receptor Grid: '/ 1033 DATA ERRCOD(25)/'175'/ , ERRMSG(25) & 1034 & /'Missing Secondary Keyword END for Receptor Grid: '/ 1035 DATA ERRCOD(26)/'180'/ , ERRMSG(26) & 1036 & /'Conflicting Secondary Keyword for Receptor Grid: '/ 1037 DATA ERRCOD(27)/'185'/ , ERRMSG(27) & 1038 & /'Missing Receptor Keywords. No Receptors Specified.'/ 1039 DATA ERRCOD(28)/'190'/ , ERRMSG(28) & 1040 & /'No Keywords for OU Path and No PERIOD/ANNUAL Aves.'/ 1041 DATA ERRCOD(29)/'195'/ , ERRMSG(29) & 1042 & /'Incompatible Option Used With SAVEFILE or INITFILE'/ 1043 DATA ERRCOD(30)/'196'/ , ERRMSG(30) & 1044 & /'Incompatible Keyword Used With GASDEPVD '/ 1045 DATA ERRCOD(31)/'197'/ , ERRMSG(31) & 1046 & /'Post-97 PM10 without MAXIFILE is incompatible with'/ 1047 DATA ERRCOD(32)/'198'/ , ERRMSG(32) & 1048 & /'TOXICS Option is Required in Order to Use Option '/ 1049 1050 DATA ERRCOD(33)/'200'/ , ERRMSG(33) & 1051 & /'Missing Parameter(s). No Options Specified For '/ 1052 DATA ERRCOD(34)/'201'/ , ERRMSG(34) & 1053 & /'Not Enough Parameters Specified For the Keyword of'/ 1054 DATA ERRCOD(35)/'202'/ , ERRMSG(35) & 1055 & /'Too Many Parameters Specified For the Keyword of '/ 1056 DATA ERRCOD(36)/'203'/ , ERRMSG(36) & 1057 & /'Invalid Parameter Specified. Troubled Parameter: '/ 1058 DATA ERRCOD(37)/'204'/ , ERRMSG(37) & 1059 & /'Option Parameters Conflict. Forced by Default to '/ 1060 DATA ERRCOD(38)/'205'/ , ERRMSG(38) & 1061 & /'No Option Parameter Setting. Forced by Default to'/ 1062 DATA ERRCOD(39)/'206'/ , ERRMSG(39) & 1063 & /'Regulatory DFAULT Overrides Non-DFAULT Option For '/ 1064 DATA ERRCOD(40)/'207'/ , ERRMSG(40) & 1065 & /'No Parameters Specified. Default Values Will Used.'/ 1066 DATA ERRCOD(41)/'208'/ , ERRMSG(41) & 1067 & /'Illegal Numerical Field Encountered in '/ 1068 DATA ERRCOD(42)/'209'/ , ERRMSG(42) & 1069 & /'Negative Value Appears For Non-negative Variable. '/ 1070 DATA ERRCOD(43)/'210'/ , ERRMSG(43) & 1071 & /'Number of Short Term Averages Exceeds Max: NAVE= '/ 1072 DATA ERRCOD(44)/'211'/ , ERRMSG(44) & 1073 & /'Duplicate Averaging Period Specified for Keyword '/ 1074 DATA ERRCOD(45)/'212'/ , ERRMSG(45) & 1075 & /'END Encountered Without (X,Y) Points Properly Set '/ 1076 DATA ERRCOD(46)/'213'/ , ERRMSG(46) & 1077 & /'ELEV Input Inconsistent With Option: Input Ignored'/ 1078 DATA ERRCOD(47)/'214'/ , ERRMSG(47) & 1079 & /'ELEV Input Inconsistent With Option: Defaults Used'/ 1080 DATA ERRCOD(48)/'215'/ , ERRMSG(48) & 1081 & /'FLAG Input Inconsistent With Option: Input Ignored'/ 1082 DATA ERRCOD(49)/'216'/ , ERRMSG(49) & 1083 & /'FLAG Input Inconsistent With Option: Defaults Used'/ 1084 DATA ERRCOD(50)/'217'/ , ERRMSG(50) & 1085 & /'More Than One Delimiter In A Field for Keyword '/ 1086 DATA ERRCOD(51)/'218'/ , ERRMSG(51) & 1087 & /'Number of (X,Y) Points Not Match With Number Of '/ 1088 DATA ERRCOD(52)/'219'/ , ERRMSG(52) & 1089 & /'Number Of Receptors Specified Exceeds Max: NREC= '/ 1090 DATA ERRCOD(53)/'220'/ , ERRMSG(53) & 1091 & /'Missing Origin (Use Default = 0,0) In GRIDPOLR '/ 1092 DATA ERRCOD(54)/'221'/ , ERRMSG(54) & 1093 & /'Missing Distance Setting In Polar Network '/ 1094 DATA ERRCOD(55)/'222'/ , ERRMSG(55) & 1095 & /'Missing Degree Or Dist Setting In Polar Network '/ 1096 DATA ERRCOD(56)/'223'/ , ERRMSG(56) & 1097 & /'Missing Distance or Degree Field in '/ 1098 DATA ERRCOD(57)/'224'/ , ERRMSG(57) & 1099 & /'Number of Receptor Networks Exceeds Max: NNET= '/ 1100 DATA ERRCOD(58)/'225'/ , ERRMSG(58) & 1101 & /'Number of X-Coords Specified Exceeds Max: IXM= '/ 1102 DATA ERRCOD(59)/'226'/ , ERRMSG(59) & 1103 & /'Number of Y-Coords Specified Exceeds Max: IYM= '/ 1104 DATA ERRCOD(60)/'227'/ , ERRMSG(60) & 1105 & /'No Receptors Were Defined on the RE Pathway. '/ 1106 DATA ERRCOD(61)/'228'/ , ERRMSG(61) & 1107 & /'Default(s) Used for Missing Parameters on Keyword '/ 1108 DATA ERRCOD(62)/'229'/ , ERRMSG(62) & 1109 & /'Too Many Parameters - Inputs Ignored on Keyword '/ 1110 DATA ERRCOD(63)/'231'/ , ERRMSG(63) & 1111 & /'Too Many Numerical Values Specified for '/ 1112 DATA ERRCOD(64)/'232'/ , ERRMSG(64) & 1113 & /'Number Of Specified Sources Exceeds Maximum: NSRC='/ 1114 DATA ERRCOD(65)/'233'/ , ERRMSG(65) & 1115 & /'Building Dimensions Specified for Non-POINT Source'/ 1116 DATA ERRCOD(66)/'234'/ , ERRMSG(66) & 1117 & /'Too Many Sectors Input for '/ 1118 DATA ERRCOD(67)/'235'/ , ERRMSG(67) & 1119 & /'Number of Source Groups Exceeds Maximum: NGRP= '/ 1120 DATA ERRCOD(68)/'236'/ , ERRMSG(68) & 1121 & /'Not Enough BUILDHGTs Specified for SourceID '/ 1122 DATA ERRCOD(69)/'237'/ , ERRMSG(69) & 1123 & /'Not Enough BUILDWIDs Specified for SourceID '/ 1124 DATA ERRCOD(70)/'239'/ , ERRMSG(70) & 1125 & /'Not Enough QFACTs Specified for SourceID '/ 1126 DATA ERRCOD(71)/'240'/ , ERRMSG(71) & 1127 & /'Inconsistent Number of Particle Categories for '/ 1128 DATA ERRCOD(72)/'241'/ , ERRMSG(72) & 1129 & /'Not Enough BUILDLENs Specified for SourceID '/ 1130 DATA ERRCOD(73)/'242'/ , ERRMSG(73) & 1131 & /'No Particle Cat. or Gas Depos. Specified for SRCID'/ 1132 DATA ERRCOD(74)/'243'/ , ERRMSG(74) & 1133 & /'Scav. Coef. may be out-of-range for SRCID '/ 1134 DATA ERRCOD(75)/'244'/ , ERRMSG(75) & 1135 & /'Too Many Particle Categories Specified for '/ 1136 DATA ERRCOD(76)/'245'/ , ERRMSG(76) & 1137 & /'No. of Particle Categories Exceeds Max: NPDMAX= '/ 1138 DATA ERRCOD(77)/'246'/ , ERRMSG(77) & 1139 & /'Not Enough XBADJs Specified for SourceID '/ 1140 DATA ERRCOD(78)/'247'/ , ERRMSG(78) & 1141 & /'Not Enough YBADJs Specified for SourceID '/ 1142 DATA ERRCOD(79)/'248'/ , ERRMSG(79) & 1143 & /'No Sources Were Defined on the SO Pathway. '/ 1144 DATA ERRCOD(80)/'250'/ , ERRMSG(80) & 1145 & /'Duplicate XPNT/DIST or YPNT/DIR Specified for GRID'/ 1146 DATA ERRCOD(81)/'252'/ , ERRMSG(81) & 1147 & /'Duplicate Receptor Network ID Specified. NETID = '/ 1148 DATA ERRCOD(82)/'254'/ , ERRMSG(82) & 1149 & /'Number of Receptor ARCs Exceeds Max: NARC= '/ 1150 DATA ERRCOD(83)/'256'/ , ERRMSG(83) & 1151 & /'EVALFILE Option Used Without EVALCART Receptors '/ 1152 DATA ERRCOD(84)/'260'/ , ERRMSG(84) & 1153 & /'Number of Emission Factors Exceeds Max: NQF= '/ 1154 DATA ERRCOD(85)/'262'/ , ERRMSG(85) & 1155 & /'First Vertex Does Not Match LOCATION for AREAPOLY '/ 1156 DATA ERRCOD(86)/'264'/ , ERRMSG(86) & 1157 & /'Too Many Vertices Specified for AREAPOLY Source '/ 1158 DATA ERRCOD(87)/'265'/ , ERRMSG(87) & 1159 & /'Not Enough Vertices Specified for AREAPOLY Source '/ 1160 DATA ERRCOD(88)/'270'/ , ERRMSG(88) & 1161 & /'Number of High Values Specified Exceeds Max: NVAL='/ 1162 DATA ERRCOD(89)/'280'/ , ERRMSG(89) & 1163 & /'Number of Max Values Specified Exceeds Max: NMAX='/ 1164 DATA ERRCOD(90)/'281'/ , ERRMSG(90) & 1165 & /'Number of OLMGROUPs Specified Exceeds Max: NOLM= '/ 1166 DATA ERRCOD(91)/'282'/ , ERRMSG(91) & 1167 & /'Following SRCID Included in Multiple OLMGROUPs: '/ 1168 DATA ERRCOD(92)/'283'/ , ERRMSG(92) & 1169 & /'Either OZONEVAL or OZONEFIL Card Needed for Option'/ 1170 DATA ERRCOD(93)/'284'/ , ERRMSG(93) & 1171 & /'Invalid POLLUTID Specified for PVMRM/OLM; Must Use'/ 1172 DATA ERRCOD(94)/'290'/ , ERRMSG(94) & 1173 & /'Number of Output Types Specified Exceeds Max:NTYP='/ 1174 DATA ERRCOD(95)/'294'/ , ERRMSG(95) & 1175 & /'PERIOD and ANNUAL averages are both selected for '/ 1176 DATA ERRCOD(96)/'295'/ , ERRMSG(96) & 1177 & /'Invalid Averaging Period Specified for SCREEN Mode'/ 1178 DATA ERRCOD(97)/'298'/ , ERRMSG(97) & 1179 & /'Error Allocating Storage for Setup Arrays! '/ 1180 DATA ERRCOD(98)/'299'/ , ERRMSG(98) & 1181 & /'Error Allocating Storage for Result Arrays! '/ 1182 1183 DATA ERRCOD(99)/'300'/ , ERRMSG(99) & 1184 & /'Specified SRCID Has Not Been Defined Yet: KEYWORD='/ 1185 DATA ERRCOD(100)/'305'/ , ERRMSG(100) & 1186 & /'Terrain Grid Does Not Cover Modeling Area, Change:'/ 1187 DATA ERRCOD(101)/'310'/ , ERRMSG(101) & 1188 & /'Attempt to Define Duplicate LOCATION Card for SRC:'/ 1189 DATA ERRCOD(102)/'313'/ , ERRMSG(102) & 1190 & /'Attempt to Define Duplicate EVENTPER card for '/ 1191 DATA ERRCOD(103)/'315'/ , ERRMSG(103) & 1192 & /'Attempt to Define Duplicate SRCPARAM Card for SRC:'/ 1193 DATA ERRCOD(104)/'319'/ , ERRMSG(104) & 1194 & /'No Sources Included in Specified Source Group: '/ 1195 DATA ERRCOD(105)/'320'/ , ERRMSG(105) & 1196 & /'Input Parameter May Be Out-of-Range for Parameter '/ 1197 DATA ERRCOD(106)/'322'/ , ERRMSG(106) & 1198 & /'Release Height Exceeds Effective Depth for OPENPIT'/ 1199 DATA ERRCOD(107)/'323'/ , ERRMSG(107) & 1200 & /'No Particle Categories Specified for OPENPIT Src. '/ 1201 DATA ERRCOD(108)/'325'/ , ERRMSG(108) & 1202 & /'Negative Exit Velocity (Set=1.0E-5) for SRCID: '/ 1203 DATA ERRCOD(109)/'330'/ , ERRMSG(109) & 1204 & /'Mass Fraction Parameters Do Not Sum to 1. for Src '/ 1205 DATA ERRCOD(110)/'332'/ , ERRMSG(110) & 1206 & /'Mass Fraction Parameter Out-of-Range for Source '/ 1207 DATA ERRCOD(111)/'334'/ , ERRMSG(111) & 1208 & /'Particle Density Out-of-Range for Source '/ 1209 DATA ERRCOD(112)/'336'/ , ERRMSG(112) & 1210 & /'Value Specified for NO2RATIO is Out-of-Range for '/ 1211 DATA ERRCOD(113)/'338'/ , ERRMSG(113) & 1212 & /'Neg Emis Rate Cannot be Used with OLM/PVMRM. Src: '/ 1213 DATA ERRCOD(114)/'340'/ , ERRMSG(114) & 1214 & /'Possible Error in PROFBASE Input: Value is < 0 '/ 1215 DATA ERRCOD(115)/'342'/ , ERRMSG(115) & 1216 & /'Src ID Mismatch in Hourly Emissions File for ID = '/ 1217 DATA ERRCOD(116)/'344'/ , ERRMSG(116) & 1218 & /'Hourly Emission Rate is Zero for KURDAT = '/ 1219 DATA ERRCOD(117)/'350'/ , ERRMSG(117) & 1220 & /'Julian Day Out Of Range at '/ 1221 DATA ERRCOD(118)/'352'/ , ERRMSG(118) & 1222 & /'Missing Field on MULTYEAR Card for Pre-1997 PM10 '/ 1223 DATA ERRCOD(119)/'353'/ , ERRMSG(119) & 1224 & /'MULTYEAR Card for PM10 Processing Applies Only for'/ 1225 DATA ERRCOD(120)/'354'/ , ERRMSG(120) & 1226 & /'High-4th-High Only Required for Post-1997 PM10 '/ 1227 DATA ERRCOD(121)/'360'/ , ERRMSG(121) & 1228 & /'2-Digit Year Specified: Valid for Range 1950-2049 '/ 1229 DATA ERRCOD(122)/'363'/ , ERRMSG(122) & 1230 & /'24HR and ANNUAL Averages Only for Post-1997 PM10 '/ 1231 DATA ERRCOD(123)/'365'/ , ERRMSG(123) & 1232 & /'Year Input is Greater Than 2147 '/ 1233 DATA ERRCOD(124)/'370'/ , ERRMSG(124) & 1234 & /'Invalid Date: 2/29 In a Non-leap Year. '/ 1235 DATA ERRCOD(125)/'380'/ , ERRMSG(125) & 1236 & /'This Input Variable is Out-of-Range: '/ 1237 DATA ERRCOD(126)/'381'/ , ERRMSG(126) & 1238 & /'Latitude in Surface File Is Not Valid: '/ 1239 DATA ERRCOD(127)/'382'/ , ERRMSG(127) & 1240 & /'Error Decoding Latitude: '/ 1241 DATA ERRCOD(128)/'383'/ , ERRMSG(128) & 1242 & /'NWETFREQ > 0, but Wet Dep/Depletion not selected '/ 1243 DATA ERRCOD(129)/'385'/ , ERRMSG(129) & 1244 & /'Averaging Period .NE. 1-Hr for TOXXFILE Option '/ 1245 DATA ERRCOD(130)/'390'/ , ERRMSG(130) & 1246 & /'Aver. Period must be .LE. 24 for EVENT Processing '/ 1247 DATA ERRCOD(131)/'391'/ , ERRMSG(131) & 1248 & /'Aspect ratio (L/W) of area source greater than 10 '/ 1249 DATA ERRCOD(132)/'392'/ , ERRMSG(132) & 1250 & /'Aspect ratio (L/W) of open pit is greater than 10 '/ 1251 DATA ERRCOD(133)/'393'/ , ERRMSG(133) & 1252 & /'Terrain Grid Value Differs >50% From Source Elev. '/ 1253 DATA ERRCOD(134)/'394'/ , ERRMSG(134) & 1254 & /'Terrain Grid Value Differs >50% From Receptor Elev'/ 1255 DATA ERRCOD(135)/'395'/ , ERRMSG(135) & 1256 & /'Met. Data Error; Incompatible Version of AERMET: '/ 1257 DATA ERRCOD(136)/'396'/ , ERRMSG(136) & 1258 & /'Met. Data Generated by Older Version of AERMET: '/ 1259 1260 DATA ERRCOD(137)/'405'/ , ERRMSG(137) & 1261 & /'Value of PHEE Exceeds 1.0 on KURDAT = '/ 1262 DATA ERRCOD(138)/'406'/ , ERRMSG(138) & 1263 & /'Increase NVMAX for Complex AREAPOLY Source '/ 1264 DATA ERRCOD(139)/'410'/ , ERRMSG(139) & 1265 & /'Wind Direction Out-of-Range. KURDAT= '/ 1266 DATA ERRCOD(140)/'413'/ , ERRMSG(140) & 1267 & /'Number of Threshold Events > 9999 for Ave Period '/ 1268 DATA ERRCOD(141)/'420'/ , ERRMSG(141) & 1269 & /'Wind Speed Out-of-Range. KURDAT= '/ 1270 DATA ERRCOD(142)/'430'/ , ERRMSG(142) & 1271 & /'Ambient Temperature Data Out-of-Range. KURDAT= '/ 1272 DATA ERRCOD(143)/'432'/ , ERRMSG(143) & 1273 & /'Friction Velocity Out-of-Range. KURDAT= '/ 1274 DATA ERRCOD(144)/'435'/ , ERRMSG(144) & 1275 & /'Surface Roughness Length Out-of-Range. KURDAT= '/ 1276 DATA ERRCOD(145)/'438'/ , ERRMSG(145) & 1277 & /'Convective Velocity Data Out-of-Range. KURDAT= '/ 1278 DATA ERRCOD(146)/'439'/ , ERRMSG(146) & 1279 & /'Monin-Obukhov Length Out-of-Range. KURDAT= '/ 1280 DATA ERRCOD(147)/'440'/ , ERRMSG(147) & 1281 & /'Calm Hour Identified in Meteorology Data File at '/ 1282 DATA ERRCOD(148)/'450'/ , ERRMSG(148) & 1283 & /'Error in Meteor. File - Record Out of Sequence at '/ 1284 DATA ERRCOD(149)/'455'/ , ERRMSG(149) & 1285 & /'Date/time Mismatch: Hourly Emission File. KURDAT ='/ 1286 DATA ERRCOD(150)/'456'/ , ERRMSG(150) & 1287 & /'Date/time Mismatch on Surface & Profile. KURDAT = '/ 1288 DATA ERRCOD(151)/'457'/ , ERRMSG(151) & 1289 & /'Date/time Mismatch on OZONEFIL Data. KURDAT = '/ 1290 DATA ERRCOD(152)/'458'/ , ERRMSG(152) & 1291 & /'Substitution made for missing ozone data. KURDAT= '/ 1292 DATA ERRCOD(153)/'459'/ , ERRMSG(153) & 1293 & /'Missing ozone data; Full conversion used. KURDAT= '/ 1294 DATA ERRCOD(154)/'460'/ , ERRMSG(154) & 1295 & /'Missing Hour Identified in Meteor. Data File at '/ 1296 DATA ERRCOD(155)/'465'/ , ERRMSG(155) & 1297 & /'Number of Profile Levels Exceeds Max: MXPLVL= '/ 1298 DATA ERRCOD(156)/'470'/ , ERRMSG(156) & 1299 & /'Mixing Height Value is < or = 0.0. KURDAT= '/ 1300 DATA ERRCOD(157)/'475'/ , ERRMSG(157) & 1301 & /'Reference height is higher than 100m. KURDAT= '/ 1302 DATA ERRCOD(158)/'480'/ , ERRMSG(158) & 1303 & /'Less Than 1 Year Found for ANNUAL Averages '/ 1304 DATA ERRCOD(159)/'485'/ , ERRMSG(159) & 1305 & /'Data Remaining After End of Year. Number of Hours='/ 1306 DATA ERRCOD(160)/'487'/ , ERRMSG(160) & 1307 & /'User Start Date is Earlier Than Start of Data File'/ 1308 1309 DATA ERRCOD(161)/'500'/ , ERRMSG(161) & 1310 & /'Fatal Error Occurs Opening the Data File of '/ 1311 DATA ERRCOD(162)/'510'/ , ERRMSG(162) & 1312 & /'Fatal Error Occurs During Reading of the File of '/ 1313 DATA ERRCOD(163)/'520'/ , ERRMSG(163) & 1314 & /'Fatal Error Occurs During Writing to the File of '/ 1315 DATA ERRCOD(164)/'530'/ , ERRMSG(164) & 1316 & /'CAUTION! Met Station ID Mismatch with SURFFILE for'/ 1317 DATA ERRCOD(165)/'540'/ , ERRMSG(165) & 1318 & /'No RECTABLE/MAXTABLE/DAYTABLE for Average Period '/ 1319 DATA ERRCOD(166)/'550'/ , ERRMSG(166) & 1320 & /'File Unit/Name Conflict for the Output Option: '/ 1321 DATA ERRCOD(167)/'560'/ , ERRMSG(167) & 1322 & /'User Specified File Unit .LE. 25 for OU Keyword: '/ 1323 DATA ERRCOD(168)/'565'/ , ERRMSG(168) & 1324 & /'Possible Conflict With Dynamically Allocated FUNIT'/ 1325 DATA ERRCOD(169)/'570'/ , ERRMSG(169) & 1326 & /'Problem Reading Temporary Event File for Event: '/ 1327 DATA ERRCOD(170)/'580'/ , ERRMSG(170) & 1328 & /'End of File Reached Trying to Read the File of '/ 1329 1330 1331 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