Warning: Undefined array key 1 in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 78

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 431

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
solvers:grid2dinterpolator [Elmer/Ice Wiki]

This is an old revision of the document!


Grid2DInterpolator

General Description

This solver interpolates data given on a regular 2D regular grid in an ASCII file (x y Value). A bilinear interpolation is used. The data in the ASCII file have to be ordered such that

  x1 y1 val11
  x2 y1 val21
  ...
  xn y1 valn1
  x1 y2 val12
  ...
  xn yn valnn 

The grid is described by giving:

  • (x0, y0) the left-bottom corner coordinate
  • (lx, ly) the x and y lengths of the covered domain
  • (Nx, Ny) the number of cells in x and y directions
  • No data are given by -9999

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. The solver gives the number of nodes outside of the DEM domain and the farthest DEM point used to evaluate the nodal values.

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, respectively.

Solver 1
  Exec Solver = Before Simulation
  Equation = "Read DEM"
  Procedure = "ElmerIceSolvers" "Grid2DInterpolator"

  Variable 1 = String "ZsDEM"
  Variable 1 data file = File "DEM_TR_surf.dat"
  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 2 = String "bedrockDEM"
  Variable 2 data file = File "DEM_TR_bed.dat"
  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

The variables bedrockDEM and ZsDEM have to exported from an other solver

Solver 3
  Equation = "Navier-Stokes"
  Exported Variable 1 = -dofs 1 "ZsDEM"
  Exported Variable 2 = -dofs 1 "bedrockDEM"
  ...
End

The variables bedrockDEM and ZsDEM can then be used to estimate the Bottom Surface and Top Surface keywords used by the StructuredMeshMapper solver (have a look here). Here, a minimum ice thickness of 1 meter is imposed.

! 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)>tx(1)+1.0) {tx(0)} else {tx(1)+1.0}"
End
solvers/grid2dinterpolator.1367004851.txt.gz · Last modified: 2013/04/26 19:34 by gag
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0