Table of Contents - The Atmospheric Studies Group at TRC

Table of Contents - The Atmospheric Studies Group at TRC Table of Contents - The Atmospheric Studies Group at TRC

29.01.2015 Views

Section 8: CALMET Model Files 8.2 Geophysical Data File (GEO.DAT) The GEO.DAT data file contains the geophysical data inputs required by the CALMET model. These inputs include land use type, elevation, surface parameters (surface roughness, length, albedo, Bowen ratio, soil heat flux parameter, and vegetation leaf area index) and anthropogenic heat flux. The land use and elevation data are entered as gridded fields. The surface parameters and anthropogenic heat flux can be entered either as gridded fields or computed from the land use data at each grid point. Default values relating each of these parameters to land use are provided in the model. A sample GEO.DAT file is shown in Table 8-5. The first line of the file identifies the file and its format to CALMET. The second line provides the number of comment lines (character strings of up to 80 characters in length) to follow that are read but are not interpreted by the model. One of these may typically contain a title to identify the data set. The next block of 5 to 6 lines contains map projection, datum, and grid information such as the number of grid cells, grid spacing, and reference coordinates. These variables define the mapping coordinates used for the modeling and are checked by CALMET for consistency and compatibility with the CALMET control file inputs. Eight sets of flags and data records follow for the land use, elevation, surface parameters, and anthropogenic heat flux data. The default CALMET land use scheme is based on the U.S. Geological Survey (USGS) land use classification system. The USGS primary land use categories are shown in Table 8-6. Two Level I USGS categories (water and wetlands) are subdivided into subcategories. Along with the default CALMET land use, the default values of the other geophysical parameters for each land use type are also shown. The default land use classification scheme contains 14 land use types. Note that a negative value of land use by CALMET is used as a flag to indicate irrigated land. Irrigated land may be assigned a different Bowen ratio than unirrigated land, and the CALPUFF dry deposition module uses the irrigated land use flag in computing the effect of moisture stress on stomatal resistance. (If the land is irrigated, it is assumed that the vegetation is not moisture stressed.) CALMET allows a more detailed breakdown of land use or a totally different classification scheme to be used by providing the option for user-defined land use categories. Currently, up to 52 user-specified land use categories are allowed. An extended 52-class land use scheme based on the USGS Level I and Level II land use categories is shown in Table 8-7. The user can specify up to "MXLU" land use categories along with new values of the other geophysical parameters for each land use type. The parameter MXLU is specified in the CALMET parameter file (PARAMS.MET). CALMET contains an option, in which temperatures over water bodies such as the ocean or large lakes are calculated by using data from only those observation stations (SEA.DAT files, usually buoys) located in it, while only land stations (SURF.DAT file) will be used to calculate temperatures over the rest of the grid. The variables JWAT1 and JWAT2 in CALMET.INP Input Group #6 specify the range of land use categories defining the water body for which this land/water temperature scheme will be implemented. A 8-55

Section 8: CALMET Model Files range is specified to allow inclusion of multiple categories, for example "bay" and "ocean," in the definition of the water body. To disable the overwater option, JWAT1 and JWAT2 are set to values greater than the highest land use category listed in the GEO.DAT file. The default values of JWAT1 and JWAT2 are both 999, indicating the overwater interpolation scheme is not applied in default mode. Because the temperature of any grid cell whose land use is included in the range defined by JWAT1 and JWAT2 will be determined by a weighting of all overwater data (SEA#.DAT files), it is recommended that smaller or distant water bodies be assigned land use categories that are distinct from those used in JWAT1 and JWAT2, to avoid use of inappropriate data in determining their surface temperatures. Thus a small reservoir will have its temperature determined by surrounding land stations, rather than by ocean buoy data. After viewing the initial temperature field that results from the CALMET run, the user may wish to "fine tune" the fields using the extended, 52-class land use system in Table 8-7 and by altering the land use assignments of particular grid cells or changing the land uses included in the JWAT1-JWAT2 range. For instance, by limiting the range to "ocean" only and then changing which near-shore cells are considered to be "bay" and which are "ocean" the user can control the appearance of the temperature field in the vicinity of the coastline. The values of IWAT1 and IWAT2 (GEO.DAT Input File) are used to determine whether the overland or overwater method will be used to produce a mixing height value for a particular grid cell. The default values of IWAT1 and IWAT2 are both 55, restricting the overwater mixing height scheme to "large" bodies of water. The user may change the values of IWAT1 and IWAT2 on a case-by-case basis to include or exclude other water bodies from being considered as overwater. For instance, the user's domain may have a bay where the mixing height should be determined using the overwater method but a series of small lakes where the overland method would be more appropriate, so the "lake" category would be excluded from the IWAT range. Alternatively, if one has a large lake that should be considered to be "overwater" and a smaller lake that should be considered to be "overland", then the land use category for the smaller lake could be changed to reflect some other category not in the IWAT range, such as forest or wetland. It is recommended that if the user creates his or her own GEO.DAT fields for roughness length, albedo, etc., they be weighted by the actual percentage of each land use in a given cell. That method is more accurate and, if one subsequently changes the dominant land use category, the variables used to calculate mixing height will still reflect the fact that there is water present in the grid cell. The surface elevation data field is entered in "user units" along with a scaling factor to convert user units to meters. The sample GEO.DAT file shown in Table 8-5 contains elevations in meters. The gridded fields are entered with the 'NXM' values on a line. NXM is the number of grid cells in the X direction. The data from left to right correspond to X=1 through NXM. The top line of a gridded field correspond to Y=NYM, the next line to Y=NYM-1, etc. All of the GEO.DAT inputs are read in FORTRAN free format. A detailed description of the GEO.DAT variables is contained in Table 8-8. 8-56

Section 8: CALMET Model Files<br />

8.2 Geophysical D<strong>at</strong>a File (GEO.DAT)<br />

<strong>The</strong> GEO.DAT d<strong>at</strong>a file contains the geophysical d<strong>at</strong>a inputs required by the CALMET model. <strong>The</strong>se<br />

inputs include land use type, elev<strong>at</strong>ion, surface parameters (surface roughness, length, albedo, Bowen<br />

r<strong>at</strong>io, soil he<strong>at</strong> flux parameter, and veget<strong>at</strong>ion leaf area index) and anthropogenic he<strong>at</strong> flux. <strong>The</strong> land use<br />

and elev<strong>at</strong>ion d<strong>at</strong>a are entered as gridded fields. <strong>The</strong> surface parameters and anthropogenic he<strong>at</strong> flux can<br />

be entered either as gridded fields or computed from the land use d<strong>at</strong>a <strong>at</strong> each grid point. Default values<br />

rel<strong>at</strong>ing each <strong>of</strong> these parameters to land use are provided in the model.<br />

A sample GEO.DAT file is shown in <strong>Table</strong> 8-5. <strong>The</strong> first line <strong>of</strong> the file identifies the file and its form<strong>at</strong><br />

to CALMET. <strong>The</strong> second line provides the number <strong>of</strong> comment lines (character strings <strong>of</strong> up to 80<br />

characters in length) to follow th<strong>at</strong> are read but are not interpreted by the model. One <strong>of</strong> these may<br />

typically contain a title to identify the d<strong>at</strong>a set. <strong>The</strong> next block <strong>of</strong> 5 to 6 lines contains map projection,<br />

d<strong>at</strong>um, and grid inform<strong>at</strong>ion such as the number <strong>of</strong> grid cells, grid spacing, and reference coordin<strong>at</strong>es.<br />

<strong>The</strong>se variables define the mapping coordin<strong>at</strong>es used for the modeling and are checked by CALMET for<br />

consistency and comp<strong>at</strong>ibility with the CALMET control file inputs. Eight sets <strong>of</strong> flags and d<strong>at</strong>a records<br />

follow for the land use, elev<strong>at</strong>ion, surface parameters, and anthropogenic he<strong>at</strong> flux d<strong>at</strong>a.<br />

<strong>The</strong> default CALMET land use scheme is based on the U.S. Geological Survey (USGS) land use<br />

classific<strong>at</strong>ion system. <strong>The</strong> USGS primary land use c<strong>at</strong>egories are shown in <strong>Table</strong> 8-6. Two Level I<br />

USGS c<strong>at</strong>egories (w<strong>at</strong>er and wetlands) are subdivided into subc<strong>at</strong>egories. Along with the default<br />

CALMET land use, the default values <strong>of</strong> the other geophysical parameters for each land use type are also<br />

shown. <strong>The</strong> default land use classific<strong>at</strong>ion scheme contains 14 land use types. Note th<strong>at</strong> a neg<strong>at</strong>ive value<br />

<strong>of</strong> land use by CALMET is used as a flag to indic<strong>at</strong>e irrig<strong>at</strong>ed land. Irrig<strong>at</strong>ed land may be assigned a<br />

different Bowen r<strong>at</strong>io than unirrig<strong>at</strong>ed land, and the CALPUFF dry deposition module uses the irrig<strong>at</strong>ed<br />

land use flag in computing the effect <strong>of</strong> moisture stress on stom<strong>at</strong>al resistance. (If the land is irrig<strong>at</strong>ed, it<br />

is assumed th<strong>at</strong> the veget<strong>at</strong>ion is not moisture stressed.)<br />

CALMET allows a more detailed breakdown <strong>of</strong> land use or a totally different classific<strong>at</strong>ion scheme to be<br />

used by providing the option for user-defined land use c<strong>at</strong>egories. Currently, up to 52 user-specified land<br />

use c<strong>at</strong>egories are allowed. An extended 52-class land use scheme based on the USGS Level I and Level<br />

II land use c<strong>at</strong>egories is shown in <strong>Table</strong> 8-7. <strong>The</strong> user can specify up to "MXLU" land use c<strong>at</strong>egories<br />

along with new values <strong>of</strong> the other geophysical parameters for each land use type. <strong>The</strong> parameter MXLU<br />

is specified in the CALMET parameter file (PARAMS.MET).<br />

CALMET contains an option, in which temper<strong>at</strong>ures over w<strong>at</strong>er bodies such as the ocean or large lakes<br />

are calcul<strong>at</strong>ed by using d<strong>at</strong>a from only those observ<strong>at</strong>ion st<strong>at</strong>ions (SEA.DAT files, usually buoys) loc<strong>at</strong>ed<br />

in it, while only land st<strong>at</strong>ions (SURF.DAT file) will be used to calcul<strong>at</strong>e temper<strong>at</strong>ures over the rest <strong>of</strong> the<br />

grid. <strong>The</strong> variables JWAT1 and JWAT2 in CALMET.INP Input <strong>Group</strong> #6 specify the range <strong>of</strong> land use<br />

c<strong>at</strong>egories defining the w<strong>at</strong>er body for which this land/w<strong>at</strong>er temper<strong>at</strong>ure scheme will be implemented. A<br />

8-55

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!