This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:grid2dinterpolator [2013/04/26 19:13] gag [General Description] |
solvers:grid2dinterpolator [2016/05/16 08:39] (current) tzwinger [General Description] |
||
---|---|---|---|
Line 4: | Line 4: | ||
==== General Description ==== | ==== General Description ==== | ||
- | This solver interpolates data given on a regular 2D regular | + | This solver interpolates data given on a regular 2D grid in an ASCII file (x y Value). A bilinear interpolation is used. By default, the data in the ASCII file have to be ordered such that |
< | < | ||
x1 y1 val11 | x1 y1 val11 | ||
Line 12: | Line 12: | ||
x1 y2 val12 | x1 y2 val12 | ||
... | ... | ||
- | xn yn valnn | + | xn ym valnm |
</ | </ | ||
+ | |||
+ | If the second column (the y-value) is changing faster, i.e., | ||
+ | < | ||
+ | x1 y1 val11 | ||
+ | x1 y2 val12 | ||
+ | ... | ||
+ | x1 ym val1m | ||
+ | ... | ||
+ | xn y1 valn1 | ||
+ | xn y2 valn2 | ||
+ | ... | ||
+ | xn ym valnm | ||
+ | </ | ||
+ | |||
+ | one can give the keyword | ||
+ | < | ||
+ | | ||
+ | </ | ||
+ | |||
+ | In any case, it is essential, that the order '' | ||
+ | |||
+ | Further, if there are those points with no-data simply missing from the file (saves space), the keyword | ||
+ | < | ||
+ | | ||
+ | </ | ||
+ | will insert either the default no-data value or the one given with the '' | ||
The grid is described by giving: | The grid is described by giving: | ||
* (x0, y0) the left-bottom corner coordinate | * (x0, y0) the left-bottom corner coordinate | ||
* (lx, ly) the x and y lengths of the covered domain | * (lx, ly) the x and y lengths of the covered domain | ||
- | * (Nx, Ny) the number of cells in x and y directions | + | * (Nx, Ny) the number of levels |
- | * No data are given by -9999 | + | * Cells with no data are identified |
+ | |||
+ | WARNING: All data with values between (noData | ||
- | WARNING: In the case a node fails in a cell of the DEM containing no data, the nodal value is assigned to the closest DEM value. This robust feature might be dangerous if too much nodes are outside of the DEM domain. | ||
==== SIF contents ==== | ==== SIF contents ==== | ||
Add this solver and execute it before the simulation. Here, it is used to read the surface DEM and the bedrock DEM given by the files DEM_TR_surf.dat and DEM_TR_bed.dat, | Add this solver and execute it before the simulation. Here, it is used to read the surface DEM and the bedrock DEM given by the files DEM_TR_surf.dat and DEM_TR_bed.dat, | ||
Line 38: | Line 65: | ||
Variable 1 Nx = Integer 268 | Variable 1 Nx = Integer 268 | ||
Variable 1 Ny = Integer 118 | Variable 1 Ny = Integer 118 | ||
+ | Variable 1 no data = Real -999.0 | ||
+ | Variable 1 no data tol = Real 0.1 | ||
Variable 2 = String " | Variable 2 = String " | ||
Line 49: | Line 78: | ||
End | End | ||
</ | </ | ||
+ | |||
+ | Note that the "no data" and "no data tol" keywords are optional for each variable. | ||
The variables '' | The variables '' | ||
Line 74: | Line 105: | ||
</ | </ | ||
+ | ==== SIF contents ==== | ||
+ | An example using the Tete Rousse surface and bedrock DEM can be found in '' |