United Statas Office of Air Quality EPA—450/4—90—007E
Environmental Planning and Standards JUNE 1990
Protection Research Triangle Park, NC 27711
Agency
AIR
E PA USER'S GUIDE FOR THE
URBAN AIRSHED MODEL
Volume V: Description and Operation of
the ROM - UAM Interface Program System
-------
TECHNICAL REPORT DATA
(Ptease read Instructions on the reverse before completing)
1. REPORT NO. 2.
EPA-450/4-90-007E
3. RECIPIENT'S ACCESSION NO.
4. TITLE ANO SUBTITLE
USER'S GUIDE FOR THE URBAN AIRSHED MODEL
Volume V: Description and Operation of the ROM-UAM
Interface Proqram System
5. REPORT DATE
June 1990
6. PERFORMING ORGANIZATION COOE
7. AUTHOR(S)
R. T. Tang, S. C. Gerry, J. S. Newsome,
a R Vsnmptpr, and R. A. Wavland
8. PERFORMING ORGANIZATION REPORT NO.
9. PERFORMING ORGANIZATION NAME ANO AOORESS
Computer Science Corporation
Research Triangle Park, N. C. 27709
10. PROGRAM ELEMENT NO.
11. CONTRACT/GRANT NO.
12. sponsor.^ VW^lfftdT0Pmect ion Agency
Office of Air Quality Planning and Standards
Research Triangle Park, N. C. 277711
13. TYPE OF REPORT ANO PERIOO COVEREO
14. SPONSORING AGENCY COOE
15. SUPPLEMENTARY NOTES
16. ABSTRACT
This document serves as a manual for the ROM-UAM Interface System that
enables the use of ROM model outputs to be used as input meteorological and
Boundary/initial condition processors.
i
: 7. KEY WOROS ANO OOCUMENT ANALYSIS
t. DESCRIPTORS
b.IOENTIFIERS/OPEN ENOED TERMS
c. COS AT I Field/Group
Ozone
Urban Airshed Model
Photochemistry
ROM-UAM Interface System
8. DISTRIBUTION STATEMENT
19. SECURITY CLASS (This Report)
21. NO. OF PAGES
20. SECURITY CLASS (Tfiis page)
22. PRICE
PA Form 2220—1 (R«v. 4—77) previous coition is obsolete
-------
EPA-450/4-90-007E
USER'S GUIDE FOR THE
URBAN AIRSHED MODEL
Volume V: Description and Operation of
the ROM - UAM Interface Program System
By
R. T. Tang
S. C. Gerry
J. S. Newsome
Allan R. Van Meter
Richard A. Wayland
Computer Science Corporation
Research Triangle Park, NC 27709
and
J. M. Godowitch
K. L. Schere
Atmospheric Sciences Modeling Division
Atmospheric Research and Exposure Assessment Laboratory
U. S. Environmental Protection Agency
Research Triangle Park, NC 27711
EPA Project Officer:
Richard D. Scheffe
OFFICE OF AIR QUALITY PLANNING AND STANDARDS
U. S. ENVIRONMENTAL PROTECTION AGENCY
RESEARCH TRIANGLE PARK, NC 27711
JUNE 1990
\
-------
NOTICE
The information in this document has been funded wholly or in part by the U.S. Environmental Protec-
tion Agency (EPA) under contract 68-01-7365 to Computer Sciences Corporation. Mention of trade names or
commercial products does not constitute endorsement or recommendation for use.
Affiliation
James M. Godowitch and Kenneth L. Schere are on assignment from the National Oceanic and Atmo-
spheric Administration (NOAA), U.S. Department of Commerce.
ii
-------
PREFACE
This user's guide for the Urban Airshed Model (UAM) is divided into five volumes as follows:
Volume I-Usefs Manual for UAM(CB-IV)
Volume n-Usefs Manual for the UAM(CB-IV) Modeling System (Preprocessors)
Volume III—User's Manual for the Diagnostic Wind Model
Volume IV-User's Manual for the Emissions Processor System -
Volume V--Description and Operation of the ROM-UAM Interface Program System . /
Volume I provides historical background on the model and describes in general the scientific basis for
the model. It describes the structure of the required unformatted (binary) files that are used directly as input
to the UAM. This volume also presents the formats of the output files and information on how to run an
actual UAM simulation. For those users that already possess a UAM modeling database or have prepared
inputs without the use of the standard UAM preprocessors, this volume should serve as a self-sufficient guide
to running the model.
Volume II describes the file formats and software for each of the standard UAM preprocessors that are
part of the UAM modeling system. The preprocessor input files are ASCII files that are generated from raw
input data (meteorological, air quality, emissions). The preprocessor input files are then read by individual
preprocessor programs to create the unformatted (binary) files that are read directly by the UAM. Included in
this volume is an example problem that illustrates how inputs were created from measurement data for an
application of the UAM in Atlanta. The preprocessors available for generating wind fields and emissions
inventories for the UAM are described separately in Volumes III and IV, respectively.
Volume III is the user's manual for the Diagnostic Wind Model (DWM). This model is a stand-alone
interpolative wind model that uses surface- and upper-level wind observations at selected sites within the
modeling domain of interest to provide hourly, gridded, three-dimensional estimates of winds using objective
techniques. It provides one means of formulating wind fields to the UAM.
Volume IV describes in detail the Emissions Processor System (EPS). This software package is used to
process anthropogenic area and point source emissions for the UAM from countywide average total hydrocar-
bon, nitrogen oxide, and carbon monoxide emissions available from national emissions inventories, such as the
National Emissions Data System or the National Acid Precipitation Assessment Program. An appendix to this
iii
-------
volume describes the Biogenic Emissions Inventory System (BEIS), which can be used to generate gridded,
speciated biogenic emissions. Software for merging the anthropogenic area, mobile, and biogenic emissions
files into UAM input format is also described in this volume.
Volume V describes the ROM-UAM interface program system, a software package that can be used to
generate UAM input files from inputs and outputs provided by the EPA Regional Oxidant Model (ROM).
-------
ACKNOWLEDGEMENTS
Since its initial conception in the early 1970% many individuals have contributed to the development of
the Urban Airshed ModeL This document reflects the latest methodology and software development and
provides a guide for new users of the modeL Based on the past efforts of the original developers of the L'AM
and the authors of the original 1978 user's manual, the first four volumes were written by the following indi-
viduals from Systems Applications, Inc.:
Volume I Ralph E. Morris, Thomas G Myers, and Jay L. Haney
Volume II Ralph E. Morris, Thomas G Myers, Edward L. Carr, Marianne C. Causley, Sharon G.
Douglas, and Jay L. Haney
Volume III Sharon G. Douglas, Robert G Kessler, and Edward L. Carr
Volume IV Marianne G Causley, Julie'L. Fieber, Michele Jimenez, and
LuAnn Gardner
Volume V, containing the ROM-UAM Interface Program Guide, as well as Appendix D in Volume IV
(Biogenics Emissions Inventory System) were written by the following individuals of Computer Sciences Cor-
poration and EPA's Atmospheric Sciences Modeling Division:
Volume V _ Ruen-Tai Tang, Susan C. Geny, Joseph S. Newsom, Allan R. Van Meter, and Richard
A. Wayland (CSC); James M. Godowitch and Kenneth L. Schere (EPA)
The U.S. Environmental Protection Agency and the South Coast Air Quality Management District
provided support for the preparation of this document. Richard D. Scheffe, Ned Meyer, Dennis Doll, and
Ellen Baldridge of the U.S. EPA's Office of Air Quality Planning and Standards contributed to this document
with their insightful technical reviews. Henry Hogo and Tom Chico of the South Coast Air Quality Manage-
ment District also reviewed the documents and provided their comments.
Others at Systems Applications, Inc. that have contributed to the continued development of the UAM
in the last few vears include Dr. Garv Whitten and Mr. Garv Moore.
v
-------
ABSTRACT
A set of computer programs has been developed in order to provide a link between select outputs of the
Regional Oxidant Model (ROM version 2.1) system and the Urban Airshed Model Carbon Bond IV [UAM
(CB-IV)] system. Thirteen individual ROM model and processor output files containing concentrations,
meteorological parameters, surface characteristic information, and biogenic emissions are needed to exercise
the seven interface programs. Output files generated by the interface codes are in compatible formats for
input to either UAM preprocessors or the model. With the exception of diffusion break heights and anthro-
pogenic emissions, the interface system generates all UAM input data from the ROM database files.
Currently, ihe ROM data are available over a domain covering the northeastern states in conjunction with the
Regional Ozone Modeling for Northeast Transport (ROMNET) program. Applications of the ROM for
other regions of the United States are anticipated.
The initial step in the interfacing process is to apply a menu-driven data retrieval program especially
developed within the Gridded Model Information Support System (GMISS). User-specified information
about the UAM domain and simulation period is employed by the retrieval program to automatically extract
the specific ROM outputs needed by the interface codes. The retrieved input data files contain the appropri-
ate spatial and temporal "windows" needed to execuie the interface programs for a particular UAM modeling
application. The ROM database and the data retrieval program reside on the EPA/NCC-IBM 3090 computer
system.
The interface programs have been written in standard FORTRAN-77 code which allows the interface
package to be readily adaptable for execution on most computer systems. Methodologies and procedures
applied to the ROM system outputs by the interface codes in order to generate compatible data files for the
UAM system are described. The interface system is also designed to allow the user to substitute or include
alternate data sets than those from the ROM database. Instructions for the preparation and application of the
individual interface programs are also given. Input control data files and example outputs for each interface
program are also provided for a test case.
vii
-------
CONTENTS
Preface ~ iii
Acknowledgments
Abstraa vii
Figures ~ ....... xiii
Tables xv
1. INTRODUCTION ..... .. 1
Z FEATURES AND LIMITATIONS OF THE INTERFACE 5
3. OVERVIEW OF THE INTERFACE SYSTEM 7
4. TECHNICAL APPROACHES .!.... 13
4.1 Attributes of the Regional and Urban Models Relevant to Interfacing 13
4.2 Treatment of Meteorological and Surface Parameters 18
4.2.1 Diffusion Break and Region Top Heights 18
4.2.2 Meteorological Scalars .. 21
4.2.3 Surface Air Temperature Field 23
4.2.4 Wind Fields 23
4.2.5 Surface Characteristics 28
4.3 Treatment of Concentrations 31
4.3.1 Initial Conditions - - 34
43.2 Lateral Boundary Conditions 35
4.3.3 Top Boundary Conditions 37
4.3.4 Summary of Concentration Interfacing 38
4.4 Treatment of Area Biogenic Emissions 38
ix
-------
CONTENTS (continued)
5. USER'S INSTRUCTIONS 40
5.1 Diffusion Break Data Processor (PDFSNBK) 41
5.1.1 Processor Function 41
5.1.2 Input/Output Components .. ™... 41
5.13 Resource Summary for a PDFSNBK Application 51
5.1.4 Example Run Stream Command File for an Interface Application 52
5.1.5 Main Program, Subroutines, Functions, and Block Data Required 52
5.1.6 I/O and Utility Library Subroutines and Functions Required 52
5.1.7 INCLUDE Files - ............... .. .. 52
5.2 Region Top Interface (IREGNTP) 53
5.2.1 Processor Function — 53
5.2.2 Input/Output Components 53
5.2.3 Resource Summary for an IREGNTP Application 65
5.2.4 Run Stream Command File for an Interface Application 66
5.2.5 Main Program, Subroutines, Functions, and Block Data Required 67
5.2.6 I/O and Utility Library Subroutines and Functions Required — 67
5.2.7 'INCLUDE Files 67
53 Temperature Interface (ITMPRTR) ........ 68
5.3.1 Processor Function - « 68
5.3.2 Input/Output Components 68
5.3.3 Resource Summary for an ITMPRTR Application 78
53.4 Run Stream Command File for an Interface Application 79
5.3.5 Main Program, Subroutines, Functions, and Block Data Required 80
5.3.6 I/O and Utility Library Subroutines and Functions Required 80
5.3.7 INCLUDE Hies ..... 80
5.4 Metscalars Interface (IMETSCL) 81
5.4.1 Processor Function ...» 81
5.4.2 Input/Output Components - «... 81
5.43 Resource Summary for an IMETSCL Application 95
5.4.4 Run Stream Command File for an Interface Application 96
5.4.5 Main Program, Subroutines, Functions, and Block Data Required 97
5.4.6 I/O and Utility Library Subroutines and Functions Required 97
5.4.7 INCLUDE Files 97
x
-------
CONTENTS (continued)
5.5 Wind Field Interface (IWIND) 98
5.5.1 Processor Function 98
5.5.2 Input/Output Components — 98
5.53 Resource Summary for an IWIND Application 112
5.5.4 Run Stream Command Hie for an IWIND Application .... 113
5.5.5 Main Program, Subroutines, Functions, and Block Data Required — 114
5.5.6 I/O and Utility Library Subroutines and Functions Required 114
5.5.7 INCLUOH Files 114
5.6 Surface Characteristics Interface (ICRETER) 115
5.6.1 Processor Function 115
5.6.2 Input/Output Components .. 115
5.6.3 Resource Summary for an ICRETER Application — 126
5.6.4 Run Stream Command File for an Interface Application 127
5.6.5 Main Program, Subroutines, Functions, and Block Data Required - 128
5.6.6 I/O and Utility Library Subroutines and Functions Required 128
5.6.7 INCLUDE Files - - - 128
5.7 Concentration Interface (ICONC) 129
5.7.1 Processor Function 129
5.7.2 Input/Output Components 129
5.7.3 Resource Summary for a UAM Application 151
5.7.4 Run Stream Command File for an Interface Application 152
5.7.5 Main Program, Subroutines, Functions, and Block Data Required 153
5.7.6 I/O and Utility Library Subroutines and Functions Required ...» 153
5.7.7 INCLUDE Files 153
5.8 Biogenic Emissions Interface (IBIOG) 154
5.8.1 Processor Function - 154
5.8.2 Input/Output Components .. 154
5.83 Resource Summary for a UAM Application 164
5.8.4 Run Stream Command File for an IBIOG Application 165
5.8.5 Main Program, Subroutines, Functions, and Block Data Required 165
5.8.6 I/O and Utility Library Subroutines and Functions Required - 165
5.8.7 INCLUDE Files ..... 165
xi
-------
CONTENTS (concluded)
REFERENCES 166
Appendix: A. Description of the Example Test Case * A-l
B. Input Control Hies for the Test Case B-l
C. Sample of ROM Input Data Files for the Test Case C-l
D. Sample Output of the Interfaces D-l
E. Conversion of Output Files (BINASC and ASCBIN) E-l
F. Magnetic Tape Listing of Programs and Data Files for the Example Test Case F-l
xii
-------
1
2
3
4
5
6
n
8
9
10
11
12
13
14
15
16
17
FIGURES
Page
Northeast regional modeling domain 3
Flow diagram showing the data retrieval and interface processing steps performed to
generate data files for the UAM preprocessors and model - 8
The gridded ROMNET region (64 columns by 52 rows); dots are situated at the ROM
grid cell corners
Wind field derived for an example UAM grid from ROM gridded wind components
1 A
Example of grid points (mid-points) of ROM cells overlaying a UAM domain. Two
ROM rows/columns extend beyond each UAM boundary 16
ROM vertical layer structure during daytime conditions 17
Time variation of the region top height (Z-f) and diffusion break height (Zog) over two
diurnal periods 20
T
Example set of ROM and UAM grid cells for the fractional area weighting method.
(ROM cells are about a factor of 4 larger than a UAM grid cell in this case.) 29
Boundary grid cells in the UAM model are the outer ceils enclosed by bold lines. ROM
grid points are shown in the lower left 36
Row diagram of the diffusion break data processor. PDFSNBK 43
Flow diagram of the IREGNTP interface program with input and output files 54
Flow diagram of the ITMPRTR interface program with input and output files 69
Flow diagram of the IMETSCL interface program with input and output files 82
Row diagram of the IWIND interface program with input and output files 99
Row diagram of the ICRETER interface program with input and output files 116
Row diagram of the ICONC interface program with input and output files 130
Row diagram of the IBIOG interface program with input and output files 155
xiii
-------
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
TABLES
Page
Summary of the UAM preprocessor programs 10
Retrieved ROM files used by interface programs ~ ~ 11
Overview of the ROM-UAM interface programs and input/output files 12
List of meteorological scalars 21
Methodology for derivation of the exposure index 23
ROM-UAM wind interfacing methodology 25
Land use categories and associated deposition factors : 30
Chemical species in the UAM (CB-FV) model 31
Vertical methodology for interfacing concentrations 33
Concentration interfacing procedures - 38
Input data files used by each interface processor - 41
DBDATA input file parameters 42
Control card variables for PDFSNBK 45
DBPACK parameters 49
PDFSNBK I/O file space requirements 51
MF166 parameter list for IREGNTP 55
PF119 parameter list for IREGNTP 56
DBDATA parameter list for IREGNTP 57
Control card variables for IREGNTP .. 59
RTPACK parameters 62
RTDATA parameters 64
IREGNTP I/O file space requirements 65
PF103 parameter list for ITMPRTR 70
Control card variables for ITMPRTR 72
TPPACK parameters 76
ITMPRTR I/O file space requirements 78
PF102 parameters for IMETSCL 84
PF117 parameters for IMETSCL - 85
PF119 parameters for IMETSCL .... 86
MF174 parameters for IMETSCL 87
DBDATA parameters for IMETSCL ... 88
RTDATA parameters for IMETSCL 88
xv
-------
33
34
35
36
37
38
39
40
41
42
4?
44
45
46
47
4S
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
90
93
95
100
101
102
102
103
104
106
108
110
112
117
118
120
124
126
131
132
133
134
136
138
140
142
145
147
149
151
156
157
159
160
162
164
TABLES (CONCLUDED)
Control card variables for IMETSCL .
MSPACK parameters
IMETSCL I/O file space requirements
MF165 parameters for IWIND
PF119 parameters for IWIND
DBDATA parameters for IWIND
RTDATA parameters for IWIND
PF115 parameters for IWIND
PF114 parameters for IWIND
WDDATA parameters (optional file)
Control card variables for rWTND
WDBIN parameters
IWIND I/O file space requirements
PF108 parameters for ICRETER
PF118 parameters for ICRETER
Control card variables for ICRETER
CRPACK parameters .... «...
ICRETER I/O file space requirements
MF165 parameters for ICONC
PF119 parameters for ICONC -
DBDATA parameters for ICONC
ROM21 parameters for ICONC
AQDATA parameters for ICONC
TCDATA parameters for ICONC
BCDATA parameters for ICONC
Control card variables for ICONC -
AQBIN parameters
TCBIN parameters ..
BCBIN parameters
ICONC I/O file space requirements
PF144 parameters for IBIOG
Anthropogenic emission parameters for IBIOG
Control card variables for IBIOG
EMBIN parameters
BIOASC parameters...»
IBIOG I/O file space requirements
XV]
-------
SECTION 1
INTRODUCTION
A set of computer programs has been developed to enable the Urban Airshed Model (UAM) system
to be driven by Regional Oxidant Model (ROM) system outputs. The computer program package has been
designated the ROM-UAM interface system. Since both the ROM and UAM are independent models with
their own separate processor systems, the interface programs serve as external links between particular ROM
output files and components of the UAM system. The principal functions of the interfaces are to interpolate
specific gridded parameters from the ROM outputs, and to generate data files in compatible formats for
input to the UAM preprocessors or the UAM program. In particular, the interface programs have been
designed for use with the ROM 2.1 system files and the UAM Carbon Bond IV [UAM (CB-IV)] system.
The impetus for the development of a ROM-UAM interface program was initiated under EPA's
Regional Ozone Modeling for Northeast Transport (ROMNET) study. One of the major objectives of
ROMNET was to apply the ROM to generate estimates of air quality concentrations for projection and
post-control scenarios in order to specify initial conditions and boundary values in urban-scale modeling to
be conducted by states in the development of State Implementation Plans. For this, the ROM 2.1 system was
applied to simulate selected episodes with emission scenarios that include a base year, a future year, and
control strategies. The procedures for the selection and number of high ozone episodes being modeled over
the northeastern United States have been described by Doll et al. f 1989).
The importance of regional and interurban transport of pollutant species in distinct urban plumes
within the ROMNET study region has been recognized from both experimental (Wolf and Lioy, 1980,
Clarke and Ching, 1983) and modeling studies (Rao, 1987). A major factor for such pollutant plumes is the
existence of emissions-rich areas associated with the corridor of large metropolitan centers along the north-
eastern coast of the United States and also the numerous large isolated urban areas within this region.
However, the specification of spatially-detailed inflow boundary conditions for urban-scale modeling has
been difficult to achieve due to a lack of spatial resolution in measurements both at the surface and aloft.
Additionally, appropriate future-year boundary conditions, which may be even more critical to urban-scale
modeling under certain emissions reduction strategies, would be difficult to quantify. Therefore, the results
1
-------
of the ROM simulations for the northeastern region (Figure 1) are to provide representative highly-resolved
boundary conditions of ozone and precursor species for use in urban model applications within this area. In
this regard, the ROM-UAM interface system is an integral component of ROMNET because it is the key
element in the transfer of the regional model results to urban model applications. Although it is to be
initially implemented in the UAM applications for urban domains within the ROMNET region, the interface
programs are expected to be suitable for the UAM applications using the ROM results from other modeled
regions.
The UAM has been identified as the preferred modeling approach for urban ozone modeling (EPA,
1986). Like the ROM, the UAM is a grid-based photochemical oxidant model that mathematically treats the
relevant physical and chemical processes important to ozone production, destruction* and removal, albeit on
a smaller spatial scale. Both models require extensive input data including concentrations for initial and
boundary conditions, and meteorological and emissions data. Consequently, the wide variety of data types
already assembled and available in the ROM system database should be exploited and the scope of the
interface effort was expanded to include more than just the ROM predicted concentrations. In all, 13
different ROM data files are applied in the interface programs, which include gridded fields of various mete-
orological parameters (e.g., winds, temperature, water vapor, etc.), land use information, and biogenic
emissions.
The ROM-UAM interface system actually consists of seven main programs. Each interface program
provides a link between certain ROM data files and a particular UAM preprocessor or the model program.
Specifically, an output file generated by an interface is in a compatible format for direct input to either a
UAM preprocessor or the modeL The ROM data files needed by the interface programs are created by the
user when exercising a data retrieval program specifically designed for this interface package. The data
retrieval program is part of the Gridded Model Information Support System (GMISS). The ROM database
and the GMISS data retrieval program are maintained on the EPA National Computer Center (NCC) IBM
3090 system. The retrieved data files generated by the user on this computer can be transferred to another
computer system if desired, where the interface programs and the UAM system will eventually be exercised.
The interface codes also exist on the IBM 3090, however, they are readily adaptable to other computer
systems since these algorithms,have been written according to American National Standard Institute (ANSI)
FORTRAN-77 language specifications.
Due te the extent of the coupling between these model systems through this interface package, the
performance of the UAM will be strongly impacted by results of the ROM simulations. The development of
the ROM system has progressed in parallel with a strong model evaluation program. The most recent
revision of the ROM (ROM 2.1), as documented in Young et al. (1989), incorporates changes based on com-
parison of the ROM 2.0 results against observed databases (Schere and Wayland, 1989). The ROM 2.1
2
-------
3
-------
model results have also been compared to ozone measurements in the ROMNET region by Pierce et aL
(1990a). The ROM evaluations have revealed the capability of the model to simulate the spatial pattern of
urban plumes and regional concentration levels. Changes in the UAM system have also occurred since the
version documented by Ames et aL (1985). Of particular significance, the photochemical mechanisms
imbedded in earlier versions of both models have been upgraded to the latest carbon bond (CB-IV)
mechanism. The details about the development of the CB-IV mechanism for these models have been
presented by Geiy et al. (1988, 1989) and appear in Volume I of the UAM(CB-IV) User's Guide (Morris
et al. 1990a). The interface programs have been specifically constructed for use with the current version of
the Urban Airshed Model, designated as UAM (CB-IV). The methodologies and user instructions for the
UAM (CB-IV) model are fully described in Volumes I and II of this User's Guide (Morris et al. 1990a,b).
For the purpose of brevity in this report, the acronyms ROM and UAM are meant to refer to the most recent
operational versions of these models.
This report is divided into sections devoted to various aspects of the interface package. Section 2 gives
a summary of some of the salient features and important limitations when applying the interface programs.
A broad overview of the interface package and its programs is provided in Section 3. A technical description
of the methods and procedures employed to derive the time dependent and spatially-varying gridded data
fields for the UAM system from the ROM outputs is given in Section 4. The detailed format specifications
of the input/output files for each interface and instructions for their execution, of particularly relevance to
computer specialists, are presented in Section 5. In addition, an example test case has also been assembled
and should be exercised after implementation of the interface codes on a particular computer. Samples of
the input/output files and run streams for an IBM system are also provided in Appendices to assist the user
in properly exercising the interface programs.
4
-------
SECTION 2
FEATURES AND LIMITATIONS OF THE INTERFACE
The ROM-UAM interface programs have been designed to utilize various ROM output files in gen-
erating input data files for the UAM preprocessor programs or the UAM model program. This coupling of
the two model systems is performed by seven individual interface programs. Notable features of the interface
package are enumerated below:
Interfaces link only components of the UAM (CB-IV) model system and certain outputs of the
ROM 2.1 system, and not earlier versions of these two models.
Four interfaces provide formatted input files for the UAM preprocessors and three interfaces
generate binary files directly for use in executing the UAM model.
Thirteen different ROM output files are applied in executing the various interface programs.
Users must create these specially ^windowed" data files for their particular UAM domain and sim-
ulation period by accessing the ROM database through a data retrieval program developed for the
Gridded Model Information Support System (GMISS) on the EPA NCC-IBM 3090 computer
system.
Interface programs have been generalized to allow the user to input specifications about the par-
ticular UAM domain and grid (i.e. origin, number oi grid ceils, and grid cell size) and the vertical
configuration of the UAM (i.e. number of lower and upper levels).
• Interfacing for initial, lateral and top boundary conditions is performed for 17 of the 23 pollutant
species which must be specified in the UAM (CB-IV) model Default values are defined for the
remaining six species.
Lateral and top boundary concentrations are resolved hourly and spatially at each UAM grid cell.
Options have been built into two interface programs to allow the user to apply non-ROM input
data files (e.g. concentrations, winds).
Quality-assured data may be incorporated into formatted packet files through an on-line editor at
a terminal.
5
-------
A utility program included with the interface package converts any binary file into an equivalent
ASCII formatted data file so that its contents may be examined prior to use in model execution.
• Horizontal wind fields generated by the wind interface are determined by matching the ROM
gridded wind components from layers 1 and 2 into the vertical levels of the UAM. The methodol-
ogies employed in the wind interface are those applied in the UAM Diagnostic Wind Model
system described in Volume III of this User's Guide (Douglas et aL 1990).
Certain limitations also exist with this version of the interface:
No interfacing of the diffusion break height is performed. The user must apply a method recom-
mended for deriving hourly values of this parameter as described in Volume II of this User's
Guide (Morris et al. 1990b). As pan of the interface package, a processor program (PDFSNBK)
produces a formatted packet file for the UAM diffusion break preprocessor program (DFSNBK)
using a user-supplied data file. A spatially-invariant diffusion break at each hour is assumed by
the interface programs.
No interfacing of anthropogenic area or point source emissions is performed. The UAM point
source preprocessor (PTSRCE) arid the Emissions Processor System (EPS) already exist for
creating emissions for these source types (Volume IV, Causiey et aL 1990). A biogenics interface
is limited to combining an existing binary anthropogenic area emissions inventory file and the
ROM gridded biogenic emissions of certain hydrocarbon and NOx species.
No interfacing is performed for the UAM preprocessors, CPREP or SPREP, the chemical param-
eters and the simulation control programs, respectively.
6
-------
SECTION 3
OVERVIEW OF THE INTERFACE SYSTEM
A general framework showing how the interface program package fits into the overall UAM model
system is depicted in Figure 2. The interface programs are executed before exercising any UAM preproces-
sor program. The important first step to be performed in the interfacing process is the creation of a set of
formatted data files which are needed as inputs to the interface programs.
The user must execute a menu-driven data retrieval program that has been specifically developed
when interfacing the UAM with the ROM system outputs. The ROM database is composed of large files of
predicted concentrations (CONC), processor output files (PF), and other ROM model files (MF). The
ROM database files and the interface data retrieval program exist exclusively on the EPA-NCC IBM 3090
computer system. In the retrieval step, the user supplies time and UAM domain-specific information to the
retrieval program which will automatically extract all the ROM parameters over the desired time period and
for the selected spatial "window*. Since the retrieved data files will be formatted, the user may decide io
transfer them to another computer system where the actual execution of the interface programs, and the
UAM preprocessor and model simulations will be conducted. Detailed instructions and procedures for
applying the retrieval program with the ROM database are described in CSC (1990).
The set of interface programs produce either formatted "packet" files which are in a compatible forma:
for direct input to particular UAM preprocessors or binary files that are ready for direct input to the UAM
model (Figure 2). For the latter, certain UAM preprocessors are bypassed with no execution required. The
interface codes have been programmed according to ANSI FORTRAN-77 full language specifications.
Therefore, although the interface programs reside on the IBM 3090 computer, these codes are easily
adaptable to other computer systems after minor revisions. The user has the ability to examine the contents
of "packet" files generated by most interface programs with an on-line editor. Furthermore, at this stage the
user may wish to supplement these files with additional quality-assured observed data before proceeding with
the execution of a UAM preprocessor. For the interface programs which generate binary files ready for the
model, an optional feature has been built into the codes that allow the user the flexibility of supplying an
alternate data set. The details in Section 5 should be consulted for the structure and format of non-ROM
/
-------
User Inputs
ROM CONC/PF/MF DATABASE
-RETRIEVED" DATA FILES
(Formatted/Transferable)
GMISS DATA RETRIEVAL PROGRAM
(EPA-NCC IBM 3090) .
-
jINTERFACE PROGRAMS;
FORMATTED "PACKET" FILES
, i ,
J UAM PREPROCESSOR PROGRAMS
I
y
^BINARY FILES!
! ,
UAM SIMULATION PROGRAM !
Figure Z Flow diagram showing the data retrieval and interface processing steps performed
to generate data files for the UAM preprocessors and modeL
8
-------
data files and the particular interface programs which have this capability. In this instance, it must be noted
that the user would be undertaking an additional effort to generate an alternate data file for an interface
program. A utility program is also part of the interface package for convening a binary file into a formatted
ASCII file version so that it can be examined before model execution.
The UAM model code requires 13 binary input files in any modeling application. Table 1 gives a
complete list of the UAM preprocessor names and the binary output file generated by each program. The
UAM preprocessor programs have been grouped into categories in Table 1 to indicate binary data files are
required for initial and boundary (lateral and top) concentrations, various meteorological parameters (e.g.,
winds, etc) and surface features, emissions, and finally, control information about the chemical species and
reactions, and the model run parameters.
A comprehensive approach was taken in interfacing ihe two model systems in order to take advantage
of the extensive variety of data sets available from the ROM model and processor network for use as inputs
to as many UAM preprocessors as possible. A complete list of data files retrieved from the ROM database
and their contents is provided in Table 2. The file names listed in Table 2 are those designated within the
ROM system and these file names will be imbedded within the extended file name of the retrieved files
generated by the user. One or more of these files are needed to exercise each interface program. With the
exception of the temperature profile data file, all of the. retrieved files contain gridded fields of the ROM
processor or model results. The structure and format of these files are documented in Section 5.
An overview of the complete ROM-UAM interface program package is provided in Table 3. It shows
the retrieved ROM files needed for each interface program. There are seven actual interface programs. The
designated name for an interface program begins with the letter "P and the remainder of its name is given by
the UAM preprocessor program name. In addition, a processor program (PDFSNBK) that reformats hourly
diffusion break data (DBDATA), supplied by the user, is part of the interface program package. The outputs
from the other interface programs are formatted "packet" files (e.g, RTPACK, etc.) for use as inputs to the
appropriate UAM preprocessor. Since these files are formatted, the user has the capability to examine their
contents and also can insert additional data before exercising the UAM preprocessors. The interface
programs for generating winds (IWIND) and for initial and boundary concentrations (ICONC) produce
binary files for use in UAM execution. Consequently, the diagnostic wind model (DWM) system and the
concentration preprocessors (i.e., AIRQUL, BNDARY, 17CONC) are not exercised when these interface
programs are applied. A combined area emissions file is created by the IBIOG interface program from a
user-supplied anthropogenic file and a ROM biogenic emissions file. The output emissions file of IBIOG is
also binary and ready for input in UAM simulations. However, interfacing of anthropogenic point or area
emissions is not performed.
9
-------
I TABLE 1. SUMMARY OF THE UAM PREPROCESSOR PROGRAMS
Name
Internal name in
binary file
Contents
CONCENTRATIONS
AIRQUL
AIRQUALITY
Initial concentration fields
BNDARY
BOUNDARY
Lateral boundary concentrations
TPCONC
TOPCONC
Top boundary concentrations
METEOROLOGY AND SURFACE CHARACTERISTICS
DFSNBK
DIFFBREAK
Diffusion break heights
REGNTP
REGIONTOP
Model region top heights
METSCL
METSCALARS
Five meteorological parameters and
photolysis rate
TMPRTR
TEMPERATURE
Surface temperature field
"DWM1
WIND
Horizontal wind component fields
CRETER
TERRAIN
Surface roughness and vegetation
fraction factor
EMISSIONS
EPS2
EMISSIONS
Anthropogenic area emissions
PTSRCE
PTSOURCE •
Major point source emissions
CONTROL
DATA
CPREP
CHEMPARAM
Species reaction rate information
SPREP
SIMCONTROL
Model simulation input information
1. DWM - Diagnostic Wind Model
2. EPS " Emissions Preprocessor System
10
-------
TABLE 2. RETRIEVED ROM FILES USED BY INTERFACE PROGRAMS
ROM system
file name
Description of file contents
ROM21
Hourly predicted concentrations of 17 species from layer 1, 2, and 3
MF165
Hourly gridded heights of layer 1
MF166
Hourly gridded heights of layer 2
MF174
Hourly gridded water vapor concentration in layer 1
PF102
Hourly vertical-interpolated temperature profiles
PF103
Hourly gridded surface air temperature
PF108
Gridded surface roughness length
PF114
Hourly gridded layer 2 horizontal wind components
PF115
Hourly gridded layer 1 horizontal wind components
PF117
Hourly gridded total sky cover fraction
PF118
Gridded fractions of eleven land use categories
PF119
Gridded terrain elevation
PF144
Hourly gridded biogenic emissions of six species
11
-------
TABLE 3. OVERVIEW OF THE ROM-UAM INTERFACE
PROGRAMS AND INPUT/OUTPUT FILES
RON or other Interface
Input files program
•PACKET' UAM
file preprocessor
Binary
file
User DBOATA
POFSNBK
-> DBPACIC
DFSNBK
0B8IN
OBOATA
MF166
PF119
IREGNTP
¦> RTOATA
-> RTPACK
REGNTP
RTBIN
NF174
PF102
PF117
PF119
RTDATA
DBOATA
IMETSCl
>> HSPAOC
user data
(optional)
HETSCL
HSBIN
PF103
ITMPRTR
-> TPPACK
user data
(optional)
TMPRTR
TPBIH
MF165
PF114
PF115
PF119
DBDATA
RTDATA
IWIND
-> UDBIN
PF108
PF118
ICRETER
•> CRPACK
User data
(optional)
CRETER
CRBIN
ROM21
MF165
PF119
DBDATA
ICONC
| AOS IN
•>| BCBIN
j TCBIH
PF144
User EMISSIONS*
IBIOG
I
BIOASC (optional)
EMBIN
Anthropogenic area emissions file
12
-------
SECTION 4
TECHNICAL APPROACHES
4.1 ATTRIBUTES OF THE REGIONAL AND URBAN MODELS RELEVANT TO INTERFACING
Although both models are Eulerian grid models, differences in the framework of their 3-dimensionaI
grids have to be addressed prior to the development of approaches for interfacing the ROM system data with
the UAM components. Fortunately, interfacing in the time dimension was straightforward since both model
systems have a 1-h time resolution in common. Output results from the ROM and processors contained in
the retrieved data files are available at hourly intervals, which is the time interval required of input data for
the (JAM preprocessor and model programs. The lime period of the ROM results in the retrieved data files
will also span a full 24-h period beginning at midnight of the day being simulated by the UAM, or two con-
secutive 24-h periods if a 2-day UAM simulation is planned.
Notable differences in the horizontal and vertical grid dimensions exist between these models that had
to be reconciled to properly interface the ROM results with the UAM system components. In the horizontal
dimension, the UAM is applied with a finer-mesh grid spacing and over a substantially smaller domain than
that of the ROM. The horizontal framework of the ROM grid is based on the latitude-longitude system.
Columns are north-south along longitude lines and rows are oriented in an east-west direction along latitude
lines with the horizontal resolution set at 1/4 degree of longitude by 1/6 degree of latitude. These specifica-
tions translate into a horizontal grid spacing of about 19 km in mid-latitudes. The ROM grid for the
ROMNET region as depicted in Figure 3 consists of 64 columns by 52 rows. Each grid point shown in
Figure 3 is situated at the lower left corner of a grid cell. In contrast, the horizontal grid framework of the
UAM is based on the Cartesian coordinate system. Horizontal grid spacing is specified by the user, and has
generally been defined to be from 2 km to 3 km. For a typical urban application, the spatial domain of the
UAM is generally on the order of 200 km, whereas the much larger ROMNET domain spans about 1,000 km
on each side. Clearly, a relatively small subregion of the ROM domain would provide sufficient overlap of a
particular UAM domain.
-------
m'ol +:
tfiO'tl M
• XQ'U A
gQO'te *¦
tfiO'tl JL
fiQ'fl *
JW'Sl *
gOQ'U M
¦¦ gpC'U JL
gOVU M
gQO'sU JL
gOO i»*» JL
gOO'Sl A
tflQ'Zl JL
fflO'U JL
gOO'SZ JL
gOO'OB JL
; ; tfiO'OS JL
• tfiO^S JL
Q00'Z8 M
¦ ffOO'ZB JL
sOOSS JL ¦
\ oOO^Q JL
ffOO'tS JL
o°0'9B JL ¦
oQQ'tB jl
q00'98 JL
cm
vn
*
«5
e
-fc
Z
2
C
££
rs
a&
2
o
££
43
<3
V
•S **
2 v
« «
V "5
C c
> *5
fi
«
hm
2
*a
14
-------
The appropriate spatial coverage of the ROM gridded parameters in a particular UAM application is
provided from information supplied by the user during the data retrieval step. Figure 4 demonstrates how
the ROM grid points overlap the UAM domain for applications with the interface programs. The data files
generated for the interface programs by the GMISS retrieval system contain parameters from all ROM grid
cells whose midpoints lie inside the UAM domain and in a surrounding buffer zone consisting of two grid
points outside each UAM boundary. The various parameters and concentrations at the ROM grid points
surrounding the UAM domain are particularly useful in the specification of boundary values for the urban
model.
Differences in the vertical structure and number of vertical levels between the models also presented a
challenge for interfacing several parameters. The height of the diffusion break, widely known as the mixing
height, is the key reference height in the UAM system that separates the lower and upper level cells. The
user specifies the number of prognostic UAM vertical levels below the diffusion break and the number of
upper levels situated between the diffusion break and the model's region top height. Lower levels are of
equal thickness, and each level expands or contracts according to the temporal behavior of the diffusion
break- The thickness of the UAM upper level(s) is controlled by the difference between the region top and
the diffusion break heights.
In contrast to the UAM vertical framework, the ROM exhibits three prognostic layers which vary
spatially and temporally. Figure 5 illustrates the configuration and vertical extent of each ROM layer during
the daytime. A shallow diagnostic surface layer (layer 0) is also depicted, however, it has not been implem-
ented in the ROM version 2.1. Consequently, predicted values from layer 0 have not been considered in any
interface approaches. Results from layer 1 have been applied as surrogates for surface values (Schere and
Wavland, 1989). During most of the daytime period, layers 1 and 2 are imbedded in the well-mixed convec-
tive boundary layer. The height of layer 1 varies from 10 to 600 m above ground and layer 2 heights are
generally from 1,000 to 1,500 m. The top model layer (layer 3) represents the synoptic scale inversion layer
capping the boundary layer or the vertical extent or a cumulus cloud layer, if present (Figure 5). At nigh:, the
contents of layer 2 represent the remnants of the bulk of the previous day's mixed layer, while the vertical
extent of the nocturnal radiative inversion is confined in layer 1. Finally, none of the ROM layers is meant to
simulate the temporal variation of the diffusion break height over the diurnal cycle.
The technical aspects and procedures incorporated into individual interface programs are described in
the remainder of jhis section.
15
-------
Fil.it t. E»i»l* '"S'w(«"•!«*"¦>n,ROM ' UAM
•|\«, ROM rows/columns extend beyond each UAM boundary.
-------
Daytime
Layer Functions
Layer 3
Layer 2
Layer 1
Inversion
or cloud
layer
Mixed
layer
Layer 0 —
\ Marina
_ _ h:Y"
Surface
Uytr
/
1.
2.
3.
\
(y.
2.
2.
\3
Downward traniport ol
stratospheric ozone
Upward traniport by
cumulus doudi
Liquid and gat phm
photochemistry
Long rang# traniport
by Iraa atmosphere
Gas phase photochemistry
Turbulence and wind shear
affects on transport and
diffusion
Deposition on mountains
Lake and marine layers
Effect on reaction vale of
uibgrkl scale segregation of
fresh and aged pollutants
Ground deposition
Spatial variation In mean
concentrations due to fine
point and area sources
Pjpint*- 5. ROM vciiicjil layer structure (luting daytime conditions.
-------
42 TREATMENT OF METEOROLOGICAL AND SURFACE PARAMETERS
4.2.1 Diffusfcon Break and Region Top Heights
The variation of vertical levels in the UAM is dictated by the diffusion break and region top heights.
In particular, the diffusion break (Zqb) is the key reference height which separates the sets of lower and
upper levels in the model and it serves as the boundary between the differing stability regimes that character-
ize these two vertical groups of levels. Hourly values of the diffusion break are needed by several interface
programs as well as a UAM preprocessor. However, as noted earlier, none of the ROM layer heights
emulates the diurnally-varying diffusion break height. Consequently, interfacing of Zdb values was not
feasible for this version of the interface programs. The user should refer to the description of DIFFBREAK
in Volume II of the User's Guide (Morris et al. 1990b) for a methodology to derive hourly diffusion break
values. Since other interface programs require a Zdb ^ development of this data file must be one of
the initial tasks to be performed.
A preprocessor program contained within the interface package has been developed to generate a
formatted "packet" file compatible with the input format specifications of the DFSNBK preprocessor. The
user should consult Section 5 for the format specifications of the input data file (DBDATA) for the interface
diffusion break processor program (PDFSNBK). A set of 24 hourly diffusion break heights is needed by
interface programs if the UAM will be simulating a single day. For a 2-day UAM application, hourly
diffusion break values for two full days must be prepared.
The region top height (Zj) defines the total extent of the UAM domain in the vertical dimension.
The magnitude and time variation of Zj is important in model applications since the thickness of upper
level(s) is determined from the difference between Zj and the diffusion break height for each hour. Addi-
tionally, Zt should be sufficiently high that elevated point source plumes remain within one of the model's
vertical levels. If a point source plume rises above Zf, its emissions are above the model domain and are not
considered in the mode! simulation.
The interface method to derive Zj uses the height of ROM layer 2. During the nocturnal period,
ROM layer 2 height was designed to represent the vertical extent of the previous day's mixed layer. In the
region top interface, a UAM domain-wide average height of layer 2 is determined from the ROM gridded
values for each hour of the simulation day. The lowest average layer 2 height value (ZR2min) is chosen as the
initial region top height at the beginning hour of the UAM simulation, which is expected be any hour prior
to sunrise on the day being modeled. The region top height is allowed to vary temporally, however; like Zdb>
no spatial variation has been imposed on Zj in this version of the interface. The hourly variation of Zj in
the interface is described in equation (1).
IS
-------
Zr(f) - z
t2mt»
*(z
DMmax ^ IZni* ~ Z Dtmt* ) t ^ DSmax Z Dtmtn )] CD
where: Z-rft) - region top height at hour t
ZdbC) = diffusion break height at hour t
^DBmm = morning minimum diffusion break height
^DBmax — afternoon maximum diffusion break height
A2 - minimum upper thickness interval = DZm • IZU
DZ„ = upper level minimum thickness criterion
IZU = number of upper levels
^R2min ^ 1,000 m (minimum criterion)
An example of the variation of Zj by applying equation (1) and the temporal variation of Zdb a**e
illustrated in Figure 6. The region top height increases gradually during the post-sunrise period and reaches
its highest value in the afternoon when Zdb reaches a maximum. The weighting function inside the brackets
in equation (1) is based on the temporal variation of Zqb anc* controls the behavior of Zj. The region top
height descends gradually during the evening and nocturnal hours, while Zdb generally decreases much more
dramatically in Figure 6. During the time period after the maximum Zqb has been reached, equation (1)
continues to specify the behavior of Zj. However, a new minimum ROM average layer 2 height, computed
from values from the next day, is substituted for ZR2min- Thus, in a single day simulation, the retrieved data
fiie of iave: 2 heights must contain the ROM results from two consecutive days, it is also evident from
Figure 6 that Zj remains above Zdb ^th this formulation because a minimum'thickness criterion (£ZU) has
also been implemented. At the time of ZoBmax* Zj is greater than Zdb by at l^e product of the upper
level minimum thickness (Z5ZU) and the number of upper levels (IZU). This requirement also ensures that
the upper levels remain above Zdb al hours of the model simulation. In the test case, the value for DZU
was set to 100 m.
The interface program for region top also requires a gridded terrain elevation file since ROM layer 2
height values have been written as altitudes above sea level. An additional input file for the region top
interface includes the hourly values of Zdb f°r a 24-h period for single day simulations. The output file
created by the region top interface (IREGNTP) is in a compatible format for direct input to the REGNTP
preprocessor program.
19
-------
2500
i ¦ i11 ri' i' i * i • ir r111 t 11"11 • iT111 'TTrn
2000
1500
1000
500
DB
) II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1 I I I I I I I I I I I I I I I
00 04 08 12 16 20 00 04 08 12 16 20 00
02 06 10 14 18 22 02 06 10 14 18 . 22
Figure 6. Time variation of the region lop height (ZT) and diffusion break height (ZDB) over two diurnal periods.
-------
4.2.2 Meteorological Sea tars
There are six parameters that must be specified on an hourly basis in an input file for the METSCL
preprocessor. Table 4 provides a list of the meteorological scalar parameters and a brief description of each
variable. No spatial variation has been built into METSCL for these parameters. Five of the six parameters
are specified or derived with retrieved data files from the ROM system. Atmospheric pressure (ATMOS-
PRESS) has not been interfaced. A default value of 1.0 atm (Le. 1 atm * 1013.25 mb for a standard atmo-
sphere) has been set for this parameter for each hour. However, a user may wish to substitute atmospheric
pressure measurements, if available, in place of this default value in the formatted packet file (MSPACK)
generated by the interface program.
The vertical temperature gradients represent layer-average values below and above the diffusion break
height. Upper-air rawinsonde temperature profiles obtained twice-daily at National Weather Service sites
have been interpolated at 50-m increments and to hourly intervals by a ROM processor (Young et al. 1989).*
During the data retrieval phase, the interpolated temperature profiles from upper-air site(s) located within
the UAM domain are provided and the user has the ability to request one or more additional sites. However,
before applying the retrieved temperature profiles in the interface, a user should examine the representative-
ness of the temperature profiles at a particular site for the meteorological conditions existing over the
. domain during the simulation period; particularly profiles from any selected site located outside the model
domain.
TABLE 4. LIST
OF METEOROLOGICAL SCALARS
Parameter
name
Internal
units
Definition
TGRADSELOW
.K/m
Vertical temperature gradient (dT/dZ)
from the surface to diffusion break height
TGRADA80VE
K/n
Vertical temperature gradient between
diffusion break and region top heights
EXPCLASS
Exposure class * integer scale indicator
of the near-surface atmospheric stability
RADFACTOR
min*1
N02 phetolysis rate constant, ki
CQNCUATER
ppa
average surface water vapor concentration
ATMOSPRESS
atm
surface atmospheric pressure
21
-------
The individual temperature gradients at 50-m intervals from the upper-air site(s) are used to compute
the hourly layer-average values below and above the diffusion break. During nocturnal hours when a
surface-based inversion layer often exists, positive values for TGRADBELOW can be expected. Although
notable spatial variations in the nocturnal low-level temperature structure have been found within large
urban areas (Godowitch et aL 19S7), values of these variables are assumed to be spatially invariant in the
current version of the METSCL. During the daytime period, values of TGRADBELOW should be close to
the adiabatic lapse rate (-0.010 K/m) or even slightly super-adiabatic, while daytime values of TGRADA-
BOVE are expected to reflect a slightly stable (d7VdZ > -0.01 K/m) layer or an inversion lapse rate (positive
d77dZ).
The exposure class (EXPCLASS) is a unitless index with values ranging between -2 and 3. It is
intended to be an indicator of the atmospheric stability near the surface due to either solar heating or radia-
tional cooling. The methodology applied to derive hourly EXPCLASS values is presented in Table 5.
EXPCLASS depends on the solar zenith angle and cloud cover. A retrieved ROM file of gridded fractional
cloud cover values interpolated from observations (Young et al. 1989) is used to compute an hourly domain-
wide average fractional cloud cover. The latitude-longitude of the middle of the UAM domain is adequate
for the solar zenith angle calculations. Table 5 reveals that positive EXPCLASS values occur during daytime
hours and negative values are restricted to the nocturnal period. The cloud cover criteria are applied to
account for the attenuation of solar insolation or reduced radiational cooling due to the presence of clouds
during the day and nocturnal hours, respectively. In the midday period when solar insolation is a maximum,
the highest possible EXPCLASS value of +3 is achieved if the cloud cover fraction is under 50%. This
methodology is identical to the scheme employed by Morris et al. (Volume II, 1990b).
The N02 photolysis rate constant (RADFACTOR) is an important parameter since it impacts the
photochemical reactions built into the carbon-bond chemical mechanism. A matrix of NO2 photolysis rate
constants (Demerjian et al., 1980) dependent on zenith angle and altitude has been incorporated into the
interface program (IMETSCL) to compute the RADFACTOR along with the date, lime, diffusion break
height and latitude-longitude position. In contrast to the other metscalar parameters, which are specified at
the beginning of each hour, the model defines RADFACTOR values at the end of each hour and performs a
linear interpolation from hourly values to individual time steps. Consequently, the RADFACTOR value is
computed with the solar zenith angle at the end of each hour. In addition, a RADFACTOR value is also
generated for the hour before the beginning hour of model simulation. During nocturnal hours,
RADFACTOR is near zero and night-time chemistry takes place. When the RADFACTOR exceeds a
threshold value of +0.011 min1, the model switches to the daytime photochemical mechanism. Clear-sky
values for RADFACTOR are currently computed by this version of the interface for use in the UAM
(CB-IV) modeL
-------
TA8LE 5. METHODOLOGY FOR DERIVATION OF
THE EXPOSURE INDEX
Solar zenith
Domain average
EXPCLASS
angle (degrees)
cloud cover (X)
(unitless)
CMOCTURNAL HOURS]
>85
150
•2
>85
>50
-1
DAYLIGHT HOURS]
<30
£50
3
<3a
>50
2
30 < 0 < 55
<50
2
30 < e < 55
>50
1
55 < 0 < 85
<50
1
55 < 9 < 85
>50
0
The concentration of water vapor (CONCWATER) in the lower atmosphere is also a metscalar
parameter needed by the model. A domain-wide average value is computed for each hour with the ROM
hourly gridded water vapor concentrationsfrom layer 1.
4*2 J Surface Air Temperature Field
Hourly surface temperatures are needed by the TMPRTR preprocessor program. A retrieved ROM
file containing hourly gridded surface air temperatures interpolated from National Weather Service sites
(Young et al. 1989) is utilized by the FTMPRTR program to generate a formatted file for input into the
TMPRTR preprocessor program. No spatial interpolation is performed in the interface program. Tne
function of the interface program is to reformat the gridded ROM values into a compatible format for
TMPRTR, which spatially interpolates temperatures to the UAM grids. The user may examine the
formatted packet file produced by the interface, and quality- assured hourly temperature data from non-
gridded sites if available, may be inserted into this file before processing it in the TMPRTR preprocessor.
4.2.4 Wind Fields
An accurate representation of the 3-dimensionaI wind flow over the domain is crucial to the model's
ability to simulate the magnitudes and spatial patterns of pollutant species. Wind fields from ROM layers 1
and 2 are used in the wind interface program (IWIND).
23
-------
In the ROM 2.1, the wind field for layer 1 is generated from observed surface data. Layer-average
wind components for layer 2 are derived from upper-air wind data, however, surface winds are also given
some weight in the determination of the gridded winds in this layer (Young et ai 1989). The wind field from
ROM layer 3 has been excluded from consideration in this interface since layer 3 generally represents the
flow above the UAM domain.
A practical methodology was developed to interface the ROM gridded wind fields into the multiple
levels of the urban model for any used-defined vertical configuration. In order to capture the important
diurnal variations that often occur in the wind structure, the gridded winds from ROM layers 1 and 2 are
applied in the wind interface. The approach designed to match the ROM layer winds into the UAM vertical
cells is outlined in Table 6. The gridded heights of ROM layer 1 are compared to the heights of the bottom
and top of each UAM level. If a UAM level is completely imbedded in ROM layer 1, then the gridded layer 1
wind components are specified for that level. For any UAM levels entirely above ROM layer 1. then ROM
layer 2 winds are applied to define transport. For the condition where a UAM level overlaps both ROM
layers, weighting factors based on the fractional amounts that the UAM level overlaps each layer are applied
to the wind components of each ROM layer to determine the wind components (Table 6).
The wind interface also applies certain methods found in the Diagnostic Wind Model (DWM) system,
one of the principal program components of the UAM system. After the wind components have been
matched into each UAM vertical level, an inverse distance-squared weighting technique (equation 2) is
applied to spatially interpolate the wind components, still at the resolution of the ROM grid points, to the
fine-mesh UAM grid points. The spatial interpolation procedure is applied to each wind component field at
each vertical level.
(u.i/),, - Xc^.O • r;2 / (2)
«• 1
where (u.v)tJ = wind components at UAM grid point i,;
= the ROM wind components at grid point n
N ¦ maximum of five surrounding ROM grid point values
r
= distance between UAM grid point and a ROM grid point
-------
TABLE 6. ROM-UAM WIND INTERFACING METHODOLOGY
Criteria:
* If Ztic i Zri ; use ROM layer 1 winds
* If ZbIc * ZR1 ; use ROM layer 2 winds
* If a UAH level overlaps both ROM layers (Zgic < Zri i Zjjc>;
determine weighting terms (U1, u2) from:
wi » (ZR1 - ZBk) / (Ztic - *bic>
U2 «
where Zri * ROM layer 1 height
Zjic ¦ top (T) of a UAM vertical level k
ZbIc s bottom CB) a a UAM vertical level k
s UR1/VR1 * m2/vR2 " W2
where ur#,vr# = ROM griddea wind comoonents in layer 4 ,2
Example configuration of models around sunrise:
sssiissauauauaat ^R2
4
3
2
1
SFC
UAH
°R2»VR2
Zqq niaasaaiaiiatsxataa 7^ ^
2t1 "R1,vri
Z3« —SFC
ROM
25
-------
Using equation (2), the wind components at UAM grid point i, j are determined from values at the
nearest surrounding ROM grid points. A default man'mnm radius of influence (RMAX) of 25 km has been
imposed for this purpose so that only the nearby ROM grid points are included in the interpolation
procedure. If RMAX was set too large, unwanted smoothing could occur in the interpolated field. Values
from up to five ROM grid points may be used in the interpolation expression in equation (2). Another con-
straint required when applying an inverse-distance weighting method is to supply a minimum distance
criterion since the distance between two grid points must always be nonzero (division by zero produces an
error on many computer systems). Therefore, a minimum distance (RMIN) ofl km is suggested between a
pair of ROM and UAM grid points.
~
The next step is to subject the interpolated wind components to a five-point filter technique, which
reduces any spatial discontinuities and dampens vertical velocities contained in the interpolated horizontal
wind field. The purpose of applying the filter is also to reduce anomalous divergence as much as possible.
The form of the five-point filter is given bv:
0.5-A'(i, ;) +0.125 ¦[*(!~ 1 . - AT(i - 1 + y - 1 )**(£.;«¦ 1)] (3)
where X is either the uorv wind component, and is the smoothed value.
Only values at the surrounding four UAM grid points are employed in this filter technique at a given grid
point (i,/). The number of times that the wind component field is subjected to the filter method is specified
by the value of NSMTH. In the test case, NSMTH was set equal to 2.
Next, an initial vertical velocity field is computed at each level from the divergence derived from the
smoothed horizontal wind component fields. Unrealistically large vertical velocities may still remain. Con-
sequently, a methuti appiied in the DWM has also been implemented in the wind interface thai progressively
diminishes vertical velocities toward zero at the region top (Douglas et al. 1990). However, the horizontal
wind component fields are not mass-cotisistent after vertical velocities have been revised in this manner.
Therefore, a final procedure is to exercise an iterative technique in order to minimize divergence which
involves slight adjustments of the horizontal wind components throughout the entire grid until a minimum
divergence criterion is reached (default minimum divergence = lx lO6 s-1). The final products of the wind
interface are gridded fields of u,v components at each UAM level. An example of the wind field at level 1 for
a UAM domain obtained from the ROM gridded winds according to these procedures is displayed in
Figure 7.
26
-------
INTERFACED WIND FIELD FOR UflM ON DGTE/HR: 80203/12 LEVEL 1
- - - S ^ i*. LLVLl- 1
0.82S£h)j
f^XlMUh v£CTOfi
Figure 7. Wind field derived for an example UAM grid from ROM gridded wind components.
2?
-------
An optional feature also exists in the interface to allow the user to input a wind field file, already
gridded for the UAM domain, which had been generated from another wind modeL The interface can accept
the wind file and create a binary wind file compatible with the UAM The user is referred to Section 5 for the
input format specifications of an alternate wind file for the interface program.
An alternative to wind interfacing is the DWM system, which is a stand-alone independent package
available to the user with the UAM (CB-IV) system [Douglas et al. 1990]. If the user elects to apply the
DWM, surface and upper-air wind data must be processed in order to exercise the computer programs asso-
ciated with this wind model The DWM system has a postprocessor program which generates a binary wind
file for the UAM.
4.2,5 Surface Characteristics
The two surface characteristics required by the CRETER preprocessor are gridded fields of surface
roughness length (ROUGHNESS) and the VEGFACTOR, a measure of the relative surface uptake capabil-
ity of a particular land use type compared to that of an alfalfa crop.
An interface program has been developed to directly apply the gridded ROM fields of surface
roughness length (ZQ) and a land use inventory available at the resolution of the ROM grid is employed to
derive grid-average vegetation factor values. Both of these surface parameters are employed to treat dry
deposition processes in the UAM and ZQ values are also applied in the derivation of vertical diffusivitv in the
model.
An area-weighting scheme was selected as a more appropriate approach than the distance-weighting
scheme for the determination of UAM gridded values for these surface parameters from gridded ROM
values. With the area-weighting technique, the fractional amounts of each UAM cell covered by different
ROM cells are determined. An algorithm based on slopes and intercepts between grid cell lines accurately
computes the fractional area of a UAM grid cell covered by any ROM cell. An example of the area-
weighting scheme is provided by a subset of the ROM and UAM grid cells in Figure 8. For UAM cell Ul,
contributions from all four ROM cells would be fractionally weighted to determine the grid-area average
value. On the other hand, grid-area average value for U4 would be specified totally by R4 since U4 is com-
pletely inside ROM cell R4. The number of ROM grid cells impacting a particular UAM cell is certainly
dependent on the horizontal grid cell size of the urban model. As grid cell size decreases, more UAM grid
cells may be completely imbedded in a single ROM grid cell since the ROM grids remains fixed at about 19
km on a side.
28
-------
%
R3
U3
U4
R4
U1
U2
R 1
R2
Figure 8. Example set of ROM and UAM grid cells for the fractional area weighting method.
(The ROM cells are about a factor of 4 larger than a UAM grid cell in this case.)
29
-------
The eleven land use categories contained in the ROM gridded inventory are presented in Table 7.
The fractional coverage of each land use in each ROM cell comprises the land use inventory data. For dry
deposition in the UAM, a deposition factor (f3) represents the relative surface uptake rate of a particular
land use category compared to an alfalfa crop. The following expression has been developed to derive an
average vegetation factor for each UAM grid cell
where: aa - T is the overall vegetative factor for ROM cell n from the fractional amounts of each
land use (Lm) times deposition factor (Jm for the land use types in Table 7.
In equation (4), the areal contribution of each ROM cell to a UAM cell's area (An JAU) is summed in order to
obtain the UAM grid-area average VEGFACTOR value.
Values of ROUGHNESS are also derived for each UAM grid by applying the area-weighting
technique to ROM gridded surface roughness values. The formatted output file of the terrain interface
(ICRETER) is in a compatible format for direct input to the CRETER preprocessor.
VECFACTOR,,-
(4)
TABLE 7. IAN0 USE CATEGORIES AND ASSOCIATED DEPOSITION FACTORS
Designation Description
Deposition factor (3 )
1 URBAN
2 AGRI
3 RANGE
4 DP
urban, little vegetation
agricultural land; aoecuate water
0.2
0.5
range Land, usually low soil moisture 0.4
7 WATER
8 BARREN
9 NFU
10 MIXED
11 ROCKY
5 EV
6 MP
deciduous forest
evergreen (coniferous) forest
nixed forest, including vet land
water bodies (fresh or salt water)
barren land, mostly desert
non-forested wetland
mixed agriculture and range land
rocky areas with low shrubs-lichens
0.4
0.3
0.3
0.03
0.2
0.3
0.5
0.3
30
-------
43 TREATMENT OF CONCENTRATIONS
The concentrations of 23 chemical species must be specified for initial conditions, lateral boundary
conditions, and top boundary conditions in model applications with the CB-TV version of the UAM. A
complete list of the chemical species and their alphanumeric names designated internally in the UAM system
are provided in Table 8. Complete details about the development of the carbon-bond chemical mechanism
have been documented by Gery et al. (1988, 1989) and its adaptation into the UAM is fully described in
Morris et al. (Volume 1,1990a).
TABLE 6. CHEMICAL SPECIES IN THE UAM (CB-IV) MOOEL
Model
nomenclature Chemical name
Interfaced
(X)
1 NO Nitric oxide
2 N02 Nitrogen dioxide
3 03 Ozone
4 OLE Olefinic carbon bond species
5 PAR Paraffinic carbon bond species
6 TOL Toluene
7 XTL Xylene
8 FORM Formaldehyde
9 ALD2 Higher molecular weight aldehydes
10 ETH Ethene
11 CRES Cresol and higher molecular
X
X
X
X
X
X
X
X
X
X
weight phenols
12 MGIY Methylglyoxal
13 OPEN Aromatic ring fragment acid
14 PNA Peroxynitric acid
15 NXOY Nitrogen species group
16 PAN Peroxyacetyl nitrate
17 CO Carbon monoxide
18 H0N0 Nitrous acid
19 H202 Hydrogen peroxide
20 HN03 Nitric acid
21 MEOK Methanol
22 ETOH Etnanol
23 I SOP Isoorene
X
X
X
X
X
X
X
NOTE: X ¦ concentrations of species interfaced from the ROM.
Default mininun value defined for six species not interfaced.
31 •
-------
Since both the regional and urban model systems include the CB-IV chemical mechanism, it is
possible to interface the ROM predicted concentrations for ail species (except for ETOH, which is not
available from the ROM simulations). However, sensitivity test simulations were undertaken to investigate
whether it is necessary to interface all chemical species. A series of model test simulations involved different
sets of species interfaced with ROM concentrations while the remaining species were prescribed by a
minimum default value (lO6 ppm). The same results for ozone were achieved by interfacing 17 species from
the ROM as with the full set, while results differed when interfacing fewer species from the ROM. Limiting
the number of species to be interfaced also helped to reduce the size of the concentration file generated by •
the retrieval program. Table 8 also indicates the 17 species interfaced from the ROM predicted concentra-
tions, which are employed in the derivation of initial, lateral boundary, and top boundary conditions for the
UAM. Although values for the other six species are not derived by the concentration interface program,
minimum default values are specified in the CHE MP ARAM file which is created by the CPREP preproces-
sor program for the model.
The most challenging aspect of interfacing concentrations was to develop a method to match the con-
centrations from the three layers of the ROM into a user-specified number of UAM vertical cells and to
design the scheme to be applicable over an entire diurnal period. The methodology outlined in Table 9 has
been incorporated into the concentration interface (1CONC) program. The approach presented in Table 9 is
versatile since it can accommodate any number of user-defined vertical levels and is applicable over the
entire diurnal cycle. It is also a realistic approach for matching of ROM layer concentrations into multiple
UAM vertical levels based on knowledge of the relationship between the diffusion break height and ROM
layer heights with time. This method is applied to obtain concentrations at the various vertical UAM levels
for initial and boundary conditions.
The key feature of the vertical interfacing methodology for lower levels is the weighting scheme which
is dependent upon the time variation of Zqb- Concentrations of species in lower levels are derived from
equation (5) in Table 9. It shows that ROM layer 1 and 2 concentrations are applied to specify lower leve;
values. There is a criterion that if ROM layer 1 height (Zrx) is greater than Zqb* ROM layer 1 concentra-
tions are exclusively employed to define UAM lower level concentrations. This condition often exists during
nocturnal or early morning hours. As the diffusion break approaches the maximum value (ZoBmax)* El
approaches zero while F2 goes to unity in equation (8) and (9), respectively. Consequently, concentrations
for lower UAM levels approach the average value of ROM layers 1 and 2. In addition, it is evident that the
same concentrations are specified for all lower levels (i.e. no vertical concentration gradient). The rationale
for the lack of a concentration gradient across the lower levels is that mixing is expected to be sufficiently
vigorous below Zqb at any h°ur that vertical gradients are quickly eliminated by the model.
32
-------
TABLE 9. VERTICAL METHODOLOGY FOR INTERFACING CONCENTRATIONS1
ROM (R) Layers (layer thickness not to scale):
ZR3
CR3 * layer 3 concentrations
ZR2
CR2 * Layer 2 concentrations
2R1
Cr] • layer 1 concentrations
_____ Q
Example Configuration of UAM Vertical Levels (thicknesses not lo scale):
Cx ¦ top concentration
_ Zj s region top height
CU<3)
Upper (U) levels between 2DB and Zj
Cjj(2) Nurber of upper levels = !ZU
Cud)
2DB
Cl(2) Lower (L) levels below the diffusion break (Zqq)
—— Number of lower levels » IZL
Cl(D
Lower level concentrations (C'Lj:
CL(k) - C„ -Fl + [(C„ + Ck)/2}'F2 for k - 1". IZU (5)
If Z„(()
-------
Vertical concentration differences have been included in the derivation of values for upper levels
according to equation (6) in Table 9. Upper level concentrations are controlled by ROM layer 2 and layer 3
values since the top concentration (Op) is dependent on ROM layer 3 values. At the beginning hour of sim-
ulation, C? equals the ROM layer 2 concentration (Cr2) a&d no vertical concentration gradient exists across
the upper levels. The rationale for this scheme follows from the specification of the height of ROM layer 2
as the initial value of Z?. In addition, layer 2 represents a rather thick residual layer of pollutants which have
been well-mixed during the previous daytime period. A strong vertical gradient may develop across the
upper levels in the UAM during the day because Op approaches ROM layer 3 concentrations. Layer 3 con-
cenlVations are generally near tropospheric background values which can be considerably lower than layer 2
concentrations in certain areas of the model domain.
The following sections describe the procedures in the concentration interface (ICONC) which derives
initial, lateral boundary, and top boundary conditions of the pollutant species for the UAM.
4J.1 Initial Conditions
The set of initial conditions represents the concentrations of all species in each cell of the model grid
at the starting hour of simulation. Model predicted concentrations are certainly impacted by initial condi-
tions for some time. However, the influence of initial conditions diminishes as a simulation progresses.
The procedure applied in the concentration interface (ICONC) for deriving initial concentrations at
each UAM grid cell begins with the use of the vertical method already described in Table 9. Once concen-
trations have been derived at each UAM vertical level, values must be spatially resolved to each UAM grid
cell by applying the inverse distance-squared interpolation technique described earlier.
where = concentration of species m at UAM grid ij
= concentration of species m at ROM grid cell n
r„ = distance between midpoints of a UAM and a ROM grid cell
In applying the spatial interpolation step prescribed by equation (10), concentrations in ROM grid
cells immediately surrounding each UAM grid cell are included in the interpolation procedure to preserve
horizontal gradients that may exist in the ROM gridded concentration field.
34
-------
Initial concentration fields are determined with the above procedure for the 17 species identified in
Table 8. The gridded arrays of initial concentrations of these pollutant species are written to a binary file for
direct use in the UAM. Thus, the UAM preprocessor for initial conditions (A1RQUL) will not be exercised
when applying the interface for concentrations. If a user wishes to examine the initial concentration file, a
binary to ASCII conversion program has been included in the interface package (Appendix E). This conver-
sion program generates a formatted ASCII data file from a binary data file so that the initial conditions may
be examined via any on-line editor or the ASCII formatted file can be listed on a line printer.
4.3.2 Lateral Boundary Conditions
One of the primary purposes for interfacing the ROM and UAM systems is to specify boundary con-
centrations for the urban model from regional model simulation results, which provides for considerably
greater spatial density of concentrations than Is available from existing ambient monitoring sites. The inflow
of pollutants along the upwind model boundary is of particular importance in UAM simulations since
pollutant boundary concentrations are advected into the heart of the domain where they become involved in
the variety of photochemical reactions with the urban emissions.
A methodology has been designed to provide for temporally-varying concentrations of each interfaced
species at each boundary cell. In the UAM, the grid cells around the outer edge of the domain at each level
constitute the group of lateral bpundary cells (Figure 9). The first step is to apply the vertical method
described earlier to obtain concentrations at each vertical UAM level for each ROM grid point location.
The next procedure consists of spatially averaging the values from three ROM grid points: the two ROM
grid points in each row (or column) exterior to the UAM domain and the ROM grid point in each row (or
column) just inside the boundary. This averaging step is illustrated with the set of ROM points shown in
Figure 9. The averaging procedure is performed with each set of these ROM grid points for each ROM row
and column around the entire perimeter of the UAM domain. These averaged values represent the
boundary values along each side at the resolution of the ROM grid. The averaging of values over the
outermost ROM grid points provides some spatial smoothing for boundary conditions. The last step is to
derive boundary concentrations at each UAM grid point. Linear interpolation is employed using the
averaged ROM values along each side of the UAM domain to derive boundary concentrations at each UAM
grid ceil. This step is repeated to determine boundary concentrations at each vertical level and the entire
procedure is also performed each hour. The lateral boundary concentrations are written to a binary file
(BCBIN) for use in UAM simulations.
35
-------
4
Figure 9. Boundary grid cells in the UAM model are the outer cells enclosed by bold lines.
ROM grid points are shown in the lower left
36
-------
4-3-3 Top Boundary Conditions
Boundary concentrations must also be defined at the top of the model domain. The procedures
installed in the ICONC interface allow top boundary values to vary both In time and space in order to take
full advantage of the ROM predicted concentration fields. Top concentrations can have a greater impact on
surface concentrations in a UAM configuration where the diffusion break height and region top height
become identical. However, this feature has been eliminated when interfacing is applied as noted in
Section 4.2.1. Nevertheless, top concentrations can still be gradually mixed into the lower levels even across
the rather shallow upper levels. Therefore, top concentrations must be properly specified.
The derivation of the top concentration (Or) begins with equation (7) in Table 9. It provides for the
hourly evolution of Cj values at each ROM grid point overlapping the UAM domain. As noted eariier, Cj
has been designed to evolve from ROM layer 2 to layer 3 concentrations during the course of the daytime
period. Then the same inverse distance-squared weighting technique described earlier is employed on CV
values at the ROM grid points to resolve concentrations at each UAM grid cell at :he top of the UAM.
These steps are repeated hourly to provide for temporally-varying top concentrations for the UAM. A
separate binary output file is generated that contains the top boundary concentrations (TCBIN).
4,3.4 Summary of Concentration Interfacing
The specification of initial, and lateral and top boundary conditions is a primary objective of the
ROM-UAM interface effort. The methodologies described in the previous sections have been designed to
provide the fullest possible temporal and spatial resolution of concentrations for these key conditions in the
UAM from ROM gridded concentrations. An overview of the steps undertaken to resolve initial, boundary,
and top conditions in the concentration interface is given in Table 10.
37
-------
TABLE 10. CONCENTRATION INTERFACING PROCEDURES
INITIAL CONDITIONS:
Perform vertical method described in Table 9 with ROM gridded concentrations at the
starting hour to derive values at each UAM level
Perform horizontal interpolation to obtain values at each UAM grid point using the
inverse distance-squared method.
LATERAL BOUNDARY CONDITIONS:
Perform vertical method with ROM gridded concentrations to derive concentrations at
each UAM level
Average the two exterior ROM grid points and the ROM grid point located immediately
inside the UAM boundary in each ROM row/column.
Perform linear interpolation to resolve UAM boundary grid values.
Iterate the above steps for each hour.
TOP CONDITIONS:
Use vertical method to determine top concentrations at each ROM "rid point.
Perform horizontal interpolation to spatially resolve concentrations to each UAM grid
Iterate to perform above steps for each hour.
4.4 TREATMENT OF AREA BIOGENIC EMISSIONS
Although an awareness of hydrocarbon emissions from naturally-occurring sources has existed for
over two decades, it was not until recently that concerted efforts were undertaken to measure emissions of
hydrocarbon species from various vegetation types (Lamb et al. 1985) and to attempt to compile a biogenic
emissions inventory (Lamb et al 1987). While uncertainties remain in the estimation of biogenic emissions,
it is also necessary to consider the best available biogenic emission rates in combination with anthropogenic
area sources of hydrocarbons and nitrogen oxides in UAM applications.
38
-------
During the development of the ROM 2.1 model, an effort was also underway to develop a new
biogenic emissions processing system to be named BEIS-Biogenics Emissions Inventory System. Briefly, the
methodology in BEIS applied broad vegetation classes and a canopy model used in estimating leaf tempera-
ture and solar intensity profiles within forest stands to derive emissions directly for isoprene, a-pinene,
monoterpene, and unknown organic species. Emission rates of these species were convened into hourly
biogenic emissions for CB-IV organic classes which include; isoprene (ISOP), paraffins (PAR), olefins
(OLE), high molecular weight aldehydes (ALD2), and nonreactive hydrocarbons (NONR). BEIS also
includes a procedure to determine NO and NO2 emission rates from grassland and soil. The allocation
scheme for these species, and the biomass and environmental factors employed in computing the biogenic
emissions for BEIS are fully described by Pierce et al. (1990b) [See also Appendix D to Volume IV of the
User's Guide].
The hourly biogenic area emissions of six species are contained in the retrieved ROM PF144 data file
and include ISOP, PAR, OLE. ALD2, NO, NO2. The tasks performed by the biogenic emissions interface
(IBIOG) are to resolve the gridded biogenic emissions to the UAM grid cells and to combine these values
with the area anthropogenic emissions file supplied by the user. The technique applied to derive UAM
gridded biogenic emissions from the ROM gridded values is the fractional-area weighting method. This is a
similar algorithm as described earlier to resolve the surface roughness and vegetative factor to the UAM grid,
except with a variation needed for its application to emissions.
The ROM gridded biogenic emissions represent emission rates over the area of each ROM grid cell.
When applying the area weighting technique, the ratio of the area of a ROM grid cell contained in a UAM
cell to the total area of the ROM cell is used to scale the ROM biogenic emission rate. This factor
is needed to preserve the emission density (QiA^). For example, in the case-of a UAM cell entirely inside a
ROM grid cell, the UAM cell biogenic emissions would be computed with the ratio of the total area of the
UAM grid cell to the ROM cell's area muitiplied by the ROM cell's biogenic emissions. In the general
application, the area of each ROM cell overlapping a given UAM grid cell is scaled by the total area of the
ROM grid cell (>!r). Then the biogenic emission rate for each species for a particular UAM grid cell is
determined by summing the scaled contributions from every ROM grid cell that overlaps a UAM cell.
The IBIOG interface combines the biogenic emissions for the six species with the corresponding area
anthropogenic emissions of these same species and generates a single binary area emissions which contains
the sum from both inventories.
39
-------
SECTION 5
USER'S INSTRUCTIONS
The interface system is a set of eight batch-oriented computer processors requiring 16 data files listed
in Table 11 and some card inputs described in the following subsections. The processor called PDFSNBK
reads in user-supplied diffusion break data instead of ROM data. The processor IBIOG combines the ROM
biogenic emissions data with user-supplied area-source emissions data. The remainder of the interface pro-
cessors require retrieved ROM data files.
This system produces five formatted data files for UAM preprocessors and five unformatted data files
directly used by UAM. 'IWIND' and 'ICONC' processors have a built-in option for user-supplied data.
Although these eight processors are independent in the interface system, the diffusion break and
regiontop data are required by 'IMETSCL' and TWIND' processors. Therefore, the user should process
diffusion break data first, followed by processing regiontop data. The remaining processors can be executed
in any order. The user needs to supply hourly observed or derived diffusion break heights for the interface
system. The user is referred to Volume II of the User's Guide (Morris et al. 1990b) for the recommended
procedures to obtain a diffusion break data file.
40
-------
TABLE 11. INPUT DATA FILES USED BY EACH INTERFACE PROCESSOR
Interlace processor
Data file pdfsnbk iregntp itmprtr imhtscl iwind icreter iconc biog
DBDATA* XX XXX
EMISSIONS** X
RTDATA*** X X
ROM21 X
MF165 X X
MF166 X
MF174 X
PF102 X
PF103 X
PF10S X
PF114 X
PF115 X
PF117 X
PF118 X
PF119 X XX X
PF144 X
* DBDATA: Diffusion break data (user-supplied).
•" EMISSIONS : Area source emission data (user-supplied).
RTDATA: Regiontopdata ocnved from 'IREGNTP' interlace processor.
(A description of the remaining data files is presented in Table Z)
41
-------
5.1 DIFFUSION BREAK DATA PROCESSOR (PDFSNBK)
5.1.1 Processor Function
This processor prepares hourly diffusion break data for use in the preprocessor 'DFSNBK' of UAM.
PDFSNBK reads in user-supplied diffusion break data and control data pack; then, it generates a formatted
packet file that is in a compatible format for input to the UAM preprocessor 'DFSNBK*. Flow of informa-
tion through the PDFSNBK interface processor is shown in Figure 10.
5.1.2 Inpqt/Output Components
5.1*2.1 Input Flies—
The data contained in input data file (DBDATA) are user-supplied diffusion break data associated
with date and time. The formats are shown below, dots indicate additional program statements. The input
parameters are listed in Table 12.
READC3,17fEN0*998) LDAT,LTIM,VHH
17 FORMAT(17,13,F10.0)
TABLE 12. DBDATA INPUT FILE PARAMETERS
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1 +
1
LDAY
IntegerM
Year/Julian date (e.g., 80203)
->
LTTM
h
Integer*4
Time, LST (e.g., 20)
3
VMH
m
Real*4
Diffusion break height AGL
42
-------
(INPUTS)
| L'SER-DERIVED HOURLY
j DIFFUSION BREAK HEIGHTS
[ (DBDATA)
CONTROL PARAMETERSI
i
PDFSXBK
PROCESSOR
DBPACK
(Formatted)
(OUTPUT)
UAM-DFSNBK
PREPROCESSOR
Note: The hourly diffusion break height file (DBDaTA), derived by the user is the primary input file. The
DBPACK outpui file :s formatted for use in the DFSNBK preDrocessor.
Figure 10. Flow diagram of the diffusion break data processor, PDFSNBIC
43
-------
5.1.2J2 Control Cards—
The control cards define the general information of the UAM model application. Thirty-five variables
included on control cards are used, in the format shown below. Dots indicate additional program statements.
Table 13 defines the control card variables. Additional information about these parameters can be found in
Volume II of the User's Guide (Morris et aL 1990b).
read<5,23> TP
READ<5,20) FIO
REAO<5,21) (I0P(I),I«1,23)
READ(5,22) ISYD,1SHM,IEYD,IEHM
READ<5]24) XLOCR,YIOCR,1ZONE1,XL0C,YL0C,XSZ,YSZ
REAO<5,25) IX,IY,IZ,IZl,1ZU,SFCH,AIH,AUH
READ<5,33) I SUB
READ<5]34) SUB ID, IROU,ICOL,ICONT
READ(5,35> SUB:0,VAR,HHD,VMIN,VHAX,NOP
(If NOP * 0, skip the following READ statements.)
REA0<5,39> PNAME,IPVAL
REA0<5,36> PNAME,PVAL
20 FORMAT(A60)
21 FORHAT<5M0,/,3M0,/,5n0,/,4I10,/,6M0)
22 FORMAT(4I10)
23 FORMAT(A10)
24 FORMAT(2F10.0,I10,/,2F10.0,/,2F10.0)
25 FORMAT(3I10,/,2I1C,F10.1,2F10.0)
33 FORMATC110)
34 FORMAT(A10,3110)
35 FORMAT(3A10,2F10.1,110)
36 FORMAT(A10,F10.2)
39 FORMAT(A1Q,110)
44
-------
TABLE 13. CONTROL CARD VARIABLES FOR PDFSNBK
Card Variable Data
number name Units type Description
1 TP Character*10 File name (Le., DIFFBREAK)
2 FID Character*60 File identifier
3 IOP(I) Integer*4 Control options
1=1 Number of species (set to 0)
1=2 Number of user-defined variables
1=3 Number of stations (i.e., ROM grid cells)
1=4 Number of subregions
1=5 Number of parameters
4 1=6 Output file number
1*7 Input print option-
0: don't print; 1: print
1=8 Output print option--
(0 or 1, as above)
5 1=9 Print unit table (0 or 1)
1=10 Print station locations (0 or 1)
1=11 Print region (0 or 1)
1=12 Print methods table (0 or 1)
1=13 Print station val ues (0 or 1)
6 1=14 Number of vertical parameters
1=15 Number of profile heights
1=16 Print vertical methods table (Oor 1)
1= 17 Print vertical profiles table (0 or 1)
7 1=18 DIFFBREAK unit number
1=19 REGIONTOP unit number
1=20 TOPCONC unit number
1=21 TEMPERATUR unit number
1=22 METSCALARS unit number
1=23 WIND unit number
(continued)
45
-------
TABLE 13. CONTROL CARD VARIABLES FOR PDFSNBK (CONCLUDED)
Card Variable Data
number name Units type Description
8
ISYD
Integer*4
ISHM
Integer*4
IEYD
Integer*4
IEHM
Integer*4
9
XLOCR
m
Real *4
YLOCR
m
Real*4
IZONEl
lnteger*4
10
XLOC
m
Real *4
YLOC
m
Real*4
11
XSZ
m
Real *4
Beginning year/Julian date (yyddd~e.g. 80203)
Beginning time (Wvnm-e.g. 0100, militaiy time)
Ending year/Julian date
Ending time (militaiy time)
Reference origin (x-coordinaie)
Reference origin (y-coordinate)
UTM zone
Origin of grid inx-direction
Origin of grid iny-direction
Cell size in x-direction
YSZ
m
Real*4
Cell size iny-direction
12
IX
IntegerM
Number of cells in x-direction
IY
Integer*4
Number of cells in v-direction
IZ
Integer"4
Number of cells in z-direction
13
IZL
Integer*4
Number of cells in lower layer
IZU
Integer*4
Number of cells in upper layer
SFCH
m
Real *4
Height of surface layer (set to zero)
ALH
m
Real*4
Minimum height of ceil in lower layers
AUH
m
Real *4
Minimum height of cell in upper layers
14
1SUB
Integer*4
Number of subregions (set equal to 1)
15
SUB ID
Charaaer'10
Subregion name
IROW
Integer*4
Beginning row number
ICOL
lnteger*4
Beginning column number
ICONT
IntegerM
Cell count; if negative, equal
to rest of model region
16
SUBID
Character* 10
Subregion name
VAR
Charaaer'10
Variable name
MHD
Charaaer'10
Method name
VMIN
m
Real *4
Minimum value
VMAX
m
Real*4
Maximum value
NOP Integer*4 Number of parameter cards that follow
17+ PNAME
IPVAL
18+ PNAME
PVAL
Charaaer'10
Integer*4
Charaaer'10
Real*4
Parameter name
Parameter value
Parameter name
Parameter value
46
-------
5.1*23 Output File
The output file (DBPACK) contains not only the diffusion break data but also general information of
UAM model application. Dots indicate additional program statements. The output parameters are listed in
Table 14.
WRITE<8,20) C*Tt
WR1TE<8,23) TP
UKITE<8,20) FID
URITE(8,21)
-------
5.2.2.3 Output Files—
'IREGNT?' generates two output files.
RTPACK-This file contains the region top data and general information of the UAM model
application. The formats are listed below; dots indicate additional program statements.
Table 5-20 defines the parameters of RTPACK.
WRITE<25,20) CNTl
WRITE<25,23) TP
URITE<25,20) FID
URtTE(25,21) <10P(I),I«1,23)
WRITE(25,22) ISYD,ISHM,IEYD,IEHN
URITE<25,23) END
URITE<25)iO) REGN
URITE<25,24) XLOCR,YLOCR,I ZONE1,X10C,YLOC,XSZ,YSZ
WRITE(25,25) IX,IY,IZ,IZl,IZU,SFCH,ALH,AUH
URITE(25,23) END
URITE(25,20) TIMTVL
URITE<25,22) IYD1,IHM1fIYD2,IHM2
URITE(25,20) SUBR
URITE<25,34) SUBID,IROU,ICOL,ICONT
WRITE(25,23) EHO
URITE(25,20) HETH
URITE<25,35) SUBID,VAR,MHD,VMIN,VMAX,NOP
(If NOP ¦ 0, skip the next two WRITE statements)
URITE<25,39) pname.ipval
URITE(25,36) PNAME,PVAL
URITE(25,23) END
URITE<25,23) CONST
URITE<25,37) SUBID,VAR,ZTOP
URtTE(25,23) END
URITEC25,23) ENOT
20 FORMAT(A60)
21 FORMAT(5l10,/,3l10,/,5110,/,4110,/,6110)
22 FORMAT(4I10)
23 FORMAT(A10)
24 FORMAT(2F10.0,I10,/,2F10.0,/,2F10.0)
25 POftMAT<3M0,/,2M0,F10.1,2F10,0)
34 FORMAT(A10,3I10)
35 FORMAT(3A10,2F10.1,110)
36 FORMAT(A10,F10.2)
37 FORMAT(A10,A10,F10.1)
39 FORMAT
-------
TABLE 20. RTPACK PARAMETERS
Card
number
Parameter
number
Parameter
name
Units
Data
type
Description
1
1
CNTL
Character*60
Control pack header
2
2
TP
Character* 10
Hie name
3
3
FID
Character*60
File identifier
4-8
4
IOP(I)
Integer*4
Control options
9
5
IS YD
Integer*4
Begimfing year/Julian date
6
ISHM
Integer* 4
Beginning time
7
IEYD
Integer* 4
Ending year/Julian date
8
IEHM
Integer*4
Ending time
10
9
END
Character* 10
Control pack terminator
11
10
REGN
Character*60
Region pack header
12
11
XLOCR
m
. Real*4
Reference origin {x-coordinate)
12
YLOCR
m
Real*4
Reference origin (y-coordinate)
13
IZONE1
Integer*4
UTM zone
13
14
XLOC
m
Real *4
Origin of grid inx-direction
15
YLOC
m
Real*4
Origin of grid iny-direction
14
16
XSZ
m
Real*4
Cell size in x-direction
17
YS2
m
Real*4
Cell size iny-direction
15
18
DC
Integer*4
Number of cells inx-direction
19
IY
Integer*4
Number of cells iny-direcxion
20
IZ
Integer*4
Number of cells in z-direction
16
21
IZL
Integer*4
Number of cells in lower layer
22
IZU
Integer"4
Number of cells in upper layer
23
SFCH
m
Real*4
Height of surface layer
24
ALH
m
Real*4
Minimum height of ceil in
lower layer
25
AUH
m
Real*4
Minimum height of cell in
upper layer
17
26
END
Character* 10
Control pack terminator
18+
27
TINTVL
Character*60
Hourly data header
19+
28
IYD1
Integer*4
Beginning year/Julian date for
hourly data
29
IHM1
Integer*4
Beginning time
30
IYD2
Integer*4
Ending year/Julian date
(continued)
62
-------
TABLE 20. RTPACK PARAMETERS (CONCLUDED)
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
31
IHM2
Integer*4
Ending time
20
32
SUBR
Character*60
Subregion pack header
21
33
SUBID
Character* 10
Subregion name
34
IROW
lntegera4
Beginning row number
35
ICOL
Integer*4
Beginning column number
36
ICONT
Integer*4
Cell count
22
37
END
Character* 10
Control pack terminator
23
38
METH
Charaaer*60
Method pack header
24
39
SUBID
Character* 10
Subregion name
40
var
Character" 10
Variable name
41
MHD
Character" 10
Method name
42
VMIN
m
Real *4
Minimum region top value
43
VMAX
m
Real*4
Maximum region top value
44
NOP
Integer*4
Number of parameter cards that followed
25 +
45
PNAME
Character* 10
Parameter name
46
IPVAL
IntegerM
Parameter value
26+
47
PNAME
Character* 10
Parameter name
48
PVAL
Real*4
Parameter value
27
49
END
Character* 10
Control pack terminator
28+
50
CONST
• Character* 10
Constant pack header
29 +
51
SUBID
Character* 10
Variable name
52
VAR
Character* 10
Subregion name
53
ZTOP
m
Real*4
Region top value
30 +
54
END
Character* 10
Control pack terminator
31 +
55
ENDT
Character* 10
Time interval pack terminator
63
-------
RTDATA-This file contains simple formal of region top data associated with date/time. Table 21
shows the parameters of RTDATA. The formats are listed below.
URITE<9,17) IY01,ISH,ZTOP
17 F08MT(I7,13,F10.1)
TABLE 21. RTDATA PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
1+ 1 rYDl Integer* 4 Year/Julian date
2 ISH Integer*4 Time
3 ZTOP m Real *4 Region top value
64
-------
5*23 Resource Summary for an IREGNTP Application
5*2.3.1 Memory Requirements-
FORTRAN source file:
1
file
13312
bytes
Object file:
1
file
12*288
bytes
Executable file:
J.
file
m*Q
bytes
3
files
35340
bytes
S232 Execution Time Requirements (Representative Values for a 48-h Scenario) ~
IBM 3090
Charged CPU time (hh:mm:ss): 00:00:01
Virtual address space: 3371 K
5.2*3.3 Space Requirements: Log and Print Files—
IREGNTP 3*584 bytes
Print Files: None
5.2.3.4 Space Requirements: Input and Output Files-
Table 22 shows the input file and output file space requirements.
TABLE 22. IREGNTP I/O FILE SPACE REQUIREMENTS
File
group
File
name
File
type
Storage
(in bytes)
Scenario
data span
Input
DBDATA
Formatted
2,048
72 hours
MF166
Formatted
198*144
72 hours
PF119
Formatted
3.584
48 hours
203,776
Output
RTPACK
Formatted
9/728
48 hours
RTDATA
Formatted
1.536
48 hours
11,264
5.23.5 Space Requirements: Tape Files-
None
65
-------
52A Run Stream Command File for an Interface Application
//* J06 CARO
//•
//*
/•ROUTE PRINT HOLD
//•
//STEP1 EXEC PGK*IREGNTP
//STEPLIB X DSN=HHAS.UAH.IMTR FACE.LOAO
//* THESE ARE THE INPUT FILES
//FT03F001 DO DSN*HMAS.UAH.INTRFACE.INPUT.DBDATA,DISP=SHR
//FT01F001 DO DSN«*HAS.UAH.INTRFACE.INPUT.PF119,01SP=SHR
//FT13F001 DO DSN*HHAS,UAH.INTRFACE.INPUT.HF166,DISP«SHR
//• THESE ARE THE OUTPUT FIL£S
//FT25F001 DD DSN».UAH.INTRFACE.OUTPUT.RTPACK,
// DISP»(NEU,CATLCfOELETE),
// SPACE«.UAH.INTRFACE.OUTPUT.RTDATA,
// DISP«(NEV,CATLG,DELETE),
// SPACE-CTRK,C5 f5)),UN IT«SYSOA f
// DCB«
-------
5.2*5 Main Program* Subroutines. Functions, and Block Data Required
5.2^.1 Main Program—
IREGNTP
5.Z3.2 Subroutines-*
None.
5.2.5.3 Functions—
None.
5.2.5.4 Block Data Files—
None.
5.2.6 I/O and Utility Library Subroutines and Functions Required
None.
5.2.7 INCLUDE Files
None.
67
-------
S3 TEMPERATURE INTERFACE (ITMPRTR)
53.1 Processor Function
This processor prepares surface air temperature data for use in.the UAM preprocessor TMPRTR\
ITMPRTR reads in ROM hourly gridded surface air temperature data from PF103 and then calculates UTM
(Universal Transverse Mercator) coordinates corresponding to the midpoint of each ROM grid cell. It also
reads in control data that contain the general information on the UAM model application. This processor
then generates a formatted packet file that is in a compatible format for input to UAM preprocessor
TMPRTRV Flow of information through the ITMPRTR interface processor is shown in Figure 12.
53.2 Input/Output Components
5.3.2.1 Input Files—
'ITMPRTR' requires a retrieved ROM MF/PF file, PF103. This data file contains hourly gridded
surface air temperature values. The formats are shown below; dots indicate additional program statements.
The input parameters are listed in Tabie 23.
READ(7,12) TlA,IMX1,TlB,IMri,T1C,ISPC
READ(7,13) T2,IDLAT1,IMLAT1,SIAT1
READ(7,13) 72,IDL0N1,1HL0N1,SL0N1
READ<7,15) T3
READ(7,300,END=999) LDATE.LTIHE
READ (7,,30, END *999) (VAL( I,J),1*1,1NX)
12 F0RHAT
-------
(INPUTS)
REGIONAL OXIDANT MODEL |
(surface air'temperatures'
PF103
I
ITMPRTR I
INTERFACEi
CONTROL PARAMETERS
I
TPPACK i
i (Formatted) '
(OUTPUT)
i
UAM-TMPRTR
! PREPROCESSOR
Note: The TPPACK output Glc is formatted for use in the TMPRTR preprocessor program.
Figure 12. Flow diagram of the rTMPRTR interface program with input and output files.
69
-------
TABLE 23.
PF103 PARAMETER LIST FOR ITMPRTR
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
T1A
Character*52
Subtitle name
2
IMX1
lnteger*4
Number of columns
3
TIB
Character*8
Subtitle name
4
IMY1
Integer*4
Number of rows
5
TIC
Character*14
Subtitle name
6
ISPC
Integer'4
Number of parameters
2
7
T2
Character*27
Subtitle name
8
IDLAT1
Integer*4
Degrees oflatitude for
southwest corner of
ROM retrieved subregion
9
IMLAT1
Integer*4
Minutes oflatitude
10
SLAT1
Real *4
Seconds oflatitude
3
11
T2
Character*27
Subtitle name
12
EDLONl
Integer*4
Degrees of longitude
•
13
IMLON1
Integer*4
Minutes of longitude
14
SLON1
Real *4
Seconds of longitude
4+
15
T3
Character* 80
Header information
5+
16
LDATE
Inieger*4
Year/Julian date
17
LTIME
h
lnteger*4
Time
6+
18
VAL(IJ)
K
Real"4
Temperature values
I: Index of columns
J: Index of rows
70
-------
533.2 Control Cards-
The control cards define the general information of UAM model application. Forty-two variables
included in twenty-three control cards are used, in the format shown below. Dots indicate additional
program statements. Table 24 defines the control card variables.
READ(5,23) TP
READ(5,20) FID
READ(5,21> (I0P
READC5.22) ISYD,ISHH,IEYD,IEHH
READ(5,24) XLOCR,YLOCR,IZONE1,XLOC,YLOC,XSZ,YSZ
READ(5,25) IX,IY,IZ,IZl,IZU,SFCH,ALH,AUH
REA£>(5*23) UNITCO
READ(5,26) FACTM,FACTA,VTMOL
REA0(S]28) HTZ,I NX,I NY
READ(S,33) I SUB
READ(5,34) SUB ID,I ROW,I COL,ICONT
READ(5,35) SUBID,VAR,MHD,VMlN,VHAX,MOP
READ(5,39) pname,ipval
READ<5,36) PNAME,PVAL
20 FORMAT(A60)
21 FORMAT(5MO,/,3I10,/,5l1O,/,4I1O,/,6l1O)
22 FORMAT(4110)
23 FORMAT(AIO)
24 FORKAT(2F10.0,!10,/,2F10.0,/,2F10.0)
25 FORMAT(3110,/,2110,F10.1,2F10.0)
26 FORMAT(3F10.1)
28 FORMAT(3X,F7.1,3110)
33 FORMAT(110)
34 FORMAT(A10,3I10)
35 FORMAT(3A10,2F10.1,110)
36 FORMAT(A10,F10.2)
39 FORMAT(A10,110)
-------
TABLE 24. CONTROL CARD VARIABLES FOR ITMPRTR
Card Variable Data
number name Units type Description
1
TP
Character* 10
File name (Le., TEMPERATURE)
2
FID
Character*60
File identifier
3
IOP(I)
IntegerM
Control options
I«1
Number of species
1=2
Number of user-defined variables
1-3
Number of stations (ROM grid points)
1=4
Number of subregions
1=5
Number of parameters
4
1=6
Output file number
1=7
Input print option
0: don't print; 1: print
1=8
Output print option
(0 or 1, as above)
5
1=9
Print unit table (0 or 1).
1=10
Print station locations (0 or 1)
1=11
Print region (0 or 1)
1=12
Print methods table (0 or 1)
1=13
Print station values (0 or 1)
6
1=14
Number of vertical parameters
1 = 15
Number of profile heights
1 = 16
Print vertical methods table (0 or 1)
1=17
Print vertical profiles table (0 or 1)
7
1=18
DIFFBREaK unit number
1 = 19
REGIONTOP unit number
1=20
TOPCONC unit number
1=21
TEMPERATUR unit number
1=22
METSCALARS unit number
1=23
WIND unit number
3
ISYD
Integer *4
Beginning year/Julian date (yyddd—e.g. 80203)
ISHM
Integer*4
Beginning time (hhmm—e.g. 0100, military time)
IEYD
Integer*4
Ending year/Julian date
EEHM
Integer*4
Ending time
9
XLOCR
m Real *4
Reference origin (x-coordinate)
(continued)
72
-------
TABLE 24. CONTROL CARD VARIABLES FOR ITMPRTR (CONCLUDED)
Card
Variable
Data
number
name
Units
type
Description
YLOCR
m
Real'4
Reference origin (y-coordinate)
IZONE1
Integer*4
UTM zone
10
XLOC
m
Real *4
Origin of grid in jxiirection
YLOC
m
Real*4
Origin of grid In ^-direction
11
xsz
m
Real*4
Cell size in x-direction
YSZ
m
Real *4
Cell size in y-direction
12
IX
lnteger*4
Number of cells in j:-direction
ry
IntegerM
Number of cells iny-direction
rz
Integer*4
Number of cells in z-direction
13
izl
lnteger*4
Number of cells in lower layer
IZU
Integer*4
Number of ceils in upper layer
SFCH
m
Reai*4
Height of surface layer
ALH
m
Real*4
Minimum height of cells in lower layer
AUH
m
Real*4
Minimum height of cells in upper layer
14
UNITCD
Character* 10
Unit name
15
• FACTM
Real*4
Multiplicative factor
facta
Real"4
Additive factor
wtmol
Real*4
Molecular weight
16
HTZ
m
Real*4
Height of station
INX
IntegerM
Number of cells injc-direction (ROM)
INY
Integer*4
Number of cells iny-direction (ROM)
17
ISUB
Integer*4
Number of subregions
18
SUBID
Character* 10
Subregion name
IROW*
IntegerM
Beginning row number
ICOL
integer*4
Beginning column number
ICONT
Integer*4
Cell count; if negative, equal
to rest of model region
19
SUBID
Character* 10
Subregion name
VAR
Character* 10
Variable name
MHD
Character*10
Method name
VMIN
m
Real*4
Minimum value
VMAX
m
Real*4
Maximum value
NOP
Integer*4
Number of parameter cards that follow
20+
PNAMB
Character* 10
Parameter name
IPVAL
Integer *4
Parameter value
21 +
PNAMB
Character* 10
Parameter name
PVAL
Real*4
Parameter value
73
-------
5-3.23 Output Files—
TTMPRTR* generates one output file, TPPACK. This file contains ROM gridded surface air temper-
ature associated with UTM coordinates and general information of UAM model application. Dots indicate
additional program statements. Table 25 shows the parameters of TPPACK.
URITE(23,20) CNTL
WRITE(23,23) TP
URITE(23,20) FID
WRITE(23,21) (IOP(I),1-1,23)
URITE(23,22) ISYO,ISHM,IEYD,IEHM
URITE(23,23) END
URITE(23,20) REGN
URITE(23,24) XlOCR,YlOCR,IZ0NE,X10C,Y10C,XS2,YSZ
URITE(23,25) IX,IY,12,121,IZU,SFCH,AIH,AUH
URITE(23,23) END
WRITE(23,20) UNIT
WRITE(23,27) TP,UNITCD,FAC7M,FACTA,WTMOL
URITE(23,23) END
URITE<23,20) STAN
URITE(23,29) !,XX(I),TT(I>,HTZ
URITE<23,23) END
VRITE(23]20) T1NTVL
WRITE(23^22) IYD1,IHM1,ITD2,IHH2
WR[TE<23,20) SUBR
WRITE(2i,34) SUB ID,IROU,I COL,ICONT
WRITE(23,23) END
WRITE(23,20) METH
URITE(23,35) SU8ID,VAR,MHO,VMIN,VMAX,NOP
WRITE(23,39) PNAHE,IPVAL
WRITE(23,36> PNAHE,PVAL
WRITE(23,23) END
URITE(23,20) STRD
WRITEC23,37) K,VAR, VAU (K)
WRITE(23,23) END
WRITE(23,23) ENOT
20 FORHAT(A60)
21 FORMAT(5!10,/,3l10,/,5n0,/,4I10,/,6I10)
22 FGRMAT<4I10)
23 FORHAT(A10)
24 FORHAT(2F1O.0,I10,/,2F10.0r/,2Fl0.0)
25 FORHAT<3I10,/,2I10,F10.1,2F10.0)
27 FO#mAT<2A10,3F10.1)
29 FORMAT(I4.3,6X,2F10.0,F10.1)
34 FORMAT(A10,3MO)
35 FORHAT(3A10,2F10.1,110)
36 FORHAT(A10,F10.2)
37 FORMAT
-------
TABLE 25. TPPACK PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
1
1
CNTL
Character*60
Control pack header
2
2
TP
Character* 10
File name
3
3
FID
Character*60
Hie identifier
4-8
4
IOP(I)
Integer*4
Control options
9
5
ISYD
Integer*4
Beginning year/Julian date
6
ISHM
Integer*4
Beginning time
7
IEYD
Integer*4
Ending year/Julian date
8
IEHM
Integer*4
Ending time
10
9
END
Character* 10
Control pack terminator
11
10
REGN
Character*60
• Region pack header
12
11
XLOCR
m
Real*4
Reference origin (x-coordinate)
12
YLOCR
m
Reai*4
Reference origin (y-coordinaie)
13
IZONE
Integer*4
UTM zone
13
14
XLOC
m
Real *4
Origin of grid inx-direction
15
YLOC
m
Real*4
Origin of grid in^v-direction
14
16
xsz
m
Real *4
Cell size in x-direction
17
YSZ
m
ReaT4
Cell size iny-direction
15
18
IX
Integer*4
Number of cells injr-direction
19
IY
Integer*4
Number of cells in v-direction
20
IZ
Integer*4
Number of cells in z-direction
16
21
IZL
Integer*4
Number of cells in lower layer
22
IZU
Integer*4
Number of cells in upper layer
23
SFCH
m
Real*4
Height of surface layer
24
ALH
m
Real*4
Minimum height of cells in lower layer
25
AUH
m
Real*4
Minimum height of cells in upper layer
17
26
END
Character* 10
Control pack terminator
18
27
UNIT
Character*60
Unit pack header
19
28
TP
Character* 10
Hie name
29
UNITCD
Character*10
Unit name
30
FACTM
Real*4
Multiplicative factor
31
FACTA
Real *4
Additive factor
32
WTMOL
Reai*4
Molecular weight
20
33
END
Character*10
Control pack terminator
21
34
STAN
Character*60
Station pack header
22
35
I
Integer*4
Index of stations
(continued)
76
-------
TABLE 25. TPPACK PARAMETERS (CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
36
XX(I)
m
Real*4
x-Iocation in UTM
37
YY(I)
m
Real*4
^-location in UTM
38
HTZ
m
Real*4
Height of station
23
39
END
Character* 10
Control pack terminator
24+
40
TINTVL
Character *60
Hourly data header
25+
41
IYD1
Integer*4
Beginning year/Julian date for hourly
42
1HM1
Integer*4
data
43
IYD2
Integer*4
Beginning time
44
IHM2
Integer*4
Ending year/Julian date
26
45
SUBR
Charaaer*60
Ending time
Subregion pack header
27
46
SUBID
Character'10
Subregion name
47
IROW
Integer*4
Beginning row number
48
ICOL
Integer*4
Beginning column number
49
ICONT
Integer*4
Ceil count
28
50
END
Character'10
Control pack terminator
29
51
METH
Charaaer*60
Method pack header
30
52
SUBID
Character* 10
Subregion name
53
VAR
Character* 10
Variable name
54
MHD-
Character* 10
Method name
55
VMIN
K
Rcal"4
Minimum temperature value
56
VMAX
K
Real*4
Maximum temperature value
57
NOP
Integer*4
Number of parameter cards that follow
31 +
5$
PNAME
Characier*10
Parameter name
59
IPVAL
integer'4
Parameter value
32+
60
PNAME
Character* 10
Parameter name
61
PVAL
Real*4
Parameter value
33
62
END
Charaaer*10
Control pack terminator
34+
63
STRD
Character *60
Station reading header
35 +
64
K
Integer*4
Index of station
65
VAR
Character* 10
Variable name
66
VAL1(K)
K
Reai*4
Temperature values
36+
67
END
Character'10
Control pack terminator
37 +
68
ENDT
Character* 10
Time interval pack terminator
77
-------
S33 Resource Summary for an ITMPRTR Application
533.1 Memory Requirements—
FORTRAN source file:
1
file
38,912
bytes
Object file:
1
file
14,848
bytes
Executable file:
1
file
nm
bytes
3
files
66,048
bytes
5.3.3.2 Execution Time Requirements (Representative Values for a 48-h Scenario)—
IBM 3090
Charged CPU time (hhrmmrss): 00:00:01
Vinual address space: 9468 K
5333 Space Requirements: Log and Print Files—
ITMPRTR 11,264 bytes
Print Files: None
533A Space Requirements: Input and Output Files-
Table 26 shows the input file and output file space requirements.
TABLE 26, ITMPRTR I/O FILE SPACE REQUIREMENTS
File File File Storage Scenario
group name type (in bytes) data span
input PF103 Formatted 132,096 48 hours
Output TPPACK Formatted 364,032 48 hours
533.5 Space Requirements: Tape Files—
None
78
-------
53.4 Run Stream Command File for an Interface Application
//* JOB CARD
/r
//*
/~ROUTE PRINT HOLD
//STEP! EXEC PQMsJTHPRTR
//STEPLIB DO DSN*MMAS.UAM.INTRFACE.LOAD
//* THIS IS THE INPUT FILE
//FT07F001 DO DSN«NHAS.UAM.INTRFACE.INPUT.PF103,DISP*SHR
//• THIS IS THE OUTPUT FILE
//FT23F001 DO DSN».UAM.INTRFACE.OUTPUT.TPPACK,
// DISP*(NEU,CATLG,DELETE),
It SPACE*(TRK,(100,10)),UNIT«STSOA ,
if DCB»(RECFM«FB,LRECLs132,BLKSI2E»2640)
//•
//* THIS IS THE CONTROL PACKET
//FT05F001 DD *
TEHPERATUR
TEMPERATURE
FILE CREATED
AT 02/07/90
0
0
224 1
10
28
0
0
0
0
0 0
0
0
0
0 0
0
0
0 0
0
80203
00
80204 2400
0.
0.
18
520000.
4460000.
8000.
8000.
31
25'
5
2
3
0.0 50.0
100.0
DEGK
1.
0.
0.
4.0
1
16
14
A 11-1
A TEHPERATURSTATINTERP 270.0 330.0 4
EXTENT 10.0
INITRAD1US 2.0
RAD IUSI NCR 1.0
MAXRADIUS 5.0
/*
//* THIS IS END OF DATA
79
-------
5.3.5 Main Program, Subroutines, Functions, and Block Data Required
5.3.5.1 Main Program—
ITMPRTR
5.3^.2 Subroutines—
UTM2
5.3.5.3 Functions—
None.
5.3.5.4 Block Data Files—
None.
5.3.6 I/O and Utility Library Subroutines and Functions Required
None.
53,7 INCLUDE Files
None.
80
-------
5.4 METSCALARS INTERFACE (IMETSCL)
5.4.1 Processor Function
This processor prepares meteorological scalar data for use in the UAM preprocessor 'METSCL'.
IMETSCL reads in ROM gridded terrain elevations, hourly gridded fractional sky coverages, hourly interpo-
lated vertical temperature profiles from upper air site(s), hourly gridded water vapor concentrations, hourly
domain-average region top heights from IREGNTP processor, and user-supplied diffusion break data.
Atmospheric pressure has been set to 1.0 atm for all hours. The hourly vertical temperature gradients for
below and above the diffusion break are calculated by averaging the individual temperature gradients at 50-m
intervals from the ROM upper air data file. Water vapor concentration is calculated by averaging ROM
hourly gridded water vapor concentrations over the entire UAM domain. The exposure class is determined
by the solar zenith angle and cloud cover. The location of the middle of the UAM domain is used for solar
zenith angle calculation and domain-wide average cloud cover is also applied in this scheme. The values used
for NO2 photolysis rate constants are based on Demerjian et al. (1980) and use the date, time, diffusion
break, and location. This processor then generates a fprmatted packet file that is in a compatible format for
input to the UAM preprocessor 'METSCL'. Flow of information through the IMETSCL interface processor
is shown in Figure 13.
5.4.2 Innut/Output Components
5.4.2.1 Input Files-
'IMETSCL' requires four retrieved ROM MF/PF files, one derived file from IREGNTP, and one
user-supplied file.
81
-------
(INPUTS)
ROM
TEMPERATURE
PROFILES
PF102
ROM
CLOUD !
FRACTIONS
PF1I7 i
ROM
TERRAIN
ELEVATION
PF i 19
ROM
WATER VAPOR
CONCENTRATIONS
MF174
{dbdataI
RTDATA
CONTROL
parameters!
IMETSCL
INTERFACE
MSPACK
(Formatted)
(OUTPUT)
I
UAM-METSCL
PREPROCESSOR
Note: The MSPaCK output fiie is formatted for use in the METSCL preprocessor program.
Figure 13. Flow diagram for the IMETSCL interface program with input and output files.
82
-------
PF102-This data Hie contains hourly vertical temperature profile data. The formats are shown below;
dots indicate additional program statements. The input parameters are listed in Table 27.
READ<12,12) TU,!MX4,T1B,!MY4,T1C,2SPC
READ<12,13) T2,IDLAT4,IMLAT4,StAT4
READ(12,13) T2,1DL0N4,IM10W4,SL0K4
REA0O2,15) T3
READ<12,111) 1ST
READ<12,300,END*996) LDAT5,LTIM5
READ(12,301,END=996) LEV(IP)
READ<12,302,EN0=996) HTCKB.IP)
REA0C12,111) 1ST
READ( 1*2,300,END*996) LDAT5.LTIM5
READ<12,301,END«996) LEV
-------
TABLE 27, PF102 PARAMETERS FOR IMETSCL
Card Parameter Parameter Data
number number name Units type Description
T1A
IMX4
TIB
IMY4
TIC
Character*52
Integer*4
Character*8
Integer*4
Character* 14
Subtitle name
Number of columns
Subtitle name
Number of rows
Subtitle name
6
7
8
9
10
ISPC
T2
IDLAT
IMLAT4
SLAT4
Integer*4 Number of parameters
Character*27 Subtitle name
Integer*4 Degrees of latitude for southwest
corner of ROM domain
Integer*4 Minutes of latitude
Real*4 Seconds of latitude
4 +
11
12
13
14
15
T2
IDLON4
IMLON4
SLON4
73
Characier*27
Integer*4
Integer*4
Real*4
Character*80
Subtitle name
Degrees of longitude
Minutes of longitude
Seconds of longitude
Header information
:>+
6+
8+
16
17
18
19
20
1ST
LDAT5
LTIM5
LEV(IP)
HT(MNT)
m
Integer*4
Integer*4
Integer*4
Integer"4
Real*4
Number of stations
Year/Julian date
Time
Number of levels-
IP : Index of stations
Height (MSL)~
M: Index of levels
N: Index of stations
9-r
10-5-
11 +
12+
21
22
23
24
25
1ST
LDAT5
LTIM5
LHV(IP)
lnteger*4
lnteger*4
lnteger*4
Integer*4
TEMP(M,N) K Real*4
Number of stations
Year/Julian date
Time
Number of levels—
IP: Index of stations
Vinuai temperature-
M &. iV are same as above
13+
14+
15+
16+
26
27
28
29
30
1ST
LDAT5
LTIM5
LEV(IP)
RM(MtN)
Integer*4
Integer*4
Integer*4
Integer*4
Real*4
Number of stations
Year/Julian date
Time
Number of levels-
IP : Index of stations
Mixing ratio fraction-
M & N are same as above
84
-------
PF117~This data file contains hourly gridded fractional cloud cover data. The formats are shown
below; dots indicate additional program statements. The input parameters are listed in Table 28.
READ<10,12> TU,IMX3,T1B,IMr3,T1C,tSPC
READ<10,13) T2,IDUT3,IMIAT3,SIAT3
REAOC10J3) T2, IDL0N3, IMLON3,SLON3
READ(10,15) T3
READ(10,300,EN0«997) LDAT3,ITIW
READ(10,30,E«>®997)
-------
PF119-This data file contains gridded terrain elevation data. The formats are shown below; dots
indicate additional program statements. The input parameters are listed in Table 29.
R£AD(1,12) T1A,IMX1,TlB,l*ri,TlC,ISPC
READ(1,13) T2,IDLAT1,1NLAT1,SLAT1
R£AD(1,13) T2,I0L0N1,IMI0M1,S10N1
READ(1,15) T3
READ(1,300,EMD=999) L0AY1,LTIM1
READ(1,30,EM>*999)
-------
MF174~This data file contains hourly gridded water vapor concentration of layer 1. The formats are
shown below; dots indicate additional program statements. The input parameters are listed in
Table 30.
READ(4,12) T1A,!HX2,T1B,1MY2,T1C,I$PC
READ<4,13) 72,IDLAT2,!MLAT2,SLAT2
READ<4,13) 72,ID10N2,IML0N2,SL0N2
READ(4,15) T3
READ(4,J00,END*996) LDAT2,LTIM2
READ(4,30,END*996) (VALCI,J),I»1,INX)
12 F0RMAT(AS2,I2,A8,I2,A14,I2)
13 FORMAT(A27,14,1X,12,1X,F5.2)
15 FORMAT(A80)
300 F0RMAT
-------
A
DBDATA-The data contained in this file are user-supplied diffusion break data associated with date
and time. The formats are shown below; dots indicate additional program statements. The input
parameters are listed in Table 31.
READ(3,17,END*998) LDAY,LT1M,VMH
17 FORMAT(17,13,F10.0)
TABLE 31. DBDATA PARAMETERS FOR IMETSCL
Card Parameter
number number
Parameter
name
Units
Data
type
Description
1 +
LDAY
LTTM
VMH
Integer*4
IntegerM
Rcal*4
Year/Julian date
Time
Diffusion break value
RTDATA-This file contains simple format of region top data associated with date/time. The formats
are shown below; dots indicate additional program statements. The input parameters are listed in
Table 32.
REA0<9,17,END»998) LDAT6,LTIM6,RT0P
17 FORMAT(I7,I3,F1Q.O)
TABLE 32. RTDATA PARAMETERS FOR IMETSCL
Card Parameter Parameter Data
number number name Units type Description
1+ 1
2
3
LDAT6 Integer*4 Year^ulian date
LTfM6 h Integer*4 Time
RTOP ra Real*4 Region top value
88
-------
5.4.2.2 Control Cards—
Twenty-nine variables included in fifteen control cards are used, in the formats shown below. It
defines The general information of UAM model application. Dots indicate additional program statements.
The control card variables are listed in Table 33.
READC5,63) ALA,A10,TZ
REAOC5,23) TP
READ(5,20) FID
READ<5,21) (I0P<1),1*1,23)
READ<5,22) isyd,ish«,ieyd,iehm
READ(5,24) XLOCR,YLOCR,:ZONE1,XLOC,YLOC,XSZ,YSZ
READ(5,25) IX,1Y,JZ,IZL,IZU,SFCH,ALH,AUH
READ(5,28) HTZ,I NX,INY,ATMOV
20 FORMAT(A60)
21 FORMAT(5l10,/,3l10,/,5M0,/,4I10,/,6I10>
22 F0RMATC4M0)
23 FORMAT(AIO)
24 FGRMAT(2F10.0,I10,/,2F10.0,/,2F10.Q)
25 F0RKAT(3l10,/,2I10,F10.1,2F10.0)
28 FORMAT(F10.1,2110,F10.4)
63 FORMAT(3F10.0)
89
-------
TABLE 33. CONTROL CARD VARIABLES FOR IMETSCL
Card Variable Data
number name Units type Description
1
ALA
*N
Real*4
Latitude of middle of UAM domain
ALO
• W
Real*4
Longitude of middle of UAM domain
TZ
Real *4
Time zone for computation of solar zenith
angle (relative to Greenwich Mean Time);
Example: 5 = EST
2
TP
Character* 10
File name (i.e. METSCALARS)
3
FID
Character*60
File identifier
4
IOP(I)
Integer*4
Control options
1=1
Number of species
1=2
Number of user-defined variables
1=3
Number of stations (ROM grid points)
1=4
Number of subregions
1=5
Number of parameters
5
1-6
Output file number
1=7
Input print option—
0: don't print; 1: print
1=8
Output print option-
(0 or 1, as above)
6
1=9
Print unit table (0 or 1)
1=10
Print station locations (0 or 1)
1=11
Print region (Oorl)
1=12
Print methods table (Oor 1)
1=13
Print station values (Oorl)
7
1=14
Number of vertical parameters
1=15
Number of profile heights
1=16
Print vertical methods table (Oor 1)
1 = 17
Print vertical profiles table (0 or 1)
8
1=18
DIFrBREAKunit number
1=19
REGIONTOP unit number
1=20
TOPCONC unit number
1=21
TEMPERATUR unit number
1=22
METSCALARS unit number
1=23
WIND unit number
i
ISYD
Integer*4
Beginning year/Juiian date (yyddd-c.g. 80203)
ISHM
Integer*4
Beginning time (hhmm—e.g. 0100,
military time)
(continued)
-------
TABLE 33, CONTROL CARD VARIABLES FOR IMETSCL (CONCLUDED)
Card
Variable
Data
number
name
Units
type
Description
IEYD
Integer*4
Ending yearflulian date
IEHM
Integer*4
Ending time
10
XLOCR
m
Real*4
Reference origin (r-coordinate)
YLOCR
m
Real*4
Reference origin (y-coordinate)
1ZONE1
Integer*4
UTMzone
11
XLOC
m
Real *4
Origin of grid inx-direction
YLOC
m
Real *4
Origin of grid in>-direcxion
12
XSZ
m
Real*4
Cell size inx-direction
YSZ
m
Real*4
Cell size iny-direction
13
IX
IntegerM
Number of cells in-r-direction
IY
Integer"4
Number of cells in v-direction
IZ
IntegerM
Number of cells in z-direction
14
IZL
IntegerM
Number of cells in lower layer
IZU
IntegerM
Number of cells in upper iayer
SFCH
m
RealM
Height of surface layer
ALH
m
RealM
Minimum height of cell in lower layer
AUH *
m
RealM
Minimum height of cell in upper layer
15
HTZ
m
RealM
Height of station
1NX
IntegerM
Number of ROM cells inx-direction
INY
IntegerM
Number of ROM cells iny-direction
ATMOV
atm
RealM
Atmospheric pressure
Example:
41.
73.
5.
METSCAIARS
METSCALARS
FILE CREATED
AT 02/07/90
0
0
224
1
10
28
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
50203
0
80204
2400
0.
0.
18
520000.
4460000.
8000.
8000.
31
25
5
2
3
0.0
50.0
100.0
4.0
16
14
1.0000
91
-------
5.4.2.3 Output FiJes—
'IMETSCL' generates one output file.
MSPACK-This file contains six meteorological scalars and general information on the UAM model
application. The formats are shown below; dots indicate additional program statements. The output
parameters in MSPACK are listed in Table 34.
WRITE(24
20)
cm
WRITE(24
23)
TP
WRITE(24
20)
FIO
WRITE{24
21)
(IOP(I),I«1,23)
WRITE(24
22)
IYDfl,IHH8,IEYD,IEHM
URITEC24
23)
END
URITEC24
20)
REGN
URITEC24
24)
XLOCR,YLOCR,IZ0NE1,
WRITE(24
25)
IXrIY,I2,I2L,I2U,SF
WRITE(24
23)
ENO
URITEC24
20)
TINTVL
URITE.C24
22)
irD8fIHM8,IYD9,IHM9
WRITE(24
23)
TITLE
WRITE(24
37)
TBEL,RTDN1
WRITE(24
37)
TABV,RTUP1
WRITE(24
47)
EXPC,EXP
WRITE(24
37)
RADF,RK1
WRITE(24
37)
CONC,H20
UR1TE(24
37)
ATMO,ATHOV
WRITE(24
23)
ENO
URITE(24
23)
ENDT
WRITE(24
20)
TINTVL
WRITE(24
22)
IYD1,IHM1,IYD2,IHM2
WRITE(24
23)
TITLE
WRITE(24
37)
T8EL,RTDM1
WRITE(24
37)
TABV,RTUP1
WRITE(24
47)
EXPC,EXP
WRITE(24
37)
RAOF,RK1
WRITE(24
37)
C0NC,H20
WRITE(24
37)
ATM0,ATMOV
WRITE(24
23)
ENO
WRITE(24
23)
ENOT
20
FORMAT(A60)
21
FORMAK5I10,/
3110,/,5I10,/,4110,,
22
FORMAT(4I10)
23
FORMAT(AIO)
25 FORMAT(3H0,/,2n0,F1Q.1,2F10.0)
37 FORMAT(A10,F10.4)
47 FORMAT(A10,F10.1)
92
-------
TABLE 34. MSPACK PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
1
1
CNTL
Character*60
Control pack header
2
2
TP
Character* 10
File name
3
3
FID
Character*60
File identifier
4-8
4
IOPfl)
Integer*4
Control options
9
5
IYD8
Integer*4
Beginning year/Julian date
6
IHM8
lnteger*4
Beginning time
7
IEYD
Integer*4
Ending year/Julian date
S
IEHM
Integer*4
Ending time
10
9
END
Character* 10
Control pack terminator
11
10
REGN
Character *60
Region pack header
m
X 4.
11
XLOCR
m
Real*4
Reference origin (^-coordinate)
12
YLOCR
m
Reai*4
Reference origin (y-coordinate)
13
IZONE1
Integer*4
UTM zone
13
14
XLOC
m
Real *4
Origin of grid inx-airection
15
YLOC
m
Real*4
Origin of grid in>-direaion
14
16
xsz
m
Real*4
Cell size in j:-direction
17
YSZ
m
Real*4
Cell size in ^-direction
15
18
IX
Integer*4
Number of cells in x-direction
19
IY
Integer*4
Number of cells iny-direction
• 20
IZ
Integer*4
Number of cells in z-direction
16
21
. IZL
Integer*4
Number of cells in lower layer
22
IZU
Integer*4
Number of cells in upper layer
23
SFCH
m
Real*4
Height of surface layer
24
ALH
m
Real *4
Minimum height of cell in lower layer
25
AUH
m
Real*4
Minimum height of cell in upper layer
17
26
END
Character* 10
Control pack terminator
18
27
TINTVL
Charaaer*60
Hourly data header
19
28
IYD8
Integer*4
Beginning year/Julian date
29
IHM8
Integer*4
Beginning time (1 h before real data
recording hour)
30
IYD9
Integer*4
Ending year/Julian date
31
IHM9
lnteger*4
Ending time
20
32
TITLE
Character* 10
Title name
21
33
TBEL
Character* 10
Parameter name
34
RTDN1
K/m
Real *4
Temperature gradient for below
diffusion break height
22
35
TABV
Character* 10
Parameter name
(continued)
93
-------
TABLE 34. MSPACK PARAMETERS (CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
36
RTUP1
K/m
Real*4
Temperature gradient for above
diffusion break height
23
37
EXPC
Character* 10
Parameter name
38
EXP
Real*4
Exposure class
24
39
RADF
Character* 10
Parameter name
40
RK1
miir1
Real*4
NO2 photolysis rate
constant
25
41
CONC
Charaaer*10
Parameter name
42
H20
ppm
Real*4
Water vapor concentration
26
43
ATMO
Character* 10
Parameter name
44
ATMOV
arm
ReaI-4
Atmospheric pressure
27
45
END
Character* 10
Control pack terminator
28
46
ENDT
Character* 10
Time interval pack terminator
29+
47
TINTVL
Character*60
Hourly data header
30+
48
IYD1
Integer*4
Beginning year/Julian date for
hourly data
49
IHM1
Integer*4
Beginning time
50
IYD2
Integer*4
Ending year/Julian date
51
IHM2
Integer*4
Ending time
31 +
52
TITLE
Character*10
Title name
32+
53
TBEL
Charaaer*10
Parameter name
54
RTDN1
K/m
Real*4
Temperature gradient for below
diffusion break height
33+
55
TABV
Character* 10
Parameter name
56
rtup:
K/m
Real *4
Temperature gradient for above
diffusion break height
34+
57
EXPC
Character* 10
Parameter name
58
EXP
Real*4
Exposure class
35+
59
RADF
Character*10
Parameter name
60
RK1
min*1
Real*4
NO2 photolysis rate constant
36 +
61
CONC
Character* 10
Parameter name
62
H20
ppm
Reai*4
Water vapor concentration
37+
63
ATMO
Character*10
Parameter name
64
ATMOV
atm
Real*4
Atmospheric pressure
38+
65
END
Character*10
Control pack terminator
39+
66
ENDT
Character*10
Time interval pack terminator
94
-------
5.4.3 Resource Summary for an IMETSCL Application
54.3.1 Memory Requirements—
FORTRAN source file:
1
file
41,472
bytes
Object file:
1
file
25,600
bytes
Executable file:
file
19,456
bytes
3
files
86,528
bytes
5.4.3.2 Execution Time Requirements (Representative Values for a 48-h Scenario)—
IBM 3090
Charged CPU time (hh:mm:ss): 00:00:02
Virtual address space: 9512 K
5.4.3.3 Space Requirements: Log and Print Files—
IMETSCL 9,216 bytes
Print Files: * None
5.4.3.4 Space Requirements: Input and Output Files-
Table 35 shows the input file and output file space requirements.
TABLE 35. IMETSCL I/O FILE SPACE REQUIREMENTS
File
group
File
name
File
type
Storage
(in bytes)
Scenario
data span
Input
PF102
Formatted
474,112
48 hours
PF117
Formatted
132,096
48 hours
PF119
Formatted
3,584
48 hours
MF174
Formatted
132,096
48 hours
DBDATA
Formatted
2,048
72 hours
RTDATA
Formatted
1436
48 hours
745,472
Output
MS PACK
Formatted
1899
48 hours
5.4.3.5 Space Requirements: Tape Files—
None
95
-------
5,4.4 Run Stream Command Flic for an Interface Application
/r J 06 CARD
/r
//*
/•ROUTE PRINT HOLD
//*
//STEP1 EXEC PGM«IMETSCL
//STEP11B DO DSN«MMAS.UAM.INTRFACE.LOAD
//* THESE ARE THE INPUT FILES
//FT03F001 DD DSN«HMAS.UAM.INTRFACE.INPUT.OBOATA,01SP*SHR
//FT04F001 DD DSN-MMAS.UAM.INTRFACE.INPUT.MF174,01SP-SHR
//FT09F001 DD DSN-MMAS.UAM.INTRFACE.INPUT.RTDATA,DISP*SHR
//FT12F001 DD DSN*MMA$.UAM.INTRFACE.INPUT.PF102,DISP'SHR
//FT01F001 DD DSN«*WAS.UAM.INTRFACE.INPUT.PF119,DISP»SHR
//FT10F001 DD DSN«MMAS.UAM.INTRFACE.INPUT.PF^17,DISP«SHR
//* THIS IS THE OUTPUT FILE
//FT24F001 DD DSN».UAM. INTRFACE.OUTPUT.MSPACIC,
// DISP>
-------
5.4.5 Main Program* Subroutines. Functions, and Block Data Required
5.4.5.1 Main Program—
IMETSCL
5.4.5.2 Subroutines—
.SOLAR
RTPHO
5.4.5.3 Functions—
None.
5.4.5.4 Block Data Files-*
None.
5.4.6 I/O and Utility Library Subroutines and Functions Required
None.
5.4.7 INCLUDE Files
None.
97
-------
SJ WIND FIELD INTERFACE (IWIND)
5.5.1 Processor Function
This processor prepares u-component and v-component wind data for the UAM model application.
IWIND reads in ROM.gridded terrain elevation, ROM hourly gridded heights of layer 1, user-supplied
diffusion break data, region top data derived from IREGNTP, and ROM hourly gridded u and v wind com-
ponents of layer 1 and 2. The methodology described in Section 4.2.4 has been applied to match the wind
data from two ROM layers into user-defined number of UAM vertical cells. An option is built in this
processor for the usdr to provide his own wind data for UAM model application. The format of an alternate
gridded wind file is described in this section. IWIND generates an unformatted file that is in a compatible
format for input to the UAM. Flow of information through IWIND interface processor is shown in
Figure 14.
5.5.2 Input/Output Components
5.5.2.1 Input Files—
'IWIND* requires four retrieved ROM MF/PF files, one user-supplied diffusion break data file, and
one derived region top data file. If the user provides his own data, the format of file should be consistent
with the format of input file WDDATA listed in this section.
MF165--This data file contains hourly gridded heights of layer 1. The formats are shown below; dots
indicate additional program statements. The input parameters are listed in Table 36.
READ(IRD2,3) IRX,IRY.NPARC
REAOCIR02.2)
READCIRD2.222) MAH£(I),IUNIT(I),IDESCCI)
REA0CIRD2,49) IDATE,ITIM,PNAHE2
READ(IRD2,51(END=999) (ZR1(I,J),I=1,IRX)
3 FORMAT(52X,I2,8X,l2,KX,l23
2 FORHATOX)
222 FORNAT<1X,A12,7X,A12,6X,A40)
49 FORHAT
-------
(INPUTS)
ROM
TERRAIN ELEVATIOi
PF119
ROM
LAYER I HEIGHTS
MFI65
i
ROM
LAYER 1 WINDS
PF115
i
ROM !
LAYER 2 WINDS ! DBDATAj RTDATAl
PFl14 | j |
J I
User-Supplied
Wind File
._(ogtjon_alJ)_
i I WIND
; INTERFACE
WDBIN
Printer
Output
(OUTPUTS)
Note: The WDBIN binary file is for use in the UAM model program.
Figure 14. Flow diagram of the IWIND interface program with input and output files.
99
-------
TABLE 36- MF165 PARAMETERS FOR IWIND
Card Parameter Parameter Data
number number name Units type Description
1
1
IRX
Integer*4
Number of columns (ROM)
2
IRY
Integer*4
Number of rows (ROM)
3
NPARC
Integer*4
Number of parameters
3+
4
NAME(I)
Character* 12
Parameter names*-
/: Index of parameter
5
IUNIT(I)
Character* 12
Unit names
6
IDESC(I)
Character*40
Parameter description
4+
7
IDATE
Integer*4
Year/3ulian date
8
ITTM
h IntegerM
Time
9
PNAME2
Character* 10
Parameter name
5+
10
2R1(U)
m Real*4
Hourly gridded height of
layer 1«
/: Index of columns
J: Index of rows
PF119--This data file contains gridded terrain elevation data. The formats are shown below; dots
indicate additional program statements. The input parameters are listed in Table 37.
READ(IR01,3) IRX,IRY,NPARB
READCIRQ 1,13) T2,IDLAT1,IMLAT1,SLAT1
READCIRD1,13) T2,IOLOW1,IML0N1,SL0H1
READCIRD1,2)
READCIRD1,222) NAMECI),IUNITCI),IDESCCI)
READCIRD1,49) IDATE,ITIH,PNAME1
READCIRD1,51)
-------
TABLE 37. PF119 PARAMETERS FOR IWIND
Card
Parameter
Parameter
Data
number
number
name
Units type
Description
1
1
IRX
Integer*4
Number of columns (ROM)
2
IRY
Integer*4
Number of rows (ROM)
3
NPARB
Integer*4
Number of parameters
2
4
T2
Character*27
Subtitle name
5
IDLAT1
Integer*4
Degrees of latitude for
the southwest corner
of the ROM domain
6
IMLAT1
Integer*4
Minutes of latitude
7
SLAT1
Real *4
Seconds of latitude
3
8
T2
Character*27
Subtitle name
9
IDLON1
Integer'4
Degrees of longitude
10
IMLON1
Integer*4
Minutes of longitude
11
SLON1
Real*4
Seconds of longitude
5
12
NAME(I)
Character* 12
Parameter names-
/: Index of parameter
13
IUNIT(I)
Character+12
Unit names
14
IDESC(I)
Character*40
Parameter description
6
15
IDATE
Integer*4
Year/Julian date
16
ITIM
h Integer*4
Time
17
PNAME1
Character* 10
Parameter name
7+
18
ELE(IJ)
m Real"4
Terrain elevation-
I: Index of columns
J: Index of rows
101
-------
DBDATA-The data contained in this file are user-supplied hourly diffusion break data associated
with date and time. The formats are shown below; dots indicate additional program statements. The
input parameters are listed in Table 38.
READ<3,17,EN0«45) LDAY(I),LTIM
-------
PF1I5—This data file contains gridded wind components for layer 1. The formats are shown below;
dots indicate additional program statements. The input parameters are listed in Table 40.
READ(IRD17,3) IRX,IRY,NPAJtD
READ(1RD17,2)
READ(iRDl7,222) NAME
-------
PF114—This data file contains gridded wind components for layer Z The formats are shown below;
dots indicate additional program statements. The input parameters are listed in Table 41.
READ(IR016,3> IRX,IRY.MPARO
READUR016,2)
READ(IR016,222) HAMECI),IUNIT(I),IDESC(I)
READ(IRD16,49,END«999) 102,IT2,PMAHEC,PNAHED
READ(IRD16,51,END«999) CUR(I,J,2),M,JRX>
READ(IRD16,51,END«999) CVR(I,J,2),'I«1,IRX)
3 F0RMAT(52X,I2,8X,I2,14X,I2)
2 FORMAT(IX)
222 FORMAT(1X,A12,7X,A12f6X,A40)
49 FORMATCI5,1X,I2,1X,2A10)
51 FORMAT(30E12.6)
TABLE 41. PF114 PARAMETERS FOR IWIND
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
IRX
Integer*4
Number of ROM columns
2
IRY
Integer*4
Number of ROM rows
3
NPARD
Integer*4
Number of parameters
3+
4
NAME(I)
Character" 12
Parameter names—
/: Index of parameter
5
IUNIT(I)
Character* 12
Unit names
6
ID ESQ I)
Character*40
Parameter description
4-f
7
ID2
Integer*4
Year/Julian date
8
IT2
Integer*4
Time
9
PNAMEC
Character* 10
Parameter name
10
PNAMED
Character* 10
Parameter name
5+
11
UR(U2)
ms-1
Real*4
Hourly gridded u-component
of layer 2-
I: Index of columns
J: Index of rows
6+
12
VR(U^)
ms"i
Real*4
Hourly gridded v-component
of layer 2
104
-------
WDDATA-This data file contains user-supplied gridded wind data. The user should set up his own
data to follow this data structure. IWtND generates a compatible file for input to the UAM applica-
tion. The formats for the user-supplied wind file are shown below; dots indicate additional program
statements. The input parameters are listed in Table 42. The formats of the output parameters are
identical to the formats of WDDATA input.
READ(33,1000) I FILE, MOTE, NSEG, MSPECS, JDATE, BEGTIH, IDATE,
& ENDTIM
REA0C33,i001> ORGX, ORGY, IZONE, UTMX, UTMY, DELTAX, DELTAY,
& NX, MY, NZ, MZLOUft, NZUPPR, HTSUR, HTLOW, HTUPP
READC33,1002) IX, IY, NXCLL, NYCLL
REAO(33!i003,END=9995 IDATE, TDATA, JDATE, TNEXT
READ(33,1004) ISEG, REF, XMX, YMX, WEST, EAST, SOUTH, CRTH
READC33,1*005) ISEG, UINDX
READ<33,1006) ((WWX(I,J,<),I»1rMX),J«1,NY)
READ(33,1005) JSEG, WINDY
READ(33,1006) <
-------
TABLE 42. WDDATA PARAMETERS (OPTIONAL FILE)
asa!=BaBsasa(ss»^nBssaass&saBSSssBSSBaBBSsssHsssssss
Card Parameter Parameter Data
number number name Units type Description
5+
21
22
23
24
25
26
27
28
29
30
Ihli.E
NOTE
NSEG
NSPECS
JDATE
HTSUR
HTLOW
HTUPP
IX
IY
NXCLL
NYCLL
IDATE
TDATA
JDATE
Integer*4 File name
Integer*4 File identifier
Integer*4 Number of segments
Integer*4 Number of species
Integer*4 Beginning year/Julian date
(yyddd-e.g. 80203)
6
BEGTIM
Real*4
Beginning lime (hhmm-Q,g. 0100,
military time)
7
IDATE
Integer*4
Ending year/Julian date
8
ENDTIM
Real*4
Ending time
9
ORGX
m
Real *4
x-coordinate (UTM units)
10
ORGY
m
Real *4
v-coordinate (UTM units)
11
IZONE
Integer*4
UTM zone
12
UTMX
m
Real *4
x-location
13
UTMY
m
Real*4
y-location
14
DELTAX
m
Real *4
Cell size inx-direction
15
DELTAY
m
Rcal*4
Cell size iny-direction
16
NX
Integer* 4
Number of cells inx-direction (UAM)
17
NY
Integer*4
Number of cells inv-direction (UAM)
18
NZ
Integer*4
Number of cells in z-direction (UAM)
19
NZLOWR
Integer'4
Number of cells between surface layer
and diffusion break
20
NZUPPR
Integer*4
Number of cells between diffusion break
Real*4
Real*4
Real*4
Integer*4
Integer*4
Integer*4
Integer*4
Imeger*4
Real *4
Integer*4
and top of region
Height of surface layer
Minimum height of cell between surface
layer and diffusion -break
Minimum height of cell between diffusion
break and top of region
x-location of segment origin with respeci to
origin of modeling region
y-location of segment origin with respeci to
origin of modeling region
Number of cells inx-direction
Number of cells in ^-direction
Beginning year/Julian date
Beginning time
Ending yeartf ulian date
(continued)
106
-------
TABLE 42. WD DATA PARAMETERS (OPTIONAL FILE CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
31
TNEXT
RealM
Ending time
6+
32
ISEG
ImegerM
Segment number
33
REF
m
RealM -
Anemometer height
34
XMX
m/h
RealM
Maximum absolute value of
%
u-component data
35
YMX
m/h
Real *4
Maximum absolute value of
v-component data
36
WEST
m/h
RealM
Average wind speed at W-boundary
37
EAST
m/h
RealM
Average wind speed at E-boundary
38
SOUTH
m/h
RealM
Average wind speed at S-boundary
39
ORTH
m/h
RealM
Average wind speed at N-boundary
; +
40
ISEG
ImegerM
Segment number
41
WINDX
RealM
Name of u-component wind
s+
42
WWX(IJJC)
m/h
RealM
u-component wind speed--
I: Index of columns
J: Index of rows
K: Index of levels
9+
43
JSEG
IntegerM
Segment number
44
WINDY
RealM
Name of v-component wind
10+
45
WWY(UK)
m/h
RealM
v-component wind speed
5*5.2,2 Control Cards-
Hie control cards define general information of the UAM model application. Twenty-five variables
included in nine control caras are used, in trie format shown below. Dots indicate additional program state-
ments. The control card variables are listed in Table 43.
READ(5,27) IUSER
READC5,101) (AOTEWCI), I®1,60)
READ(5,23) I80Y,TBEG,1EDY,TEND
READ(5,24) 0RGX,0RGY,IZ0NE1,UTHXOR,UTMYORrDElTAX,DELTAY
READ<5,25) NX, NY,NZUAM,NZL0WR,NZUPPR,HTSUR,HT10U,HTUPP
READ(5,26) IPR1,IPR2,IOBR,NSMTH
101 FORMAT(60A1)
23 FORMAT(2(I10,F10.0))
24 FORMAT(2MO.O,HO,/,2F10.0,/,2F10.0)
25 FORMAT(3t10,/,2I10,Fl0.1,2Fl0.0)
26 F0RMAT(4I10)
27 FORMAT(HO)
107
-------
TABLE 43. CONTROL CARD VARIABLES FOR IWIND
Card Parameter Variable Data
number number name Units type Description
1
1
IUSER
Integer*4
Index of user-suppplied data—
0: ROM data
> (h user-supplied data
2
2
AOTEW
Character*60
Hie name of initial conditions
3
3
IBDY
*
Integer*4
Beginning year/Julian date
4
TBEG
Real *4
Beginning time
5
IEDY
Integer*4
Ending year/Julian date
6
TEND
Real*4
Ending time
4
7
ORGX
m
Real*4
Reference origin (r-coordinate)
8
ORGY
m
Real*4
Reference origin (y-coordinate)
9
IZONE1
Integer*4
LTM zone
5
10
LTMXOR
m
Real *4
Origin of grid in .r-direction
11
UTMYOR
m
Real"4
Origin of grid iny-direction
6
12
DELTAX
m
Real *4
Cell size in .r-direcuon
13
DELTAY
m
Real*4
Cell size in y-direction
7
14
NX
Integer*4
Number of cells inx-direction (UAM)
15
NY
Integer*4
Number of cells iny-direction (UAM)
16
NZUAM
Integer*4
Number of cells inz-direction (UAM)
8
17
NZLOWR
Integer*4
Number of cells in lower layer
18
NZUPPR
Integer*4
Number of cells in upper layer
19
HTSUR
• m
Real*4
Height of surface layer
20
HTLOW
m
Real*4
Minimum height of cell in lower layer
21
HTUPP
m
Real *4
Minimum height of cell in upper layer
9
IPR1
lnteger*4
Print option for ROM parameters-- -
0 : Don't print; 1 : Print
23
IPR2
Integer*4
Print option for UAM level winds-
0: Don't print; 1: Print
24
IOBR
Integer*4
O'Brien procedure option-
0: Don't apply this option
1: Apply this option
25
NSMTH
Integer*4
Smoothing scheme option-
0: Don't apply this option
> 0: Apply option NSMTH times
108
-------
Example:
INTERFACE ROM WINDS TEST
RUN FOR
21-22 JULY 1980(3/07/90)
80203
0.
80204
24.
0.
0.
18
520000.
4460000.
8000.
8000.
31
25
5
2
3
0.0
50.0 100.0
1
0
1
2
5.5.2.3 Output Files—
'IWIND* generates one output file.
WppiN—This file contains the w-component and v-component wind data and general information for
the UAM model application. The formats are shown below; dots indicate additional program state-
ments. The output parameters are listed in Table 44.
URITECIWR28) AUFHDR,AOTEU,NSEG,NSPECSfIBDY,TBEG,IEDY,TEND
URITEUUR28) 0RGX,0RGY,IZ0NE1 ,UTMXOR,UTMYOR,DELTAX,DELTAY,
& NX,NY,NZUAM,NZLOUR,MZUPPR,HTSUR,HTLOU,HTUPP
URITEOUR28) IX,IY,NX,NY
WRITECIUR28) IYJD1,RTIM1,IYJD2,RTIM2
URITECIWR28) ISEG,REF,XMX,YMX,WBND,E9ND,S8ND,RNBND
WRITECIWR28) ISEG.AUINDX, <(UM(I,J,K), 1=1,HX), J=1,NY>
WRITE(IWR28) JSEG,AUINDY, <
-------
TABLE 44. WDBIN PARAMETERS
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
AWFHDR
Character* 10
Hie name
2
AOTEW
Character*60
File identifier
3
NSEG
Integer*4
Number of segments
4
NSPECS
Integer*4
Number of species
5
IBDY
Integer*4
Beginning year/Julian date
6
TBEG
Rea!*4
Beginning time
7
IEDY
Integer* 4
Ending year/Julian date
8
TEND
Real*4
Ending time
2
9
ORGX
m
Real *4
x-coordinate
10
ORGY
m
Real *4
v-coordinate
11
IZONE1
Integer*4
UTM zone
12
UTMXOR
m
Real *4
x-location
13
UTMYOR
m
Rear 4
y-location
14
DELTAX
m
Real"4
Cell size in x-direction
15
DELTAY
m
Real*4
Cell size iny-dircction
16
NX
Integer*4
Number of cells in x-direction
17
NY
Integer*4
Number of cells iny-direction
18
NZUAM
Integer*4
Number of cells in z-direction
19
NZLOWR
Integer*4
Number of cells between surface
layer and diffusion break
20
NZUPPR
Integer*4
Number of cells between diffusion
break and top of region
21
HTSUR
m
Reai*4
Height of surface layer
22
HTLOW
m
Real*4
Minimum height of cell between surface
layer and diffusion break
23
HTUPP
m
Real *4
Minimum height of ceil between diffusion
break and top of region
3
24
IX
IntegerM
x-location of segment origin with
respect to origin of modeling region
25
IY
Integer*4
y-locaiion of segment origin with
respect to origin of modeling region
26
NX
Integer*4
Number of ceils in x-direction
27
NY-
Integer*4
Number of cells iny-direction
4+
28
IYJD1
Integer*4
Beginning year/Julian date
for hourly data
29
RTIM1
Real*4
Beginning time
30
IYJD2
Integer*4
Ending year/Julian date
(continued)
110
-------
TABLE 44. WD BIN PARAMETERS (CONCLUDED)
« , ,i , ,in ,¦! i ¦ ¦¦ ¦¦ ¦¦
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
31
RTIM2
Real *4
Ending time
5+
32
ISEG
Integer*4
Segment number
33
REF
m
Real'4
Anemometer height
34
XMX
m/h
Real*4
Maximum absolute value of u-component
35
YMX
m/h
Real *4
data
Maximum absolute value of v-component
data
36
WBND
m/h
Real *4
Average wind at W-boundary
37
EBND
m/h
Real *4
Average wind at E-boundary
38
SBND
m/h
Real'4
Average wind at S-boundary
39
RNBND
m/h
Reala4
Average wind at N-boundary
6 +
40
ISEG
Integer*4
Segment number
41
AWINDX
Character* 10
Name of u-component wind
42
UM(UK)
m/h
Real*4
u-component wind speed—
1: Index of columns
/: Index of rows
K: Index of levels
7 +
43
JSEG
Integer*4
Segment number
44
awindy
Character'lO
Name of v-component wind
45
VM(U,K)
m/h
Real*4
v-component wind speed
111
-------
5-53 Resource Summary for ail IWIND Application
5.53.1 Memory Requirements—
FORTRAN source file:
1
file
74,240
bytes
Object file:
1
file
41,472
bytes
Executable file:
-1
file
30.20?
bytes
3
files
145,920
bytes
5.53.2 Execution Time Requirements (Representative Values for a 48-h Scenario)—
IBM 3090
Charged CPU time (hh:mm:ss): 00:10:21
Vinual address space: 10136 K
5.533 Space Requirements: Log and Print Files—
IWIND 296,448 bytes
Print Files: None
5.53.4 Space Requirements: Input and Output Files-
Table 45 shows the input file and output file space requirements.
TABLE 45. IWIND I/O FILE SPACE REQUIREMENTS
File
group
File
name
File
type
Storage
(in bytes)
Scenario
data span
Input
DBDATA
Formatted
2,048
72 hours
RTDATA
Formatted
1,536
48 hours
MF165
Formatted
197,632
72 hours
PF119
Formatted
3,584
48 hours
.
PF114
Formatted
263,168
48 hours
PF115
Formatted
263.168
48 hours
731,136
WDDATA
Formatted
(optional)
Output
WDBIN
Unformatted
269,132
48 hours
5.53.5 Space Requirements: Tape Files—
None.
112
-------
5.5,4 Run Stream Command File for an IWIND Application
//* JOB CARD
//*
//*
/•ROUTE PRINT HOLD
//STEP1 EXEC PGH«IUIND
//STEPLIB DO DSN*HHAS.UAM.INTRFACEAOAD.
//* THESE ARE THE INPUT FILE
//FT16F001 DO DSIWWAS.UAM.INTRFACE.INPUT.PF114,DISP«SHR
//FT17F001 DO DSN=MMAS.UAM.1NTRFACE.INPUT.PF115,DISP*SHR
//FT02F001 DO DSN=MMAS.UAM.INTRFACE.INPUT.MF165,DISP*SHR
//FT01F001 DD DSN=MMAS.UAM.1NTRFACE.INPUT.PF119,DISP*SHR
//FT03F001 DD DSN=MMAS.UAM.INTRFACS.INPUT.DBOATA,DISP»SHR
//FT09F001 DD DSN*MMA$.UAM.INTRFACE.INPUT.RTDATA,DI$Pa$HR
//• THIS IS THE OUTPUT FILE
//FT28F001 DD DSN*.UAM.INTRFACE.OUTPUT.UD!N,
// DI$P=(NEU,CATLG,CAUC),
// SPACE=(TRJC, (100,10)) ,UNlT*SYSOA,
// DCB=(RECFH=VB,LRECL*5000,BLKSIZE=5004)
//*
//* THIS IS THE CONTROL CARD
//FT05F001 DD *
0
INTERFACE ROM WINDS TEST RUN FOR 21-22 JULY 1980(2/12/90)
80203 0. 80204 24.
0. 0, 18
520000. 4460000.
8000. 8000.
31
2
25
3
0
5
0.0
50.0 100.0
2
/*
//* THIS IS ENO OF DATA
//
113
-------
Sii Main Program. Subroutines, Functions* and Block Dam Required
5.5.5.1 Main Program—
IWIND
5.5*5.2 Subroutines-
UTM2
GRIDROM
GRIDPRT
XMIT
SMTH
CELLHT
HORINTERP
DIVCEL
SPD
RTHETA
FMINF
MINIM
USER WD
5.5.5.3 Functions—
None.
5.5.5.4 Block Data Files—
None.
5*5.6 I/O and Utility Library Subroutines and Functions Required
None.
5.5.7 INCLUDE Files
None.
114
-------
5.6 SURFACE CHARACTERISTICS INTERFACE (ICRETER)
5.6.1 Processor Function
This processor prepares surface characteristic parameters for use in the UAM preprocessor
TERAINV ICRETER reads in ROM gridded effective surface roughness and gridded land use fractions for
each land use category. The surface roughness length (ROUGHNESS) and the vegetation factor (VEG-
FACTOR) are the two surface characteristics interpolated to the UAM grid cells using an area-weighting
method described earlier. This processor also reads in control data that define the general information of
the UAM model application and then generates a formatted packet file that is in a compatible format for
input to the UAM preprocessor 'CRETER*. Flow of information through the ICRETER interface
processor is shown in Figure 15.
5.6.2 Input/Output Components
5.6.2.1 Input Files—
'ICRETER' requires two retrieved ROM MF/PF files.
PF108-This data file contains the ROM gridded effective surface roughness length values. The
formats are shown below; dots indicate additional program statements. The input parameters are
listed in Table 46.
READ<15,12) T1A,IMX2,T1B,IMY2,71C,ISPC
R6ADC15,13) T2,IDIAT2,IMIAT2,SLAT2
REA0O5,13) T2, 1D10N2, IML0N2,SL0N2
READC15,15) T3
READC15,300,END*999) LDAT2,UIH2
READ(15,39,END*999) (A
-------
REGIONAL OXIDANT MODEL
SURFACE ROUGHNESS LENGTHS
PF108
(INPUTS)
REGIONAL OXIDANT MODEL
LAND USE INVENTORY
PF1 18
ICRETER ;
INTERFACE i
I
I CRPACK |
I (Formatted) |
CONTROL PARAMETERSI
(OUTPUT)
IAM-CRETER
PREPROCESSOR
Note: The output file (CRPACK) is in a compatible format for use in the CRETER preprocessor program.
Figure 15. Flow diagram of the ICRETER interface program with input and output files.
116
-------
TABLE 46. PF108 PARAMETERS FOR ICRETER
Card Parameter Parameter Data
number number name Units type Description
1
1
T1A
Characier*52
Subtitle name
2
IMX2
Integer*4
Number of columns
3
TIB
Character*8
Subtitle name
4
IMY1
Integer*4
Number of rows
5
TIC
Character* 14
Subtitle name
6
ISPC
Integer*4
Number of parameters
2
7
T2
Characier"27
Subtitle name
8
IDLAT2
Integer*4
Degrees of latitude for the southwest
corner of the ROM domain
9
rMLAT2
Integer*4
Minutes of latitude
10
SLAT2
Real*4
Seconds of latitude
3
11
T2
Character'27
Subtitle name
12
IDLON2
Integer*4
Degrees of longitude for the southwest
corner of the ROM domain
13
rMLON2
Integer*4
Minutes of longitude
SLON2
Reai*4
Seconds of longitude
4 +
15
T3
Charaaer*80
Header information
5
16
LDAT2
Integer*4
Year/Julian date
17
LTTM2
h Integer*4
Time
6+
18
A(IJ)
m Real*4
Surface roughness length—
/: Index of columns
J: Index of rows
117
-------
PF118—This data file contains the ROM gridded land use fractions for each land use category. The
formats are shown below; dots indicate additional program statements. The input parameters are
listed in Table 47.
READ<14,12) TU,IMXl,TlB,IMn,TlC,ISPC
READ(14,13) T2,I0LAT1,IMLAT1,SLAT1
READ<14,13) T2,IDI0N1,IML0N1,$L0N1
REA0(14,15) T3
READ(14r300(END*999) LDAT2,LTIN2
READ,I*1,INX)
12 FORMAT(A52,I2,A8,12,A14,12)
13 F0RMAT(A27,I4,1X,I2,1X,F5.2)
15 FORMAT(AflO)
300 F0RMATCI5,1X,12)
39 FORHAT(30E12.6)
TABLE 47. PF118 PARAMETERS FOR ICRETER
Card
Parameter
Parameter
Data
number
number
name
Units type
Description
1
1
TlA
Character* 52
Subtitle name
2
tMXl
Integer*4
Number of columns
3
TIB
Character'8
Subtitle name
4
IMY1
IntegerM
Number of rows
5
TIC
Character* 14
Subtitle name
6
ISPC
Integer*4
Number of parameters
2
7
T2
Character*27
Subtitle name
S
IDLAT1
IntegerM
Degrees of latitude for the southwest
corner of the ROM domain
9
tMLATl
Integer*4
Minutes of latitude
10
SLAT1
Real*4
Seconds of latitude
j
11
T2
Character'27
Subtitle name
12
IDLON1
IntegerM
Degrees of longitude for the southwest
corner of the ROM domain
13
IMLON1
IntegerM
Minutes of longitude
14
SLON1
Real *4
Seconds of longitude
4+
15
T3
Character*80
Header information
5
16
LDAT2
Integer*4
Year/Julian date
17
LTIM2
h Integer*4
"Time
6+
18
VEG(IJ,K)
RealM
Land use fractional amounts-
/: Index of columns
J: Index of rows
K: Index of land use category
-------
5.6.2*2 Control Cards—
Fifty-two variables included on control cards are used in the format shown below,
tional program statements. The control card variables are listed in Table 48.
Dots indicate addi-
READ(5,10) NX,*T,Xl,n,HX,Mr,DX,Dr
READ(5,1000) IDOLAT,IMHLAT,SECLAT, IDDLON,IMMLON,SECLON
READ(5,23) TP
READ(5,20) FID
REAItf5,21)
33 FORMATC10)
34 F0RMAT
-------
TABLE 48. CONTROL CARD VARIABLES FOR ICRETER
^————^— ¦¦
Card Card Variable Data
number number name Units type Description
1
1
NX
Integer*4
Number of cells in ^-direction (ROM)
2
NY
Integer*4
Number of cells iny-direction (ROM)
3
XI
m
Real *4
jr-location of origin
4
Y1
m
Real*4
^location of origin
5
MX
Integer*4
Number of cells in ^-direction (UAM)
6
• MY
Integer*4
Number of cells iny-direction (UAM)
7
DX
m
Real *4
UAM cell size in j:-direction
$
DY
m
Real*4
UAM cell size iny-direction
2
9
IDDLAT
Integer*4
Degrees of latitude for origin
of windowed ROM domain
10
IMMLAT
Integer*4
Minutes of latitude
11
SECLAT
Real*8
Seconds of latitude
12
IDDLON
Integer*4
Degrees oflongitude for origin
of windowed ROM domain
13
IMMLON
Integer*4
Minutes oflongitude
14
SECLON
Real*8
Seconds of longitude
3
15
TP
Character* 10
File name
4
16
FID
Character*60
File identifier
5
17
IOP(I) -
Integer*4
Control options
1=1 .
Number of species
1=2
Number of user-defined variables
1=3
Number of stations (ROM grid points)
I»4
Number of subregions
1=5
Number of parameters
6
1=6
Output file number
1=7
Input print option—
0: don't print; 1: print
1=8
Output print option-
(0 or 1, as above)
7
1=9
Print unit table (0 or 1)
1=10
Print station locations (0 or 1)
1 = 11
Print region (Oor 1)
1=12
Print methods table (0 or 1)
1=13
Print station values (0 or 1)
8
1=14
Number of vertical parameters
1=15
Number of profile heights •
1 = 16
Print vertical methods table (0 or 1)
1=17
Print vertical profiles table (0 or 1)
9
1=18
DIFFBREAK unit number
1-19
REGIONTOP unit number
1=20
TOPCONC unit number
(continued)
120
-------
TABLE 48. CONTROL CARD VARIABLES FOR ICRETER (CONTINUED)
Card Card Variable Data
number number name Units type Description
1-21
TEMPERATUR unit number
1=22
METSCALARS unit number
.
1-23
WIND unit number
10
18
ISYD
IntegerM
Beginning year/Julian date
(yyddd-e. g. 80203)
19
ISHM
IntegerM
Beginning time (hhmm-e.g. 0100,
military time)
20
IEYD
IntegerM
Ending ye$r/Julian date
21
IEHM
IntegerM
Ending time
11
22
XLOCR
m
RealM
Reference origin (x-coordinate)
23
YLOCR
m
RealM
Reference origin (y-coordinate)
24
I20NE1
IntegerM
UTM zone
12
25
XLOC
m
RealM
Origin oi grid in x-direction
26
YLOC
m
RealM
Origin of grid in y-direction
13
27
XS2
m
RealM
Cell size inx-direction
28
YS2
m
RealM
Cell size in y-direction
14
29
DC
IntegerM
Number of cells inx-direction
30
IY
Integer*4
Number of cells in y-direction
31
12
Integer*4
Number of cells in z-direction
32
I2L
Integer*4
Number of cells in lower layer
33
I2U
IntegerM
Number of cells in upper layer
34
SFCH
m
RealM
Height of surface layer
15
35
ALH
m
RealM
Minimum height of cell in lower layer
36
AUH
m
RealM
Minimum height of cell in upper layer
16
37
HT2
m
RealM
Height of station
38
INX
IntegerM
Number of cells inx-direction (ROM)
39
INY
IntegerM
Number of cells iny-direction (ROM)
17
40
ISUB
IntegerM
Number of subregions
IS
41
SUBrD
Character* 10
Subregion name
42
IROW
Integer*4
Beginning row number
43 •
ICOL
IntegerM
Beginning column number
44
ICONT
IntegerM
Cell count; if negative, equal
to rest of model region
19
45
SUBrD
Characier'lO
Subregion name
46
VAR1
Character* 10
Roughness name
47
MHD
Character* 10
Method name
48
VMIN
RealM
Minimum value
49
VMAX
RealM
Maximum value
50
NOP
IntegerM
Number of parameter, cards that followed
20+
51
PNAME
Character*10
Parameter name
52
IPVAL
IntegerM
Parameter value
(continued)
121
-------
TABLE 48. CONTROL CARD VARIABLES FOR ICRETER (CONCLUDED)
Card Card Variable Data
number number name Units type Description
21 +
53
PNAME
Character*10
Parameter name
54
PVAL
Real*4
Parameter value
22
55
SUB ID
Character* 10
Subregion name
56
VAR2
Character* 10
Vegetative factor name
57
MHD
Character'lO
Method name
58
VMIN
Real* 4
, Minimum value
59
VMAX
Real*4
Maximum value
60
NOP
Integer*4
Number of parameter cards that follow
23+
61
PNAME
Character* 10
Parameter name
62
IPVAL
Integer*4
Parameter value
24+
63
PNAME
Character* TO
Parameter name
64
PVAL
Real*4
Parameter value
Example:
16 14 520000. 4460000. 31 25 8000. 8000.
40 00 00.000 075 15 00.000
TERRAIN
SURFACE
ROUGHNESS AND
VEGETATIVE FACTORS
0
0
224
1
20
11
1
1
0
0
0
1
0
0
0
0
0
0
0
0
0
0
30001
0
80366
2400
0.
0.
18 .
520000.
4460000.
8000.
8000.
31
25
5
2
3
0.0
50.0
100.0
4.0
1
16
14
i
1
1
-1
ROUGHNESS GRID
VALUE
0.0000
20.
VEGFACTOR GRID
VALUE
0.0000
2.
-------
5.6.2*3 Output Files—
'ICRETER' generates one output file, CRPACK. This file contains a gridded surface roughness and
vegetative faaor for each UAM cell. The formats are shown below; dots indicate additional program state-
ments. The output parameters are listed in Table 49.
URITEC26,20) CNTL
WRITEC26,23) TP
WR!TE<26,20) FID
UR1TE(26,21) (I0P(I),I*1,23)
URITE(26,22) ISYD,ISHM,IEYD,IEHM
URITE(26,23) END
URITEC26i20) REGN
WRITE(26,24) XLOCR,YLOCR,IZONEl,XLOC,YLOC,XSZ,rSZ
UR1TE(26,25) IX,IY,IZ,IZL,IZU,SFCH,ALH,AUH
URITE(26,23) END
WRITE(26*,20) SU8R
WRITE(26,34) SUB ID,1 ROW,1 COL,IC3NT
URITE(26,23) END
URITE(26,20) METH
WRITE<26*35) SUBID,VAR1,MHD,VMIN,VMAX,NOP
WRITE(26,49) PNAME,IPVAL
WRITE<26,36) PNAME,PVAL
WRITE(26,35) SUBID,VAR2,MHD,VMIN,VMAX,N0P
URITEC26,49) PNAME,IPVAL
URITE(26,36) PNAME,PVAL
URITE(26,23) END
WRITE(26,23) CRIDV
WRITE(26,37) SUB ID,VAR1,1,J,VAL(I, J)
URITE<26,37) SUBID,VAR2,I,J,VAV(I,J)
WR[TEC26*23) END
URITE(26,23) ENDT
20 FORMAT(A60)
21 FORMAT<5I10,/,3I10,/,5I10,/,4I10,/,6I10>
22 FORMAT(4I10)
23 FORMATCA10)
24 FORMATC2F10.0,I10,/,2F10.0,/,2F10.0)
25 FORMAT(3110,/,2110,F10.1,2F10.0)
34 FORMAT(A10,3I10)
35 FQRMAT(3A10,2F10.1,110)
36 FORMAT(A10,F10.2)
37 FORMAT(A10,A10,2I10,F10.4)
49 FORMAT(A10,110)
123
-------
TABLE 49. CRPACK PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
1
1
CNTL
Character*60
Control pack header
2
2
TP
Character* 10
Hie name
3
3
FID
Character* 60
File identifier
4-8
4
IOP(I)
Integer*4
Control options
9
5
ISYD
Integer*4
Beginning yeartJulian date
6
ISHM
Integer*4
Beginning time
7
IEYD
Integer*4
Ending year/Julian date
8
IEHM
Integer* 4
Ending time
10
9
END
Character* 10
Control pack terminator
11
10
REGN
Character"60
Region pack header
12
11
XLOCR
m
Real *4
Reference origin (x-coordinate)
12
YLOCR
m
Real *4
Reference origin (y-coordinate)
13
IZONE1
Integer*4
UTM zone
13
14
XLOC
m
Real *4
Origin of grid in x-direction -
15
YLOC
m
Real*4
Origin of grid in y-direction
14
16
xsz
m
Real*4
Cell size in .r-direction
17
YSZ
m
Real*4
Cell size iny-direction
15
18
DC
integer*4
Number of cells in x-direction
19
IY
IntegerM
Number of cells iny-direction
20
12
lnteger*4
Number of cells in z-direction
16
21
IZL
Integer*4
Number of cells in lower layer
22
IZU
IntegerM
Number of cells in upper layer
23
SFCH
m
Real*4
Height of surface layer
24
ALH
m
Real*4
Minimum height of cell in lower layer
25
AUH '
m
Real*4
Minimum height of cell in upper layer
17
26
END
Character* 10
Control pack terminator
18
27
SUBR
Character* 60
Subregion pack header
19
28
SUBID
Character* 10
Subregion name
29
IROW
Integer*4
Beginning row number
30
ICOL
lnteger*4
Beginning column number
31
ICONT
Integer *4
Cell count
20
32 '
END
Character* 10
Control pack terminator
21
33
METH
Character*60
Method pack header
22
34
SUBID
Charaaer*10
Subregion name
35
VARl
Charaaer*10
Roughness name
(continued)
124
-------
TABLE 49. CRPACK PARAMETERS (CONCLUDED)
Card
Parameter
Parameter
Data
number
number
name
Units type
Description
36
MHD
Character* 10
Method name
37
VMIN
Real *4
Minimum value
38
VMAX
RealM
Maximum value
39
NOP
Integer*4
Number of parameter cards that follow
23+
40
PNAME
Character*10
Parameter name
41
IPVAL
Integer*4
Parameter value
24+
42
PNAME
Character* 10
Parameter name
43
PVAL
Real*4
Parameter name
25
44
SUBID
Character*10
Subregion name
45
VAR2
Character* 10
Vegetative factor name
46
MHD
Character* 10
Method name
47
VMIN
Real *4
Minimum value
48
VMAX
Real*4
Maximum value
49
NOP
Integer*4
Number of parameter cards that follow
26+'
50
PNAME
Character* 10
Parameter name
51
IPVAL
Integer*4
Parameter value
27+
52
PNAME
Character* 10
Parameter name
53
PVAL .
Real*4
Parameter value
28
54
END
Character* 10
Control pack terminator
29
55
GRIDV
Character* 10
Station reading header
30+
56
SUBID
Character* 10
Subregion name
57
VAR1
Character* 10
Roughness name
58
I
Integer*4
Index ofx-coordinate
59
J
Integer*4
Index ofy-coordinate
60
VAL(U)
m Real*4
Roughness value-
/: Index of columns
J: Index of rows
31 +
61
SUBID
Character'10
Subregion name
62
VAR2
Character* 10
Vegetative factor name
63
I
Integer*4
Index ofx-coordinate
64
J
Integer *4
Index ofy-coordinate
65
VAV(IJ)
Real*4
Vegetative factor value
I Sl J are same as above
32
66
END
Character* 10
Control pack terminator
33
67
ENDT
Character* 10
Time interval pack terminator
125
-------
5.6.3 Resource Summary for an ICRETER Application
5.63,1 Memory Requirements—
FORTRAN source file:
1
file
52,736
bytes
Object file:
1
file
23,552
bytes
Executable file:
.I
file
1W44
bytes
3
files
95,232
bytes
5,632 Execution Time Requirements (Representative Values for a 48-h Scenario)-
Charged CPU time (hh:mm:ss):
Virtual address space:
5.6.33 Space Requirements: Log and Print Files-
ICRETER
Print Files:
IBM3Q9Q
00:0(h01
9620 K
4,096 bytes
None
5.6.3.4 Space Requirements: Input and Output Files-
Table 50 shows the input file and output file space requirements.
TABLE 50. ICRETER I/O FILE SPACE REQUIREMENTS
File
group
File
name
File
type
Storage
(in bytes)
Scenario
data span
Input
PF108
PF118
Formatted
Formatted
3J584
31.232
34,816
48 hours
48 hours
Output
CRPACK
Formatted
81,920
48 hours
5.6.3.5 Space Requirements: Tape Files—
None
126
-------
5,6.4 Run Stream Command File for an Interface Application
//* J06 CARD
//*
//*
/•ROUTE PRINT HOLD
//STEP1 EXEC PGM»ICRETER
//STEPLIB 00 DSN=MMAS.UAM.1NTRFACE.LOAD
//* THESE ARE THE INPUT FILES
//FT14F001 00 DSN«MMAS.UAM.INTRFACE.INPUT.PF118,D!SP3$HR
//FT15F001 DO DSNsHMAS.UAM.INTRFACE.INPUT.PF108,DISP*SHR
//~ THIS IS THE OUTPUT FILE
//FT26F001 00 DSN«.UAM.INTRFACE.OUTPUT.CRPACK,
// DI5P=(NEW,CATLG,DELETE),
// SPACE=CTRK,<5,5)),UNIT*SYSDA,
// DC8"(RECFH-FB,tRECt«132,BLK$IZEa264Q)
//•
//• THIS IS THE CONTROL CARD
//FT05F001 00 •
16 14 520000. 4460000. 31 25 8000. 8000.
40 00 00.000 075 15 00.000
TERRAIN
SURFACE ROUGHNESS AND VEGETATIVE FACTORS
0
0
224
1
20
11
1
1
0
0
0
1
0
0
0
0
0
0
0
0
0
0
80203
0
80204
2400
0:
0.
18
520000. 4460000.
8000. 8000.
31
25
5
2
3
0.0
50.0
100.0
4.0
16
14
1
1
1
•1
ROUGHNESS GRID VALUE 0.0000 20. 0
VEGFACTOR GRID VALUE 0.0000 2. 0
/• THIS IS END OF DATA
/
127
-------
5.6.5 Main Program, Subroutines* Functions, and Block Data Required
5.6.5.1 Main Program—
ICRETER
5.6.5.2 Subroutines—
UTM2
QUAD
5.6iJ Functions—
None.
5.6.5.4 Block Data Files-
None.
5.6.6 I/O and Utility Library Subroutines and Functions Required
None.
5.6.7 INCLUDE Files
None
128
-------
5.7 CONCENTRATION INTERFACE (ICONC)
5.7.1 Processor Function
This processor prepares initial conditions, lateral boundary conditions, and top boundary conditions
of 17 species for the UAM model application. ICONC reads in ROM gridded terrain elevation, ROM
hourly gridded heights of layer 1, user-supplied diffusion break data, and ROM 3-h running average
predicted concentrations. The methodology described in Section 43 has been applied to match the concen-
trations from the three layers of ROM into a user-defined number of UAM vertical cells.
The user should be aware that: (1) this vertical interfacing methodology does not specify vertical
gradients below the diffusion break height, (2) concentrations among the lower layers are determined by
ROM layer 1 and 2 concentrations, and (3) concentrations above the diffusion break height are specified by
ROM layer 2 and layer 3 concentrations. After applying the vertical interfacing methodology, a horizontal
interpolation scheme (inverse distance-squared weighting) is incorporated into this processor for initial con-
ditions and top boundary conditions. For lateral boundary conditions, two extra ROM columns/rows of data
and one interior ROM data point are averaged. The same vertical method applied for initial conditions is
used for each UAM level and then a linear interpolation method is performed to derive UAM boundary
values at each grid cell. An option is built in this processor for the user to provide concentration files for
initial conditions, lateral conditions, and top boundary conditions. This processor then generates an unfor-
matted file that is in a format compatible for input to the UAM model application. Flow of information
through the ICONC interface processor is shown in Figure 16.
5.7.2 Input/Output Components
5.7.2.1 Input Files-
'ICONC' requires two retrieved ROM MF/PF files (MF165 and PF119), one user-supplied data file
(DBDATA), and a ROM concentration data file (ROM21). If the user provides data files, the format of the
file should be consistent with the formats of input file AQDATA, TCDATA, and BCD ATA listed below.
MF165—This data file contains hourly gridded heights of layer 1. The formats are shown below; dots
indicate additional program statements. The input parameters are listed in Table 51.
129
-------
(INPUTS)
I ROM
] TERRAIN ELEVATION
i PFII9.
ROM
LAYER 1 HEIGHTS
MF165
ROM
CONC FILE
ROM21
DBDATA
t
CONTROL
ARAMETERSl
ALTERNATE
CONCENTRATION FILES
(for initial conditions, _[ogtiona[)
lateral conditions, and
top boundary conditions)
ICONC i
INTERFACE ;
I
(OUTPUTS)
! AQB1N !
/Unformatted)!
I
BCBIN |
(Unformatted)!
TCBIN |
Unformatted)!
Note The output files are unformaued for direct use in the UAM program.
Figure 16. Flow diagram of the ICONC interface program with input and output files.
130
-------
1
2
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
READ(2,12) T1A,IMX2,T1B,IHT2,T1C,ISPC
READ(2,13) T2,10LAT2,IHIAT2,SIAT2
READC2,13) T2,IDION2,IH10N2,SL0N2
READ(2,15) T3
READ<2,300,E»«996) LDATE,LTIHE
READ(2jio,EMDa996) (Z1(I,J),1*1,IMX)
12 F0RHAT
300 FORMAT(t5,1X,I2)
30 FORMATC30E12.6)
TABLE 51. MT165 PARAMETERS FOR ICONC
Parameter
name
Units
Data
type
Description
T1A
IMX2
TIB
IMY2
TIC
Character*52
Integer*4
Character'8
Integer*4
Character*14
Subtitle name
Number of columns
Subtitle name
Number of rows
Subtitle name
ISPC
T2
IDLAT2
IMLAT2
SLAT2
IntegerM
Character*27
Integer*4
Integer*4
RealM
Number of parameters
Subtitle name
Degrees of latitude for
the southwest corner
of the ROM domain
Minutes of latitude
Seconds of latitude
T2
IDLON2
IMLON2
SLON2
T3
Ciiaracter"27
IntegerM
Integer*4
RealM
Character'80
Subtitle name
Degrees of longitude for
the southwest corner
of the ROM domain
Minutes of longitude
Seconds of longitude
Header information
LDATE
LTIME
Zl(IJ)
IntegerM
h IntegerM
m RealM
Year/Julian date
Time
Hourly gridded height of layer 1-
/: Index of columns
J: Index of rows
131
-------
PFII9—This data file contains gridded terrain elevations. The formats are shown below; dots indicate
additional program statements. The input parameters are listed in Table 52.
READ<1,12) T1A,IMX1,T1B,IMY1,T1C,ISPC
READ(1,13> T2,IDL>T1,IML>T1,$LAT1
READC1,13) T2,1DL0N1,IML0N1,SL0N1
READ<1,15) T3
READ(11300,EHD*999) L0AT2,LTIM2
READ(i*30fEMD»999) (SFCCI,J),I»1,MX)
12 FORMAT(A52,12,AS,I2,A14,12)
13 FORMAT(A27,14,1X,I2,1X,F5.2)
15 FORMAT(A60)
300 F0RMAT(15,1X,12)
30 FORMAT(30E12.6)
TABLE 52. PF119 PARAMETERS FOR ICONC
Card
Parameter
Parameter
Data
number
number
name
Units type
Description
1
1
T1A
Character*52
Subtitle name
2
IMX1
lnteger*4
Number of columns
j
TIB
Character'8
Subtitle name
4
IMY1
Integer*4
Number of rows
5
TIC
Character* 14
Subtitle name
6
ISPC
Integer*4
Number of parameters
2
7
T2
Character'27
Subtitle name
8
IDLAT1
Integer*4
Degrees of latitude for
the southwest corner
of the ROM domain
9
IMLAT1
Integer*4
Minutes of latitude
10
SLAT1
Real"4
Seconds of latitude
3
11
T2
Character'27
Subtitle name
12
IDLON1
Integer*4
Degrees of longitude for
the southwest corner
of the ROM domain
13
IMLON1
Integer'4
Minutes of longitude
14
SLON1
Real "4
Seconds of longitude
4+
15
T3
Character'80
Header information
5 +
16
LDAT2
lnteger'4
Year/Julian date
17
LTIM2
h Integer*4
Time
6+
18
SFC(IJ)
m Real "4
Hourly gridded terrain elevation-
I: Index of columns
J: Index of rows
-------
DBDATA-'The data contained in this file are user-supplied diffusion break data associated with date
and time. The formats are shown below; dots indicate additional program statements. The input
parameters are listed in Table 53.
READ(3,17jEND*998) LOAYCI>#LTIH(I>
17 F0RNAT(i7,l3,F1Q.1)
~
TABLE 53. DBDATA PARAMETERS FOR ICONC
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1 +
1
LDAY(I)
Integer*4
Year/Julian date
2
LTIM(I)
h
Integer"4
Time
3
VMH(I)
m
Real*4
Diffusion break value--
I: Index of hours
133
-------
ROM21~This data file contains modified ROM gridded 3-h running average concentration data. The
formats are shown below; dots indicate additional program statements. The input parameters are
listed in Table 54.
READ<8,400) MCOL, NROU, NlEV, NSPC
READ&SOO) <,J=1,4), M,NSPC)
READ(8,339) ((LVMAHE(J,I>,J=1,4), M.MLEV)
• READ(8*441) TEXT1, TEXT2
READ(8,1000) LDATE, ITIHE
READ(8,900X (C1
-------
AODATA-This data Hie contains user-supplied gridded concentration data for initial conditions. The
user should set up his own data to follow this data structure. ICONC generates a compatible file for
input to the UAM application. The formats for the user-supplied data file are shown below, dots
Indicate additional program statements. The input parameters are listed in Table 55. The formats of
output parameters are identical with formats of input, but it is a binary unformatted file.
READ (17,2100) IFlLE,NOTE,NSEG,NSPECS,IDATE,BEGTlM,JDATE,EN0TIM
READ*(17,2001> ORGX, ORGY, IZQNE, UTMX, UTMY, DELTAX, DELTAY,
& NX, MY, NZ, NZLOUR, NZUPPR, HTSUR, HTLOU, HTUPP
READ*(17,1002) IX, IY, NXCLL, NYCLL
READ*(17,1003) ((MSP£C(H,L),M«1, 10),L-1 ,MSPECS)
READ*(17,1005,END«999) IBGDAT.BEGTIM,IENDAT,ENDTIM
READ (17,1006) ISEG, (MSPEC(M,l),M * 1,10)
READ (17,1007) ((EMOB(I,J,1),1=1,NX),J=1,HY)
1002 F0RHAT(4I5)
1003 FORMATC10A1)
1005 FORMAT(5X,2(I1C,F10.2))
1006 FORMAT(!A,10A1)
1007 F0RMAT(9E14.7)
2001 FORMAT(F10.1,1X,F10.1,1X,I3,F10.1,1X,?10.1,1X,2F6.0,5I4,3F7.0)
2100 FORHAT(10A1,60A1,/,12,1X,12,1X,16,F6.0,16,F6.0)
135
-------
TABLE 55. AQDATA PARAMETERS FOR ICONC
Card Parameter Parameter Data
number number name Units type Description
1
IFILE
Character* 10
File name
2
NOTE
Character*60
File identifier
3
NSEG
IntegerM
Number of segments
4
NSPECS
IntegerM
Number of species
5
IDATE
IntegerM
Beginning year/Julian date
(yyddd~s.g. 80203)
6
BEGTIM
RealM
Beginning time 0100,
7
JDATE
IntegerM
military time)
Ending year/Julian date
8
ENDTIM
RealM
Ending time
9
ORGX
m
Real *4
r-coordinate (UTM units)
10
ORGY
m
RealM
y-coordinate (UTM units)
11
IZONE
IntegerM
UTM zone
12
UTMX
m
RealM
x-location
13
UTMY
m
RealM
^•location
14
DELTAX
m
RealM
Cell size inx-direction
15
DELTAY
m
Rear 4
Cell size iny-direction
16
NX
IntegerM
Number of cells inx-direction
17
NY
IntegerM
Number of cells iny-direction
18
NZ
Integer* 4
Number of cells in z-direction
19
NZLOWR
IntegerM
Number of cells between surface
layer and diffusion break
20
NZUPPER
IntegerM
Number of cells between diffusion
break and top of region
21
22
23
24
25
HTSUR
HTLOW
HTUPP
EX
IY
m Real*4
m RealM
m RealM
IntegerM
Integer*4
(continued)
Height of surface layer
Minimum height of cell between
surface layer and diffusion break
Minimum height of cell between
diffusion break and top of region
x-location of segment origin
with respect to origin of
modeling region
y-location of segment origin
with respect to origin of
modeling region
136
-------
TABLE 55. AQDATA PARAMETERS FOR ICONC (CONCLUDED)
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
26
NXCLL
Imeger*4
Number of cells injc-direction
27
NYCLL
Integer*4
Number of cells in ^-direction
5
28
MS?EC(IJ)
Character* 10
Species names
I: Index of name
/: Index of species
6+
29
BGDAT
Integer*4
Beginning year/Julian date
for hourly data
30
BEGTIM
Real*4
Beginning time
31
IENDAT
Imeger*4
Ending year/Julian date
32
ENDTIM
Real*4
Ending time
7 +
33
ISEG
Integer*4
Segment number
34
MSPEC(U)
Character* 10
Species names
8+
35
EMOB(U,L)
ppm
Real*4
Concentrations for initial condition
/: Index of columns
J: Index of rows
L : Index of species
TCDATA-This data file contains user-supplied gridded concentration data for top boundary
condition. The user should set up his own data to follow ihis data structure. ICONC generates a
compatible file for input to the UAM application. The formats for the user-supplied data file are
shown below, dots indicate additional program statements. The input parameters are listed in
Table 56. The formats of output parameters are identical with formais of input, but it is a binary
unformatted file.
READ (18,2100) IMLE,N0TE,NSEG,N5PECS,IDATE,BEC7iH,JDATE,ENDTiH
READ (18,2001) ORGX, ORGY, IZCNE, UTHX, UTMY, DELTAX, DELTAT,
$ NX, NT, HZ, NZLOUR, NZUPPR, HTSUR, HTLOW, HTUPP
READ (18,1002) IX, IT, NXCLL, NTCLL
READ (18,1003) <(HSPEC(H,L),M*1,10),L*1.NSPECS)
READ (18,1005,END«999) IBG0AT.BEGTIM,IENDAT,ENDTIM
READ (18,1006) ISEG, (MSPEC(M,L),M«1,10)
READ (18,1007) ((EHOS(I,J),1=1,NX)fJ=1,NY)
FORMAT(1OA1,60A1,/,12,1X,12,1X,I6,F6.0,16,F6.0)
FORMAT(F10.1,1X,F10.1,1X,I3,F10.1,1X,F10.1I1X,2F6.0,5K,3F7.0)
FORMAT(AI5)
FORMAT(10A1)
FORMAT(5X,2(I10,F10.2))
FORMAT(14,10A1)
FORMAT(9EH.7)
2100
2001
1002
1003
1005
1006
1007
137
-------
TABLE 56. TCDATA PARAMETERS FOR ICONC
Card Parameter Parameter Data
number number name Units type Description
1
lhli-E
CharacterMO
File name
2
NOTE
Character*60
File identifier
3
NSEG
Integer*4
Number of segments
4
NSPECS
Integer"4
Number of species
5
IDATE
Integer*4
Beginning year/JuIian date
6
BEGTIM
Real *4
Beginning time
7
JDATE
lnteger*4
Ending year/Julian date
8
ENDTIM
Real *4
Ending time
9
ORGX
m
Real*4
^'Coordinate (TJTM units)
10
ORGY
m
Real*4
^-coordinate (TJTM units)
11
IZONE
Integers
UTM zone
12
UTMX
m
Real*4
x-Iocation
13
UTMY
m
Real*4
^-location
14
DELTAX
m
Real*4
Cell size in-c-direction
15
DELTAY
m
Real*4
Cell size in>>-direction
16
NX
Integer*4
Number of cells in ^-direction
17
NY
Integer*4
Number of cells iny-direction
18
NZ
Integer*4
Number of cells in z-direction
19
NZLOWR
Integer*4
Number of cells between
surface layer and
diffusion break
20
NZ'UPPR
Inieger"4
Number of cells between
diffusion break and top
of region
21
HTSUR
m
Real*4
Height of surface layer
22
HTLOW
m
Real*4
Minimum height of cell
between surface layer and
diffusion break
23
HTUPP
m
Real *4
Minimum height of cell
24
25
IX
IY
Inieger*4
Integer*4
between diffusion break
and top of region
^-location of segment
origin with respect to
origin of modeling region
^-location of segment
origin with respect to
origin of modeling region
(continued)
138
-------
TABLE 56. TCDATA PARAMETERS FOR ICONC (CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
6+
7+
8 +
26
27
28
29
30
31
32
33
34
35
NXCLL
NYCLL
MSPECtfJ)
IBGDAT
BEGTIM
IENDAT
ENDTIM
ISEG
MSPEC
EMOB0J)
Integer*4
Integer*4
Character* 10
Integer*4
Real"4
Integer* 4
Reai*4
Integer*4
Integer*4
ppm Reai*4
Number of cells inx-direction
Number of cells iny-direction
Species names
I: Index of name
J: Index of species
Beginning year/Jiilian date
for hourly data
Beginning time
Ending year/Julian date
Ending time
Segment number
Species names
Concentrations at top of region
/: Index of columns
J: Index of rows
BCt)ATA--This data file contains user-supplied gridded concentration data for lateral boundary
condition. The.user should set up his own data to follow this data structure. ICONC generates a
compatible file for input to the UAM application. The formats for the user-supplied data file are
shown below, dots indicate additional program statements. The input parameters are listed in
Table 57. The formats of output parameters are identical with formats of input, but it is a binary
unformatted file.
READ(19,2100)tFIL£,NOTE,NSEG,MSPECS, IDATE,9EGTIM, JDATE,EN0TlM
READ(19,2001) ORGX, ORGY, IZONE, UTMX, UTMY, DELTAX, DEITAY,
& MX, MY, HZ, NZLOVR, NZUPPR, HTSUR, HTLOW, HTUPP
READ<19,1002) (IX(l),IY(I), NXClL(I), NYCLL(I),1=1,NSEG)
READ(19,1003) ((KSPEC(I,J),1=1,10),J = 1,NSPECS)
READ(19.1021) ISEGH,IEDGE,NCELLS,CCIIOCC1I,J),11 = 1,4),J=1,MCELIS)
READ<19,1007,END=999) IBGDAT,BEGTIM,JEMOAT,ENDTIM
READ<19,1009,iostat*iostat)ISEGNHf(SPNAMECI1,<),11=1,10),NEDGNO,
& (CBCONCCII,JJ),11*1,NZ),JJ=1,NCEL)
1002 FORMATC415)
1003 FORMAK10A1)
1007 FORMAT(SX,2CI10IF10.2))
1009 FORMAT(I10f10A1,I10/C9E14.7))
1021 FORMAT<3110/(9114))
2001 FORMATCF10.1,1X,F10.1,1X,I3,F10.1,1X,F10.1,1X,2F6.0,5I4,3F7.0)
2100 FORMAT(10A1,60A1,/,t2,lX,I2,1X,I6,F6.0,I6,F6.0)
139
-------
TABLE 57. BCDATA PARAMETERS FOR ICONC
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
1FILE
Character*10
Hie name
2
NOTE
Character*60
File identifier
2
3
NSEG
lnteger*4
Number of segments
4
NSPECS
lnteger*4
Number of species
5
IDATE
Integer*4
Beginning year/Julian date
6
BEGTIM
Real *4
Beginning time
7
JDATE
lnteger*4
Ending year^ulian date
8
ENDTTM
Real*4
Ending time
3
9
ORGX
m
Real *4
x-coordinate (UTM units)
10
ORGY
m
Real *4
y-coordinate (UTM units)
11
IZONE
Integer"4
UTM zone
12
UTMX
m
Real*4
jr-location
13
UTMY
m
Real*4
^location
14
DELTAX
m
Real *4
Cell size in x-direction
15
DELTAY
m
Real*4
Cell size iny-direction
16
NX
Integer*4
Number of cells in x-direction
17
NY
Integer*4
Number of cells in y-direction
IS
NZ
Integer*4
Number of cells in z-direction
19
NZLOWR
Integer*4
Number of cells between
surface layer and
diffusion break
20
NZUPPR
Integer*4
Number of cells between
diffusion break and top of
region
21
HTSUR
m
RealN
Height of surface layer
22
HTLOW
m
Real*4
Minimum height of cell
between surface layer and
diffusion break
23
HTUPP
m
Real*4
Minimum height of cell
between diffusion break
and top of region
4
24
IX
Integer*4
x-location of segment
origin with respect to
origin of modeling region
25
rY
Integer*4
y-Iocation of segment
origin with res pea to
origin of modeling region
(continued)
140
-------
TABLE 57. BCDATA PARAMETERS FOR ICONC (CONCLUDED)
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
26
NXCLL
IntegerM
Number of cells injr-direction
27
NYCLL
IntegerM
Number of cells iny-direction
5
28
MSPECflJ)
Character* 10
Species names
/: Index of name
% /: Index of species
6+
29
ISEGM
IntegerM
Segment number
30
IEDGE
IntegerM
Edge number
31
NCELLS
IntegerM
Number of cells on edge
32
ILOC(U)
Integer*4
Index of each grid
I: Index indicator
J: Grid indicator
7 -r
33
IBGDAT
IntegerM
Beginning year/Julian date
for hourly data
34
BEGTIM
RealM
Beginning time
35
rENDAT
Integer*4
Ending year/Julian date
36
ENDTIM
Real*4
Ending time
8+
37
ISEGNM
Integer*4
Segment number
38
SPNAME(U)
Character* 10
Species names
39
NEDGNO
Integer*4
Edge number
40
BCONC(U)
ppm
Real'4
Concentrations for lateral boundary
*
/: Index of columns
J: Index of rows
141
-------
5.7XI Control Cards-
Thirty-six variables included in sixteen control cards are used, in the format shown below. Dots
indicate additional program statements. The control card variables are listed in Table 58.
READ(5r29) IUSE*
READ&ii) TP
READ(5,23) TP2
READ<5,23) TP3
READ<5,20) FID
READ(5,20) FID2
READ(5,20) FID3
REAO(5,22) ISYD,ISHM,IEYD,IEHM
REAO<5,24) XLOCR,YLOCR,I ZONE1,X10C,Y10C,XSZ,YSZ
READ(5,25) IX,|Y,IZ,IZL,IZU,SFCH,AIH,AUH
READCS^M) IMX.INY
READC5,26) ROUS,3D INC,RDMAX
READ(S,33) ISUB
20 FORMAT(A60)
22 FORMATC4I10)
23 FORMATCA10)
24 FORMAT(2F10.0,I10,/,2F10.0,/,2F10.0)
25 FORMAT(3I10,/,2I10,F10.1,2F1Q.O)
26 FORMATC3F10.1)
28 FORMATC2I10)
29 FORMAT(HO)
33 FORMAT(HO)
TABLE 58. CONTROL CARD VARIABLES FOR ICONC
Card
Variable
Data
number
name
Units type
Description
1
IUSER
Integer*4
Index of user-suppplied data-
0: ROM data
>0: user-supplied data
2
TP
Character* 10
File name of initial conditions
(i.e., AIRQUALITY)
3
TP2
Ciiaracter'10
File name of top boundary conditions
(i.e., TOPCONC)
4
TP3
Character* 10
File name of lateral boundary conditions
(i.e., BOUNDARY)
5
FID
Character*60
File identifier of initial conditions
(continued)
142
-------
TABLE 58. CONTROL CARD VARIABLES FOR ICONC (CONCLUDED)
Card
Variable
Data
number
name
Units
type
Description
6
FTD2
Character*60
File identifier of top boundary conditions
7
FTD3
Character*60
File identifier of lateral boundary conditions
8
ISYD
IntegerM
Beginning year/Julian date
ISHM
IntegerM
Beginning time
EEYD
IntegerM
Ending yearfJulian date
IEHM
IntegerM
Ending time
9
XLOCR
m
Real*4
Reference origin (x-coordinate)
YLOCR
m
RealM
Reference origin (y-coordinate)
IZONE1
IntegerM
UTM zone
10
XLOC
m
RealM
Origin of grid in x-direction
YLOC
m
RealM
Origin of grid in ^-direction
11
xsz
m
Real *4
Cell size in jc-direction
YSZ
m
RealM
Cell size iny-direction
12
IX
Integer*4
Number of cells in jc-direction (UAM)
rY
IntegerM
Number of cells in^-direction (CAM)
IZ
IntegerM
Number of cells in z-direction (UAM)
13
IZL
Integer*4
Number of cells in lower layer
IZU
Inieger*4
Number of cells in upper layer
SFCH
m
Real*4
Height of surface layer
ALH
m
RealM
Minimum height of cell in lower layer
A UH
m
RealM
Minimum height of cell in upper layer
14
INX
IntegerM
Number of cells in jc-direction (ROM)
1NY
Integer*4
Number of cells in v-direction (ROM)
15
RDUS
km
RealM
Initial radius for influence of
spatial interpolation
RD1NC
km
RealM
Radius increment for influence
of spatial interpolation
RDMAX
km
RealM
Maximum radius for influence
of spatial interpolation
16
ISUB
IntegerM
Number of subregions
143
-------
Example:
0
A1R0UALITY
TOPCONC
BOUNDARY
AIRQUALITY FILE CREATED AT 03/08/90 (UPPER VERTICAL VARIATION)
TOP CONCENTRATION FILE CREATED AT 03/08/90 (UPPER VERTICAL VARIATION}
BOUNOARY FILE CREATED AT 03/08/90 (UPPER VERTICVAL VARIATION)
80203
0
80204
2400
0.
0.
18
520000.
4460000.
8000.
8000.
31
25
5
2
3
0.0
50.0
16
14
20.0
10.0
50.0
1
5.7.23 Output Files—
'ICONC generates three output files.
AOBIN—This file contains the initial conditions of chemical concentrations and general information
for the UAM model application. The formats are shown below; dots indicate additional program
statements. The output parameters are listed in Table 59.
URITE(20) TP,nD,ISUB,NSPC,ISY0,SHM,ISY0fSHM+1.
URITEC20) XLOCR,YLOCR,I ZONE 1,XL0C,YL0C,XSZ,YSZ,
& IX,IY,IZ,IZL,IZU,SFCH,ALH,AUH
URITEC20) 1X1,IY1,IX,IY
WRITEC20) <(SPC(J,I),J=1,10),1=1,NSPC)
URITEC20) IYD1,HM1,iTD1,HM2
WRITE(20)"ISUB,(SPC(J,K),J=1,10),((C3(I,J,M,K),1-1,IX),J»1,IY)
144
-------
TABLE 59. AQBIN PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
1
TP
Character* 10
File name (i.e., AIRQUALITY)
2
FID
Character*60
File identifier
3
ISUB
Integer*4
Number of segments (Le., 1)
4
NSPC
Integer*4
Number of chemical species
5
ISYD
Integer*4
Beginning year/Julian date
6
SHM
Real *4
Beginning time
7
XLOCR
m
Real *4
x-coordinate (UTM units)
8
YLOCR
m
Real "4
^-coordinate (UTM units)
9
IZONE1
Integer*4
UTM zone
10
XLOC
m
Real*4
x-location
11
YLOC
m
Real*4
y-location
12 -
XSZ
m
Real'4
Cell size in j:-direction
13
YSZ
m
RealM
Cell size iny-direction
14
DC
Integer*4
Number of cells in x-direction
15
IY
Integer*4
Number of cells iny-direction
16
IZ
Integer*4
Number of cells in z-direction
17
IZL
Integer*4
Number of cells between surface layer
and diffusion break
18
IZU
IntegerM
Number of cells between diffusion break
and top of region
19
SFCH
m
Real*4
Height of surface layer
20
ALH
m
RealM
Minimum height of cell between surface layer
and diffusion break
21
AUH
m
Real*4
Minimum height of cell between diffusion
break and top of region
22
1X1
Integer*4
x-location of segment origin with respect
to origin of modeling region
23
IY1
IntegerM
y-location of segment origin with respect
to origin of modeling region
24
DC
IntegerM
Number of cells in ^-direction
25
IY
IntegerM
Number of cells iny-direction
(continued)
145
-------
TABLE 59. AQBIN PARAMETERS (CONCLUDED)
Card Parameter Parameter Data
number number name Units type
Description
26 SPC(J,I)
6+
27
28
29
30
IYD1
HM1
HM2
ISUB
Character* 10
Integer*4
Real *4
Real*4
Integer*4
Species names—
/: Index of name
J: Index of species
Beginning year/Julian date for hourly data
Beginning time
Ending time
Segment number
31 SPC(JJC) Character*10 Species names
32 C3(IJ,M,K) ppm Real*4 Concentrations for initial conditions-
I: Index of columns
J: Index of rows
M: Index of levels
K: Index of species
TCBIN-This file contains the top boundary conditions of concentrations and general information for
the UAM model application. The formats are shown below; dots indicate additional program state-
ments. The output parameters are listed in Table 60.
URITE(21) TP2,FtD2,1SU8,NSPC,ISYD,SHM,IEYD,EHM
URITE(21) XLOCR,YIOCR,IZONE1,XLOC,YIOC,XSZ,YSZ,
$ IX,IY,12,IZl,IZU,SFCH,ALH,AUH
URITE(21) 1X1,IY1,IX,IY
WRITEC21) ((SPC
WRITEC21) IYD1,HM1,IY02,HM2
WRITEC21) ISUB,,J-1,1Q>,aC3(I,J,IZ-M,K),I=1,IX),J = 1,ir)
-------
TABLE 60. TCBIN PARAMETERS
Card Parameter Parameter Data
number number name Units type Description
5+
1
TP2
Character*10
File name
2
FID2
Character*60
File identifier
3
ISUB
Integer*4
Number of segments
4
NSPC
IntegerM
Number of species
5
ISYD
IntegerM
Beginning year/Julian date
6
SHM
Real*4
•
Beginning time
7
IEYD
IntegerM
Ending year/Julian date
8
EHM
RealM
Ending time
9
XLOCR
m
RealM
x-coordinate (UTM units)
10
YLOCR
m
Real *4
y-coordinate (UTM units)
11
IZONE1
IntegerM
UTM zone
12
XLOC
m
Real"4
x-location
13
YLOC
m
Real *4
y-locaiion
14
xsz
m
Real'4
Cell size inx-direction
15
YSZ
m
RealM
Cell size in y-direction
16
IX
Integer*4
Number of cells in .r-direction
17
ry
Integer*4
Number of cells in y-direction
18
IZ
Integer*4
Number of cells inz-direaion
19
IZL
Integer*4
Number of cells between surface layer
and diffusion break
20
IZU
Integer*4
Number of cells between diffusion break
21
22
23
24
25
26
27
28
29
30
SFCH
ALH
AUH
1X1
rYl
IX
IY
SPC(J,I)
IYD1
HM1
m Real*4
m Real*4
m Real *4
Integer*4
Integer*4
Integer*4
Integer*4
Character* 10
lnteger*4
Real*4
and top of region
Height of surface layer
Minimum height of cell between surface
layer and diffusion break
Minimum height of cell between diffusion
break and lop of region
x-Iocation of segment origin with respect
to origin of modeling region
^-location of segment origin with respect
to origin of modeling region
Number of cells in x-direction
Number of cells in v-direction
Species names-
I: Index of name
J: Index of species
Beginning year/Julian date for hourly data
Beginning time
(continued)
147
-------
TABLE 60. TCBIN PARAMETERS (CONCLUDED)
SBsatsssaasassaassaassssssssssssacsBss&ss
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
31
IYD2
Integer*4
Ending year/Julian date
32
HM2
Real *4
Ending time
6+
33
ISUB
ppm
Real *4
Concentrations for initial conditions-
34
SPC(J,K)
I: Index of columns
35
C3(IJMJ£)
/: Index of rows
M: Index of levels #
K: Index of species
BCBTN~This file contains the lateral boundary conditions of concentrations and general information
for the UAM model application. The formats are shown below; dots indicate additional program
statements. The output parameters are listed in Table 61.
WRITEC22) TP3,FID3,ISU8,NSPC,!SYD,SHM,IEYD,EKM
URITE(22) XL0CR,Y10CR,I ZONE1,X10C,YL0C,XSZ,YSZ,
I IX,IY,IZ,IZL,I2U,SFCH,ALH,AUH
URITE(22) 1X1,IY1#IX,IT
WRITEC22) ((SPC(J,I),J*1,10),1*1,NSPC)
WRITE(22)"ISUB,M,ICELL(H),((ILOC(I,J),1*1,4),J*1,ICELL(M))
WRITEC22) IYD1,HM1,IYD1,HM2
WRITEC22) ISUB,(SPC(J,O,J=1,10),Ml,((BDC(I,J),I«1,IZ),J*1,ir)
WRITE(22)*ISUB, (SPC( J,)(), J = 1,10),M2,((BDCCI,J),1=1,12),J*1,IT)
WRITE(22)* ISUB,(SPC(J,K),J=1,10),M3,((80CCI,J),1=1,12),J=1,IY)
WRITE(22)'lSUB,(SPCCJ,K),J-M0),M4,«BDC(t,J),I = 1 ,IZ)tJ = 1,IY)
14$
-------
TABLE 61. BCBIN PARAMETERS
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
TP3
Character* 10
File name
2
FID3
Character*60
File identifier
3
ISUB
Integer*4
Number of segments
4
NSPC
Integer*4
Namber of species
5
ISYD
Integer*4
Beginning year/Julian date
6
SHM
Real*4
Beginning time
7
IEYD
Integer*4
Ending year^Juiian date
8
EHM
Real *4
Ending time
2
9
XLOCR
m
Real *4
x-coordinate (UTM units)
10
YLOCR
m
Real*4
y-coordinate (UTM units)
11
IZONE1
Integer'4
UTM zone
12
XLOC
m
Real *4
x-location
13
YLOC
m
Real *4
v-location
14
xsz
m
Real'4
Cell size in x-direction
15
YSZ
m
Real'4
Cell size in y-direction
16
IX
IntegerM
Number of cells in x-direction
17
rY
Integer*4 *
Number of cells iny-direction
18
IZ
Integer*4
Number of cells in z-direction
19
IZL
Integer*4
Number of cells between diffusion break
and surface layer
20
IZU
Integer*4
Number of cells between diffusion break
and top of region
21
SFCH
m
Real "4
Height of surface layer
22
ALH
m
Real*4
Minimum height of cell between surface
23
AUH
m
Real"4
layer and diffusion break
3
24
1X1
Integer*4
Minimum height of cell between diffusion
25
nri
IntegerM
break and top of region
x-Iocation of segment origin with respect
to origin of modeling region
,v-location of segment origin with respect
to origin of modeling region
26
IX
Integer*4
Number of cells in x-direction
27
ry
Integer*4
Number of cells in ^-direction
4
28
SPC(J,K)
Character* 10
Species names-
J: Index of name
K: Index of species
5+
29
ISUB
Integer*4
Segment number
30
M
Integer*4
Edge number
(continued)
149
-------
TABLE 61. BCBIN PARAMETERS (CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
6+
7+
31
32
33
34
35
36
37
38
39
ICELL(M)
ILOCftJ)
IYD1
HM1
IYD2
HM2
ISUB
SPC(J,K)
Ml
40 BDC(IJ)
Integer*4
Integer*4
Integer*4
Real *4
Integer*4
RealM
IntegerM
Character* 10
IntegerM
RealM
Number of ceils on edge
Index of location—
1: Index of indicator
J: Index of grids
Beginning year/Julian date for hourly data
Beginning time
Ending year/Julian date
Ending time
Segment number
Species name
Index of boundarv-
1 : West
2: East
3: South
4: North
Concentrations for lateral
boundary condition-
I: Index of columns
J: Index of rows
8+
41
ISUB
Integer*4
Segment number
42
SPC(J JC)
Character* 10
Species name
43
M2
Integer*4
Edge number
44
BDC(IJ)
Real*4
Concentration for lateral boundary
9+
45
ISUB
Integer*4
Segment number
46
SPC(J,K)
Character* 10
Species name
47
M3
integer*4
Edge number
48
BDC(IJ)
Real*4
Concentration for lateral boundary
10+
49
ISUB
Integer*4
Segment number
50
SPC(JJC)
Charaaer*10
Species name
51
M4
Integer*4
Edge number
52
BDC(IJ)
Real*4
Concentration for lateral boundary
150
-------
5.73 Resource Summary for a (JAM Application
5.7.3.1 Memory Requirements—
FORTRAN source file:
1
file
61,440
bytes
Object file:
1
file
37^88
bytes
Executable file:
-I
file
3U3?
bytes
3
files
130,560
bytes
5,132 Execution Time Requirements (Representative Values for a 48*h Scenario)*
Charged CPU time (hh:mmss):
Virtual address space:
5.7.3-3 Space Requirements: Log and Print Files—
ICONC
Print Files:
IBM 3090
00:01:31
10136 K
19,456 bytes
None
5.7.3.4 Space Requirements: Input and Output Files-
Table 62 shows the input file and output file space requirements.
TABLE 62. ICONC I/O FILE SPACE REQUIREMENTS
File
group
File
name
File
type
Storage
(in bytes)
Scenario
data span
Input
DBDATA
MF165
PF119
ROM21
Formatted
Formatted
Formatted
Formatted
2,048
197,632
3,584
7,067.648
7,270,912
72 hours
72 hours
48 hours
48 hours
Output
AQBIN
BCBIN
TCBIN
Unformatted
Unformatted
Unformatted
269312
2,001,920
2,574336
4,845,568
48 hours
48 hours
48 hours
5.7Space Requirements: Tape Files—
None
151
-------
5,7,4 Run Stream Command File for an Interface Application
//• joe card
//*
//*
/~ROUTE PRINT HOLD
//STEP1 EXEC PGK-ICONC
//STEPLIB DO 0SN«MHAS.UAH.INTRFACE.LOAD
//* THESE ARE THE INPUT FILES
//FT01F001 DD DSN«MMAS.UAM.INTRFACE.INPUT.PF119,DISP»SHR
//FT02F001 DD DSN«#MAS.UAM.INTRFACE.INPUT.MF165,DISP«SHR
//FT03F001 DD DSW-MMAS.UAM.INTRFACE.INPUT.DBDATA,D1SP«SHR
//FT08F001 DD DSN-HMAS.UAM.INTRFACE.INPUT.R0N21,DISP«SHR
//* THESE ARE THE OUTPUT FILES
//FT20F001 DD DSN«.UAK.INTRFACE.OUTPUT.AQ8IN,
// DISP»(NEW,CATLG,DELETE),
// SPACE«CTRIC,{5,5)),UNITsSYSDA,
// DC8*(RECFMsVB,LRECL*5000,BLKSl2Es5004)
//FT21F001 DD DSN*.UAM.INTRFACE.OUTPUT.TCBIN,
// DISP»(NEW,CATLG,DELETE),
// SPACE«(TRIC, (5,5)),UMIT=SYSOA,
// DC8*(RECFM=VB,LRECL=5000,BLKSIZEs5004)
//FT22F001 DD DSN».UAM.INTRFACE.OUTPUT.BC8IN,
// DISP=(NEU,CATLG,DELETE),
// $PACE=(TRK,(5,5)),UNlTsSYSDA,
// DC8=(RECFM=VB,LRECL=5000,BLICSIZE«5004)
//* THIS IS THE CONTROL CARD
//FT05F001 DD *
0
AIRQUALITY
TOPCONC
BOUNDARY
AIRQUALITY FILE CREATED AT 02/07/90 (UPPER VERTICAL VARIATION)
TOP CONCENTRATION FILE CREATED AT 02/07/90 (UPPER VERTICAL VARIATION)
BOUNDARY FILE CREATED AT 02/07/90 (UPPER VERTICVAL VARIATION)
80203 0 80204 2400
0. 0. 18
520000. 4460000.
8000. 8000.
31
25
5
2
3
0.0
16
14
20.0
10.0
50.0
1
/*
//« THIS IS
END OF DATA
//
152
-------
5.7-5 Mnin Program, Subroutines, Functions, and Block Data Required
5.7.5.1 Main Program—
ICONC
5.7 J £ Subroutines—
USERAQ
USERBC
USERTC.
TSFM
STGNV
UTM2
5.7.5.3 Functions—
None.
5.7.5.4 Block Data Files—
None.
5.7.6 I/O and Utility Library Subroutines and Functions Required
None.
5.7.7 INCLUDE Files
None.
153
-------
SA BIOGENIC EMISSIONS INTERFACE (IBIOG)
5.8.1 Processor Function
This preprocessor prepares a combined biogenic-anthropogenic emissions data file for use in the
UAM. IBIOG reads in ROM biogenics emissions from PF144 for six species: olefin, paraffin, isoprene,
aldyhyde, NO, and NO2. IBIOG computes UTM coordinates for ROM grid points, applies an area-
weighting method to derive biogenic emissions values for the UAM grid cells, and merges the UAM gridded
biogenics with the area source emissions file. The combined file is an unformatted output file for use in the
modeL An optional ASCII version of just the biogenic emissions data can also be written. Flow of informa-
tion through the IBIOG interface processor is shown in Figure 17.
5.8.2 Input/Output Components
5.8.2.1 Input FHes—
'IBIOG' requires one retrieved ROM MF/PF file and one user-supplied data file.
PF144--PF144 iS the windowed ROM biogenics file. The formats are shown below; dots indicate
additional program statements. The input parameters are listed in Table 63.
READ(19,1103, IOSTAT = IQST) PFCOL,PFROU,NPARH
1103 FORMAT <52X,I2,8X,12,14X,IS)
READ(19,1104, IOSTAT = IOST) DUMMY
1104 FORMAT(A4)
READ<19,1101,IOSTAT«IOST)LDATE, LTIHE, (PNAMECI>,1=1,NPARM)
1101 FORMAT(I5,1X,I2,1X,7(1X,A10))
READ(19,1102,IOSTAT"IOST)
& (A(I,H,l),,PFCOl)
1102 F0RMAT(30(E12.6))
154
-------
(INPUTS)
REGIONAL QXIDANT MODEL !
BIOGENIC EMISSIONS |
PF144
USER'S
ANTHROPOGENIC AREA EMISSIONS
(From UAM - EPS)
iCONTROL PARAMETERS
IBIOG
' INTERFACE
(Optional)
I
BIOASC !
(Formatted)
EMBIN j
(Unformatted)
(OUTPUTS)
Note: The area anthropogenic emissions (unformatted) file is generated by the user from the UAM Emissions
Processor System (EPS). The EMBIN (unformatted) File contains both biogenic and anthropogenic area
emissions for use in the UAM. An optional formatted biogenics file (BIOASC) may also be generated if
desired.
Figure 17. Flow diagram of the rBIOG interface program with input and output files.
155
-------
TABLE 63. PF144 PARAMETERS FOR IBIOG
Card
number
Parameter
number
Parameter
name
Units
Data
type
Description
4+
9
10
PFCOL
PFROW
NPARM
DUMMY
LDATE ,
LTTME
PNAME(I)
A(U,1)
A(IJ.2)
A (U3)
moles/h
moies/h
Integer*2
lnteger*2
Integer*2
Character*4
Integer*5
Integer*2
Character* 10
moles/h Real*4
Real1
Real'
Number of columns
Number of rows
Number of parameters
Dummy variable
Date
Time
Parameter name-
I: Index of parameters
Olefin biogenic emission rates-
I: Index of columns
J: Index of rows
Paraffin biogenic emission rates
Isoprene biogenic emission rates
11 A(IJ,4) moles/h Real*4
12 A(IJt5) moles/h Real*4
13 A(IJ,6) moles/h Real *4
Aldehydes biogenic emission rates
Nitric oxide biogenic emission rates
Nitrogen dioxide biogenic emission rates
EMISSION-EMISSION is the binary area source emissions file generated by exercising the UAM
Emissions Processor System (Causley et al. 1990). The formats are shown below; dots indicate addi-
tional program statements. The input parameters are listed in Table 64.
READ (18,IQSTAT = I0ST) IN FIIE,MOTE,NSEG,NSPECS,IDATE,
& BEGTIM,JDATE,ENDTIM
READ (18) ORGX,ORGY,IZ0NE,UTMX,U7MY,DEITAX,DELTAY,NXA,NYA,
& HZ rNZL0UR,NZUPPR,HTSUR,HT10U,HTUPP
READ (18) IX,IY,NXCLL,NYCIL
READ (18) ((HSPECd,J),I«1,10),J»1,NSPECS)
READ (18) IBGDAT, 8EGTIM, IENDAT, ENDTIH
READ (18) ISEG,
-------
TABLE 64. ANTHROPOGENIC EMISSION PARAMETERS FOR IBIOG
Card Parameter Parameter Data
number number name Units type Description
1
1
INFILE
Integer*4
Hie name
2
NOTE
Integer*4
FQe identifier
3
NSEG
Integer*4
Number of segments
4
NSPECS
Integer*4
Number of species
5
IDATE
Integer*4
Beginning date (yyddd—c.g. 80203)
6
BEGTIM
Real *4
Beginning time (hhmm~e.g. 0100,
military time)
7
JDATE
Integer*4
Ending date
8
ENDTIM
Reai*4
Ending time
2
9
ORGX
m
Real*4
Reference origin (z-coordinate)
10
ORGY
m
Real*4
Reference origin (y-coordinate)
11
IZONE
IntegerM
UTM zone
12
UTMX
Real*4
Origin of grid in j:-direction
13
UTMY
Rea]*4
Origin of grid iny-direction
14
DELTAX
m
Real*4
Cell size in x-direction
15
DELTAY
m
Real*4
Cell size iny-direction
16
NXA
Integer*4
Number of cells in x-direction
17
NYA
Integer*4
Number of cells iny-direction
18
N2
Integer*4
Number of cells in z-direction
19
NZLOWR
Integer*4
Number of cells in lower layer
20
NZUPPR
Integer*4
Number of cells in upper layer
21
HTSUR
m
Real*4
Height of surface layer
22
HTLOW
m
Real*4
Minimum height of cell in lower layer
23
HTUPP
m
Real*4
Minimum height of cell in upper layer
3
24
IX
Integer*4
Dummy variable
25
rY
Integer*4
Dummy variable
26
NXCLL
Integer*4
Number of cells in x-direction
27
NYCLL
Integer*4
Number of cells iny-direction
4
28
MSPEC(IJ)
Integer*4
Species in the header-
I: Index of species name
J: Index of species
5
29
IBGDAT
Integer*4
Beginning date
30
BEGTIM
Real*4
Beginning hour
(continued)
157
-------
TABLE 64. ANTHROPOGENIC EMISSION PARAMETERS FOR IBIOG
(CONCLUDED)
Card Parameter Parameter Data
number number name Units type Description
31
32
33
34
35
IENDAT
ENDTIM
ISEG
MSPEC(M,L)
Integer*4
Real*4
Integer*4
Integer*4
EMOB(ltJ,L) moles/h Real*4
Ending date
Ending hour
Segment number
Species in hour—
M - 1,10
L = l.NSPECS
Emission rate-
I: Index of columns
J: Index of rows
L: Index of species
5.8.2.2 Control Cards-
Control card information is changed each time in the command file. The formats are shown below;
dots indicate additional program statements. The control card variables are listed in Table 65.
READ(5,10) BDATE,BKOUR,EDATE,SHOUR
10 FORMAT(1X,I5,1X,I3,1X,I5,1X,I3)
REAOC5,15) MX,NT,X1,Y1,HX,MY,DX,DY
15 FORMAT(1X,2I5,2F10.0,2I7,2F10.0)
READ<5,20)IPT,IDDLAT, IMML.AT,SECLAT, IDDLON,
& XIMHL0N,SEC10N
20 FORMAT(A4,1X,12/IX,J2,1X,F6.3,1X,I3,1X,12,1X,F6.3)
READ(5,25)ASCURITE
25 FORKAT(1X,12)
158
-------
TABLE 65. CONTROL CARD VARIABLES FOR IBIOG
Card
Variable
Data
number
name
Units
type
Description
1
BDATE
Integer*5
Beginning date
BHOUR
Integer*3
Beginning hour
EDATE
Integer*5
Ending date
EHOUR
Integer*3
Ending hour
2
NX
Integer*4
Number of cells inx-direction (ROM)
NY
lnteger*4
Number of cells in ^-direction (ROM)
-XI
m
Real*4
Origin of grid inx-direction
Yl
m
Real *4
Origin of grid in ^-direction
MX
Integer*4
Number of cells in x-direction (UAM)
MY
Integer'4
Number of cells in y-direction (UAM)
DX
m
Real*4
UAM cell size inx-direction
DY
m
Real*4
UAM cell size in ^-direction
3
IPT
Character^
Station ID
IDDLAT
Integer*4
Degrees of latitude
IMMLAT
integer*4
Minutes of latitude
SECLAT
Real*4
Seconds of latitude
IDDLON
Integer*4
Degrees of longitude
XIMMLON
Integer"4
Minutes of longitude
SECLON
Real*4
Seconds of longitude
4
ASCWRITE
Integer*2
Index of output
Example:
80202 0 80204 24
16 14 520000. 4460000. 31 25 8000. 8000.
PLL1 40 00 00.000 075 15 00.00
1
159
-------
5.8.2-3 Output Files—
'IBIOG' generates two output files-one binary file and one optional ASCII file.
EMBIN-EMBIN contains the merged biogenic and area hourly gridded emissions, including: olefin,
paraffin, isoprene, aldehydes, NO, and NO2. The formats are shown below; dots indicate additional
program statements. The output parameters are listed in Table 66.
WRITE (29) INFILE, MOTE, NSEG, NSPECS+1, BDATE, RBHOUR, EDATE,
I REHOUR
WRITE (29) OftGX, ORGY, IZONE, UTMX, UTMY, DELTAX, DELTA?,
& NXA, HXA, HZ, NZLOUR, NZUPPR, HTSUR, HTLOU, HTUPP
WRITE (29) IX, IY, NXCLL, NYCLL
WRITE (29) ((HSP€C(H,L),H«1,10),L«1,NSPECS+1)
WRITE (29) IBGDAT,BEGT1M,IENOAT,ENDTIM
WRITE (29) ISEGr(MSPEC(M,K),H*1,10),((AREABIO(I,J,)C),
& I « 1, MX), J » 1f MY)
TABLE 66. EMBIN PARAMETERS
Card
Parameter
Parameter
Data
number
number
name
Units
type
Description
1
1
INFILE
Integer*4
File name (i.e., EMISSIONS)
2
NOTE
Integer*4
File identifier
3
NSEG
Integer*4
Number of segments
4
NSPECS+1
Integer*4
Number of species
5
BDATE
Integer*4
Beginning date tyyddd-t.g. 80203)
6
RBHOUR
Real*4
Beginning time (hhmm-c.g. 0100,
military time)
7
EDATE
Integer*4
Ending date
8
REHOUR
Real *4
Ending time
2
9
ORGX
m
Real*4
Reference origin (r-coordinate)
10
ORGY
m
Real* 4
Reference origin (y-coordinate)
(continued)
160
-------
TABLE 66. EMBIN PARAMETERS (CONCLUDED)
Card Parameter
number number
Parameter
name
Units
Data
type
Description
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
IZONE
UTMX
UTMY
DELTAX
DELTAY
NXA
NYA
NZ
N2LOWR
NZUPPR
HTSUR
HTLOW
HTUPP
rx
IY
NXCLL
NYCLL
MSPECS(U)
rBGDAT
BEGTIM
IntegerM UTM zone
m RealM Origin of grid in x-direction
m Real*4 Origin of grid in y-direction
m Real *4 Cell size in x-direction
m RealM Cell size iny-direction
IntegerM Number of cells in x-direction
IntegerM Number of cells in y-direction
IntegerM Number of cells in z-direction
IntegerM Number of cells in lower layer
IntegerM Number of cells in upper layer
m RealM Height of surface layer
m Real *4 Minimum height of cell in lower layer
m RealM Minimum height of cell in upper layer
IntegerM Dummy variable
IntegerM Dummy variable
IntegerM Number of cells in x-direction
IntegerM Number of cells iny-direction
IntegerM Species in header
IntegerM Beginning date
RealM Beginning time
31 IENDAT
32 ENDTIM
33 ISEG
34 MSPECS(U)
35 AREABIO(UK)
• Integer*4
RealM
Integer*4
IntegerM
moles/h Real*4
Ending date
Ending time
Segment number
Species each hour
Emission rate-
I: Index of columns
J : Index of rows
K: Index of species
161
-------
BIOASC-BIOASC contains onJy the biogenic emissions data. This file is optional. The formats are
shown below; dots indicate additional program statements. The output parameters are listed in
Table 67.
WRITE (9,60) INFllE,KOTE,NSEG,BIOSPEC$2,BOATE,RBHOUR,EDATE,
& REHOUR
WRITE (9,70) ORGX, ORGY, IZOME, UTMX, UTMY, DELTAX, DELTAY,
$ NXA, NYA, NZ, N2L0UR, NZUPPR, HTSUR, HTLOW, HTIPP
WRITE (9,80) IX, IT, NXCll, MYCLl
WRITE (9,90) ((8I0MAME2(N,O,H»1,10),l»1,8I0$PEC$2)
60 FORNAT(10A1,60A1,/,I2,1X,I2,1X,I6,F6.0,I6,F6.0)
70 FORMAT(F10.1,1X,F10.1,1X,I3,F10.1,1X,F10.1,1X,2F6.0,5I4,3F7.0)
80 FORHAT(4l5)
90 FORMAT(10A1)
WRITE*(9,140) I8GDAT,BEGTIM,IENDAT,ENDTIM
HO FORMAT(5X,2(I10,F10.2))
WRITE (9,380)*ISEG, (BIOMAME2(M,IO,M"1,10)
380 F0RMAT(I4,10A1)
WRITE (9,390) ((BIOGENO,J,K),1=1,MX5,J«1,MY)
390 F0RMAT(9E14.7)
TABLE 67. BIOASC PARAMETERS
Card Parameter Parameter Data
number number name1 Units type Description
1
1
INFILE
Integer*4
File name
2
NOTE
Integer*4
File identifier
3
NSEG
Integer*4
Number of segments
4
BIOSPECS2
Integer*4
Number of species
5
BDATE
Integer*4
Beginning date (yyddd-c.g. 80203)
6
RBHOUR
Real'4
Beginning time (hhmm-c.g. 0100,
military time)
7
EDATE
Integer*4
Ending date
8
REHOUR
Real*4
Ending time
2
9
ORGX
m Real *4
Reference origin (x-coordinate)
10*
ORGY
m Real*4
Reference origin (y-coordinate)
(continued)
162
-------
TABLE 67. BIOASC PARAMETERS (CONCLUDED)
Card Parameter Parameter Data
number number name1 Units type Description
11
IZONE
Integer*4
UTM zone
12
UTMX
m
Real*4
Origin of grid in jc-direction
13
UTMY
m
Real *4
Origin of grid in ^-direction
14
DELTAX
m
Real* 4
Cell size in je-direction
15
DELTAY
m
Real *4
Cell size iny-
-------
5.83 Resource Summary for a UAM Application
5.83.1 Memory Requirements—
FORTRAN source file:
1
file
59392
bytes
Object file:
1
file
33,280
bytes
Executable file:
-I
file
27.13$
bytes
3
files
119,808
bytes
5.8*3 J Execution Time Requirements (Representative Values for a 48-h Scenario)'
Charged CPU time (hh:mm:ss):
Virtual address space:
5.8.3.3 Space Requirements: Log and Print Files—
IBIOG
Print Files:
IBM 3090
00:00-33
3371 K
43,008 bytes
None
5.8.3.4 Space Requirements: Input and Output Files-
Table 68 shows the input file and output file space requirements
TABLE 68. IBIOG I/O FILE SPACE REQUIREMENTS
File name
File type
Storage
(in bytes)
Scenario
data span
INPUT:
PF144
EMISSION
Formatted
Unformatted
918,016
1.515.008
2,433,024
48 hours
48 hours
BINARY OUTPUT:
EMBIN
Unformatted
1,666,560
48 hours
ASCn OUTPUT:
BIOASC
Formatted
5,832,192
48 hours
5.8J.5 Space Requirements: Tape Files—
None
164
-------
5.S.4 Run Stream Command File for an IBIOG Application
//*J08 CARD
//*
//•
/•ROUTE PRINT HOLD
//STEP1 EXEC PGtt*IBIOG
//STEPLIB DD DSN=MMAS.UAM.lNTRFACE.LOAD
//* THESE ARE THE INPUT FILE
//FT18F001 DO DSN«HHAS.UAM.INTRFACE.INPUT.EMISSION,DISP=SHR
//FT19F001 DO DSNs»4AS.UAM.INTRFACE.INPUT.PF144,DISP=SHR
//* THESE ARE THE OUTPUT FILES
//FT09F001 OD DSN«.UAM.INTRFACE.OUTPUT.BIOASC,
// DISP*(NEV,CATLG,DELETE),
// SPACE«(TRK,(100,10)),UNIT«SYSOA,
// DC8«.UAM.INTRFACE.OUTPUT.EMBIN,
// DISP*(NEU,CATLG,DELETE),
// SPACE«< TRK,C100,10)),UNIT»SYSOA,
// DC8*(RECFM«VB,LRECLs5000,BLKSIZE=5004)
//• THIS IS THE CONTROL CARD
//FT05F001 DO *
80203 0 80204 24
16 14 520000. 4460000.
PLL1 40 00 00.000 075 15 00.000
1
r
It* THIS IS END OF DATA
//
31 25 3000. 8000.
NOTE: =USERID AND ACCOUNT
5.8.5 Main Program* Subroutines, Functions* and Block Data Required
5.8.5.1 Main Program—
IBIOG
5.8.5.2 Subroutines—
BIOUTM6
QUAD
UTMCON
5.8.5.3 Functions—
None.
5.8.5.4 Block Data Files—
None.
5.8.6 I/O and Utility Library Subroutines and Functions Required
None.
5.8.7 INCLUDE Files
None.
165
-------
REFERENCES
Ames, J., T.C. Myers, L.E. Reid, D.C Whimey, S.H. Golding, S.R. Hayes and SD. Reynolds. 1985. SAI
Airshed Model Operations Manuals, Volume I: User's Manual. EPA/600/8-85/007A (NTIS PB85-
191567). U.S. Environmental Protection Agency, Research Triangle Park, NC
Causley, M.C., J.L. Fieber, M. Jimenez, and L. Gardner. 1990. User's Guide for the Urban Airshed Model,
Volume IV: User's Guide for the Emissions Processor System. Prepared for the Office of .Air
Quality Planning and Standards, U.S. Environmental Protection Agency, Research Triangle Park,
NC.
Clarke, J.F. and J.K.S. Ching. 1983. Aircraft observations of regional transport of ozone in the northeastern
United States. Atmos. Environ., 17: 1703-1712.
Computer Sciences Corporation (CSC). 1990. Gridded Model Information Support System (GMISS):
UAM subsystem design. EPA Office of Air Quality Planning and Standards, Source Receptor
Analysis Branch, Research Triangle Park, NC.
Demerjian, K.L., K.L. Schere, and J.T. Peterson. 1980. Theoretical estimates of actinic (spherically inte-
grated) flux and photoiytic rate constants of atmospheric species in the lower troposphere. In:
Advances in Environmental Science and Technolo$>y Vol. 10, J. Pitts and R. Metcalf, Eds., Wiley Publ.,
New York. pp. 369-459.
Doll, D.C., T.E Pierce, and N.C. Possiel. 1989. Regional Ozone Modeling in the Northeastern U.S.:
- Selection of Meteorological Episodes. Sixth Joint Conference on Applications of Air Pollution
Meteorology, Jan. 30-Feb. 3,1989, Anaheim, CA. Preprints, American Meteorological Society,
Boston, MA., pg. 40-43.
Douglas, S.G., R.C. Kessler, and EL Carr. 1990. User's Guide for the Urban Airshed Model, Volume III:
User's Manual for the Diagnostic Wind Model (Version 1.1). Prepared for the Office of Air Quality
Planning and Standards, U.S. Environmental Protection Agency, Research Triangle Park, NC.
166
-------
EPA. 1986. Guideline on Air Quality Models (Revised). EPA-450/2-78-027R. U.S. Environmental Protec-
tion Agency. Research Triangle Park, NC
Geiy, M.W., G.Z. Whitten, and J.P. Killus. 1988. Development and Testing of the CBM-IV for Urban and
Regional Modeling. EPA/600/3-88-012, U.S. Environmental Protection Agency, Research Triangle
Park, NC
Geiy, M.W., G.Z. Whitten, J.P. Killus, and M.C Dodge. 1989. A Photochemical Kinetics Mechanism for
Urban and Regional Scale Computer Modeling. J. Geophys. Res., 94(10): 12925-12956.
Godowitch, J.M., J.K.S. Ching, and J.F. Garke. 1987. Spatial Variation of the Evolution and Structure of
the Urban Boundary Layer. Bound. Layer-MeteoroL, 38: 249-272.
Lamb, B., H. Westberg, and G. Allwine. 1985. Biogenic Hydrocarbon Emissions from Deciduous and
Coniferous Trees in the United States. J. Geophys. Res.y 90:2380-2390.
Lamb, B., A Guenther, D. Gay, and H. Westberg. 1987. A National Inventory of Biogenic Hydrocarbon
Emissions. Atmos. Environ., 21:1695-1705.
Morris, R.E., T.C. Myers, and J.L. Haney. 1990a. User's Guide for the Urban Airshed Model, Volume I:
User's Manual for the UAM(CB-P/). Prepared for the Office of Air Quality Planning and Standards,
U.S. Environmental Protection Agency, Research Triangle Park, NC.
Morris, R.E., T.C. Mvers, E.L. Carr, M.C. Causley, S.G. Douglas, and J.L Haney. 1990b. User's Guide for
the Urban Airshed Model, Volume II: User's Manual for the UAM(CB-P/) Modeling System (Pre-
processors). Prepared for the Office of Air Quality Planning and Standards, U.S. Environmental Pro-
tection Agency, Research Triangle Park, NC.
Pierce, T.E.. K.L. Schere, D.C. Doll, and W.E. Heilman. 1990a. Evaluation of the Regional Oxidant Model
(Version 2.1) Using Ambient and Diagnostic Simulations. Atmospheric Research and Exposure
Assessment Laboratory, U.S. Environmental Protection Agency.
Pierce, T.E., B.K. Lamb, and AR. Van Meter. 1990b. Development of a Biogenics Emissions Inventory-
System for Regional Scale Air Pollution Models. Paper No. 90-94.3 presented at the 83rd Air and
Waste Management Association Annual Meeting, June 24-29,1990, Pittsburgh, PA.
Rao, S.T. 1987. Application of the Urban Airshed Model to the New York Metropolitan Area. EPA/450/4-
87-011. NT1S PB87-201422. U.S. Environmental Protection Agency, Research Triangle Park, NC.
167
-------
Schere, K.L. and RA Wayland. 1989. EPA Regional Oxidant Model (ROM2.0): Evaluation on 1980
NEROS Data Bases. EPA/600/3-89/057. NTIS PB89-200828. Atmospheric Research and Exposure
Assessment Laboratory, U.S. Environmental Protection Agency, Research Triangle Park, NC
Wolf, G.T. and PJ. Lioy. 1980. Development of an Ozone River Associated with Synoptic Scale Episodes
In the Eastern United States. Environ. Set and Techn14(10): 1257-1260.
Young, J.O., M Aissa, T.L. Boehm, CJ. Coats, J.R. Eichinger, S J. Roselle, AR Van Meter, R A- Wayland,
and T.R Pierce. 1989. Development of the Regional Oxidant Model Version 2.1. EPA/600/3-89/044.
(NTIS PB89-194252) Atmospheric Research and Exposure Assessment Laboratory, U.S. Environ-
mental Protection Agency, Research Triangle Park, NC
168
-------
APPENDIX A
DESCRIPTION OF THE EXAMPLE TEST CASE
A test case of input/output files has been assembled for a user to exercise when initially implementing
the interface program package with the UAM system. The test case is the two-day of July 21-22, 1980. The
example UAM domain is the New York metropolitan area. A set of 13 ROM retrieved data files contains all
parameters needed for exercising the interfaces for the 48-h period.
Appendix B presents the contents of the control data files input to the interface programs and
Appendix C contains a partial listing of each ROM data file used in this test case. In addition, Appendix D
presents the printer output produced by each interface and the initial contents of each output data file
generated by the interface programs for comparison with the results obtained by the user.
Appendix E describes the utility programs to convert a binary data file to ASCII formatted file
(BINASC) or convert a formatted file into a binary file (ASCBIN). The BINASC program can be applied to
convert the binary files generated by the concentration, wind, and biogenic emissions interface programs to
compare user results with those provided in Appendix D for these particular files. Finally, Appendix F gives
a complete listing of the contents of a magnetic tape that contains the input/output files for this test case
along with the interface program source codes.
A-l
-------
APPENDIX C
SAMPLE OF THE ROM INPUT DATA FILES FOR THE TEST CASE
C-l
-------
1. KOM2I
000207900014441800080202000000230000360000000000NEROSXXX 75.125
AL02C0 MTHLETH F0RMH202HN02HN03IS0PN0 N02 03 OLE PAN PAR TOL
1 2 3
R0H2.1
* ROM 2 J RUN ON IBH 3090
* BC FROM FORMULATION FV1072A, NEWIC IS FV1073A
* ROH/CONC FILE FV1073A PROOUCED
AVERAGED EXTRACTIONS. NUMBER OF TIHESTEPS AVERAGED: 2
OUTPUTS DERIVED FROM ADDING OTHER SPECIES OR LEVELS
THIS CONCENTRATION FILE HAS BEEN CONVERTED TO ASCII FORMAT
40.083 71.37b 42.250-0.25000 0.16667 16 14 3 17
L
** GRIMES 06-06-89
* BENCHMARK TEST OF R0M2 1 CODE UlTH BUM = 10
* BMATRIX FROM FORMULATION FV1073A
THIS IS AN EXTRACTION
TIME SHIFTED TO BEGIN TIME
80203000000 3600 1
0.222582E-020.494094E-020.3602786-020.307475E-020.261280e-020.l570296-020.9932466-030.750864E-030.646239E-030.485347E-030.310672E-030.18318lE-03
0.1546066-030.1618866-030.1418716-030.1198716-03
0.1844076-020.1845526-020.1353426-020.8954986-030.7507776-030.7326536-030.6764706-030.5772256-030.4637126-030.377411E-030.3167556-030.266311E-03
0.2101836-030.1533496-030.1165316-030.9364446-04
0.2823396-030.846461E-030.698181E-030.463226E-030.4085976-030.4151896-030.4504176-030.4357456 030.4049346-030.3827216-030.3472206-030.3117796-03
0.281868E-030.244681E- 030.18605 7E-030.140866E-03
0.2654626+000.4384596*000.4025456+000.2379426*000.2006386*000.1594416+000.1420706+000.1332346+000.1354176+000.145185E+000.1461086+000.136010E+00
0.1244366+000.1206106+000.117051E+000.112343E+00
0.233856E+000.254743E+000.209077E+000.1518906+000.1310456+000.1340326+000.1355956+000.1332826*000.130210E+000.128313E+000.1300206+000.1285626+00
0.122178E+000.116109E+000.114374E+000.112938E+00
0.103085E+000.160922E+000.145131E+000.123292E+000.1175596+000.118001E+000.121020E+000.120844E+G00.120389E+000.122380E+000.123676E+000.126417E+00
0.127160E+000.123475E+000.118569E+000.115093E+00
0.1000006-150.1000006-150.1000006*150.1000006-150.1000006-150.999999E-160.999999E-160.999999E-160.9999996-160.9999991-160.9999996-160.9999996-16
0.9999996-160.9999996-160.999999E-160.999999E-16
0.1000006-150.1000006-150.1000006-150.100000E-150.1000006*150.1000006-150.1000006-150.100000E-150.100000E-150.100000E-150.1000006-150.1000006-15
0.10OOOOE *150.100000E-150.100000E *150.1OOOOOE-15
0.10OOOOE-150.1OOOOOE-150.1000006-150.1OOOOOE-150.1000006-150.1OOOOOE-150.1OOOOOE-150.1000006-150.1000006-150.100000E-150.100000E-150.1000006-15
0.1OOOOOE-150.1OOOOOE-150.1OOOOOE-150.1OOOOOE-15
0.164683E-020.2971216-020.4348476-020.1929896-020.162764E-020.132350E-020.1110826-020.8983356-u30.6046106-030.231768E-030.436297E-040.9543726-06
0.866671E-050.980672E-050.475024E-050.337862E-05
0.1344866-020.1281796*020.9356936-030.7218456-030.7232156-030.7265646-030.638008E-030.497715E OiO.3263736-030.2028396-030.1054466-030.522091E-04
0.3474776-040.251123E-040.150660E-040.737244E-05
0.1275276-030.5418496-030.4322966-030.2382666-030.2239446-030.2544716-030.3178766-030.3086136 GJO.270448E-030.234669E-030.1789406-030.1211396-03
0.896614E-040.850496E-040.655945E-040.443403E-04
0.4876396-020.7080156-020.5574496-020.4954606-020.4572166-020.3028306-020.2281106-020.1827476-020.1239826-020.6535946-030.3552426-030.2135856-03
0.1861966-030.166796E-030.163260E-030.163888E-03
0.466940E-020.471600E-020.402145E-020.325404E-020.5035426-020.302365E-020.278461E-020.238644E 020.1991146-020.1717126-020.1503876-020.1342886-02
0.1208616-020.108191E-020.100116E-020.973808E-03
0.195197E-020.2938226-020.2690006-020.2269396-020.2183206-020.2199866-020.2263206-020.218668E-U20.2058056-020.1913836-020.1772156-020.159917E-02
0.1457476-020.1344526-020.1192226-020.109265602
0.5 326796-020.5734 506 - 020.605 7916 - 020.5 5 34456-020.504 2 74E-020.399067E-020.321740E-020.2663746-u20.1885756*020.1179286-020.8170176-030.4882346-03
-------
8. IBIOG
80203 0 80203 24
16 14 520000. 4460000. 31 25 8000. 8000.
Ml 40 00 00.000 075 15 00.000
1
8EGIN
DATE» 80203 8EGIN
TIME*
0.00
ENO
DATE* 80203
ENO
TIME*
1.00
BEGIN
DATE* 80203 BEGIN
TIME*
1.00
Em
DATE* 80203
ENO
TIME*
2.00
BEGIN
DATE" 80203
BEGIN
TIME*
2.00
ENO
DATE* 80203
ENO
TIME-
3.00
BEGIN
DATE* 80203
BEGIN
TIME*
3.00
ENO
DATE* 80203
ENO
TIME*
4.00
BEGIN
DATEa 80203
BEGIN
TIME*
4.00
ENO
DATE* 80203
ENO
TIME*
5.00
BEGIN
DATE* 80203
BEGIN
TIME*
5.00 ENO
DATE* 80203
ENO
TIME*
6.00
BEGIN
DATE* 80203
BEGIN
TIME*
6.00
ENO
DATE* 80203
ENO
TIME*
7.00
BEGIN
DATE* 80203
BEGIN
TIME*
7.00
ENO
DATE* 80203
ENO
TIME*
8.00
BEGIN
DATE* 80203
BEGIN
TIME*
8.00
ENO
DATE* 80203
ENO
TIME*
9.00
BEGIN
DATE* 80203
8EGIN
TIME*
9.00
ENO
DATE* 80203
ENO
TIME-
10.00
BEGIN
DATE* 80203
BEGIN
TIME*
10.00
ENO
OATE* 80203
ENO
TIME*
11.00
BEGIN
DATE* 80203
8EGIN
TIME*
11.00
ENO
DATE* 80203
ENO
TIME*
12.00
BEGIN
DATE* 80203
8EGIN
TIME*
12.00
ENO
DATE* 80203
ENO
TIME*
13.00
BEGIN
DATE* 80203
8EGIN
TIME-
13.00
ENO
DATE* 80203
ENO
TIME*
14.00
BEGIN
DATE* 80203
BEGIN
TIME*
14.00
ENO
DATE* 80203
ENO
TIME*
15.00
BEGIN
DATE* 30203
3EGIN
TIME*
15.30
END
OATE* 30203
END
TIME*
<6.CO
3EGIN
DATE' 30203
SEGIN
"!ME*
:6.00
END
OATE* 80203
END
TIME*
17.00
8EGIN
DATE* 30203
BEGIN
TIME*
*7.00
ENO
DATE* 80203
ENO
TIME*
18.00
BEGIN
DATE* 30203
BEGIN
TIME*
18.00
ENO
OATE* 80203
ENO
TIME*
19.00
BEGIN
DATE* 80203
8EGIN
TIME*
19.00
ENO
DATE* 80203
ENO
TIME*
20.00
BEGIN
DATE* 80203
BEGIN
TIME*
20.00
ENO
OATE* 80203
END
TIME*
21.00
BEGIN
DATE* 80203
BEGIN
TIME*
21.00
ENO
OATE* 80203
ENO
TIME*
22.00
8EGIN
DATE* 80203
8EGIN
TIME*
22.00
ENO
DATE* 80203
ENO
TIME*
23.00
BEGIN
DATE* 30203
BEGIN
TIME*
23.00
ENO
OATE* 80203
ENO
TIME*
24.00
. •.
.Normal Completion
Of UAMBIO
D-ll
-------
IL OUTPUT DATA FILES
L DBPACK
CONTROL
DIFFBREAK
HEIGHT Of THE MIXED LAYER
0.
4460000.
8000.
25
224
1
0
0
0
80205
18
5
0.0
30203
0 0
23 1
0 0
0 0
0 0
80203 0
END
REGION
Q.
520000.
8000.
31
2
ENO
TIME INTERVAL
30203
:U8REGI0N
A 1 1
END
METHOO
A DIFFBREAK constant
ENO
CONSTANTS
A DIFFBREAK 285.0
END
ENDTIME
TIME INTERVAL
30203 100 80203
CONSTANTS
A DIFFBREAK 285.0
END
ENDTIHE
TIME INTERVAL
30203 200 80203
CONSTANTS
A DIFFBREAK 285.0
ENO
ENDTIME
TIME INTERVAL
30203 300 80203
CONSTANTS
A DIFFBREAK 285.0
END
ENDTIME
TIME INTERVAL
30203 400 80203
CONSTANTS
A DIFFBREAK 285.0
END
ENDTIME
TIME INTERVAL
30203 500 80203
CONSTANTS
A DIFFBREAK 300'. 0
ENO
ENDTIME
TIME INTERVAL
80203 600 30203
CONSTANTS
A OIFFBREAK 360.0
END
ENDTIME
1
0
0
0
100
200
300
400
500
600
700
10
0
0
50. 100.
100
-1
0.0 3000.0
D-12
-------
2. RTPACK
CONTROL
REGIONTOP
REGIONTOP FILE CREATED
AT 01/05/90
0 0
224
1
10
28 1
1
0 0
0
0
0
0 0
0
0
0 0
0
0
0
80203 0
80204
2400
EMO
REGION
0. 0.
18
520000. 4460000.
8000. 8000.
31 22
5
2 3
0.0
50.
100.
END
TIME INTERVAL
80203 0
80203
100
5UBREGICN
A 1
1
-1
END
METHOD
A REGIONTOP CONSTANT
200.0
3000.0
END
CONSTANTS
A REGIONTOP
1213.6
END
ENOTIHE
TIME INTERVAL
80203 100
80203
200
CONSTANTS
A REGIONTOP
1213.6
END
ENOTIHE
TIME INTERVAL
80203 200
80203
300
CONSTANTS
A
END
ENDTIME
TIME INTERVAL
80203
CONSTANTS
A REGIONTOP
END
EHDTIME
REGIONTOP
300
1213.6
80203
1213.6
400
D-13
-------
X TPPACK
CONTROL
TENPERATUR
TEMPERATURE PILE CREATED AT 02/07/90
0
28
0
0
0
60203
END
REGION
0.
520000.
8000.
31
2
ENO
UNITS
TEMPERATURDEGK
ENO
STATIONS
0.
4460000.
8000.
25
3
224
0
0
0
0
80204
18
5
0.0
1.0
1
0
0
0
2400
50.
0.0
001
489343.
4436804.
4.0
002
489369.
4455302.
4.0
003
489395.
4473801.
4.0
304
489421.
4492301.
4.0
005
489447.
'4510801.
4.0
006
489474.
4529302.
4.0
007
489500.
4547803.
4.0
008
489527.
4566305.
4.0
009
489554.
4584808.
4.0
010
489580.
4603311.
4.0
10
0
0
100.
0.0
TINE INTERVAL
80203 0 80203
SUBREGION
A 1 1
ENO
NETHOO
A TENPERATURSTATINTERP
EXTENT 10.00
INITRAOIUS 2.00
RAO IUSI NCR 1.00
NAXRAOIUS 5.00
END
STATION REAOtNGS
001
7ENPERATUR
303.55
002
TENPERATUR
303.70
003
TEMPERATUR
303.83
004
TENPERATUR
303.91
005
TEMPERATUR
303.92
006
TEMPERATUR
303.35
007
TENPERATUR
303.70
008
TENPERATUR
303.47
009
TENPERATUR
303.17
010
TENPERATUR
302.82
100
-1
270.0 330.0 4
ENOTINE*'
TINE INTERVAL
80203 100
80203
STATION READINGS
001
TEMPERATUR
303.13
002
TENPERATUR
303.33
003
TENPERATUR
303.46
004
TENPERATUR
303.55
005
TENPERATUR
303.58
006
TENPERATUR
303.53
007
TENPERATUR
303.39
008
TENPERATUR
303.17
009
TEMPERATUR
302.86
010
TENPERATUR
302.48
D-14
-------
fc CRPACK
CONTROL
TERRAIN
SURFACE ROUGHNESS ANO
VEGETATIVE
FACTORS
0 0
224
1
20
11 1
1
0 0
0
1
0
0 0
0
0
0 0
0
0
0
80203 0
80204
2400
ENO •
REGION
0. 0.
18
520000. 4460000.
8000. 8000.
31 25
5
2 3
0.0
50.
100.
ENO
SUBREGION
A 1
1
•1
5.ND
METHOO
A ROUGHNESS GRIO
VALUE
0.0
20.0
A VEGFACTOR GRIO
VALUE
0.0
2.0
ENO
GRI0VALUE5
A ROUGHNESS
1
0.4413
A ROUGHNESS
2
I
0.4625
A ROUGHNESS
3
0.5002
A ROUGHNESS
4
0.6595
A ROUGHNESS
5
0.6582
A ROUGHNESS
6
0.6260
A ROUGHNESS
7
0.6006
A ROUGHNESS
a
0.6023
A ROUGHNESS
9
0.1564
A ROUGHNESS
10
0.0954
A VEGFACTOR
1
0.4234
A VEGFACTOR
2
0.4227
A VEGFACTOR
3
0.4185
A VEGFACTOR
4
0.4014
A VEGFACTOR
5
0.4021
A VEGFACTOR
6
0.3455
A VEGFACTOR
7
0.2970
A VEGFACTOR
8
0.2983
A VEGFACTOR
9
0.0773
A VEGFACTOR
10
0.0470
D-17
-------
APPENDIX E
CONVERSION OF OUTPUT FILES (BINASC AND ASCBIN)
1. BINASC (BINARY TO ASCII: UNFORMATTED DATA TO FORMATTED DATA)
I. PROCESSOR FUNCTION
This processor converts binary UAM model input and output data to ASCII formatted data for inves-
tigation of pre- and post-model application. BINASC reads in UAM preprocessors' or ROM-UAM
interface processors' output data files and control cards. The average concentration output data file
from UAM model can aiso be converted to ASCII formatted data.
II. I/O COMPONENTS
A. Input Files
The input files consist of any input data file for UAM model application and UAM average con-
centration output data file.
3. Control Cards
Two control cards are used in this processor. The control card variables are listed in Table E-l.
READ (*,100) £ FILE
READ (*,100) IPATH3
100 FORMAT (A)
E-l
-------
TABLE E-l. BINASC CONTROL CARD VARIABLES
Variable Data
name Units type Description
IFTLE Character* 10 File name
IPATH3 Character* 80 Index of unit for log file
Example: The run stream is listed below.
//« J08 CARD
//*
/r
/•ROUTE PRINT HOLD
//STEP1 EXEC PGH=8INASC
//STEPLIB DO DSN«.LOAD,D[SP=SHR
//* THIS IS "HE INPUT ?!LE
//FT07F001 3D DSN=, 0 T SP-SHR
//* THIS IS THE OUTPUT FILES
//FT09F001 DO DSN=
-------
2. ASCBIN (ASOI TO BINARY: FORMATTED DATA TO UNFORMATTED DATA)
L PROCESSOR FUNCTION
This processor converts ASCII formatted UAM model input and output data back to binary data for
model application. ASCBIN reads in UAM preprocessors' or ROM-UAM interface processors'
output data files (ASCII formatted) and control cards. The average concentration output data file
from the UAM model can also be converted to binary unformatted data for postprocessor application.
n. I/O COMPONENTS
A- Input Files
The input files consist of any input data file for the UAM model application and the UAM
average concentration output data file.
B. Control Cards
Two control cards are used in the processor. The control card variables are listed in Table E-2.
READ (*,100) IFILE
READ 100) IPATH3
100 FORMAT (A)
TABLE E-2. ASCBIN CONTROL CARD VARIABLES
Variable
Data
name
Units type
Description
IFILE
Character* 10
File name
IPATH3
Character*80
Index of unit for log file
E-3
-------
Example: The run stream is listed below.
tr Jos CARD
//*
//*
/•ROUTE PRINT HOLD
//STEP1 EXEC PGM-ASCBIN
//STEPLIB DO 0SN«« LOAD,DISP«SHR
//* THIS IS THE INPUT FILE
//FT07F001 DO 0SN»,DISP«SHR
//* THIS IS THE OUTPUT FILES
//FT09F001 DO 0SN»
-------
APPENDIX F
MAGNETIC TAPE LISTING OF PROGRAMS AND DATA FILES
FOR THE EXAMPLE TEST CASE
FT le
nunber
file name
1. 131OG
INPUTS-
EMISSION
PF144
IBIOG
INTERFACE
RUNSTREAH- RUN8G
LOGFILE- BIOLOG
OUTPUT FILES
6IOASC
EMSFN (ASCII VERSION)
2. ICONC
8
9
10
11
12
13
14
15
16
17
INPUTS- PF119
MF165
OBOATA
R0M21
INTERFACE- ICONC
RUNSTREAH- RUNCC
LOGFILE- CONCLOG
OUTPUT FILES- AOBIN (ASCII VERSION)
TC8IN (ASCII VERSION)
BCBIN (ASCII .VERSION)
3. ICRETER
18
19
20
21
22
23
INPUTS- PF118
PF108
INTERFACE- ICRETER
RUNSTREAH- RUNCR
LOGFILE- TERLOG
OUTPUT FILES- CRPAC1C
4. IHETSCL
24
25
26
27
28
29
30
31
32
33
INPUTS- OBOATA
MF174
RTDATA
PF102
PF119
PF117
INTERFACE- IHETSCL
RUNSTREAH- RUNHS
LOGFILE- NETLOG
OUTPUT FILES- MSPACX
F-I