Appendix 3 CALMET/CALPUFF Input Files

Transcription

1 Innovative Development & Design Engineers Ltd. Engineering Services 3D Modelling Anderson House, Tel: (709) Power s Court, Cell: (709) St. John s, NL A1A 1B6 Web: Appendix 3 CALMET/CALPUFF Input Files 31 January 2018

2 1 CALMET.INP 2.2 Hour Start and End Times with Seconds 2 Example CALMET simulation (not intended as a guide to option selections) 3 MM4 data, 5 surface met stations, 1 overwater station, - for CALMET V upper air met and 16 precip stations- New parameter: ISTEPPGS Run title (3 lines) CALMET MODEL CONTROL FILE INPUT GROUP: 0 -- Input and Output File Names Subgroup (a) Default Name Type File Name GEO.DAT input! GEODAT=GEO_SUMMER.DAT! 20 SURF.DAT input * SRFDAT= * 21 CLOUD.DAT input * CLDDAT= * 22 PRECIP.DAT input * PRCDAT= * 23 WT.DAT input * WTDAT= * CALMET.LST output! METLST=CALMET.LST! 26 CALMET.DAT output! METDAT=CALMET.DAT! 27 PACOUT.DAT output * PACDAT= * All file names will be converted to lower case if LCFILES = T 30 Otherwise, if LCFILES = F, file names will be converted to UPPER CASE 31 T = lower case! LCFILES = F! 32 F = UPPER CASE NUMBER OF UPPER AIR & OVERWATER STATIONS: Number of upper air stations (NUSTA) No default! NUSTA = 0! 37 Number of overwater met stations 38 (NOWSTA) No default! NOWSTA = 0! NUMBER OF PROGNOSTIC and IGF-CALMET FILEs: Number of MM4/MM5/3D.DAT files 43 (NM3D) No default! NM3D = 7! Number of IGF-CALMET.DAT files 46 (NIGF) No default! NIGF = 0! 47 48!END! Subgroup (b) Upper air files (one per station) Default Name Type File Name UP1.DAT input 1 * UPDAT=UPPWM.DAT* *END* 57 UP2.DAT input 2 * UPDAT=UPALBR.DAT* *END* 58 UP3.DAT input 3 * UPDAT=UPCHH.DAT* *END* Subgroup (c) Overwater station files (one per station) Default Name Type File Name SEA1.DAT input 1 * SEADAT=4007.DAT* *END* Subgroup (d)

3 70 MM4/MM5/3D.DAT files (consecutive or overlapping) Default Name Type File Name MM51.DAT input 1! M3DDAT=WRF_ _St_Johns_29.DAT!!END! 75 MM52.DAT input 2! M3DDAT=WRF_ _St_Johns_30.DAT!!END! 76 MM53.DAT input 3! M3DDAT=WRF_ _St_Johns_31.DAT!!END! 77 MM54.DAT input 4! M3DDAT=WRF_ _St_Johns_32.DAT!!END! 78 MM55.DAT input 5! M3DDAT=WRF_ _St_Johns_33.DAT!!END! 79 MM56.DAT input 6! M3DDAT=WRF_ _St_Johns_34.DAT!!END! 80 MM56.DAT input 7! M3DDAT=WRF_ _St_Johns_35.DAT!!END! Subgroup (e) IGF-CALMET.DAT files (consecutive or overlapping) Default Name Type File Name IGFn.DAT input 1 * IGFDAT=CALMET0.DAT * *END* Subgroup (f) Other file names Default Name Type File Name DIAG.DAT input * DIADAT= * 98 PROG.DAT input * PRGDAT= * TEST.PRT output * TSTPRT= * 101 TEST.OUT output * TSTOUT= * 102 TEST.KIN output * TSTKIN= * 103 TEST.FRD output * TSTFRD= * 104 TEST.SLP output * TSTSLP= * 105 DCST.GRD output * DCSTGD= * NOTES: (1) File/path names can be up to 70 characters in length 109 (2) Subgroups (a) and (f) must have ONE 'END' (surrounded by 110 delimiters) at the end of the group 111 (3) Subgroups (b) through (e) are included ONLY if the corresponding 112 number of files (NUSTA, NOWSTA, NM3D, NIGF) is not 0, and each must have 113 an 'END' (surround by delimiters) at the end of EACH LINE !END! INPUT GROUP: 1 -- General run control parameters Starting date: Year (IBYR) -- No default! IBYR = 2016! 124 Month (IBMO) -- No default! IBMO = 5! 125 Day (IBDY) -- No default! IBDY = 16! 126 Starting time: Hour (IBHR) -- No default! IBHR = 0! 127 Second (IBSEC) -- No default! IBSEC = 0! Ending date: Year (IEYR) -- No default! IEYR = 2016! 130 Month (IEMO) -- No default! IEMO = 11! 131 Day (IEDY) -- No default! IEDY = 1! 132 Ending time: Hour (IEHR) -- No default! IEHR = 0! 133 Second (IESEC) -- No default! IESEC = 0! UTC time zone (ABTZ) -- No default! ABTZ= UTC-0300! 136 (character*8) 137 PST = UTC-0800, MST = UTC-0700, GMT = UTC CST = UTC-0600, EST = UTC-0500

4 Length of modeling time-step (seconds) 141 Must divide evenly into 3600 (1 hour) 142 (NSECDT) Default:3600! NSECDT = 3600! 143 Units: seconds Run type (IRTYPE) -- Default: 1! IRTYPE= 1! = Computes wind fields only = Computes wind fields and micrometeorological variables 149 (u*, w*, L, zi, etc.) 150 (IRTYPE must be 1 to run CALPUFF or CALGRID) Compute special data fields required 153 by CALGRID (i.e., 3-D fields of W wind 154 components and temperature) 155 in additional to regular Default: T! LCALGRD = T! 156 fields? (LCALGRD) 157 (LCALGRD must be T to run CALGRID) Flag to stop run after 160 SETUP phase (ITEST) Default: 2! ITEST= 2! 161 (Used to allow checking 162 of the model inputs, files, etc.) 163 ITEST = 1 - STOPS program after SETUP phase 164 ITEST = 2 - Continues with execution of 165 COMPUTATIONAL phase after SETUP Test options specified to see if 169 they conform to regulatory 170 values? (MREG) No Default! MREG = 0! = NO checks are made = Technical options must conform to USEPA guidance 174 IMIXH -1 Maul-Carson convective mixing height 175 over land; OCD mixing height overwater 176 ICOARE 0 OCD deltat method for overwater fluxes 177 THRESHL 0.0 Threshold buoyancy flux over land needed 178 to sustain convective mixing height growth 179 ISURFT > 0 in OBS mode (pick one representative station) in NOOBS mode (itprog=2) (average all 181 surface prognostic temperatures to get 182 a single representative sf. temp) 183 IUPT > 0 in OBS mode (pick one representative station) in NOOBS mode (ITPROG>0) (average all surface 185 prognostic temperatures to get a single 186 representative sf. temp) 187 IZICRLX 0 Do NOT use convective mixing height relaxation 188 to equilibrium value !END! INPUT GROUP: 2 -- Map Projection and Grid control parameters Projection for all (X,Y): Map projection 201 (PMAP) Default: UTM! PMAP = UTM! UTM : Universal Transverse Mercator 204 TTM : Tangential Transverse Mercator 205 LCC : Lambert Conformal Conic 206 PS : Polar Stereographic 207 EM : Equatorial Mercator

5 208 LAZA : Lambert Azimuthal Equal Area False Easting and Northing (km) at the projection origin 211 (Used only if PMAP= TTM, LCC, or LAZA) 212 (FEAST) Default=0.0! FEAST = 0.000! 213 (FNORTH) Default=0.0! FNORTH = 0.000! UTM zone (1 to 60) 216 (Used only if PMAP=UTM) 217 (IUTMZN) No Default! IUTMZN = 22! Hemisphere for UTM projection? 220 (Used only if PMAP=UTM) 221 (UTMHEM) Default: N! UTMHEM = N! 222 N : Northern hemisphere projection 223 S : Southern hemisphere projection Latitude and Longitude (decimal degrees) of projection origin 226 (Used only if PMAP= TTM, LCC, PS, EM, or LAZA) 227 (RLAT0) No Default! RLAT0 = 0N! 228 (RLON0) No Default! RLON0 = 0E! TTM : RLON0 identifies central (true N/S) meridian of projection 231 RLAT0 selected for convenience 232 LCC : RLON0 identifies central (true N/S) meridian of projection 233 RLAT0 selected for convenience 234 PS : RLON0 identifies central (grid N/S) meridian of projection 235 RLAT0 selected for convenience 236 EM : RLON0 identifies central meridian of projection 237 RLAT0 is REPLACED by 0.0N (Equator) 238 LAZA: RLON0 identifies longitude of tangent-point of mapping plane 239 RLAT0 identifies latitude of tangent-point of mapping plane Matching parallel(s) of latitude (decimal degrees) for projection 242 (Used only if PMAP= LCC or PS) 243 (XLAT1) No Default! XLAT1 = 0N! 244 (XLAT2) No Default! XLAT2 = 0N! LCC : Projection cone slices through Earth's surface at XLAT1 and XLAT2 247 PS : Projection plane slices through Earth at XLAT1 248 (XLAT2 is not used) Note: Latitudes and longitudes should be positive, and include a 252 letter N,S,E, or W indicating north or south latitude, and 253 east or west longitude. For example, N Latitude = 35.9N E Longitude = 118.7E Datum-region The Datum-Region for the coordinates is identified by a character 262 string. Many mapping products currently available use the model of the 263 Earth known as the World Geodetic System 1984 (WGS-84). Other local 264 models may be in use, and their selection in CALMET will make its output 265 consistent with local mapping products. The list of Datum-Regions with 266 official transformation parameters is provided by the National Imagery and 267 Mapping Agency (NIMA) NIMA Datum - Regions(Examples) WGS-84 WGS-84 Reference Ellipsoid and Geoid, Global coverage (WGS84) 272 NAS-C NORTH AMERICAN 1927 Clarke 1866 Spheroid, MEAN FOR CONUS (NAD27) 273 NAR-C NORTH AMERICAN 1983 GRS 80 Spheroid, MEAN FOR CONUS (NAD83) 274 NWS-84 NWS 6370KM Radius, Sphere 275 ESR-S ESRI REFERENCE 6371KM Radius, Sphere 276

6 277 Datum-region for output coordinates 278 (DATUM) Default: WGS-84! DATUM = WGS-84! Horizontal grid definition: Rectangular grid defined for projection PMAP, 285 with X the Easting and Y the Northing coordinate No. X grid cells (NX) No default! NX = 80! 288 No. Y grid cells (NY) No default! NY = 80! Grid spacing (DGRIDKM) No default! DGRIDKM = 0.25! 291 Units: km Reference grid coordinate of 294 SOUTHWEST corner of grid cell (1,1) X coordinate (XORIGKM) No default! XORIGKM = ! 297 Y coordinate (YORIGKM) No default! YORIGKM = ! 298 Units: km Vertical grid definition: No. of vertical layers (NZ) No default! NZ = 9! Cell face heights in arbitrary 307 vertical grid (ZFACE(NZ+1)) No defaults 308 Units: m 309! ZFACE = 0.,20.,40.,80.,160.,300.,600.,1000.,1500.,2200.! !END! INPUT GROUP: 3 -- Output Options DISK OUTPUT OPTION Save met. fields in an unformatted 323 output file? (LSAVE) Default: T! LSAVE = T! 324 (F = Do not save, T = Save) Type of unformatted output file: 327 (IFORMO) Default: 1! IFORMO = 1! = CALPUFF/CALGRID type file (CALMET.DAT) = MESOPUFF-II type file (PACOUT.DAT) LINE PRINTER OUTPUT OPTIONS: Print met. fields? (LPRINT) Default: F! LPRINT = F! 336 (F = Do not print, T = Print) 337 (NOTE: parameters below control which 338 met. variables are printed) Print interval 341 (IPRINF) in hours Default: 1! IPRINF = 1! 342 (Meteorological fields are printed 343 every 1 hours)

7 346 Specify which layers of U, V wind component 347 to print (IUVOUT(NZ)) -- NOTE: NZ values must be entered 348 (0=Do not print, 1=Print) 349 (used only if LPRINT=T) Defaults: NZ*0 350! IUVOUT = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0! Specify which levels of the W wind component to print 355 (NOTE: W defined at TOP cell face values) 356 (IWOUT(NZ)) -- NOTE: NZ values must be entered 357 (0=Do not print, 1=Print) 358 (used only if LPRINT=T & LCALGRD=T) Defaults: NZ*0 361! IWOUT = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0! Specify which levels of the 3-D temperature field to print 365 (ITOUT(NZ)) -- NOTE: NZ values must be entered 366 (0=Do not print, 1=Print) 367 (used only if LPRINT=T & LCALGRD=T) Defaults: NZ*0 370! ITOUT = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0! Specify which meteorological fields 373 to print 374 (used only if LPRINT=T) Defaults: 0 (all variables) Variable Print? 379 (0 = do not print, = print) ! STABILITY = 0! - PGT stability class 384! USTAR = 0! - Friction velocity 385! MONIN = 0! - Monin-Obukhov length 386! MIXHT = 0! - Mixing height 387! WSTAR = 0! - Convective velocity scale 388! PRECIP = 0! - Precipitation rate 389! SENSHEAT = 0! - Sensible heat flux 390! CONVZI = 0! - Convective mixing ht Testing and debug print options for micrometeorological module Print input meteorological data and 396 internal variables (LDB) Default: F! LDB = F! 397 (F = Do not print, T = print) 398 (NOTE: this option produces large amounts of output) First time step for which debug data 401 are printed (NN1) Default: 1! NN1 = 1! Last time step for which debug data 404 are printed (NN2) Default: 1! NN2 = 2! Print distance to land 407 internal variables (LDBCST) Default: F! LDBCST = F! 408 (F = Do not print, T = print) 409 (Output in.grd file DCST.GRD, defined in input group 0) Testing and debug print options for wind field module 412 (all of the following print options control output to 413 wind field module's output files: TEST.PRT, TEST.OUT, 414 TEST.KIN, TEST.FRD, and TEST.SLP)

8 Control variable for writing the test/debug 417 wind fields to disk files (IOUTD) 418 (0=Do not write, 1=write) Default: 0! IOUTD = 0! Number of levels, starting at the surface, 421 to print (NZPRN2) Default: 1! NZPRN2 = 1! Print the INTERPOLATED wind components? 424 (IPR0) (0=no, 1=yes) Default: 0! IPR0 = 0! Print the TERRAIN ADJUSTED surface wind 427 components? 428 (IPR1) (0=no, 1=yes) Default: 0! IPR1 = 0! Print the SMOOTHED wind components and 431 the INITIAL DIVERGENCE fields? 432 (IPR2) (0=no, 1=yes) Default: 0! IPR2 = 0! Print the FINAL wind speed and direction 435 fields? 436 (IPR3) (0=no, 1=yes) Default: 0! IPR3 = 0! Print the FINAL DIVERGENCE fields? 439 (IPR4) (0=no, 1=yes) Default: 0! IPR4 = 0! Print the winds after KINEMATIC effects 442 are added? 443 (IPR5) (0=no, 1=yes) Default: 0! IPR5 = 0! Print the winds after the FROUDE NUMBER 446 adjustment is made? 447 (IPR6) (0=no, 1=yes) Default: 0! IPR6 = 0! Print the winds after SLOPE FLOWS 450 are added? 451 (IPR7) (0=no, 1=yes) Default: 0! IPR7 = 0! Print the FINAL wind field components? 454 (IPR8) (0=no, 1=yes) Default: 0! IPR8 = 0! !END! INPUT GROUP: 4 -- Meteorological data options NO OBSERVATION MODE (NOOBS) Default: 0! NOOBS = 2! = Use surface, overwater, and upper air stations = Use surface and overwater stations (no upper air observations) 467 Use MM4/MM5/3D for upper air data = No surface, overwater, or upper air observations 469 Use MM4/MM5/3D for surface, overwater, and upper air data NUMBER OF SURFACE & PRECIP. METEOROLOGICAL STATIONS Number of surface stations (NSSTA) No default! NSSTA = 0! Number of precipitation stations 476 (NPSTA=-1: flag for use of MM5/3D precip data) 477 (NPSTA) No default! NPSTA = -1! CLOUD DATA OPTIONS 480 Output option - output a CLOUD.DAT file (yes or no) 481 0=no, 1=yes 482 (ICLDOUT) Default:999! ICLDOUT = 0! 483

9 484 Method to compute cloud fields: 485 (MCLOUD) Default: 999! MCLOUD = 4! 486 MCLOUD = 1 - Clouds data generated from surface observations 487 MCLOUD = 2 - Gridded CLOUD.DAT read from CLOUD.DAT file (no output 488 is possible since already exist) 489 MCLOUD = 3 - Gridded cloud cover from Prognostic Rel. Humidity 490 at 850mb (Teixera) 491 MCLOUD = 4 - Gridded cloud cover from Prognostic Rel. Humidity 492 at all levels (MM5toGrads algorithm) FILE FORMATS Surface meteorological data file format 497 (IFORMS) Default: 2! IFORMS = 2! 498 (1 = unformatted (e.g., SMERGE output)) 499 (2 = formatted (free-formatted user input)) Precipitation data file format 502 (IFORMP) Default: 2! IFORMP = 2! 503 (1 = unformatted (e.g., PMERGE output)) 504 (2 = formatted (free-formatted user input)) Cloud data file format 507 (IFORMC) Default: 2! IFORMC = 2! 508 (1 = unformatted - CALMET unformatted output) 509 (2 = formatted - free-formatted CALMET output or user input) !END! INPUT GROUP: 5 -- Wind Field Options and Parameters WIND FIELD MODEL OPTIONS 521 Model selection variable (IWFCOD) Default: 1! IWFCOD = 1! = Objective analysis only = Diagnostic wind module Compute Froude number adjustment 526 effects? (IFRADJ) Default: 1! IFRADJ = 1! 527 (0 = NO, 1 = YES) Compute kinematic effects? (IKINE) Default: 0! IKINE = 0! 530 (0 = NO, 1 = YES) Use O'Brien procedure for adjustment 533 of the vertical velocity? (IOBR) Default: 0! IOBR = 0! 534 (0 = NO, 1 = YES) Compute slope flow effects? (ISLOPE) Default: 1! ISLOPE = 1! 537 (0 = NO, 1 = YES) Extrapolate surface wind observations 540 to upper layers? (IEXTRP) Default: -4! IEXTRP = 1! 541 (1 = no extrapolation is done, = power law extrapolation used, = user input multiplicative factors 544 for layers 2 - NZ used (see FEXTRP array) = similarity theory used 546-1, -2, -3, -4 = same as above except layer 1 data 547 at upper air stations are ignored Extrapolate surface winds even 550 if calm? (ICALM) Default: 0! ICALM = 0! 551 (0 = NO, 1 = YES) 552

10 553 Layer-dependent biases modifying the weights of 554 surface and upper air stations (BIAS(NZ)) 555-1<=BIAS<=1 556 Negative BIAS reduces the weight of upper air stations 557 (e.g. BIAS=-0.1 reduces the weight of upper air stations 558 by 10%; BIAS= -1, reduces their weight by 100 %) 559 Positive BIAS reduces the weight of surface stations 560 (e.g. BIAS= 0.2 reduces the weight of surface stations 561 by 20%; BIAS=1 reduces their weight by 100%) 562 Zero BIAS leaves weights unchanged (1/R**2 interpolation) 563 Default: NZ*0 564! BIAS = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0! Minimum distance from nearest upper air station 567 to surface station for which extrapolation 568 of surface winds at surface station will be allowed 569 (RMIN2: Set to -1 for IEXTRP = 4 or other situations 570 where all surface stations should be extrapolated) 571 Default: 4.! RMIN2 = 4.0! Use gridded prognostic wind field model 574 output fields as input to the diagnostic 575 wind field model (IPROG) Default: 0! IPROG = 14! 576 (0 = No, [IWFCOD = 0 or 1] = Yes, use CSUMM prog. winds as Step 1 field, [IWFCOD = 0] = Yes, use CSUMM prog. winds as initial guess field [IWFCOD = 1] = Yes, use winds from MM4.DAT file as Step 1 field [IWFCOD = 0] = Yes, use winds from MM4.DAT file as initial guess field [IWFCOD = 1] = Yes, use winds from MM4.DAT file as observations [IWFCOD = 1] = Yes, use winds from MM5/3D.DAT file as Step 1 field [IWFCOD = 0] = Yes, use winds from MM5/3D.DAT file as initial guess field [IWFCOD = 1] = Yes, use winds from MM5/3D.DAT file as observations [IWFCOD = 1] Timestep (seconds) of the prognostic 587 model input data (ISTEPPGS) Default: 3600! ISTEPPGS = 3600! Use coarse CALMET fields as initial guess fields (IGFMET) 590 (overwrites IGF based on prognostic wind fields if any) 591 Default: 0! IGFMET = 0! RADIUS OF INFLUENCE PARAMETERS Use varying radius of influence Default: F! LVARY = F! 596 (if no stations are found within RMAX1,RMAX2, 597 or RMAX3, then the closest station will be used) Maximum radius of influence over land 600 in the surface layer (RMAX1) No default! RMAX1 = 20.! 601 Units: km 602 Maximum radius of influence over land 603 aloft (RMAX2) No default! RMAX2 = 20.! 604 Units: km 605 Maximum radius of influence over water 606 (RMAX3) No default! RMAX3 = 20.! 607 Units: km OTHER WIND FIELD INPUT PARAMETERS Minimum radius of influence used in 613 the wind field interpolation (RMIN) Default: 0.1! RMIN = 1.! 614 Units: km 615 Radius of influence of terrain 616 features (TERRAD) No default! TERRAD = 20.! Units: km 619 Relative weighting of the first 620 guess field and observations in the

11 621 SURFACE layer (R1) No default! R1 = 20.! 622 (R1 is the distance from an Units: km 623 observational station at which the 624 observation and first guess field are 625 equally weighted) Relative weighting of the first 628 guess field and observations in the 629 layers ALOFT (R2) No default! R2 = 20.! 630 (R2 is applied in the upper layers Units: km 631 in the same manner as R1 is used in 632 the surface layer) Relative weighting parameter of the 635 prognostic wind field data (RPROG) No default! RPROG = 0.! 636 (Used only if IPROG = 1) Units: km Maximum acceptable divergence in the 640 divergence minimization procedure 641 (DIVLIM) Default: 5.E-6! DIVLIM= 5.0E-06! Maximum number of iterations in the 644 divergence min. procedure (NITER) Default: 50! NITER = 50! Number of passes in the smoothing 647 procedure (NSMTH(NZ)) 648 NOTE: NZ values must be entered 649 Default: 2,(mxnz-1)*4! NSMTH = 650 2, 4, 4, 4, 4, 4, 4, 4, 4, 4! Maximum number of stations used in 653 each layer for the interpolation of 654 data to a grid point (NINTR2(NZ)) 655 NOTE: NZ values must be entered Default: 99.! NINTR2 = , 99, 99, 99, 99, 99, 99, 99, 99, 99! Critical Froude number (CRITFN) Default: 1.0! CRITFN = 1.! Empirical factor controlling the 661 influence of kinematic effects 662 (ALPHA) Default: 0.1! ALPHA = 0.1! Multiplicative scaling factor for 665 extrapolation of surface observations 666 to upper layers (FEXTR2(NZ)) Default: NZ* ! FEXTR2 = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.! 668 (Used only if IEXTRP = 3 or -3) BARRIER INFORMATION Number of barriers to interpolation 674 of the wind fields (NBAR) Default: 0! NBAR = 0! Level (1 to NZ) up to which barriers 677 apply (KBAR) Default: NZ! KBAR = 12! THE FOLLOWING 4 VARIABLES ARE INCLUDED 680 ONLY IF NBAR > NOTE: NBAR values must be entered No defaults 682 for each variable Units: km X coordinate of BEGINNING 685 of each barrier (XBBAR(NBAR))! XBBAR = 0.! 686 Y coordinate of BEGINNING 687 of each barrier (YBBAR(NBAR))! YBBAR = 0.! X coordinate of ENDING

12 690 of each barrier (XEBAR(NBAR))! XEBAR = 0.! 691 Y coordinate of ENDING 692 of each barrier (YEBAR(NBAR))! YEBAR = 0.! DIAGNOSTIC MODULE DATA INPUT OPTIONS Surface temperature (IDIOPT1) Default: 0! IDIOPT1 = 0! = Compute internally from 699 hourly surface observations or prognostic fields = Read preprocessed values from 701 a data file (DIAG.DAT) Surface met. station to use for 704 the surface temperature (ISURFT) Default: -1! ISURFT = -1! 705 (Must be a value from 1 to NSSTA 706 or -1 to use 2-D spatially varying 707 surface temperatures). 708 or -2 to use a domain-average prognostic 709 lapse rate (only with ITPROG=2) 710 (Used only if IDIOPT1 = 0) Temperature lapse rate used in the Default: 0! IDIOPT2 = 0! 714 computation of terrain-induced 715 circulations (IDIOPT2) = Compute internally from (at least) twice-daily 717 upper air observations or prognostic fields = Read hourly preprocessed values 719 from a data file (DIAG.DAT) Upper air station to use for 722 the domain-scale lapse rate (IUPT) Default: -1! IUPT = -1! 723 (Must be a value from 1 to NUSTA 724 or -1 to use 2-D spatially varying lapse rate) 725 or -2 to use a domain-average prognostic 726 lapse rate (only with ITPROG>0) 727 (Used only if IDIOPT2 = 0) Depth through which the domain-scale 731 lapse rate is computed (ZUPT) Default: 200.! ZUPT = 200.! 732 (Used only if IDIOPT2 = 0) Units: meters Initial Guess Field Winds 736 (IDIOPT3) Default: 0! IDIOPT3 = 0! = Compute internally from 738 observations or prognostic wind fields = Read hourly preprocessed domain-average wind values 740 from a data file (DIAG.DAT) Upper air station to use for 743 the initial guess winds (IUPWND) Default: -1! IUPWND = -1! 744 (Must be a value from -1 to NUSTA, with indicating 3-D initial guess fields, 746 and IUPWND>1 domain-scaled (i.e. constant) IGF 747 (Used only if IDIOPT3 = 0 and noobs=0) Bottom and top of layer through 751 which the domain-scale winds 752 are computed 753 (ZUPWND(1), ZUPWND(2)) Defaults: 1., 1000.! ZUPWND= 1., 1000.! 754 (Used only if IDIOPT3 = 0, NOOBS>0 and IUPWND>0) Units: meters Observed surface wind components 758 for wind field module (IDIOPT4) Default: 0! IDIOPT4 = 0!

13 759 0 = Read WS, WD from a surface 760 data file (SURF.DAT) = Read hourly preprocessed U, V from 762 a data file (DIAG.DAT) Observed upper air wind components 765 for wind field module (IDIOPT5) Default: 0! IDIOPT5 = 0! = Read WS, WD from an upper 767 air data file (UP1.DAT, UP2.DAT, etc.) = Read hourly preprocessed U, V from 769 a data file (DIAG.DAT) LAKE BREEZE INFORMATION Use Lake Breeze Module (LLBREZE) 774 Default: F! LLBREZE = F! Number of lake breeze regions (NBOX)! NBOX = 0! X Grid line 1 defining the region of interest 779! XG1 = 0.! 780 X Grid line 2 defining the region of interest 781! XG2 = 0.! 782 Y Grid line 1 defining the region of interest 783! YG1 = 0.! 784 Y Grid line 2 defining the region of interest 785! YG2 = 0.! X Point defining the coastline (Straight line) 788 (XBCST) (KM) Default: none! XBCST = 0.! Y Point defining the coastline (Straight line) 791 (YBCST) (KM) Default: none! YBCST = 0.! X Point defining the coastline (Straight line) 794 (XECST) (KM) Default: none! XECST = 0.! Y Point defining the coastline (Straight line) 797 (YECST) (KM) Default: none! YECST = 0.! Number of stations in the region Default: none! NLB = 0! 801 (Surface stations + upper air stations) Station ID's in the region (METBXID(NLB)) 804 (Surface stations first, then upper air stations) 805! METBXID = 0! !END! INPUT GROUP: 6 -- Mixing Height, Temperature and Precipitation Parameters EMPIRICAL MIXING HEIGHT CONSTANTS Neutral, mechanical equation 818 (CONSTB) Default: 1.41! CONSTB = 1.41! 819 Convective mixing ht. equation 820 (CONSTE) Default: 0.15! CONSTE = 0.15! 821 Stable mixing ht. equation 822 (CONSTN) Default: 2400.! CONSTN = 2400.! 823 Overwater mixing ht. equation 824 (CONSTW) Default: 0.16! CONSTW = 0.16! 825 Absolute value of Coriolis 826 parameter (FCORIOL) Default: 1.E-4! FCORIOL = 1.0E-04! 827 Units: (1/s)

14 SPATIAL AVERAGING OF MIXING HEIGHTS Conduct spatial averaging 832 (IAVEZI) (0=no, 1=yes) Default: 1! IAVEZI = 1! Max. search radius in averaging 835 process (MNMDAV) Default: 1! MNMDAV = 10! 836 Units: Grid 837 cells 838 Half-angle of upwind looking cone 839 for averaging (HAFANG) Default: 30.! HAFANG = 30.! 840 Units: deg. 841 Layer of winds used in upwind 842 averaging (ILEVZI) Default: 1! ILEVZI = 1! 843 (must be between 1 and NZ) CONVECTIVE MIXING HEIGHT OPTIONS: 847 Method to compute the convective 848 mixing height(imihxh) Default: 1! IMIXH = 1! 849 1: Maul-Carson for land and water cells 850-1: Maul-Carson for land cells only OCD mixing height overwater 852 2: Batchvarova and Gryning for land and water cells 853-2: Batchvarova and Gryning for land cells only 854 OCD mixing height overwater Threshold buoyancy flux required to 857 sustain convective mixing height growth 858 overland (THRESHL) Default: 0.0! THRESHL = 0.0! 859 (expressed as a heat flux units: W/m3 860 per meter of boundary layer) Threshold buoyancy flux required to 864 sustain convective mixing height growth 865 overwater (THRESHW) Default: 0.05! THRESHW = 0.05! 866 (expressed as a heat flux units: W/m3 867 per meter of boundary layer) Flag to allow relaxation of convective mixing height 871 to equilibrium value when 0<QH<THRESHL (overland) 872 or 0<QH<THRESHW (overwater) 873 (IZICRLX) Default: 1! IZICRLX = 1! : do NOT use convective mixing height relaxation 875 to equilibrium value (treatment identical to CALMET v5.8) : use convective mixing height relaxation 877 to equilibrium value Relaxation time of convective mixing height to 880 equilibrium value when 0<QH<THRESHL (overland) 881 or 0<QH<THRESHW (overwater) 882 (Used only if IZICRLX = 1 and TZICRLX must be >= 1.) 883 (TZICRLX) Default: 800.! TZICRLX = 800.! 884 Units: seconds Option for overwater lapse rates used 887 in convective mixing height growth 888 (ITWPROG) Default: 0! ITWPROG = 0! : use SEA.DAT lapse rates and deltat (or assume neutral 890 conditions if missing) : use prognostic lapse rates (only if IPROG>2) 892 and SEA.DAT deltat (or neutral if missing) : use prognostic lapse rates and prognostic delta T 894 (only if iprog>12 and 3D.DAT version# 2.0 or higher) Land Use category ocean in 3D.DAT datasets

15 897 (ILUOC3D) Default: 16! ILUOC3D = 16! 898 Note: if 3D.DAT from MM5 version 3.0, iluoc3d = if MM4.DAT, typically iluoc3d = OTHER MIXING HEIGHT VARIABLES Minimum potential temperature lapse 905 rate in the stable layer above the 906 current convective mixing ht. Default: 0.001! DPTMIN = 0.001! 907 (DPTMIN) Units: deg. K/m 908 Depth of layer above current conv. 909 mixing height through which lapse Default: 200.! DZZI = 200.! 910 rate is computed (DZZI) Units: meters Minimum overland mixing height Default: 50.! ZIMIN = 50.! 913 (ZIMIN) Units: meters 914 Maximum overland mixing height Default: 3000.! ZIMAX = 2500.! 915 (ZIMAX) Units: meters 916 Minimum overwater mixing height Default: 50.! ZIMINW = 50.! 917 (ZIMINW) -- (Not used if observed Units: meters 918 overwater mixing hts. are used) 919 Maximum overwater mixing height Default: 3000.! ZIMAXW = 2500.! 920 (ZIMAXW) -- (Not used if observed Units: meters 921 overwater mixing hts. are used) OVERWATER SURFACE FLUXES METHOD and PARAMETERS 925 (ICOARE) Default: 10! ICOARE = 10! 926 0: original deltat method (OCD) : COARE with no wave parameterization (jwave=0, Charnock) : COARE with wave option jwave=1 (Oost et al.) 929 and default wave properties : COARE with wave option jwave=1 (Oost et al.) 931 and observed wave properties (must be in SEA.DAT files) : COARE with wave option 2 (Taylor and Yelland) 933 and default wave properties : COARE with wave option 2 (Taylor and Yelland) 935 and observed wave properties (must be in SEA.DAT files) Note: When ICOARE=0, similarity wind profile stability PSI functions 939 based on Van Ulden and Holtslag (1985) are substituted for 940 later formulations used with the COARE module, and temperatures 941 used for surface layer parameters are obtained from either the 942 nearest surface station temperature or prognostic model 2D 943 temperatures (if ITPROG=2) Coastal/Shallow water length scale (DSHELF) 948 (for modified z0 in shallow water) 949 ( COARE fluxes only) 950 Default : 0.! DSHELF = 0.! 951 units: km COARE warm layer computation (IWARM)! IWARM = 0! 954 1: on - 0: off (must be off if SST measured with 955 IR radiometer) Default: COARE cool skin layer computation (ICOOL)! ICOOL = 0! 958 1: on - 0: off (must be off if SST measured with 959 IR radiometer) Default: RELATIVE HUMIDITY PARAMETERS D relative humidity from observations or 964 from prognostic data? (IRHPROG) Default:0! IRHPROG = 1! 965

16 966 0 = Use RH from SURF.DAT file 967 (only if NOOBS = 0,1) = Use prognostic RH 969 (only if NOOBS = 0,1,2) TEMPERATURE PARAMETERS D temperature from observations or 974 from prognostic data? (ITPROG) Default:0! ITPROG = 2! = Use Surface and upper air stations 977 (only if NOOBS = 0) = Use Surface stations (no upper air observations) 979 Use MM5/3D for upper air data 980 (only if NOOBS = 0,1) = No surface or upper air observations 982 Use MM5/3D for surface and upper air data 983 (only if NOOBS = 0,1,2) Interpolation type 986 (1 = 1/R ; 2 = 1/R**2) Default:1! IRAD = 1! Radius of influence for temperature 989 interpolation (TRADKM) Default: 500.! TRADKM = 500.! 990 Units: km Maximum Number of stations to include 993 in temperature interpolation (NUMTS) Default: 5! NUMTS = 5! Conduct spatial averaging of temp- 996 eratures (IAVET) (0=no, 1=yes) Default: 1! IAVET = 1! 997 (will use mixing ht MNMDAV,HAFANG 998 so make sure they are correct) Default temperature gradient Default: ! TGDEFB = ! 1001 below the mixing height over Units: K/m 1002 water (TGDEFB) Default temperature gradient Default: ! TGDEFA = ! 1005 above the mixing height over Units: K/m 1006 water (TGDEFA) Beginning (JWAT1) and ending (JWAT2) 1009 land use categories for temperature! JWAT1 = 999! 1010 interpolation over water -- Make! JWAT2 = 999! 1011 bigger than largest land use to disable PRECIP INTERPOLATION PARAMETERS Method of interpolation (NFLAGP) Default: 2! NFLAGP = 2! 1016 (1=1/R,2=1/R**2,3=EXP/R**2) 1017 Radius of Influence (SIGMAP) Default: 100.0! SIGMAP = 100.! 1018 (0.0 => use half dist. btwn Units: km 1019 nearest stns w & w/out 1020 precip when NFLAGP = 3) 1021 Minimum Precip. Rate Cutoff (CUTP) Default: 0.01! CUTP = 0.01! 1022 (values < CUTP = 0.0 mm/hr) Units: mm/hr 1023!END! INPUT GROUP: 7 -- Surface meteorological station parameters SURFACE STATION VARIABLES 1032 (One record per station -- 0 records in all)

17 Name ID X coord. Y coord. Time Anem (km) (km) zone Ht.(m) Four character string for station name 1042 (MUST START IN COLUMN 9) Six digit integer for station ID !END! INPUT GROUP: 8 -- Upper air meteorological station parameters UPPER AIR STATION VARIABLES 1056 (One record per station -- 3 records in all) Name ID X coord. Y coord. Time zone 1060 (km) (km) Four character string for station name 1065 (MUST START IN COLUMN 9) Five digit integer for station ID !END! INPUT GROUP: 9 -- Precipitation station parameters PRECIPITATION STATION VARIABLES 1079 (One record per station records in all) 1080 (NOT INCLUDED IF NPSTA = 0) Name Station X coord. Y coord Code (km) (km) Four character string for station name 1091 (MUST START IN COLUMN 9) Six digit station code composed of state 1095 code (first 2 digits) and station ID (last digits) !END! 1099

18 1 CALPUFF.INP 7.0 Groups 0f,0g added; new emission scaling 2 CALPUFF Demonstration Run 3 (Not intended as a guide for configuring options) Run title (3 lines) MODEL: Version TNG 7 SOFTWARE: CALApps v (Beta) - September 22, CALPUFF MODEL CONTROL FILE INPUT GROUP: 0 -- Input and Output File Names Default Name Type File Name CALMET.DAT input! METDAT =! 21 or 22 ISCMET.DAT input * ISCDAT = * 23 or 24 PLMMET.DAT input * PLMDAT = * 25 or 26 PROFILE.DAT input * PRFDAT = * 27 SURFACE.DAT input * SFCDAT = * 28 RESTARTB.DAT input * RSTARTB= * CALPUFF.LST output! PUFLST =CALPUFF.LST! 31 CONC.DAT output! CONDAT =CALPUFF.CON! 32 DFLX.DAT output * DFDAT = * 33 WFLX.DAT output * WFDAT = * VISB.DAT output! VISDAT =CALPUFF.VIS! 36 TK2D.DAT output * T2DDAT = * 37 RHO2D.DAT output * RHODAT = * 38 RESTARTE.DAT output * RSTARTE= * Other Files OZONE.DAT input * OZDAT = * 43 VD.DAT input * VDDAT = * 44 CHEM.DAT input * CHEMDAT= * 45 AUX input! AUXEXT =aux! 46 (Extension added to METDAT filename(s) for files 47 with auxiliary 2D and 3D data) 48 H2O2.DAT input * H2O2DAT= * 49 NH3Z.DAT input * NH3ZDAT= * 50 HILL.DAT input * HILDAT= * 51 HILLRCT.DAT input * RCTDAT= * 52 COASTLN.DAT input * CSTDAT= * 53 FLUXBDY.DAT input * BDYDAT= * 54 BCON.DAT input * BCNDAT= * 55 DEBUG.DAT output * DEBUG = * 56 MASSFLX.DAT output * FLXDAT= * 57 MASSBAL.DAT output * BALDAT= * 58 FOG.DAT output * FOGDAT= * 59 RISE.DAT output * RISDAT= * 60 PFTRAK.DAT output * TRKDAT= * All file names will be converted to lower case if LCFILES = T 63 Otherwise, if LCFILES = F, file names will be converted to UPPER CASE 64 T = lower case! LCFILES = F! 65 F = UPPER CASE 66 NOTE: (1) file/path names can be up to 132 characters in length Provision for multiple CALMET Domains and files

19 Number of CALMET.DAT Domains (NMETDOM) 73 Default: 1! NMETDOM = 1! Number of CALMET.DAT files (NMETDAT) 76 (Total for ALL Domains) 77 Default: 1! NMETDAT = 12! Variable point/area/volume/flare emissions input files Number of POINT source files (PTEMARB.DAT) 83 with time-varying data (NPTDAT) 84 Default: 0! NPTDAT = 0! Number of BUOYANT AREA source files (BAEMARB.DAT) 87 with time-varying data (NARDAT) 88 Default: 0! NARDAT = 0! Number of VOLUME source files (VOLEMARB.DAT) 91 with time-varying data (NVOLDAT) 92 Default: 0! NVOLDAT = 0! Number of FLARE source files (FLEMARB.DAT) 95 with time-varying data (NFLDAT) 96 Default: 0! NFLDAT = 0! Number of ROAD source files (RDEMARB.DAT) 99 with time-varying data (NRDDAT) 100 Default: 0! NRDDAT = 0! Number of BUOYANT LINE source files (LNEMARB.DAT) 103 with time-varying data (NLNDAT) 104 Default: 0! NLNDAT = 0! 105 Note: Only 1 BUOYANT LINE source file is allowed !END! Subgroup (0a) Provide a name for each CALMET domain if NMETDOM > Enter NMETDOM lines. 117 a,b 118 Default Name Domain Name none * DOMAIN1= * *END* 121 none * DOMAIN2= * *END* 122 none * DOMAIN3= * *END* The following CALMET.DAT filenames are processed in sequence 126 if NMETDAT > Enter NMETDAT lines, 1 line for each file name a,c,d 131 Default Name Type File Name none input! METDAT1= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_WINTER_14.DAT!!END! 134 none input! METDAT2= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SPRING_14.DAT!!END! 135 none input! METDAT3= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SUMMER_14.DAT!!END!

20 136 none input! METDAT4= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_AUTUMN_14.DAT!!END! 137 none input! METDAT5= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_WINTER_15.DAT!!END! 138 none input! METDAT6= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SPRING_15.DAT!!END! 139 none input! METDAT7= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SUMMER_15.DAT!!END! 140 none input! METDAT8= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_AUTUMN_15.DAT!!END! 141 none input! METDAT9= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_WINTER_16.DAT!!END! 142 none input! METDAT10= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SPRING_16.DAT!!END! 143 none input! METDAT11= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_SUMMER_16.DAT!!END! 144 none input! METDAT12= C:\Users\User\Desktop\CALPUFF\CALMET Outputs\CALMET_AUTUMN_16.DAT!!END! a 148 The name for each CALMET domain and each CALMET.DAT file is treated 149 as a separate input subgroup and therefore must end with an input 150 group terminator. 151 b 152 Use DOMAIN1= to assign the name for the outermost CALMET domain. 153 Use DOMAIN2= to assign the name for the next inner CALMET domain. 154 Use DOMAIN3= to assign the name for the next inner CALMET domain, etc When inner domains with equal resolution (grid-cell size) 157 overlap, the data from the FIRST such domain in the list will 158 be used if all other criteria for choosing the controlling 159 grid domain are inconclusive c 162 Use METDAT1= to assign the file names for the outermost CALMET domain. 163 Use METDAT2= to assign the file names for the next inner CALMET domain. 164 Use METDAT3= to assign the file names for the next inner CALMET domain, etc. 165 d 166 The filenames for each domain must be provided in sequential order Subgroup (0b) PTEMARB.DAT files POINT Source File Names 174 The following PTEMARB.DAT filenames are processed if NPTDAT>0 175 A total of NPTDAT lines is expected with one file name assigned per line 176 Each line is treated as an input group and must terminate with END 177 (surrounded by delimiters) 178 (Each file contains emissions parameters for the entire period modeled 179 for 1 or more sources) Default Name Type File Name none input * PTDAT= * *END* Subgroup (0c) BAEMARB.DAT files BUOYANT AREA Source File Names 191 The following BAEMARB.DAT filenames are processed if NARDAT>0 192 A total of NARDAT lines is expected with one file name assigned per line 193 Each line is treated as an input group and must terminate with END 194 (surrounded by delimiters) 195 (Each file contains emissions parameters for the entire period modeled

21 196 for 1 or more sources) Default Name Type File Name none input * ARDAT= * *END* Subgroup (0d) VOLEMARB.DAT files VOLUME Source File Names 208 The following VOLEMARB.DAT filenames are processed if NVOLDAT>0 209 A total of NVOLDAT lines is expected with one file name assigned per line 210 Each line is treated as an input group and must terminate with END 211 (surrounded by delimiters) 212 (Each file contains emissions parameters for the entire period modeled 213 for 1 or more sources) Default Name Type File Name none input * VOLDAT= * *END* Subgroup (0e) FLEMARB.DAT files FLARE Source File Names 225 The following FLEMARB.DAT filenames are processed if NFLDAT>0 226 A total of NFLDAT lines is expected with one file name assigned per line 227 Each line is treated as an input group and must terminate with END 228 (surrounded by delimiters) 229 (Each file contains emissions parameters for the entire period modeled 230 for 1 or more sources) Default Name Type File Name none input * FLDAT= * *END* Subgroup (0f) RDEMARB.DAT files ROAD Source File Names 242 The following RDEMARB.DAT filenames are processed if NRDDAT>0 243 A total of NRDDAT lines is expected with one file name assigned per line 244 Each line is treated as an input group and must terminate with END 245 (surrounded by delimiters) 246 (Each file contains emissions parameters for the entire period modeled 247 for 1 or more sources) Default Name Type File Name none input * RDDAT= * *END* Subgroup (0g) LNEMARB.DAT file BUOYANT LINE Source File Name (not more than 1) 259 The following LNEMARB.DAT filename is processed if NLNDAT>0 260 The assignment is treated as an input group and must terminate with END 261 (surrounded by delimiters) Default Name Type File Name

22 265 LNEMARB.DAT input * LNDAT= * *END* INPUT GROUP: 1 -- General run control parameters Option to run all periods found 274 in the met. file (METRUN) Default: 0! METRUN = 0! METRUN = 0 - Run period explicitly defined below 277 METRUN = 1 - Run all periods in met. file Starting date: Year (IBYR) -- No default! IBYR = 2014! 280 Month (IBMO) -- No default! IBMO = 1! 281 Day (IBDY) -- No default! IBDY = 1! 282 Starting time: Hour (IBHR) -- No default! IBHR = 0! 283 Minute (IBMIN) -- No default! IBMIN = 0! 284 Second (IBSEC) -- No default! IBSEC = 0! Ending date: Year (IEYR) -- No default! IEYR = 2017! 287 Month (IEMO) -- No default! IEMO = 1! 288 Day (IEDY) -- No default! IEDY = 1! 289 Ending time: Hour (IEHR) -- No default! IEHR = 0! 290 Minute (IEMIN) -- No default! IEMIN = 0! 291 Second (IESEC) -- No default! IESEC = 0! (These are only used if METRUN = 0) Base time zone: (ABTZ) -- No default! ABTZ= UTC-0300! 296 (character*8) 297 The modeling domain may span multiple time zones. ABTZ defines the 298 base time zone used for the entire simulation. This must match the 299 base time zone of the meteorological data. 300 Examples: 301 Greenwich Mean Time (GMT) = UTC EST = UTC CST = UTC MST = UTC PST = UTC Los Angeles, USA = UTC New York, USA = UTC Santiago, Chile = UTC UK = UTC Western Europe = UTC Rome, Italy = UTC Cape Town, S.Africa = UTC Sydney, Australia = UTC Length of modeling time-step (seconds) 316 Equal to update period in the primary 317 meteorological data files, or an 318 integer fraction of it (1/2, 1/3...) 319 Must be no larger than 1 hour 320 (NSECDT) Default:3600! NSECDT = 3600! 321 Units: seconds Number of chemical species (NSPEC) 324 Default: 5! NSPEC = 1! Number of chemical species 327 to be emitted (NSE) Default: 3! NSE = 1! Flag to stop run after 330 SETUP phase (ITEST) Default: 2! ITEST = 2! 331 (Used to allow checking 332 of the model inputs, files, etc.) 333 ITEST = 1 - STOPS program after SETUP phase

23 334 ITEST = 2 - Continues with execution of program 335 after SETUP Restart Configuration: Control flag (MRESTART) Default: 0! MRESTART = 0! = Do not read or write a restart file = Read a restart file at the beginning of 343 the run = Write a restart file during run = Read a restart file at beginning of run 346 and write a restart file during run Number of periods in Restart 349 output cycle (NRESPD) Default: 0! NRESPD = 0! = File written only at last period 352 >0 = File updated every NRESPD periods Meteorological Data Format (METFM) 355 Default: 1! METFM = 1! METFM = 1 - CALMET binary file (CALMET.MET) 358 METFM = 2 - ISC ASCII file (ISCMET.MET) 359 METFM = 3 - AUSPLUME ASCII file (PLMMET.MET) 360 METFM = 4 - CTDM plus tower file (PROFILE.DAT) and 361 surface parameters file (SURFACE.DAT) 362 METFM = 5 - AERMET tower file (PROFILE.DAT) and 363 surface parameters file (SURFACE.DAT) Meteorological Profile Data Format (MPRFFM) 366 (used only for METFM = 1, 2, 3) 367 Default: 1! MPRFFM = 1! MPRFFM = 1 - CTDM plus tower file (PROFILE.DAT) 370 MPRFFM = 2 - AERMET tower file (PROFILE.DAT) Sigma-y is adjusted by the factor (AVET/PGTIME)**0.2 to either 374 decrease it if the averaging time selected is less than the base 375 averaging time, or increase it if the averaging time is greater. 376 The base averaging time is denoted as PGTIME due to historical 377 reasons as this adjustment was originally applied to the PG sigma 378 option. It is now applied to all dispersion options. 379 The factor is applied to the ambient turbulence sigma-v (m/s) and 380 does not alter buoyancy enhancement or far-field Heffter growth Averaging Time (minutes) (AVET) 383 Default: 60.0! AVET = 60.! 384 Base Averaging Time (minutes) (PGTIME) 385 Default: 60.0! PGTIME = 60.! Output units for binary concentration and flux files 389 written in Dataset v2.2 or later formats 390 (IOUTU) Default: 1! IOUTU = 2! = mass - g/m3 (conc) or g/m2/s (dep) = odour - odour_units (conc) = radiation - Bq/m3 (conc) or Bq/m2/s (dep) !END! INPUT GROUP: 2 -- Technical options