FD2BASIS

Basic Information

Name of Program

FD2BASIS

Version-Date

April 1998

Description-Date

April 1998

Catchwords

finite differences methods
grid generation
distribution of land and water

Short Description of Functionality

The program FD2BASIS generates a FIDISOR/FIDIRB grid with correct distribution of land and water areas. The latter is constructed from digitized boundary structure lines like the high water boundary ( simulations of normal tidal situations) or dike line (for storm surge simulations).

Input-Files

1. input control file (type fd2basis.dat)
2. one or more files with digitized boundary structures (type digi.gkk)

Output-Files

1. basic FIDISOR/FIDIRB grid (type topo.bin)
2. file randpkte.dat with boundary cell indices and the intersecting polygon part (is used internally for boundary fit).
3. buffer files liste_1.dat and liste_2.dat
4. printer file fd2basis.sdr
5. (optional) trace file fd2basis.trc

Methodology

First a matrix IPOLY is allocated with given grid dimension which contains 0 overall. Then all cells which are touched by the given boundary structures are set to the values 99 or 98 (fix point or interpolated point respectively). Now land is initialized with -99, beginning with the landpoints given by the user. If a cell is reached, which has a value different from 0 the search continues in another main direction or at the next landpoint. All initialized land cells become starting points itself. For purpose of memory saving parts of the starting points are buffered on disk.

Last step is the correction of the boundary points. If the area of land inside the cell is greater than that of water, the cell is set to land. The area is calculated from the intersection polygon, which was stored in randpkte.dat.

Remarks:

• The structures can be stored as parts and there is no need of sorting before using this program. The only constraint is that they must in total represent closed structures (islands) or their beginning and ending shall both lie outside the grid area. Holes in a boundary can lead to unwanted results!
• A structure which contains only one point is ignored. Double points following each other are used like one point.
• A manual postprocessing can be neccessary for thin areas of water or land (land connections, small channels) if resolution of the grid is not high enough to resolve theese structures.

Program(s) to run before this Program

digitizer program (PC)

Program(s) to run after this Program

FD2ADDTOPO, FDGITTER05

Additional Information

Language

Fortran90

Additional software

-

Original Version

G. Seiß

Maintenance

G. Seiß

Documentation/Literature

$PROGHOME/examples/fd2basis/*

---

back to Program Descriptions

---

Overview