SHARE
TWEET

Untitled

a guest Dec 16th, 2018 81 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. CALPUFF Demonstration Run
  2.  
  3.  
  4. ---------------- Run title (3 lines) ------------------------------------------
  5.  
  6.                     CALPUFF MODEL CONTROL FILE
  7.                     --------------------------
  8.  
  9. -------------------------------------------------------------------------------
  10.  
  11. INPUT GROUP: 0 -- Input and Output File Names
  12.  
  13. --------------
  14. Default Name  Type          File Name
  15. ------------  ----          ---------
  16. CALMET.DAT    input    ! METDAT =calmet.met   !
  17.     or
  18. ISCMET.DAT    input    * ISCDAT =             *
  19.     or
  20. PLMMET.DAT    input    * PLMDAT =             *
  21.     or
  22. PROFILE.DAT   input    * PRFDAT =             *
  23. SURFACE.DAT   input    * SFCDAT =             *
  24. RESTARTB.DAT  input    * RSTARTB=             *
  25. --------------------------------------------------------------------------------
  26. CALPUFF.LST   output   ! PUFLST =caca.LST  !
  27. CONC.DAT      output   ! CONDAT =caca.CON     !
  28. DFLX.DAT      output   ! DFDAT  =caca.DRY     !
  29. WFLX.DAT      output   ! WFDAT  =caca.WET     !
  30.  
  31. VISB.DAT      output   ! VISDAT =caca.VIS     !
  32. TK2D.DAT      output   * T2DDAT =             *
  33. RHO2D.DAT     output   * RHODAT =             *
  34. RESTARTE.DAT  output   ! RSTARTE= rstart.dat  !
  35. --------------------------------------------------------------------------------
  36. Emission Files
  37. --------------
  38. PTEMARB.DAT   input    ! PTDAT  =    PTEMARB.DAT          !
  39. VOLEMARB.DAT  input    * VOLDAT =             *
  40. BAEMARB.DAT   input    * ARDAT  =             *
  41. LNEMARB.DAT   input    * LNDAT  =             *
  42. --------------------------------------------------------------------------------
  43. Other Files
  44. -----------
  45. OZONE.DAT     input    ! OZDAT  =             !
  46. VD.DAT        input    * VDDAT  =             *
  47. CHEM.DAT      input    * CHEMDAT=             *
  48. H2O2.DAT      input    * H2O2DAT=             *
  49. HILL.DAT      input    * HILDAT=             *
  50. HILLRCT.DAT   input    * RCTDAT=             *
  51. COASTLN.DAT   input    * CSTDAT=             *
  52. FLUXBDY.DAT   input    * BDYDAT=             *
  53. BCON.DAT      input    * BCNDAT=             *
  54. DEBUG.DAT     output   * DEBUG =             *
  55. MASSFLX.DAT   output   * FLXDAT=             *
  56. MASSBAL.DAT   output   * BALDAT=             *
  57. FOG.DAT       output   * FOGDAT=             *
  58. --------------------------------------------------------------------------------
  59. All file names will be converted to lower case if LCFILES = T
  60. Otherwise, if LCFILES = F, file names will be converted to UPPER CASE
  61.          T = lower case      ! LCFILES = F !
  62.          F = UPPER CASE
  63. NOTE: (1) file/path names can be up to 70 characters in length
  64.  
  65.  
  66. Provision for multiple input files
  67. ----------------------------------
  68.  
  69.      Number of CALMET.DAT files for run (NMETDAT)
  70.                                      Default: 1       ! NMETDAT =   1   !
  71.  
  72.      Number of PTEMARB.DAT files for run (NPTDAT)
  73.                                      Default: 0       ! NPTDAT =  0  !
  74.  
  75.      Number of BAEMARB.DAT files for run (NARDAT)
  76.                                      Default: 0       ! NARDAT =  0  !
  77.  
  78.      Number of VOLEMARB.DAT files for run (NVOLDAT)
  79.                                      Default: 0       ! NVOLDAT =  0  !
  80.  
  81. !END!
  82.  
  83. -------------
  84. Subgroup (0a)
  85. -------------
  86.  
  87.   The following CALMET.DAT filenames are processed in sequence if NMETDAT>1
  88.  
  89. Default Name  Type          File Name
  90. ------------  ----          ---------
  91.  none         input    * METDAT=     *   *END*
  92.  
  93.  
  94. --------------------------------------------------------------------------------
  95.  
  96. INPUT GROUP: 1 -- General run control parameters
  97. --------------
  98.  
  99.     Option to run all periods found
  100.     in the met. file     (METRUN)   Default: 0       ! METRUN =   0  !
  101.  
  102.          METRUN = 0 - Run period explicitly defined below
  103.          METRUN = 1 - Run all periods in met. file
  104.  
  105.      Starting date:   Year (IBYR) -- No default       ! IBYR =  2017  !
  106.      (used only if   Month (IBMO) -- No default       ! IBMO =  1  !
  107.       METRUN = 0)      Day (IBDY) -- No default       ! IBDY =  1  !
  108.                       Hour (IBHR) -- No default       ! IBHR =  0  !
  109.  
  110.      Note: IBHR is the time at the END of the first hour of the simulation
  111.            (IBHR=1, the first hour of a day, runs from 00:00 to 01:00)
  112.  
  113.      Base time zone        (XBTZ) -- No default       ! XBTZ = 3.0  !
  114.      The zone is the number of hours that must be
  115.      ADDED to the time to obtain UTC (or GMT)
  116.      Examples: PST = 8., MST = 7.
  117.                CST = 6., EST = 5.
  118.  
  119.      Length of run (hours) (IRLG) -- No default       ! IRLG =  145  !
  120.  
  121.      Number of chemical species (NSPEC)
  122.                                      Default: 5       ! NSPEC =  7   !
  123.  
  124.      Number of chemical species
  125.      to be emitted  (NSE)            Default: 3       ! NSE =  4   !
  126.  
  127.      Flag to stop run after
  128.      SETUP phase (ITEST)             Default: 2       ! ITEST =  2   !
  129.      (Used to allow checking
  130.      of the model inputs, files, etc.)
  131.            ITEST = 1 - STOPS program after SETUP phase
  132.            ITEST = 2 - Continues with execution of program
  133.                        after SETUP
  134.  
  135.      Restart Configuration:
  136.  
  137.         Control flag (MRESTART)      Default: 0       ! MRESTART =  0   !
  138.  
  139.            0 = Do not read or write a restart file
  140.            1 = Read a restart file at the beginning of
  141.                the run
  142.            2 = Write a restart file during run
  143.            3 = Read a restart file at beginning of run
  144.                and write a restart file during run
  145.  
  146.         Number of periods in Restart
  147.         output cycle (NRESPD)        Default: 0       ! NRESPD =  0   !
  148.  
  149.            0 = File written only at last period
  150.           >0 = File updated every NRESPD periods
  151.  
  152.      Meteorological Data Format (METFM)
  153.                                      Default: 1       ! METFM =  1   !
  154.  
  155.            METFM = 1 - CALMET binary file (CALMET.MET)
  156.            METFM = 2 - ISC ASCII file (ISCMET.MET)
  157.            METFM = 3 - AUSPLUME ASCII file (PLMMET.MET)
  158.            METFM = 4 - CTDM plus tower file (PROFILE.DAT) and
  159.                        surface parameters file (SURFACE.DAT)
  160.            METFM = 5 - AERMET tower file (PROFILE.DAT) and
  161.                        surface parameters file (SURFACE.DAT)
  162.  
  163.      Meteorological Profile Data Format (MPRFFM)
  164.             (used only for METFM = 1, 2, 3)
  165.                                      Default: 1       ! MPRFFM =  1   !
  166.  
  167.            MPRFFM = 1 - CTDM plus tower file (PROFILE.DAT)
  168.            MPRFFM = 2 - AERMET tower file (PROFILE.DAT)
  169.  
  170.      PG sigma-y is adjusted by the factor (AVET/PGTIME)**0.2
  171.      Averaging Time (minutes) (AVET)
  172.                                      Default: 60.0    ! AVET = 60. !
  173.      PG Averaging Time (minutes) (PGTIME)
  174.                                      Default: 60.0    ! PGTIME = 60. !
  175.  
  176.  
  177. !END!
  178.  
  179.  
  180. -------------------------------------------------------------------------------
  181.  
  182. INPUT GROUP: 2 -- Technical options
  183. --------------
  184.  
  185.  
  186.      Vertical distribution used in the
  187.      near field (MGAUSS)                   Default: 1     ! MGAUSS =  1   !
  188.         0 = uniform
  189.         1 = Gaussian
  190.  
  191.      Terrain adjustment method
  192.      (MCTADJ)                              Default: 3     ! MCTADJ =  3   !
  193.         0 = no adjustment
  194.         1 = ISC-type of terrain adjustment
  195.         2 = simple, CALPUFF-type of terrain
  196.             adjustment
  197.         3 = partial plume path adjustment
  198.  
  199.      Subgrid-scale complex terrain
  200.      flag (MCTSG)                          Default: 0     ! MCTSG =  0   !
  201.         0 = not modeled
  202.         1 = modeled
  203.  
  204.      Near-field puffs modeled as
  205.      elongated slugs? (MSLUG)              Default: 0     ! MSLUG =  0   !
  206.         0 = no
  207.         1 = yes (slug model used)
  208.  
  209.      Transitional plume rise modeled?
  210.      (MTRANS)                              Default: 1     ! MTRANS =  1   !
  211.         0 = no  (i.e., final rise only)
  212.         1 = yes (i.e., transitional rise computed)
  213.  
  214.      Stack tip downwash? (MTIP)            Default: 1     ! MTIP =  1  !
  215.         0 = no  (i.e., no stack tip downwash)
  216.         1 = yes (i.e., use stack tip downwash)
  217.  
  218.      Method used to simulate building
  219.      downwash? (MBDW)                      Default: 1     ! MBDW =   1  !
  220.         1 = ISC method
  221.         2 = PRIME method
  222.  
  223.      Vertical wind shear modeled above
  224.      stack top? (MSHEAR)                   Default: 0     ! MSHEAR =  0  !
  225.         0 = no  (i.e., vertical wind shear not modeled)
  226.         1 = yes (i.e., vertical wind shear modeled)
  227.  
  228.      Puff splitting allowed? (MSPLIT)      Default: 0     ! MSPLIT =  0  !
  229.         0 = no (i.e., puffs not split)
  230.         1 = yes (i.e., puffs are split)
  231.  
  232.      Chemical mechanism flag (MCHEM)       Default: 1     ! MCHEM =  3   !
  233.         0 = chemical transformation not
  234.             modeled
  235.         1 = transformation rates computed
  236.             internally (MESOPUFF II scheme)
  237.         2 = user-specified transformation
  238.             rates used
  239.         3 = transformation rates computed
  240.             internally (RIVAD/ARM3 scheme)
  241.         4 = secondary organic aerosol formation
  242.             computed (MESOPUFF II scheme for OH)
  243.  
  244.      Aqueous phase transformation flag (MAQCHEM)
  245.      (Used only if MCHEM = 1, or 3)        Default: 0     ! MAQCHEM =  0   !
  246.         0 = aqueous phase transformation
  247.             not modeled
  248.         1 = transformation rates adjusted
  249.             for aqueous phase reactions
  250.  
  251.      Wet removal modeled ? (MWET)          Default: 1     ! MWET =  1   !
  252.         0 = no
  253.         1 = yes
  254.  
  255.      Dry deposition modeled ? (MDRY)       Default: 1     ! MDRY =  1   !
  256.         0 = no
  257.         1 = yes
  258.         (dry deposition method specified
  259.          for each species in Input Group 3)
  260.  
  261.  
  262.      Gravitational settling (plume tilt)
  263.      modeled ? (MTILT)                     Default: 0     ! MTILT =  0   !
  264.         0 = no
  265.         1 = yes
  266.         (puff center falls at the gravitational
  267.          settling velocity for 1 particle species)
  268.  
  269.      Restrictions:
  270.          - MDRY  = 1
  271.          - NSPEC = 1  (must be particle species as well)
  272.          - sg    = 0  GEOMETRIC STANDARD DEVIATION in Group 8 is
  273.                       set to zero for a single particle diameter
  274.  
  275.      Method used to compute dispersion
  276.      coefficients (MDISP)                  Default: 3     ! MDISP =  3   !
  277.  
  278.         1 = dispersion coefficients computed from measured values
  279.             of turbulence, sigma v, sigma w
  280.         2 = dispersion coefficients from internally calculated
  281.             sigma v, sigma w using micrometeorological variables
  282.             (u*, w*, L, etc.)
  283.         3 = PG dispersion coefficients for RURAL areas (computed using
  284.             the ISCST multi-segment approximation) and MP coefficients in
  285.             urban areas
  286.         4 = same as 3 except PG coefficients computed using
  287.             the MESOPUFF II eqns.
  288.         5 = CTDM sigmas used for stable and neutral conditions.
  289.             For unstable conditions, sigmas are computed as in
  290.             MDISP = 3, described above.  MDISP = 5 assumes that
  291.             measured values are read
  292.  
  293.      Sigma-v/sigma-theta, sigma-w measurements used? (MTURBVW)
  294.      (Used only if MDISP = 1 or 5)         Default: 3     ! MTURBVW =  3  !
  295.         1 = use sigma-v or sigma-theta measurements
  296.             from PROFILE.DAT to compute sigma-y
  297.             (valid for METFM = 1, 2, 3, 4, 5)
  298.         2 = use sigma-w measurements
  299.             from PROFILE.DAT to compute sigma-z
  300.             (valid for METFM = 1, 2, 3, 4, 5)
  301.         3 = use both sigma-(v/theta) and sigma-w
  302.             from PROFILE.DAT to compute sigma-y and sigma-z
  303.             (valid for METFM = 1, 2, 3, 4, 5)
  304.         4 = use sigma-theta measurements
  305.             from PLMMET.DAT to compute sigma-y
  306.             (valid only if METFM = 3)
  307.  
  308.      Back-up method used to compute dispersion
  309.      when measured turbulence data are
  310.      missing (MDISP2)                      Default: 3     ! MDISP2 =  3  !
  311.      (used only if MDISP = 1 or 5)
  312.         2 = dispersion coefficients from internally calculated
  313.             sigma v, sigma w using micrometeorological variables
  314.             (u*, w*, L, etc.)
  315.         3 = PG dispersion coefficients for RURAL areas (computed using
  316.             the ISCST multi-segment approximation) and MP coefficients in
  317.             urban areas
  318.         4 = same as 3 except PG coefficients computed using
  319.             the MESOPUFF II eqns.
  320.  
  321.      [DIAGNOSTIC FEATURE]
  322.      Method used for Lagrangian timescale for Sigma-y
  323.      (used only if MDISP=1,2 or MDISP2=1,2)
  324.      (MTAULY)                              Default: 0     ! MTAULY =  0  !
  325.         0 = Draxler default 617.284 (s)
  326.         1 = Computed as Lag. Length / (.75 q) -- after SCIPUFF
  327.        10 < Direct user input (s)             -- e.g., 306.9
  328.  
  329.  
  330.      [DIAGNOSTIC FEATURE]
  331.      Method used for Advective-Decay timescale for Turbulence
  332.      (used only if MDISP=2 or MDISP2=2)
  333.      (MTAUADV)                             Default: 0     ! MTAUADV =  0  !
  334.         0 = No turbulence advection
  335.         1 = Computed (OPTION NOT IMPLEMENTED)
  336.        10 < Direct user input (s)   -- e.g., 300
  337.  
  338.  
  339.      Method used to compute turbulence sigma-v &
  340.      sigma-w using micrometeorological variables
  341.      (Used only if MDISP = 2 or MDISP2 = 2)
  342.      (MCTURB)                              Default: 1     ! MCTURB =  1  !
  343.         1 = Standard CALPUFF subroutines
  344.         2 = AERMOD subroutines
  345.  
  346.      PG sigma-y,z adj. for roughness?      Default: 0     ! MROUGH =  0  !
  347.      (MROUGH)
  348.         0 = no
  349.         1 = yes
  350.  
  351.      Partial plume penetration of          Default: 1     ! MPARTL =  1  !
  352.      elevated inversion?
  353.      (MPARTL)
  354.         0 = no
  355.         1 = yes
  356.  
  357.      Strength of temperature inversion     Default: 0     ! MTINV =  0  !
  358.      provided in PROFILE.DAT extended records?
  359.      (MTINV)
  360.         0 = no (computed from measured/default gradients)
  361.         1 = yes
  362.  
  363.      PDF used for dispersion under convective conditions?
  364.                                            Default: 0     ! MPDF =  0  !
  365.      (MPDF)
  366.         0 = no
  367.         1 = yes
  368.  
  369.      Sub-Grid TIBL module used for shore line?
  370.                                            Default: 0     ! MSGTIBL = 0  !
  371.      (MSGTIBL)
  372.         0 = no
  373.         1 = yes
  374.  
  375.      Boundary conditions (concentration) modeled?
  376.                                            Default: 0     ! MBCON = 0  !
  377.      (MBCON)
  378.         0 = no
  379.         1 = yes, using formatted BCON.DAT file
  380.         2 = yes, using unformatted CONC.DAT file
  381.  
  382.      Note:  MBCON > 0 requires that the last species modeled
  383.             be 'BCON'.  Mass is placed in species BCON when
  384.             generating boundary condition puffs so that clean
  385.             air entering the modeling domain can be simulated
  386.             in the same way as polluted air.  Specify zero
  387.             emission of species BCON for all regular sources.
  388.  
  389.      Individual source contributions saved?
  390.                                            Default: 0     ! MSOURCE = 0  !
  391.      (MSOURCE)
  392.         0 = no
  393.         1 = yes
  394.  
  395.  
  396.      Analyses of fogging and icing impacts due to emissions from
  397.      arrays of mechanically-forced cooling towers can be performed
  398.      using CALPUFF in conjunction with a cooling tower emissions
  399.      processor (CTEMISS) and its associated postprocessors.  Hourly
  400.      emissions of water vapor and temperature from each cooling tower
  401.      cell are computed for the current cell configuration and ambient
  402.      conditions by CTEMISS. CALPUFF models the dispersion of these
  403.      emissions and provides cloud information in a specialized format
  404.      for further analysis. Output to FOG.DAT is provided in either
  405.      'plume mode' or 'receptor mode' format.
  406.  
  407.      Configure for FOG Model output?
  408.                                            Default: 0     ! MFOG =  0   !
  409.      (MFOG)
  410.         0 = no
  411.         1 = yes  - report results in PLUME Mode format
  412.         2 = yes  - report results in RECEPTOR Mode format
  413.  
  414.  
  415.      Test options specified to see if
  416.      they conform to regulatory
  417.      values? (MREG)                        Default: 1     ! MREG =  0   !
  418.  
  419.         0 = NO checks are made
  420.         1 = Technical options must conform to USEPA
  421.             Long Range Transport (LRT) guidance
  422.                        METFM    1 or 2
  423.                        AVET     60. (min)
  424.                        PGTIME   60. (min)
  425.                        MGAUSS   1
  426.                        MCTADJ   3
  427.                        MTRANS   1
  428.                        MTIP     1
  429.                        MCHEM    1 or 3 (if modeling SOx, NOx)
  430.                        MWET     1
  431.                        MDRY     1
  432.                        MDISP    2 or 3
  433.                        MPDF     0 if MDISP=3
  434.                                 1 if MDISP=2
  435.                        MROUGH   0
  436.                        MPARTL   1
  437.                        SYTDEP   550. (m)
  438.                        MHFTSZ   0
  439.                        SVMIN    0.5 (m/s)
  440.  
  441.  
  442. !END!
  443.  
  444.  
  445. -------------------------------------------------------------------------------
  446.  
  447. INPUT GROUP: 3a, 3b -- Species list
  448. -------------------
  449.  
  450. ------------
  451. Subgroup (3a)
  452. ------------
  453.  
  454.   The following species are modeled:
  455.  
  456. z
  457.  
  458.                                                        Dry                OUTPUT GROUP
  459.     SPECIES          MODELED          EMITTED       DEPOSITED                NUMBER
  460.      NAME         (0=NO, 1=YES)    (0=NO, 1=YES)    (0=NO,                 (0=NONE,
  461.    (Limit: 12                                        1=COMPUTED-GAS        1=1st CGRUP,
  462.     Characters                                       2=COMPUTED-PARTICLE   2=2nd CGRUP,
  463.     in length)                                       3=USER-SPECIFIED)     3= etc.)
  464.  
  465. !          SO2  =         1,               1,           1,                 0   !
  466. !          SO4  =         1,               0,           2,                 0   !
  467. !           NO  =         1,               1,           1,                 0   !
  468. !          NO2  =         1,               1,           1,                 0   !
  469. !         HNO3  =         1,               0,           1,                 0   !
  470. !          NO3  =         1,               0,           2,                 0   !
  471. !         PM10  =         1,               1,           2,                 0   !
  472.  
  473. !END!
  474.  
  475.   Note:  The last species in (3a) must be 'BCON' when using the
  476.          boundary condition option (MBCON > 0).  Species BCON should
  477.          typically be modeled as inert (no chem transformation or
  478.          removal).
  479.  
  480.  
  481. -------------
  482. Subgroup (3b)
  483. -------------
  484.   The following names are used for Species-Groups in which results
  485.   for certain species are combined (added) prior to output.  The
  486.   CGRUP name will be used as the species name in output files.
  487.   Use this feature to model specific particle-size distributions
  488.   by treating each size-range as a separate species.
  489.   Order must be consistent with 3(a) above.
  490.  
  491.  
  492.  
  493. -------------------------------------------------------------------------------
  494.  
  495.  
  496. INPUT GROUP: 4 -- Map Projection and Grid control parameters
  497. --------------
  498.  
  499.      Projection for all (X,Y):
  500.      -------------------------
  501.  
  502.      Map projection
  503.      (PMAP)                     Default: UTM    ! PMAP = LCC  !
  504.  
  505.          UTM :  Universal Transverse Mercator
  506.          TTM :  Tangential Transverse Mercator
  507.          LCC :  Lambert Conformal Conic
  508.           PS :  Polar Stereographic
  509.           EM :  Equatorial Mercator
  510.         LAZA :  Lambert Azimuthal Equal Area
  511.  
  512.      False Easting and Northing (km) at the projection origin
  513.      (Used only if PMAP= TTM, LCC, or LAZA)
  514.      (FEAST)                    Default=0.0     ! FEAST  = 254.763  !
  515.      (FNORTH)                   Default=0.0     ! FNORTH = 6341.324  !
  516.  
  517.      UTM zone (1 to 60)
  518.      (Used only if PMAP=UTM)
  519.      (IUTMZN)                   No Default      ! IUTMZN =  19   !
  520.  
  521.      Hemisphere for UTM projection?
  522.      (Used only if PMAP=UTM)
  523.      (UTMHEM)                   Default: N      ! UTMHEM = S  !
  524.          N   :  Northern hemisphere projection
  525.          S   :  Southern hemisphere projection
  526.  
  527.      Latitude and Longitude (decimal degrees) of projection origin
  528.      (Used only if PMAP= TTM, LCC, PS, EM, or LAZA)
  529.      (RLAT0)                    No Default      ! RLAT0 = 33.039S  !
  530.      (RLON0)                    No Default      ! RLON0 = 71.626W  !
  531.  
  532.          TTM :  RLON0 identifies central (true N/S) meridian of projection
  533.                 RLAT0 selected for convenience
  534.          LCC :  RLON0 identifies central (true N/S) meridian of projection
  535.                 RLAT0 selected for convenience
  536.          PS  :  RLON0 identifies central (grid N/S) meridian of projection
  537.                 RLAT0 selected for convenience
  538.          EM  :  RLON0 identifies central meridian of projection
  539.                 RLAT0 is REPLACED by 0.0N (Equator)
  540.          LAZA:  RLON0 identifies longitude of tangent-point of mapping plane
  541.                 RLAT0 identifies latitude of tangent-point of mapping plane
  542.  
  543.      Matching parallel(s) of latitude (decimal degrees) for projection
  544.      (Used only if PMAP= LCC or PS)
  545.      (XLAT1)                    No Default      ! XLAT1 = 33.039S  !
  546.      (XLAT2)                    No Default      ! XLAT2 = 33.039S  !
  547.  
  548.          LCC :  Projection cone slices through Earth's surface at XLAT1 and XLAT2
  549.          PS  :  Projection plane slices through Earth at XLAT1
  550.                 (XLAT2 is not used)
  551.  
  552.      ----------
  553.      Note:  Latitudes and longitudes should be positive, and include a
  554.             letter N,S,E, or W indicating north or south latitude, and
  555.             east or west longitude.  For example,
  556.             35.9  N Latitude  =  35.9N
  557.             118.7 E Longitude = 118.7E
  558.  
  559.  
  560.      Datum-region
  561.      ------------
  562.  
  563.      The Datum-Region for the coordinates is identified by a character
  564.      string.  Many mapping products currently available use the model of the
  565.      Earth known as the World Geodetic System 1984 (WGS-84).  Other local
  566.      models may be in use, and their selection in CALMET will make its output
  567.      consistent with local mapping products.  The list of Datum-Regions with
  568.      official transformation parameters is provided by the National Imagery and
  569.      Mapping Agency (NIMA).
  570.  
  571.      NIMA Datum - Regions(Examples)
  572.      ------------------------------------------------------------------------------
  573.      WGS-84    WGS-84 Reference Ellipsoid and Geoid, Global coverage (WGS84)
  574.      NAS-C     NORTH AMERICAN 1927 Clarke 1866 Spheroid, MEAN FOR CONUS (NAD27)
  575.      NAR-C     NORTH AMERICAN 1983 GRS 80 Spheroid, MEAN FOR CONUS (NAD83)
  576.      NWS-84    NWS 6370KM Radius, Sphere
  577.      ESR-S     ESRI REFERENCE 6371KM Radius, Sphere
  578.  
  579.      Datum-region for output coordinates
  580.      (DATUM)                    Default: WGS-84    ! DATUM = NWS-84  !
  581.  
  582.  
  583. METEOROLOGICAL Grid:
  584.  
  585.      Rectangular grid defined for projection PMAP,
  586.      with X the Easting and Y the Northing coordinate
  587.  
  588.             No. X grid cells (NX)      No default     ! NX =  29   !
  589.             No. Y grid cells (NY)      No default     ! NY =  29   !
  590.          No. vertical layers (NZ)      No default     ! NZ =  10   !
  591.  
  592.            Grid spacing (DGRIDKM)      No default     ! DGRIDKM = 3.0 !
  593.                                        Units: km
  594.  
  595.                 Cell face heights
  596.                     (ZFACE(nz+1))      No defaults
  597.                                        Units: m
  598.    ! ZFACE = .0, 20.0, 40.0, 80.0, 160.0, 320.0, 640.0, 1200.0, 2000.0, 3000.0,
  599.               4000.0 !
  600.  
  601.             Reference Coordinates
  602.            of SOUTHWEST corner of
  603.                  grid cell(1, 1):
  604.  
  605.             X coordinate (XORIGKM)     No default     ! XORIGKM = -43.5 !
  606.             Y coordinate (YORIGKM)     No default     ! YORIGKM = -43.5 !
  607.                                       Units: km
  608.  
  609.  
  610. COMPUTATIONAL Grid:
  611.  
  612.      The computational grid is identical to or a subset of the MET. grid.
  613.      The lower left (LL) corner of the computational grid is at grid point
  614.      (IBCOMP, JBCOMP) of the MET. grid.  The upper right (UR) corner of the
  615.      computational grid is at grid point (IECOMP, JECOMP) of the MET. grid.
  616.      The grid spacing of the computational grid is the same as the MET. grid.
  617.  
  618.         X index of LL corner (IBCOMP)      No default     ! IBCOMP =  1   !
  619.                   (1 <= IBCOMP <= NX)
  620.  
  621.         Y index of LL corner (JBCOMP)      No default     ! JBCOMP =  1   !
  622.                   (1 <= JBCOMP <= NY)
  623.  
  624.  
  625.         X index of UR corner (IECOMP)      No default     ! IECOMP =  29   !
  626.                   (1 <= IECOMP <= NX)
  627.  
  628.         Y index of UR corner (JECOMP)      No default     ! JECOMP =  29   !
  629.                   (1 <= JECOMP <= NY)
  630.  
  631.  
  632.  
  633. SAMPLING Grid (GRIDDED RECEPTORS):
  634.  
  635.      The lower left (LL) corner of the sampling grid is at grid point
  636.      (IBSAMP, JBSAMP) of the MET. grid.  The upper right (UR) corner of the
  637.      sampling grid is at grid point (IESAMP, JESAMP) of the MET. grid.
  638.      The sampling grid must be identical to or a subset of the computational
  639.      grid.  It may be a nested grid inside the computational grid.
  640.      The grid spacing of the sampling grid is DGRIDKM/MESHDN.
  641.  
  642.         Logical flag indicating if gridded
  643.         receptors are used (LSAMP)         Default: T     ! LSAMP = T !
  644.         (T=yes, F=no)
  645.  
  646.         X index of LL corner (IBSAMP)      No default     ! IBSAMP =  1   !
  647.          (IBCOMP <= IBSAMP <= IECOMP)
  648.  
  649.         Y index of LL corner (JBSAMP)      No default     ! JBSAMP =  1   !
  650.          (JBCOMP <= JBSAMP <= JECOMP)
  651.  
  652.  
  653.         X index of UR corner (IESAMP)      No default     ! IESAMP =  29   !
  654.          (IBCOMP <= IESAMP <= IECOMP)
  655.  
  656.         Y index of UR corner (JESAMP)      No default     ! JESAMP =  29   !
  657.          (JBCOMP <= JESAMP <= JECOMP)
  658.  
  659.  
  660.        Nesting factor of the sampling
  661.         grid (MESHDN)                      Default: 1     ! MESHDN =  1  !
  662.         (MESHDN is an integer >= 1)
  663.  
  664. !END!
  665.  
  666.  
  667. -------------------------------------------------------------------------------
  668.  
  669.  
  670. INPUT GROUP: 5 -- Output Options
  671. --------------
  672.                                              *                          *
  673.      FILE                       DEFAULT VALUE             VALUE THIS RUN
  674.      ----                       -------------             --------------
  675.  
  676.    Concentrations (ICON)              1                   !  ICON =  1   !
  677.    Dry Fluxes (IDRY)                  1                   !  IDRY =  1   !
  678.    Wet Fluxes (IWET)                  1                   !  IWET =  1   !
  679.    2D Temperature (IT2D)              0                   !  IT2D =  0   !
  680.    2D Density (IRHO)                  0                   !  IRHO =  0   !
  681.    Relative Humidity (IVIS)           1                   !  IVIS =  1   !
  682.     (relative humidity file is
  683.      required for visibility
  684.      analysis)
  685.    Use data compression option in output file?
  686.    (LCOMPRS)                           Default: T         ! LCOMPRS = T !
  687.  
  688.    *
  689.     0 = Do not create file, 1 = create file
  690.  
  691.  
  692.     QA PLOT FILE OUTPUT OPTION:
  693.  
  694.        Create a standard series of output files (e.g.
  695.        locations of sources, receptors, grids ...)
  696.        suitable for plotting?
  697.        (IQAPLOT)                       Default: 1         !  IQAPLOT =  1   !
  698.          0 = no
  699.          1 = yes
  700.  
  701.     DIAGNOSTIC MASS FLUX OUTPUT OPTIONS:
  702.  
  703.        Mass flux across specified boundaries
  704.        for selected species reported?
  705.        (IMFLX)                         Default: 0         ! IMFLX =  0  !
  706.          0 = no
  707.          1 = yes (FLUXBDY.DAT and MASSFLX.DAT filenames
  708.                   are specified in Input Group 0)
  709.  
  710.        Mass balance for each species
  711.        reported?
  712.        (IMBAL)                         Default: 0         ! IMBAL =  0  !
  713.          0 = no
  714.          1 = yes (MASSBAL.DAT filename is
  715.               specified in Input Group 0)
  716.  
  717.  
  718.     LINE PRINTER OUTPUT OPTIONS:
  719.  
  720.        Print concentrations (ICPRT)    Default: 0         ! ICPRT =  1   !
  721.        Print dry fluxes (IDPRT)        Default: 0         ! IDPRT =  0   !
  722.        Print wet fluxes (IWPRT)        Default: 0         ! IWPRT =  0   !
  723.        (0 = Do not print, 1 = Print)
  724.  
  725.        Concentration print interval
  726.        (ICFRQ) in timesteps            Default: 1         ! ICFRQ =  1   !
  727.        Dry flux print interval
  728.        (IDFRQ) in timesteps            Default: 1         ! IDFRQ =  1   !
  729.        Wet flux print interval
  730.        (IWFRQ) in timesteps            Default: 1         ! IWFRQ =  1   !
  731.  
  732.        Units for Line Printer Output
  733.        (IPRTU)                         Default: 1         ! IPRTU =  3   !
  734.                        for            for
  735.                   Concentration    Deposition
  736.            1 =       g/m**3         g/m**2/s
  737.            2 =      mg/m**3        mg/m**2/s
  738.            3 =      ug/m**3        ug/m**2/s
  739.            4 =      ng/m**3        ng/m**2/s
  740.            5 =     Odour Units
  741.  
  742.        Messages tracking progress of run
  743.        written to the screen ?
  744.        (IMESG)                         Default: 2         ! IMESG =  2   !
  745.          0 = no
  746.          1 = yes (advection step, puff ID)
  747.          2 = yes (YYYYJJJHH, # old puffs, # emitted puffs)
  748.  
  749.  
  750.      SPECIES (or GROUP for combined species) LIST FOR OUTPUT OPTIONS
  751.  
  752.                  ---- CONCENTRATIONS ----   ------ DRY FLUXES ------   ------ WET FLUXES ------   -- MASS FLUX --
  753.    SPECIES
  754.    /GROUP        PRINTED?  SAVED ON DISK?   PRINTED?  SAVED ON DISK?   PRINTED?  SAVED ON DISK?   SAVED ON DISK?
  755.    -------       ------------------------   ------------------------   ------------------------   ---------------
  756. !          SO2 =     1,           1,           0,           1,           0,           1,           0   !
  757. !          SO4 =     0,           1,           0,           1,           0,           1,           0   !
  758. !           NO =     0,           1,           0,           1,           0,           1,           0   !
  759. !          NO2 =     0,           1,           0,           1,           0,           1,           0   !
  760. !         HNO3 =     0,           1,           0,           1,           0,           1,           0   !
  761. !          NO3 =     0,           1,           0,           1,           0,           1,           0   !
  762. !         PM10 =     0,           1,           0,           1,           0,           1,           0   !
  763.  
  764.   Note:  Species BCON (for MBCON > 0) does not need to be saved on disk.
  765.  
  766.  
  767.      OPTIONS FOR PRINTING "DEBUG" QUANTITIES (much output)  
  768.  
  769.        Logical for debug output
  770.        (LDEBUG)                                 Default: F     ! LDEBUG = F !
  771.  
  772.        First puff to track
  773.        (IPFDEB)                                 Default: 1     ! IPFDEB =  1  !
  774.  
  775.        Number of puffs to track
  776.        (NPFDEB)                                 Default: 1     ! NPFDEB =  1  !
  777.  
  778.        Met. period to start output
  779.        (NN1)                                    Default: 1     ! NN1 =  1   !
  780.  
  781.        Met. period to end output
  782.        (NN2)                                    Default: 10    ! NN2 =  10  !
  783.  
  784. !END!
  785.  
  786.  
  787. -------------------------------------------------------------------------------
  788.  
  789.  
  790. INPUT GROUP: 6a, 6b, & 6c -- Subgrid scale complex terrain inputs
  791. -------------------------
  792.  
  793. ---------------
  794. Subgroup (6a)
  795. ---------------
  796.        Number of terrain features (NHILL)       Default: 0     ! NHILL =  0   !
  797.  
  798.        Number of special complex terrain
  799.        receptors  (NCTREC)                      Default: 0     ! NCTREC =  0   !
  800.  
  801.        Terrain and CTSG Receptor data for
  802.        CTSG hills input in CTDM format ?
  803.        (MHILL)                                  No Default     ! MHILL =  2   !
  804.        1 = Hill and Receptor data created
  805.            by CTDM processors & read from
  806.            HILL.DAT and HILLRCT.DAT files
  807.        2 = Hill data created by OPTHILL &
  808.            input below in Subgroup (6b);
  809.            Receptor data in Subgroup (6c)
  810.  
  811.        Factor to convert horizontal dimensions  Default: 1.0   ! XHILL2M = 1.0 !
  812.        to meters (MHILL=1)
  813.  
  814.        Factor to convert vertical dimensions    Default: 1.0   ! ZHILL2M = 1.0 !
  815.        to meters (MHILL=1)
  816.  
  817.        X-origin of CTDM system relative to      No Default     ! XCTDMKM = 0 !
  818.        CALPUFF coordinate system, in Kilometers (MHILL=1)
  819.  
  820.        Y-origin of CTDM system relative to      No Default     ! YCTDMKM = 0 !
  821.        CALPUFF coordinate system, in Kilometers (MHILL=1)
  822.  
  823. ! END !
  824.  
  825. ---------------
  826. Subgroup (6b)
  827. ---------------
  828.  
  829.                       1 **
  830.      HILL information
  831.  
  832.  
  833. HILL           XC        YC       THETAH  ZGRID  RELIEF    EXPO 1    EXPO 2   SCALE 1    SCALE 2    AMAX1     AMAX2
  834.  NO.          (km)      (km)      (deg.)   (m)     (m)      (m)       (m)       (m)        (m)       (m)       (m)
  835. ----          ----      ----      ------  -----  ------    ------    ------   -------    -------    -----     -----
  836.  
  837. ---------------
  838. Subgroup (6c)
  839. ---------------
  840.  
  841.     COMPLEX TERRAIN RECEPTOR INFORMATION
  842.  
  843.                       XRCT         YRCT        ZRCT          XHH
  844.                       (km)         (km)         (m)
  845.                      ------        -----      ------         ----
  846.  
  847.  
  848. -------------------
  849. 1
  850.      Description of Complex Terrain Variables:
  851.           XC, YC  = Coordinates of center of hill
  852.           THETAH  = Orientation of major axis of hill (clockwise from
  853.                     North)
  854.           ZGRID   = Height of the  0  of the grid above mean sea
  855.                     level
  856.           RELIEF  = Height of the crest of the hill above the grid elevation
  857.           EXPO 1  = Hill-shape exponent for the major axis
  858.           EXPO 2  = Hill-shape exponent for the major axis
  859.           SCALE 1 = Horizontal length scale along the major axis
  860.           SCALE 2 = Horizontal length scale along the minor axis
  861.           AMAX    = Maximum allowed axis length for the major axis
  862.           BMAX    = Maximum allowed axis length for the major axis
  863.  
  864.           XRCT, YRCT = Coordinates of the complex terrain receptors
  865.           ZRCT    = Height of the ground (MSL) at the complex terrain
  866.                     Receptor
  867.           XHH     = Hill number associated with each complex terrain receptor
  868.                     (NOTE: MUST BE ENTERED AS A REAL NUMBER)
  869.  
  870.    **
  871.      NOTE: DATA for each hill and CTSG receptor are treated as a separate
  872.            input subgroup and therefore must end with an input group terminator.
  873.  
  874. -------------------------------------------------------------------------------
  875.  
  876.  
  877. INPUT GROUP: 7 -- Chemical parameters for dry deposition of gases
  878. --------------
  879.  
  880.       SPECIES     DIFFUSIVITY      ALPHA STAR      REACTIVITY    MESOPHYLL RESISTANCE     HENRY'S LAW COEFFICIENT
  881.        NAME        (cm**2/s)                                            (s/cm)                (dimensionless)
  882.       -------     -----------      ----------      ----------    --------------------     -----------------------
  883.  
  884. !          SO2 =      .1509,        1000.0,           8.0,                .0,                   .04 !
  885. !           NO =      .1345,           1.0,           2.0,              25.0,                  18.0 !
  886. !          NO2 =      .1656,           1.0,           8.0,               5.0,                   3.5 !
  887. !         HNO3 =      .1628,           1.0,          18.0,                .0,              .0000001 !
  888.  
  889. !END!
  890.  
  891.  
  892. -------------------------------------------------------------------------------
  893.  
  894.  
  895. INPUT GROUP: 8 -- Size parameters for dry deposition of particles
  896. --------------
  897.  
  898.      For SINGLE SPECIES, the mean and standard deviation are used to
  899.      compute a deposition velocity for NINT (see group 9) size-ranges,
  900.      and these are then averaged to obtain a mean deposition velocity.
  901.  
  902.      For GROUPED SPECIES, the size distribution should be explicitly
  903.      specified (by the 'species' in the group), and the standard deviation
  904.      for each should be entered as 0.  The model will then use the
  905.      deposition velocity for the stated mean diameter.
  906.  
  907.       SPECIES      GEOMETRIC MASS MEAN        GEOMETRIC STANDARD
  908.        NAME             DIAMETER                   DEVIATION
  909.                         (microns)                  (microns)
  910.       -------      -------------------        ------------------
  911. !          SO4 =           .48,                     2.0   !
  912. !          NO3 =           .48,                     2.0   !
  913. !         PM10 =           .48,                     2.0   !
  914.  
  915. !END!
  916.  
  917.  
  918. -------------------------------------------------------------------------------
  919.  
  920.  
  921. INPUT GROUP: 9 -- Miscellaneous dry deposition parameters
  922. --------------
  923.  
  924.      Reference cuticle resistance (s/cm)
  925.      (RCUTR)                           Default: 30    !  RCUTR = 30.0 !
  926.      Reference ground resistance  (s/cm)
  927.      (RGR)                             Default: 10    !    RGR = 5.0 !
  928.      Reference pollutant reactivity
  929.      (REACTR)                          Default: 8     ! REACTR = 8.0 !
  930.  
  931.      Number of particle-size intervals used to
  932.      evaluate effective particle deposition velocity
  933.      (NINT)                            Default: 9     !   NINT =  9  !
  934.  
  935.      Vegetation state in unirrigated areas
  936.      (IVEG)                            Default: 1     !   IVEG =  1   !
  937.         IVEG=1 for active and unstressed vegetation
  938.         IVEG=2 for active and stressed vegetation
  939.         IVEG=3 for inactive vegetation
  940.  
  941. !END!
  942.  
  943.  
  944. -------------------------------------------------------------------------------
  945.  
  946.  
  947. INPUT GROUP: 10 -- Wet Deposition Parameters
  948. ---------------
  949.  
  950.                                                          
  951.                       Scavenging Coefficient -- Units: (sec)**(-1)
  952.  
  953.        Pollutant      Liquid Precip.       Frozen Precip.
  954.        ---------      --------------       --------------
  955. !          SO2 =         3.0E-05,              0.0E00 !
  956. !          SO4 =         1.0E-04,             3.0E-05 !
  957. !         HNO3 =         6.0E-05,              0.0E00 !
  958. !          NO3 =         1.0E-04,             3.0E-05 !
  959. !         PM10 =         1.0E-04,             3.0E-05 !
  960.  
  961. !END!
  962.  
  963.  
  964. -------------------------------------------------------------------------------
  965.  
  966.  
  967. INPUT GROUP: 11 -- Chemistry Parameters
  968. ---------------
  969.  
  970.      Ozone data input option (MOZ)     Default: 1            ! MOZ =  0   !
  971.      (Used only if MCHEM = 1, 3, or 4)
  972.         0 = use a monthly background ozone value
  973.         1 = read hourly ozone concentrations from
  974.             the OZONE.DAT data file
  975.  
  976.      Monthly ozone concentrations
  977.      (Used only if MCHEM = 1, 3, or 4 and
  978.       MOZ = 0 or MOZ = 1 and all hourly O3 data missing)
  979.      (BCKO3) in ppb                    Default: 12*80.
  980.      !  BCKO3 = 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00, 40.00 !
  981.  
  982.      Monthly ammonia concentrations
  983.      (Used only if MCHEM = 1, or 3)
  984.      (BCKNH3) in ppb                   Default: 12*10.      
  985.      !  BCKNH3 = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00 !
  986.  
  987.      Nighttime SO2 loss rate (RNITE1)
  988.      in percent/hour                   Default: 0.2          ! RNITE1 = .2 !
  989.  
  990.      Nighttime NOx loss rate (RNITE2)
  991.      in percent/hour                   Default: 2.0          ! RNITE2 = 2.0 !
  992.  
  993.      Nighttime HNO3 formation rate (RNITE3)
  994.      in percent/hour                   Default: 2.0          ! RNITE3 = 2.0 !
  995.  
  996.      H2O2 data input option (MH2O2)    Default: 1            ! MH2O2 =  1   !
  997.      (Used only if MAQCHEM = 1)
  998.         0 = use a monthly background H2O2 value
  999.         1 = read hourly H2O2 concentrations from
  1000.             the H2O2.DAT data file
  1001.  
  1002.      Monthly H2O2 concentrations
  1003.      (Used only if MQACHEM = 1 and
  1004.       MH2O2 = 0 or MH2O2 = 1 and all hourly H2O2 data missing)
  1005.      (BCKH2O2) in ppb                  Default: 12*1.        
  1006.      !  BCKH2O2 = 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00 !
  1007.  
  1008.  
  1009.  --- Data for SECONDARY ORGANIC AEROSOL (SOA) Option
  1010.      (used only if MCHEM = 4)
  1011.  
  1012.      The SOA module uses monthly values of:
  1013.           Fine particulate concentration in ug/m^3 (BCKPMF)
  1014.           Organic fraction of fine particulate     (OFRAC)
  1015.           VOC / NOX ratio (after reaction)         (VCNX)
  1016.      to characterize the air mass when computing
  1017.      the formation of SOA from VOC emissions.
  1018.      Typical values for several distinct air mass types are:
  1019.  
  1020.         Month    1    2    3    4    5    6    7    8    9   10   11   12
  1021.                 Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov  Dec
  1022.  
  1023.      Clean Continental
  1024.         BCKPMF   1.   1.   1.   1.   1.   1.   1.   1.   1.   1.   1.   1.
  1025.         OFRAC  .15  .15  .20  .20  .20  .20  .20  .20  .20  .20  .20  .15
  1026.         VCNX    50.  50.  50.  50.  50.  50.  50.  50.  50.  50.  50.  50.
  1027.  
  1028.      Clean Marine (surface)
  1029.         BCKPMF  .5   .5   .5   .5   .5   .5   .5   .5   .5   .5   .5   .5
  1030.         OFRAC  .25  .25  .30  .30  .30  .30  .30  .30  .30  .30  .30  .25
  1031.         VCNX    50.  50.  50.  50.  50.  50.  50.  50.  50.  50.  50.  50.
  1032.  
  1033.      Urban - low biogenic (controls present)
  1034.         BCKPMF  30.  30.  30.  30.  30.  30.  30.  30.  30.  30.  30.  30.
  1035.         OFRAC  .20  .20  .25  .25  .25  .25  .25  .25  .20  .20  .20  .20
  1036.         VCNX     4.   4.   4.   4.   4.   4.   4.   4.   4.   4.   4.   4.
  1037.  
  1038.      Urban - high biogenic (controls present)
  1039.         BCKPMF  60.  60.  60.  60.  60.  60.  60.  60.  60.  60.  60.  60.
  1040.         OFRAC  .25  .25  .30  .30  .30  .55  .55  .55  .35  .35  .35  .25
  1041.         VCNX    15.  15.  15.  15.  15.  15.  15.  15.  15.  15.  15.  15.
  1042.  
  1043.      Regional Plume
  1044.         BCKPMF  20.  20.  20.  20.  20.  20.  20.  20.  20.  20.  20.  20.
  1045.         OFRAC  .20  .20  .25  .35  .25  .40  .40  .40  .30  .30  .30  .20
  1046.         VCNX    15.  15.  15.  15.  15.  15.  15.  15.  15.  15.  15.  15.
  1047.  
  1048.      Urban - no controls present
  1049.         BCKPMF 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100.
  1050.         OFRAC  .30  .30  .35  .35  .35  .55  .55  .55  .35  .35  .35  .30
  1051.         VCNX     2.   2.   2.   2.   2.   2.   2.   2.   2.   2.   2.   2.
  1052.  
  1053.      Default: Clean Continental
  1054.      !  BCKPMF = 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00 !
  1055.      !  OFRAC  = 0.15, 0.15, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.15 !
  1056.      !  VCNX   = 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00, 50.00 !
  1057.  
  1058.  
  1059. !END!
  1060.  
  1061.  
  1062. -------------------------------------------------------------------------------
  1063.  
  1064.  
  1065. INPUT GROUP: 12 -- Misc. Dispersion and Computational Parameters
  1066. ---------------
  1067.  
  1068.      Horizontal size of puff (m) beyond which
  1069.      time-dependent dispersion equations (Heffter)
  1070.      are used to determine sigma-y and
  1071.      sigma-z (SYTDEP)                           Default: 550.   ! SYTDEP = 5.5E02 !
  1072.  
  1073.      Switch for using Heffter equation for sigma z          
  1074.      as above (0 = Not use Heffter; 1 = use Heffter
  1075.      (MHFTSZ)                                   Default: 0      ! MHFTSZ =  0   !
  1076.  
  1077.      Stability class used to determine plume
  1078.      growth rates for puffs above the boundary
  1079.      layer (JSUP)                               Default: 5      ! JSUP =  5   !
  1080.  
  1081.      Vertical dispersion constant for stable
  1082.      conditions (k1 in Eqn. 2.7-3)  (CONK1)     Default: 0.01   ! CONK1 = .01 !
  1083.  
  1084.      Vertical dispersion constant for neutral/
  1085.      unstable conditions (k2 in Eqn. 2.7-4)
  1086.      (CONK2)                                    Default: 0.1    ! CONK2 = .1 !
  1087.  
  1088.      Factor for determining Transition-point from
  1089.      Schulman-Scire to Huber-Snyder Building Downwash
  1090.      scheme (SS used for Hs < Hb + TBD * HL)
  1091.      (TBD)                                      Default: 0.5    ! TBD = .5 !
  1092.         TBD < 0   ==> always use Huber-Snyder
  1093.         TBD = 1.5 ==> always use Schulman-Scire
  1094.         TBD = 0.5 ==> ISC Transition-point
  1095.  
  1096.      Range of land use categories for which
  1097.      urban dispersion is assumed
  1098.      (IURB1, IURB2)                             Default: 10     ! IURB1 =  1  !
  1099.                                                          19     ! IURB2 =  1  !
  1100.  
  1101.      Site characterization parameters for single-point Met data files ---------
  1102.      (needed for METFM = 2,3,4,5)
  1103.  
  1104.         Land use category for modeling domain
  1105.         (ILANDUIN)                              Default: 20     ! ILANDUIN =  20  !
  1106.  
  1107.         Roughness length (m) for modeling domain
  1108.         (Z0IN)                                  Default: 0.25   ! Z0IN = .25 !
  1109.  
  1110.         Leaf area index for modeling domain
  1111.         (XLAIIN)                                Default: 3.0    ! XLAIIN = 3.0 !
  1112.  
  1113.         Elevation above sea level (m)
  1114.         (ELEVIN)                                Default: 0.0    ! ELEVIN = .0 !
  1115.  
  1116.         Latitude (degrees) for met location
  1117.         (XLATIN)                                Default: -999.  ! XLATIN = .0 !
  1118.  
  1119.         Longitude (degrees) for met location
  1120.         (XLONIN)                                Default: -999.  ! XLONIN = .0 !
  1121.  
  1122.      Specialized information for interpreting single-point Met data files -----
  1123.  
  1124.         Anemometer height (m) (Used only if METFM = 2,3)
  1125.         (ANEMHT)                                Default: 10.    ! ANEMHT = 10.0 !
  1126.  
  1127.         Form of lateral turbulance data in PROFILE.DAT file
  1128.         (Used only if METFM = 4,5 or MTURBVW = 1 or 3)
  1129.         (ISIGMAV)                               Default: 1      ! ISIGMAV =  1  !
  1130.             0 = read sigma-theta
  1131.             1 = read sigma-v
  1132.  
  1133.         Choice of mixing heights (Used only if METFM = 4)
  1134.         (IMIXCTDM)                              Default: 0      ! IMIXCTDM =  0  !
  1135.             0 = read PREDICTED mixing heights
  1136.             1 = read OBSERVED mixing heights
  1137.  
  1138.      Maximum length of a slug (met. grid units)
  1139.      (XMXLEN)                                   Default: 1.0    ! XMXLEN = 1.0 !
  1140.  
  1141.      Maximum travel distance of a puff/slug (in
  1142.      grid units) during one sampling step
  1143.      (XSAMLEN)                                  Default: 1.0    ! XSAMLEN = 1.0 !
  1144.  
  1145.      Maximum Number of slugs/puffs release from
  1146.      one source during one time step            
  1147.      (MXNEW)                                    Default: 99     ! MXNEW =  99   !
  1148.  
  1149.      Maximum Number of sampling steps for    
  1150.      one puff/slug during one time step            
  1151.      (MXSAM)                                    Default: 99     ! MXSAM =  99   !
  1152.  
  1153.      Number of iterations used when computing
  1154.      the transport wind for a sampling step
  1155.      that includes gradual rise (for CALMET
  1156.      and PROFILE winds)
  1157.      (NCOUNT)                                   Default: 2      ! NCOUNT =  2   !
  1158.  
  1159.      Minimum sigma y for a new puff/slug (m)      
  1160.      (SYMIN)                                    Default: 1.0    ! SYMIN = 1.0  !
  1161.  
  1162.      Minimum sigma z for a new puff/slug (m)    
  1163.      (SZMIN)                                    Default: 1.0    ! SZMIN = 1.0  !
  1164.  
  1165.      Default minimum turbulence velocities sigma-v and sigma-w
  1166.      for each stability class over land and over water (m/s)
  1167.      (SVMIN(12) and SWMIN(12))
  1168.  
  1169.                      ----------  LAND  ----------       ---------  WATER  ----------
  1170.         Stab Class :  A    B    C    D    E    F         A    B    C    D    E    F
  1171.                      ---  ---  ---  ---  ---  ---       ---  ---  ---  ---  ---  ---
  1172.      Default SVMIN : .50, .50, .50, .50, .50, .50,      .37, .37, .37, .37, .37, .37
  1173.      Default SWMIN : .20, .12, .08, .06, .03, .016,     .20, .12, .08, .06, .03, .016
  1174.  
  1175.            ! SVMIN = 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.370, 0.370, 0.370, 0.370, 0.370, 0.370!
  1176.            ! SWMIN = 0.200, 0.120, 0.080, 0.060, 0.030, 0.016, 0.200, 0.120, 0.080, 0.060, 0.030, 0.016!
  1177.  
  1178.      Divergence criterion for dw/dz across puff
  1179.      used to initiate adjustment for horizontal
  1180.      convergence (1/s)
  1181.      Partial adjustment starts at CDIV(1), and
  1182.      full adjustment is reached at CDIV(2)
  1183.      (CDIV(2))                                  Default: 0.0,0.0  ! CDIV = .0, .0 !
  1184.  
  1185.      Minimum wind speed (m/s) allowed for
  1186.      non-calm conditions. Also used as minimum
  1187.      speed returned when using power-law
  1188.      extrapolation toward surface
  1189.      (WSCALM)                                   Default: 0.5    ! WSCALM = .5 !
  1190.  
  1191.      Maximum mixing height (m)                      
  1192.      (XMAXZI)                                   Default: 3000.  ! XMAXZI = 3000.0 !
  1193.  
  1194.      Minimum mixing height (m)                    
  1195.      (XMINZI)                                   Default: 50.    ! XMINZI = 20.0 !
  1196.  
  1197.      Default wind speed classes --
  1198.      5 upper bounds (m/s) are entered;
  1199.      the 6th class has no upper limit
  1200.      (WSCAT(5))                      Default   :
  1201.                                      ISC RURAL : 1.54, 3.09, 5.14, 8.23, 10.8 (10.8+)
  1202.  
  1203.                               Wind Speed Class :  1     2     3     4     5  
  1204.                                                  ---   ---   ---   ---   ---
  1205.                                        ! WSCAT = 1.54, 3.09, 5.14, 8.23, 10.80 !
  1206.  
  1207.      Default wind speed profile power-law
  1208.      exponents for stabilities 1-6
  1209.      (PLX0(6))                       Default   : ISC RURAL values
  1210.                                      ISC RURAL : .07, .07, .10, .15, .35, .55
  1211.                                      ISC URBAN : .15, .15, .20, .25, .30, .30
  1212.  
  1213.                                Stability Class :  A     B     C     D     E     F
  1214.                                                  ---   ---   ---   ---   ---   ---
  1215.                                         ! PLX0 = 0.07, 0.07, 0.10, 0.15, 0.35, 0.55 !
  1216.  
  1217.      Default potential temperature gradient
  1218.      for stable classes E, F (degK/m)
  1219.      (PTG0(2))                       Default: 0.020, 0.035
  1220.                                         ! PTG0 = 0.020,   0.035 !
  1221.  
  1222.      Default plume path coefficients for
  1223.      each stability class (used when option
  1224.      for partial plume height terrain adjustment
  1225.      is selected -- MCTADJ=3)
  1226.      (PPC(6))                  Stability Class :  A     B     C     D     E     F
  1227.                                   Default  PPC : .50,  .50,  .50,  .50,  .35,  .35
  1228.                                                  ---   ---   ---   ---   ---   ---
  1229.                                         !  PPC = 0.50, 0.50, 0.50, 0.50, 0.35, 0.35 !
  1230.  
  1231.      Slug-to-puff transition criterion factor
  1232.      equal to sigma-y/length of slug
  1233.      (SL2PF)                               Default: 10.        ! SL2PF = 10.0 !
  1234.  
  1235.      Puff-splitting control variables ------------------------
  1236.  
  1237.        VERTICAL SPLIT
  1238.        --------------
  1239.  
  1240.        Number of puffs that result every time a puff
  1241.        is split - nsplit=2 means that 1 puff splits
  1242.        into 2
  1243.        (NSPLIT)                            Default:   3        ! NSPLIT =  3  !
  1244.  
  1245.        Time(s) of a day when split puffs are eligible to
  1246.        be split once again; this is typically set once
  1247.        per day, around sunset before nocturnal shear develops.
  1248.        24 values: 0 is midnight (00:00) and 23 is 11 PM (23:00)
  1249.        0=do not re-split    1=eligible for re-split
  1250.        (IRESPLIT(24))                      Default:  Hour 17 = 1
  1251.        !  IRESPLIT = 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0 !
  1252.  
  1253.        Split is allowed only if last hour's mixing
  1254.        height (m) exceeds a minimum value
  1255.        (ZISPLIT)                           Default: 100.       ! ZISPLIT = 100.0 !
  1256.  
  1257.        Split is allowed only if ratio of last hour's
  1258.        mixing ht to the maximum mixing ht experienced
  1259.        by the puff is less than a maximum value (this
  1260.        postpones a split until a nocturnal layer develops)
  1261.        (ROLDMAX)                           Default: 0.25       ! ROLDMAX = 0.25 !
  1262.  
  1263.  
  1264.        HORIZONTAL SPLIT
  1265.        ----------------
  1266.  
  1267.        Number of puffs that result every time a puff
  1268.        is split - nsplith=5 means that 1 puff splits
  1269.        into 5
  1270.        (NSPLITH)                           Default:   5        ! NSPLITH =  5  !
  1271.  
  1272.        Minimum sigma-y (Grid Cells Units) of puff
  1273.        before it may be split
  1274.        (SYSPLITH)                          Default:  1.0       ! SYSPLITH = 1.0 !
  1275.  
  1276.        Minimum puff elongation rate (SYSPLITH/hr) due to
  1277.        wind shear, before it may be split
  1278.        (SHSPLITH)                          Default:  2.        ! SHSPLITH = 2.0 !
  1279.  
  1280.        Minimum concentration (g/m^3) of each
  1281.        species in puff before it may be split
  1282.        Enter array of NSPEC values; if a single value is
  1283.        entered, it will be used for ALL species
  1284.        (CNSPLITH)                          Default:  1.0E-07   ! CNSPLITH = 1.0E-07 !
  1285.  
  1286.      Integration control variables ------------------------
  1287.  
  1288.        Fractional convergence criterion for numerical SLUG
  1289.        sampling integration
  1290.        (EPSSLUG)                           Default:   1.0e-04  ! EPSSLUG = 1.0E-04 !
  1291.  
  1292.        Fractional convergence criterion for numerical AREA
  1293.        source integration
  1294.        (EPSAREA)                           Default:   1.0e-06  ! EPSAREA = 1.0E-06 !
  1295.  
  1296.        Trajectory step-length (m) used for numerical rise
  1297.        integration
  1298.        (DSRISE)                            Default:   1.0      ! DSRISE = 1.0 !
  1299.  
  1300.        Boundary Condition (BC) Puff control variables ------------------------
  1301.  
  1302.        Minimum height (m) to which BC puffs are mixed as they are emitted
  1303.        (MBCON=2 ONLY).  Actual height is reset to the current mixing height
  1304.        at the release point if greater than this minimum.
  1305.        (HTMINBC)                           Default:   500.     ! HTMINBC = 500.0 !
  1306.  
  1307.        Search radius (km) about a receptor for sampling nearest BC puff.
  1308.        BC puffs are typically emitted with a spacing of one grid cell
  1309.        length, so the search radius should be greater than DGRIDKM.
  1310.        (RSAMPBC)                           Default:   10.      ! RSAMPBC = 10.0 !
  1311.  
  1312.        Near-Surface depletion adjustment to concentration profile used when
  1313.        sampling BC puffs?
  1314.        (MDEPBC)                            Default:   1        ! MDEPBC =  1  !
  1315.           0 = Concentration is NOT adjusted for depletion
  1316.           1 = Adjust Concentration for depletion
  1317.  
  1318. !END!
  1319.  
  1320.  
  1321. -------------------------------------------------------------------------------
  1322.  
  1323.  
  1324. INPUT GROUPS: 13a, 13b, 13c, 13d -- Point source parameters
  1325. --------------------------------
  1326.  
  1327. ---------------
  1328. Subgroup (13a)
  1329. ---------------
  1330.  
  1331.      Number of point sources with
  1332.      parameters provided below      (NPT1)  No default  !  NPT1 =  1  !
  1333.  
  1334.      Units used for point source
  1335.      emissions below                (IPTU)  Default: 1  !  IPTU =   1  !
  1336.            1 =        g/s
  1337.            2 =       kg/hr
  1338.            3 =       lb/hr
  1339.            4 =     tons/yr
  1340.            5 =     Odour Unit * m**3/s  (vol. flux of odour compound)
  1341.            6 =     Odour Unit * m**3/min
  1342.            7 =     metric tons/yr
  1343.  
  1344.      Number of source-species
  1345.      combinations with variable
  1346.      emissions scaling factors
  1347.      provided below in (13d)        (NSPT1) Default: 0  !  NSPT1 =  0  !
  1348.  
  1349.      Number of point sources with
  1350.      variable emission parameters
  1351.      provided in external file      (NPT2)  No default  !  NPT2 =  0  !
  1352.  
  1353.      (If NPT2 > 0, these point
  1354.      source emissions are read from
  1355.      the file: PTEMARB.DAT)
  1356.  
  1357. !END!
  1358.  
  1359. ---------------
  1360. Subgroup (13b)
  1361. ---------------
  1362.                                       a
  1363.           POINT SOURCE: CONSTANT DATA
  1364.           -----------------------------
  1365.                                                                               b          c
  1366.   Source       X         Y       Stack    Base     Stack    Exit  Exit    Bldg.  Emission
  1367.    No.     Coordinate Coordinate Height Elevation Diameter  Vel.  Temp.   Dwash   Rates
  1368.               (km)      (km)       (m)      (m)       (m)  (m/s) (deg. K)        
  1369.   ------   ---------- ---------- ------  ------   -------- ----- -------- ----- --------
  1370.    1 ! SRCNAM = STK1 !
  1371.    1 ! X =        1.0,   1.0,    40.0, 160.0,       3.0,   5.0,  355.0,   .0, 1.0E01,  0.0E00,  4.0E00,  
  1372.        1.0E00,  0.0E00,  0.0E00,  1.0E01 !
  1373.    1 ! ZPLTFM  =       .0 !
  1374.    1 ! FMFAC  =      1.0 !   !END!
  1375.  
  1376. --------
  1377.  
  1378.     a
  1379.      Data for each source are treated as a separate input subgroup
  1380.      and therefore must end with an input group terminator.
  1381.  
  1382.      SRCNAM  is a 12-character name for a source
  1383.              (No default)
  1384.      X       is an array holding the source data listed by the column headings
  1385.              (No default)
  1386.      SIGYZI  is an array holding the initial sigma-y and sigma-z (m)
  1387.              (Default: 0.,0.)
  1388.      FMFAC   is a vertical momentum flux factor (0. or 1.0) used to represent
  1389.              the effect of rain-caps or other physical configurations that
  1390.              reduce momentum rise associated with the actual exit velocity.
  1391.              (Default: 1.0  -- full momentum used)
  1392.      ZPLTFM  is the platform height (m) for sources influenced by an isolated
  1393.              structure that has a significant open area between the surface
  1394.              and the bulk of the structure, such as an offshore oil platform.
  1395.              The Base Elevation is that of the surface (ground or ocean),
  1396.              and the Stack Height is the release height above the Base (not
  1397.              above the platform).  Building heights entered in Subgroup 13c
  1398.              must be those of the buildings on the platform, measured from
  1399.              the platform deck.  ZPLTFM is used only with MBDW=1 (ISC
  1400.              downwash method) for sources with building downwash.
  1401.              (Default: 0.0)
  1402.  
  1403.     b
  1404.      0. = No building downwash modeled
  1405.      1. = Downwash modeled for buildings resting on the surface
  1406.      2. = Downwash modeled for buildings raised above the surface (ZPLTFM > 0.)
  1407.      NOTE: must be entered as a REAL number (i.e., with decimal point)
  1408.  
  1409.     c
  1410.      An emission rate must be entered for every pollutant modeled.
  1411.      Enter emission rate of zero for secondary pollutants that are
  1412.      modeled, but not emitted.  Units are specified by IPTU
  1413.      (e.g. 1 for g/s).
  1414.  
  1415. ---------------
  1416. Subgroup (13c)
  1417. ---------------
  1418.  
  1419.            BUILDING DIMENSION DATA FOR SOURCES SUBJECT TO DOWNWASH
  1420.            -------------------------------------------------------
  1421. Source                                                                     a
  1422.  No.       Effective building height, width, length and X/Y offset (in meters)
  1423.            every 10 degrees.  LENGTH, XBADJ, and YBADJ are only needed for
  1424.            MBDW=2 (PRIME downwash option)
  1425. ------     --------------------------------------------------------------------
  1426.  
  1427.  
  1428. --------
  1429.  
  1430.     a
  1431.      Building height, width, length, and X/Y offset from the source are treated
  1432.      as a separate input subgroup for each source and therefore must end with
  1433.      an input group terminator.  The X/Y offset is the position, relative to the
  1434.      stack, of the center of the upwind face of the projected building, with the
  1435.      x-axis pointing along the flow direction.
  1436.  
  1437. ---------------
  1438. Subgroup (13d)
  1439. ---------------
  1440.                                                 a
  1441.           POINT SOURCE: VARIABLE EMISSIONS DATA
  1442.           ---------------------------------------
  1443.  
  1444.      Use this subgroup to describe temporal variations in the emission
  1445.      rates given in 13b.  Factors entered multiply the rates in 13b.
  1446.      Skip sources here that have constant emissions.  For more elaborate
  1447.      variation in source parameters, use PTEMARB.DAT and NPT2 > 0.
  1448.  
  1449.      IVARY determines the type of variation, and is source-specific:
  1450.      (IVARY)                                Default: 0
  1451.            0 =       Constant
  1452.            1 =       Diurnal cycle (24 scaling factors: hours 1-24)
  1453.            2 =       Monthly cycle (12 scaling factors: months 1-12)
  1454.            3 =       Hour & Season (4 groups of 24 hourly scaling factors,
  1455.                                     where first group is DEC-JAN-FEB)
  1456.            4 =       Speed & Stab. (6 groups of 6 scaling factors, where
  1457.                                     first group is Stability Class A,
  1458.                                     and the speed classes have upper
  1459.                                     bounds (m/s) defined in Group 12
  1460.            5 =       Temperature   (12 scaling factors, where temperature
  1461.                                     classes have upper bounds (C) of:
  1462.                                     0, 5, 10, 15, 20, 25, 30, 35, 40,
  1463.                                     45, 50, 50+)
  1464.  
  1465.  
  1466.  
  1467. --------
  1468.     a
  1469.      Data for each species are treated as a separate input subgroup
  1470.      and therefore must end with an input group terminator.
  1471.  
  1472.  
  1473. -------------------------------------------------------------------------------
  1474.  
  1475.  
  1476. INPUT GROUPS: 14a, 14b, 14c, 14d -- Area source parameters
  1477. --------------------------------
  1478.  
  1479. ---------------
  1480. Subgroup (14a)
  1481. ---------------
  1482.  
  1483.      Number of polygon area sources with
  1484.      parameters specified below (NAR1)       No default  !  NAR1 =  0   !
  1485.  
  1486.      Units used for area source
  1487.      emissions below            (IARU)       Default: 1  !  IARU =   1  !
  1488.            1 =        g/m**2/s
  1489.            2 =       kg/m**2/hr
  1490.            3 =       lb/m**2/hr
  1491.            4 =     tons/m**2/yr
  1492.            5 =     Odour Unit * m/s  (vol. flux/m**2 of odour compound)
  1493.            6 =     Odour Unit * m/min
  1494.            7 =     metric tons/m**2/yr
  1495.  
  1496.      Number of source-species
  1497.      combinations with variable
  1498.      emissions scaling factors
  1499.      provided below in (14d)        (NSAR1) Default: 0  !  NSAR1 =  0  !
  1500.  
  1501.      Number of buoyant polygon area sources
  1502.      with variable location and emission
  1503.      parameters (NAR2)                      No default  !  NAR2 =  0   !
  1504.      (If NAR2 > 0, ALL parameter data for
  1505.      these sources are read from the file: BAEMARB.DAT)
  1506.  
  1507. !END!
  1508.  
  1509. ---------------
  1510. Subgroup (14b)
  1511. ---------------
  1512.                                      a
  1513.           AREA SOURCE: CONSTANT DATA
  1514.           ----------------------------
  1515.                                                          b
  1516. Source           Effect.    Base      Initial    Emission
  1517.  No.             Height   Elevation   Sigma z     Rates
  1518.                    (m)       (m)        (m)      
  1519. -------          ------    ------     --------   ---------
  1520.  
  1521.  
  1522. --------
  1523.     a
  1524.      Data for each source are treated as a separate input subgroup
  1525.      and therefore must end with an input group terminator.
  1526.     b
  1527.      An emission rate must be entered for every pollutant modeled.
  1528.      Enter emission rate of zero for secondary pollutants that are
  1529.      modeled, but not emitted.  Units are specified by IARU
  1530.      (e.g. 1 for g/m**2/s).
  1531.  
  1532. ---------------
  1533. Subgroup (14c)
  1534. ---------------
  1535.  
  1536.            COORDINATES (km) FOR EACH VERTEX(4) OF EACH POLYGON
  1537.            --------------------------------------------------------
  1538. Source                                                               a
  1539.  No.       Ordered list of X followed by list of Y, grouped by source
  1540. ------     ------------------------------------------------------------
  1541.  
  1542.  
  1543. --------
  1544.     a
  1545.      Data for each source are treated as a separate input subgroup
  1546.      and therefore must end with an input group terminator.
  1547.  
  1548.  
  1549. ---------------
  1550. Subgroup (14d)
  1551. ---------------
  1552.                                                a
  1553.           AREA SOURCE: VARIABLE EMISSIONS DATA
  1554.           --------------------------------------
  1555.  
  1556.      Use this subgroup to describe temporal variations in the emission
  1557.      rates given in 14b.  Factors entered multiply the rates in 14b.
  1558.      Skip sources here that have constant emissions.  For more elaborate
  1559.      variation in source parameters, use BAEMARB.DAT and NAR2 > 0.
  1560.  
  1561.      IVARY determines the type of variation, and is source-specific:
  1562.      (IVARY)                                Default: 0
  1563.            0 =       Constant
  1564.            1 =       Diurnal cycle (24 scaling factors: hours 1-24)
  1565.            2 =       Monthly cycle (12 scaling factors: months 1-12)
  1566.            3 =       Hour & Season (4 groups of 24 hourly scaling factors,
  1567.                                     where first group is DEC-JAN-FEB)
  1568.            4 =       Speed & Stab. (6 groups of 6 scaling factors, where
  1569.                                     first group is Stability Class A,
  1570.                                     and the speed classes have upper
  1571.                                     bounds (m/s) defined in Group 12
  1572.            5 =       Temperature   (12 scaling factors, where temperature
  1573.                                     classes have upper bounds (C) of:
  1574.                                     0, 5, 10, 15, 20, 25, 30, 35, 40,
  1575.                                     45, 50, 50+)
  1576.  
  1577.  
  1578.  
  1579. --------
  1580.     a
  1581.      Data for each species are treated as a separate input subgroup
  1582.      and therefore must end with an input group terminator.
  1583.  
  1584.  
  1585. -------------------------------------------------------------------------------
  1586.  
  1587. INPUT GROUPS: 15a, 15b, 15c -- Line source parameters
  1588. ---------------------------
  1589.  
  1590. ---------------
  1591. Subgroup (15a)
  1592. ---------------
  1593.  
  1594.      Number of buoyant line sources
  1595.      with variable location and emission
  1596.      parameters (NLN2)                              No default  !  NLN2 =  0   !
  1597.  
  1598.      (If NLN2 > 0, ALL parameter data for
  1599.       these sources are read from the file: LNEMARB.DAT)
  1600.  
  1601.      Number of buoyant line sources (NLINES)        No default   ! NLINES =  0  !
  1602.  
  1603.      Units used for line source
  1604.      emissions below                (ILNU)          Default: 1  !  ILNU =   1  !
  1605.            1 =        g/s
  1606.            2 =       kg/hr
  1607.            3 =       lb/hr
  1608.            4 =     tons/yr
  1609.            5 =     Odour Unit * m**3/s  (vol. flux of odour compound)
  1610.            6 =     Odour Unit * m**3/min
  1611.            7 =     metric tons/yr
  1612.  
  1613.      Number of source-species
  1614.      combinations with variable
  1615.      emissions scaling factors
  1616.      provided below in (15c)        (NSLN1) Default: 0  !  NSLN1 =  0  !
  1617.  
  1618.      Maximum number of segments used to model
  1619.      each line (MXNSEG)                             Default: 7   ! MXNSEG =  7  !
  1620.  
  1621.      The following variables are required only if NLINES > 0.  They are
  1622.      used in the buoyant line source plume rise calculations.
  1623.  
  1624.         Number of distances at which                Default: 6   ! NLRISE =  6  !
  1625.         transitional rise is computed
  1626.  
  1627.         Average building length (XL)                No default   ! XL = .0 !
  1628.                                                     (in meters)
  1629.  
  1630.         Average building height (HBL)               No default   ! HBL = .0 !
  1631.                                                     (in meters)
  1632.  
  1633.         Average building width (WBL)                No default   ! WBL = .0 !
  1634.                                                     (in meters)
  1635.  
  1636.         Average line source width (WML)             No default   ! WML = .0 !
  1637.                                                     (in meters)
  1638.  
  1639.         Average separation between buildings (DXL)  No default   ! DXL = .0 !
  1640.                                                     (in meters)
  1641.  
  1642.         Average buoyancy parameter (FPRIMEL)        No default   ! FPRIMEL = .0 !
  1643.                                                     (in m**4/s**3)
  1644.  
  1645. !END!
  1646.  
  1647. ---------------
  1648. Subgroup (15b)
  1649. ---------------
  1650.  
  1651.           BUOYANT LINE SOURCE: CONSTANT DATA
  1652.           ----------------------------------
  1653.                                                                                           a
  1654. Source     Beg. X      Beg. Y      End. X    End. Y     Release    Base        Emission
  1655.  No.     Coordinate  Coordinate  Coordinate Coordinate  Height    Elevation      Rates
  1656.             (km)        (km)        (km)       (km)       (m)       (m)          
  1657. ------   ----------  ----------  ---------  ----------  -------   ---------    ---------
  1658.  
  1659. --------
  1660.  
  1661.     a
  1662.      Data for each source are treated as a separate input subgroup
  1663.      and therefore must end with an input group terminator.
  1664.  
  1665.     b
  1666.      An emission rate must be entered for every pollutant modeled.
  1667.      Enter emission rate of zero for secondary pollutants that are
  1668.      modeled, but not emitted.  Units are specified by ILNTU
  1669.      (e.g. 1 for g/s).
  1670.  
  1671. ---------------
  1672. Subgroup (15c)
  1673. ---------------
  1674.                                                        a
  1675.           BUOYANT LINE SOURCE: VARIABLE EMISSIONS DATA
  1676.           ----------------------------------------------
  1677.  
  1678.      Use this subgroup to describe temporal variations in the emission
  1679.      rates given in 15b.  Factors entered multiply the rates in 15b.
  1680.      Skip sources here that have constant emissions.
  1681.  
  1682.      IVARY determines the type of variation, and is source-specific:
  1683.      (IVARY)                                Default: 0
  1684.            0 =       Constant
  1685.            1 =       Diurnal cycle (24 scaling factors: hours 1-24)
  1686.            2 =       Monthly cycle (12 scaling factors: months 1-12)
  1687.            3 =       Hour & Season (4 groups of 24 hourly scaling factors,
  1688.                                     where first group is DEC-JAN-FEB)
  1689.            4 =       Speed & Stab. (6 groups of 6 scaling factors, where
  1690.                                     first group is Stability Class A,
  1691.                                     and the speed classes have upper
  1692.                                     bounds (m/s) defined in Group 12
  1693.            5 =       Temperature   (12 scaling factors, where temperature
  1694.                                     classes have upper bounds (C) of:
  1695.                                     0, 5, 10, 15, 20, 25, 30, 35, 40,
  1696.                                     45, 50, 50+)
  1697.  
  1698.  
  1699.  
  1700. --------
  1701.     a
  1702.      Data for each species are treated as a separate input subgroup
  1703.      and therefore must end with an input group terminator.
  1704.  
  1705.  
  1706. -------------------------------------------------------------------------------
  1707.  
  1708.  
  1709. INPUT GROUPS: 16a, 16b, 16c -- Volume source parameters
  1710. ---------------------------
  1711.  
  1712. ---------------
  1713. Subgroup (16a)
  1714. ---------------
  1715.  
  1716.      Number of volume sources with
  1717.      parameters provided in 16b,c (NVL1)     No default  !  NVL1 =  0   !
  1718.  
  1719.      Units used for volume source
  1720.      emissions below in 16b       (IVLU)     Default: 1  !  IVLU =   1  !
  1721.            1 =        g/s
  1722.            2 =       kg/hr
  1723.            3 =       lb/hr
  1724.            4 =     tons/yr
  1725.            5 =     Odour Unit * m**3/s  (vol. flux of odour compound)
  1726.            6 =     Odour Unit * m**3/min
  1727.            7 =     metric tons/yr
  1728.  
  1729.      Number of source-species
  1730.      combinations with variable
  1731.      emissions scaling factors
  1732.      provided below in (16c)      (NSVL1)    Default: 0  !  NSVL1 =  0  !
  1733.  
  1734.      Number of volume sources with
  1735.      variable location and emission
  1736.      parameters                   (NVL2)     No default  !  NVL2 =   0   !
  1737.  
  1738.      (If NVL2 > 0, ALL parameter data for
  1739.       these sources are read from the VOLEMARB.DAT file(s) )
  1740.  
  1741. !END!
  1742.  
  1743. ---------------
  1744. Subgroup (16b)
  1745. ---------------
  1746.                                         a
  1747.            VOLUME SOURCE: CONSTANT DATA
  1748.            ------------------------------
  1749.                                                                                b
  1750.          X           Y        Effect.    Base     Initial    Initial    Emission
  1751.      Coordinate  Coordinate   Height   Elevation  Sigma y    Sigma z     Rates
  1752.         (km)       (km)         (m)       (m)        (m)       (m)      
  1753.      ----------  ----------   ------    ------    --------   --------   --------
  1754.  
  1755.  
  1756. --------
  1757.     a
  1758.      Data for each source are treated as a separate input subgroup
  1759.      and therefore must end with an input group terminator.
  1760.  
  1761.     b
  1762.      An emission rate must be entered for every pollutant modeled.
  1763.      Enter emission rate of zero for secondary pollutants that are
  1764.      modeled, but not emitted.  Units are specified by IVLU
  1765.      (e.g. 1 for g/s).
  1766.  
  1767. ---------------
  1768. Subgroup (16c)
  1769. ---------------
  1770.                                                  a
  1771.           VOLUME SOURCE: VARIABLE EMISSIONS DATA
  1772.           ----------------------------------------
  1773.  
  1774.      Use this subgroup to describe temporal variations in the emission
  1775.      rates given in 16b.  Factors entered multiply the rates in 16b.
  1776.      Skip sources here that have constant emissions.  For more elaborate
  1777.      variation in source parameters, use VOLEMARB.DAT and NVL2 > 0.
  1778.  
  1779.      IVARY determines the type of variation, and is source-specific:
  1780.      (IVARY)                                Default: 0
  1781.            0 =       Constant
  1782.            1 =       Diurnal cycle (24 scaling factors: hours 1-24)
  1783.            2 =       Monthly cycle (12 scaling factors: months 1-12)
  1784.            3 =       Hour & Season (4 groups of 24 hourly scaling factors,
  1785.                                     where first group is DEC-JAN-FEB)
  1786.            4 =       Speed & Stab. (6 groups of 6 scaling factors, where
  1787.                                     first group is Stability Class A,
  1788.                                     and the speed classes have upper
  1789.                                     bounds (m/s) defined in Group 12
  1790.            5 =       Temperature   (12 scaling factors, where temperature
  1791.                                     classes have upper bounds (C) of:
  1792.                                     0, 5, 10, 15, 20, 25, 30, 35, 40,
  1793.                                     45, 50, 50+)
  1794.  
  1795.  
  1796.  
  1797. --------
  1798.     a
  1799.      Data for each species are treated as a separate input subgroup
  1800.      and therefore must end with an input group terminator.
  1801.  
  1802.  
  1803. -------------------------------------------------------------------------------
  1804.  
  1805. INPUT GROUPS: 17a & 17b -- Non-gridded (discrete) receptor information
  1806. -----------------------
  1807.  
  1808. ---------------
  1809. Subgroup (17a)
  1810. ---------------
  1811.  
  1812.      Number of non-gridded receptors (NREC)  No default  !  NREC =  0   !
  1813.  
  1814. !END!
  1815.  
  1816. ---------------
  1817. Subgroup (17b)
  1818. ---------------
  1819.                                                a
  1820.            NON-GRIDDED (DISCRETE) RECEPTOR DATA
  1821.            ------------------------------------
  1822.  
  1823.                    X            Y          Ground        Height   b
  1824. Receptor       Coordinate   Coordinate    Elevation   Above Ground
  1825.   No.             (km)         (km)          (m)           (m)
  1826. --------       ----------   ----------    ---------   ------------
  1827.  
  1828.  
  1829.  
  1830. -------------
  1831.     a
  1832.      Data for each receptor are treated as a separate input subgroup
  1833.      and therefore must end with an input group terminator.
  1834.  
  1835.     b
  1836.      Receptor height above ground is optional.  If no value is entered,
  1837.      the receptor is placed on the ground.
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top