Model Domain and Grid  

Model_Domain_and_Grid_dfs2.jpg

     

Model_Domain_and_Grid_shp.jpg

     

 

Model Domain and Grid

dialogue Type:

Special; Integer Grid Codes

EUM Data Units

Integer Grid Code

Valid Values

-1e-35, 1, 2

    

One of the first steps in your model development is to define the model area. On a catchment scale, the model boundary is typically a topographic divide, a groundwater divide or some combination of the two. In general, there are no constraints on the definition of the model boundaries. However, the model boundaries should be chosen carefully, keeping in mind the boundary condi­tions that will be used for both the surface water and groundwater compo­nents.

Any non-gridded data, such as point .shp files and xyz ASCII files are interpo­lated to the grid defined in this dialogue.

The grid defined in this dialogue is the primary grid. Other .dfs2 gridded data does not have to use the same grid. However, if another .dfs2 file uses a dif­ferent grid then Real data is interpolated. If the two grids are coincident, that is the cells are the same size and the grids line up, then the data is bilinearly interpolated to the Model Grid. If the grids are not coincident, then the data is treated as if it were point values (i.e. the same as XYZ or .shp data). Integer Grid Code must use coincident grids, as it is impossible to interpolated inte­ger values.

The model domain and grid can be rotated any angle measured clockwise from the x-axis.

Using a dfs2 file

If you define your model domain using a dfs2 grid file, then you must define the cell values as follows:

·         Grid cells outside of the model domain must be assigned a delete value - usually -1e-35.

·         Grid cells inside the model domain must be assigned a value of 1.

·         Grid cells on the model boundary must be assigned a value of 2.

This distinction between interior grid cells and boundary cells is to facilitate the definition of boundary conditions. For example, drainage flow can be routed to external boundaries but not to internal boundaries.

The catchment definition is displayed in the greyed out text boxes but is not editable, since the catchment definition is part of the .dfs2 file format. If you want to change the cell size, origin, number of cells etc., you must change the .dfs2 file itself. For more information on editing and setting up the Model Domain and Grid see Conceptual geologic model for the finite difference approach (V1 p. 59).

Using a shp file

It is much easier to define your Model Domain and Grid via an ArcView .shp file (i.e a grid independent polygon). In this case, the definition of integer code values is taken care of automatically. Further, the definition of the grid (num­ber of rows and columns, cell size and origin) can be easily adjusted.

Map projections

A projection is the way the earth’s curved surface is mapped onto a flat plane. The projection causes distortions as you move away from the center of the projection. One of the most common projection systems is the UTM system, which divides the earth into zones to minimize the distortion within a particu­lar zone. Every place on earth has a UTM zone.

The earth is not a sphere and the actual curvature varies locally across the globe. Thus, the other part of the projection is the asssumed shape of the earth, which is defined as the Ellipsoid type. The ellipsoid used for the UTM system is WGS 84.

Thus, unless the projection and the ellipsoid are consistent between two maps, they will not overlap correctly. This is typically the reason if your model results do not appear to be located in the right place when you try to map them to Google Earth.

Map projections in MIKE SHE

Up until the 2009 Release, the only map projection allowed in MIKE SHE was NON-UTM. This is essentially the same as local grid coordinates. From the 2009 Release onwards, MIKE SHE uses the same map projection utilities as the rest of the MIKE Zero products, which is based on a installed library of standard world map projections.  

You must define the map projection when you create a dfs2 file and this pro­jection becomes part of the properties of the file. It can be changed in the Properties dialogue in the Grid Editor. Likewise, the map projection is a prop­erty of all shape files and can be changed in your GIS program.

note.jpg 

Note: In the Model Domain and Grid, you must define the map projection. All of the rest of the dfs2 and shape files used in the model setup must use the same projection.