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 18:26] gag |
solvers:grid2dinterpolator [2016/05/16 08:39] (current) tzwinger [General Description] |
||
---|---|---|---|
Line 1: | Line 1: | ||
===== Grid2DInterpolator ===== | ===== Grid2DInterpolator ===== | ||
- | This solver interpolates data given on a regular 2D regular | + | |
+ | ==== General Description ==== | ||
+ | |||
+ | 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 | ||
x2 y1 val21 | x2 y1 val21 | ||
Line 8: | Line 12: | ||
x1 y2 val12 | x1 y2 val12 | ||
... | ... | ||
- | xn yn valnn | + | xn ym valnm |
+ | </ | ||
- | The grid is describe | + | 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 | ||
* (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 by value (within a given tolerance) |
+ | |||
+ | WARNING: All data with values between (noData - tol) and (noData + tol) will be ignored (where noData and tol are the no data value and tolerance, see below). | ||
+ | |||
+ | ==== 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, | ||
+ | < | ||
+ | Solver 1 | ||
+ | Exec Solver = Before Simulation | ||
+ | Equation = "Read DEM" | ||
+ | Procedure = " | ||
+ | |||
+ | Variable 1 = String " | ||
+ | Variable 1 data file = File " | ||
+ | Variable 1 x0 = REal 947700.0d0 | ||
+ | Variable 1 y0 = REal 2104850.0d0 | ||
+ | Variable 1 lx = REal 800.0 | ||
+ | Variable 1 ly = REal 350.0 | ||
+ | Variable 1 Nx = Integer 268 | ||
+ | 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 data file = File " | ||
+ | Variable 2 x0 = REal 947700.0d0 | ||
+ | Variable 2 y0 = REal 2104850.0d0 | ||
+ | Variable 2 lx = REal 600.0 | ||
+ | Variable 2 ly = REal 350.0 | ||
+ | Variable 2 Nx = Integer 301 | ||
+ | Variable 2 Ny = Integer 176 | ||
+ | End | ||
+ | </ | ||
+ | |||
+ | Note that the "no data" and "no data tol" keywords are optional for each variable. | ||
+ | |||
+ | The variables '' | ||
+ | < | ||
+ | Solver 3 | ||
+ | Equation = " | ||
+ | Exported Variable 1 = -dofs 1 " | ||
+ | Exported Variable 2 = -dofs 1 " | ||
+ | ... | ||
+ | End | ||
+ | </ | ||
+ | The variables '' | ||
+ | < | ||
+ | ! Bedrock BC | ||
+ | Boundary Condition 2 | ||
+ | Bottom Surface = Equals bedrockDEM | ||
+ | ... | ||
+ | End | ||
+ | |||
+ | ! Upper Surface | ||
+ | Boundary Condition 3 | ||
+ | Top Surface = Variable ZsDEM, bedrockDEM | ||
+ | Real MATC "if (tx(0)> | ||
+ | End | ||
+ | </ | ||
+ | |||
+ | ==== SIF contents ==== | ||
+ | An example using the Tete Rousse surface and bedrock DEM can be found in '' |